Emergent behavior of Cucker-Smale model with normalized weights and distributed time delays

TL;DR

This paper analyzes a flocking model with normalized weights and distributed delays, establishing conditions for velocity alignment, approximation by a Vlasov equation, and global solutions with long-term behavior.

Contribution

It introduces a novel analysis of a Cucker-Smale model with distributed delays and normalized weights, including new conditions for flocking and a rigorous continuum limit.

Findings

Velocity alignment conditions derived

Approximation by delayed Vlasov equation validated

Global existence and asymptotic behavior established

Abstract

We study a Cucker-Smale-type flocking model with distributed time delay where individuals interact with each other through normalized communication weights. Based on a Lyapunov functional approach, we provide sufficient conditions for the velocity alignment behavior. We then show that as the number of individuals N tends to infinity, the N-particle system can be well approximated by a delayed Vlasov alignment equation. Furthermore, we also establish the global existence of measure-valued solutions for the delayed Vlasov alignment equation and its large-time asymptotic behavior.

Equations

$${\tau }^{\prime }(t)\le 0\text{mbox}and\tau (t)\ge {\tau }_{\ast },\text{mbox}fort\ge 0,$$

$${\tau }_{\ast }\le \tau (t)\le {\tau }_0\text{mbox}fort\ge 0.$$

$$\begin{array}{rl}\frac{d{x}_i(t)}{dt}& =v_i(t),i=1,\cdots ,N,t>0,\\ \frac{d{v}_i(t)}{dt}& =\frac{1}{h(t)}\sum _{k=1}^N{\int }_{t-\tau (t)}^t\alpha (t-s)\varphi (x_k(s),x_i(t))(v_k(s)-v_i(t))ds,\end{array}$$

$$\phi(x_{k}(s),x_{i}(t))=\left\{\begin{array}[]{ll}\displaystyle\frac{\psi(|x_{k}(s)-x_{i}(t)|)}{\sum_{j\neq i}\psi(|x_{j}(s)-x_{i}(t)|)}&\textrm{if $k\neq i$,}\\[11.38109pt] 0&\textrm{if $k=i$,}\end{array}\right.$$

$${\int }_0^{{\tau }_{\ast }}\alpha (s)ds>0,$$

$$h(t):={\int }_0^{\tau (t)}\alpha (s)ds,t\ge 0.$$

$$x_i(s)=:x_i^0(s),v_i(s)=:v_i^0(s),i=1,\cdots ,N,s\in [-{\tau }_0,0],$$

$${\partial }_t{f}_t+v\cdot {\nabla }_x{f}_t+{\nabla }_v\cdot \left(\frac{1}{h(t)}{\int }_{t-\tau (t)}^t\alpha (t-s)F[f_s]ds{f}_t\right)=0,$$

$$F[f_s](x,v):=\frac{{\int }_{{\mathbb{R}}^d\times {\mathbb{R}}^d}\psi (|x-y|)(w-v)f_s(y,w)dydw}{{\int }_{{\mathbb{R}}^d\times {\mathbb{R}}^d}\psi (|x-y|)f_s(y,w)dydw}.$$

$$d_X(t):=\underset{1\le i,j\le N}{\max }|x_i(t)-x_j(t)|\text{mbox}and{d}_V(t):=\underset{1\le i,j\le N}{\max }|v_i(t)-v_j(t)|.$$

$$\underset{t\ge 0}{\sup }{d}_X(t)<\infty \text{mbox}and\underset{t\to \infty }{\lim }{d}_V(t)=0.$$

$$R_v:=\underset{s\in [-{\tau }_0,0]}{\max }\underset{1\le i\le N}{\max }|v_i^0(s)|.$$

$$h(0)d_V(0)+{\int }_0^{{\tau }_0}\alpha (s)\left({\int }_{-s}^0{d}_V(z)dz\right)ds<{\beta }_N{\int }_0^{{\tau }_{\ast }}\alpha (s){\int }_{d_X(-s)+R_v{\tau }_0}^{\infty }\psi (z)dzds,$$

$$\underset{t\ge 0}{\sup }{d}_X(t)<\infty$$

$$d_V(t)\le \underset{s\in [-{\tau }_0,0]}{\max }{d}_V(s)e^{-\gamma t}\text{mbox}fort\ge 0,$$

$$\frac{d{x}_i(t)}{dt}$$

$$\frac{d{v}_i(t)}{dt}$$

$$d_V(0)+{\int }_{-\tau }^0{d}_V(z)dz<{\beta }_N{\int }_{d_X(-\tau )+R_v\tau }^{\infty }\psi (z)dz.$$

$$\underset{1\le i\le N}{\max }|v_i(t)|\le R_v\text{mbox}fort\ge -{\tau }_0.$$

$$S^{ϵ}:=\{t>0:\underset{1\le i\le N}{\max }|v_i(s)|<R_v+ϵ\forall s\in [0,t)\}.$$

$$\underset{1\le i\le N}{\max }|v_i(T^{ϵ})|=R_v+ϵ.$$

$$\begin{array}[]{l}\displaystyle{\frac{d|v_{i}(t)|^{2}}{dt}\leq\frac{2}{h(t)}\sum_{k=1}^{N}\int_{t-\tau(t)}^{t}\alpha(t-s)\phi(x_{k}(s),x_{i}(t))(|v_{k}(s)|\,|v_{i}(t)|-|v_{i}(t)|^{2})\,ds}\\ \displaystyle{\hskip 42.67912pt\leq\frac{2}{h(t)}\sum_{k=1}^{N}\int_{t-\tau(t)}^{t}\alpha(t-s)\phi(x_{k}(s),x_{i}(t))\max_{1\leq k\leq N}|v_{k}(s)|\,|v_{i}(t)|\,ds-2|v_{i}(t)|^{2}}\,.\end{array}$$

$$\underset{1\le k\le N}{\max }|v_k(s)|\le R_v+ϵ\text{mbox}and\sum _{k=1}^N\varphi (x_k(s),x_i(t))=1,$$

$$\frac{d|v_i(t){|}^2}{dt}\le 2[(R_v+ϵ)|v_i(t)|-|v_i(t){|}^2],$$

$$\frac{d|v_i(t)|}{dt}\le (R_v+ϵ)-|v_i(t)|.$$

$$\underset{t\to {T^{ϵ}}^{-}}{\lim }\underset{1\le i\le N}{\max }|v_i(t)|\le e^{-T^{ϵ}}\left(\underset{1\le i\le N}{\max }|v_i(0)|-R_v-ϵ\right)+R_v+ϵ<R_v+ϵ.$$

$$D^{+}F(t):=\underset{h\to 0+}{\mathrm{limsup}}\frac{F(t+h)-F(t)}{h}.$$

$$\left|D^{+}{d}_X(t)\right|$$

$$D^{+}{d}_V(t)$$

$$\bigcup _{\sigma \in \mathbb{N}}\overline{{\mathcal{I}}_{\sigma }}=[0,\infty )$$

$$d_V(t)=|v_{i(\sigma )}(t)-v_{j(\sigma )}(t)|\text{mbox}fort\in {\mathcal{I}}_{\sigma }.$$

Full text

1. Introduction

In the last years the study of collective behavior of multi-agent systems has attracted the interest of many researchers in different scientific fields, such as biology, physics, control theory, social sciences, economics. The celebrated Cucker-Smale model has been proposed and analyzed in [21, 22] to describe situations in which different agents, e.g. animals groups, reach a consensus (flocking), namely they align and move as a flock, based on a simple rule: each individual adjusts its velocity taking into account other agents’ velocities.

In the original papers a symmetric interaction potential is considered. Then, the case of non-symmetric interactions has been studied by Motsch and Tadmor [31]. Several generalizations and variants have been introduced to cover various applications’ fields, e.g. more general interaction rates and singular potentials [8, 10, 19, 27, 30, 32], cone-vision constraints [40], presence of leadership [17, 37], noise terms [20, 23, 25], crowds dynamics [18, 29], infinite-dimensional models [1, 2, 11, 26, 28, 39], control problems [3, 5, 16, 33]. We refer to [6, 12] for recent surveys on the Cucker-Smale type flocking models and its variants.

It is natural to introduce a time delay in the model, as a reaction time or a time to receive environmental information. The presence of a time delay makes the problem more difficult to deal with. Indeed, the time delay destroys some symmetry features of the model which are crucial in the proof of convergence to consensus. For this reason, in spite of a great amount of literature on Cucker-Smale models, only a few papers are available concerning Cucker-Smale model with time delay [13, 14, 15, 23, 36]. Cucker-Smale models with delay effects are also studied in [34, 35] when a hierarchical structure is present, namely the agents are ordered in a specific order depending on which other agents they are leaders of or led by.

Here we consider a distributed delay term, i.e. we assume that the agent $i,$ $i=1,\dots,N,$ changes its velocity depending on the information received from other agents on a time interval $[t-\tau(t),t].$ Moreover, we assume normalized communication weights (cf [13]). Let us consider a finite number $N\in\mathbb{N}$ of autonomous individuals located in $\mathbb{R}^{d}$, $d\geq 1$. Let $x_{i}(t)$ and $v_{i}(t)$ be the position and velocity of $i$th individual. Then, in the current work, we will deal with a Cucker-Smale model with distributed time delays. Namely, let $\tau:[0,+\infty)\rightarrow(0,+\infty),$ be the time delay function belonging to $W^{1,\infty}(0,+\infty)$. Throughout this paper, we assume that the time delay function $\tau$ satisfies

[TABLE]

for some positive constant $\tau_{*}.$ Then, denoting $\tau_{0}:=\tau(0),$ we have

[TABLE]

It is clear that the constant time delay $\tau(t)\equiv\bar{\tau}>0$ satisfies the conditions above.

Our main system is given by

[TABLE]

where $\phi(x_{k}(s),x_{i}(t))$ are the normalized communication weights given by

[TABLE]

with the influence function $\psi:[0,\infty)\to(0,\infty).$ Throughout this paper, we assume that the influence function $\psi$ is bounded, positive, nonincreasing and Lipschitz continuous on $[0,\infty)$, with $\psi(0)=1$. Moreover, $\alpha:[0,\tau_{0}]\to[0,\infty)$ is a weight function satisfying

[TABLE]

and

[TABLE]

We consider the system subject to the initial datum

[TABLE]

i.e., we prescribe the initial position and velocity trajectories $x^{0}_{i},v^{0}_{i}\in\mathcal{C}([-\tau_{0},0];\mathbb{R}^{d})$.

For the particle system (1.3), we will first discuss the asymptotic behavior of solutions in Section 2. Motivated from [13, 26, 31], we derive a system of dissipative differential inequalities, see Lemma 2.2, and construct a Lyapunov functional. This together with using Halanay inequality enables us to show the asymptotic velocity alignment behavior of solutions under suitable conditions on the initial data.

We next derive, analogously to [13] where the case of a single pointwise time delay is considered, a delayed Vlasov alignment equation from the particle system (1.3) by sending the number of particles $N$ to infinity:

[TABLE]

where $f_{t}=f_{t}(x,v)$ is the one-particle distribution function on the phase space $\mathbb{R}^{d}\times\mathbb{R}^{d}$ and the velocity alignment force $F$ is given by

[TABLE]

We show the global-in-time existence and stability of measure-valued solutions to (1.6) by employing the Monge-Kantorowich-Rubinstein distance. As a consequence of the stability estimate, we discuss a mean-field limit providing a quantitative error estimate between the empirical measure associated to the particle system (1.3) and the measure-valued solution to (1.6). We then extend the estimate of large behavior of solutions for the particle system (1.3) to the one for the delayed Vlasov alignment equation (1.6). For this, we use the fact that the estimate of large-time behavior of solutions to the particle system (1.3) is independent of the number of particles. By combining this and the mean-field limit estimate, we show that the diameter of velocity-support of solutions of (1.6) converges to zero as time goes to infinity. Those results will be proved in Section 3.

2. Asymptotic flocking behavior of the particle system (1.3)

We start with presenting a notion of flocking behavior for the system (1.3), and for this we introduce the spatial and, respectively, velocity diameters as follows:

[TABLE]

Definition 2.1**.**

We say that the system with particle positions $x_{i}(t)$ and velocities $v_{i}(t)$, $i=1,\dots,N$ and $t\geq 0$, exhibits asymptotic flocking if the spatial and velocity diameters satisfy

[TABLE]

We then state our main result in this section on the asymptotic flocking behavior of the system (1.3).

Theorem 2.1**.**

*Assume $N>2.$ Suppose that the initial data $x_{i}^{0},v_{i}^{0}$ are continuous on the time interval $[-\tau_{0},0]$ and denote *

[TABLE]

*Moreover, denoted $\beta_{N}=\frac{N-2}{N-1}$, assume that *

[TABLE]

where $d_{X}$ and $d_{V}$ denote, respectively, the spatial and velocity diameters defined in (2.1). Then the solution of the system (1.3)–(1.5) is global in time and satisfies

[TABLE]

and

[TABLE]

for a suitable positive constant $\gamma$ independent of $t$. The constant $\gamma$ can be chosen also independent of $N.$

Remark 2.1**.**

If the influence function $\psi$ is not integrable, i.e., it has a heavy tail, then the right hand side of (2.3) becomes infinite, and thus the assumption (2.3) is satisfied for all initial data and time delay $\tau(t)$ satisfying (1.1). This is reminiscent of the unconditional flocking condition for the Cucker-Smale type models, see [11, 12, 13]. On the other hand, if the influence function $\psi$ is integrable and the weight function $\alpha(s)$ is given by $\alpha(s)=\delta_{\tau}(s)$, then the system (1.3) with a constant time delay $\tau(t)\equiv\tau$ becomes the one with discrete time delay:

[TABLE]

Note the above system is studied in [13]. For this system, the assumption (2.3) reduces to

[TABLE]

Since $d_{V}(z)\geq 0$, the left hand side of the above inequality increases as the size of time delay $\tau>0$ increases. On the other hand, the integrand on the right hand side decreases as $\tau$ increases. This asserts that small size of the time delay provides a larger set of initial data. **

Remark 2.2**.**

Observe that our theorem above gives a flocking result when the number of agents $N$ is greater than two. The result for $N=2$ is trivial for our normalized model.

Remark 2.3**.**

Note that $0<\beta_{N}\uparrow 1$ as $N\to\infty.$ Then, this implies that for any $\beta\in(0,1)$ there exists $N_{*}\in\mathbb{N}$ such that $\beta_{N}\geq\beta$ for $N\geq N_{*}$.**

For the proof of Theorem 2.1, we will need several auxiliary results. Inspired by [13], we first show the uniform-in-time bound estimate of the maximum speed of the system (1.3).

Lemma 2.1**.**

Let $R_{v}>0$ be given by (2.2) and let $(\mathbf{x},\mathbf{v})$ be a local-in-time $\mathcal{C}^{1}$-solution of the system (1.3)–(1.5). Then the solution is global in time and satisfies

[TABLE]

Proof.

Let us fix $\epsilon>0$ and set

[TABLE]

By continuity, $S^{\epsilon}\neq\emptyset.$ Let us denote $T^{\epsilon}:=\sup S^{\epsilon}>0.$ We want to show that $T^{\epsilon}=+\infty.$ Arguing by contradiction, let us suppose $T^{\epsilon}$ finite. This gives, by continuity,

[TABLE]

From (1.3)–(1.5), for $t<T^{\epsilon}$ and $i=1,\dots,N,$ we have

[TABLE]

Note that

[TABLE]

for $s\in[t-\tau(t),t]$ with $t<T^{\epsilon}$. Thus we obtain

[TABLE]

which gives

[TABLE]

From (2.5) we obtain

[TABLE]

This is in contradiction with (2.4). Therefore, $T^{\epsilon}=+\infty\,.$ Being $\epsilon>0$ arbitrary, the claim is proved. ∎

In the lemma below, motivated from [13, 38] we derive the differential inequalities for $d_{X}$ and $d_{V}$. We notice that the diameter functions $d_{X}$ and $d_{V}$ are not $\mathcal{C}^{1}$ in general. Thus we introduce the upper Dini derivative to consider the time derivative of these functions: for a given function $F=F(t)$, the upper Dini derivative of $F$ at $t$ is defined by

[TABLE]

Note that the Dini derivative coincides with the usual derivative when the function is differentiable at $t$.

Lemma 2.2**.**

Let $R_{v}>0$ be given by (2.2) and let $(\mathbf{x},\mathbf{v})$ be the global $\mathcal{C}^{1}$-solution of the system (1.3)–(1.5) constructed in Lemma 2.1. Then, for almost all $t>0$, we have

[TABLE]

Proof.

The first inequality is by now standard, Then, we concentrate on the second one. Due to the continuity of the velocity trajectories $v_{i}(t),$ there is an at most countable system of open, mutually disjoint intervals $\{\mathcal{I}_{\sigma}\}_{\sigma\in\mathbb{N}}$ such that

[TABLE]

and thus for each ${\sigma\in\mathbb{N}}$ there exist indices $i(\sigma)$, $j(\sigma)$ such that

[TABLE]

Then, by using the simplified notation $i:=i(\sigma)$, $j:=j(\sigma)$ (we may assume that $i\neq j$), we have for every $t\in\mathcal{I}_{\sigma}$,

[TABLE]

Set

[TABLE]

Note that, from the definition (1.4) of $\phi(x_{k}(s),x_{i}(t)),$ it results $\bar{\phi}_{ij}(s,t)<1$ for all $i\neq j,$ $1\leq i,j\leq N$ and for all $t>0,\ s\in[t-\tau(t),t].$ Then, we get

[TABLE]

where

[TABLE]

Observe that $a_{ij}^{k}(s,t)\geq 0$ and $\sum_{k=1}^{N}a_{ij}^{k}(s,t)=1.$ Thus, if we consider the convex hull of a finite velocity point set $\{v_{1}(t),\dots,v_{N}(t)\}$ and denote it by $\Omega(t),$ then

[TABLE]

This gives

[TABLE]

which, used in (2.7), implies

[TABLE]

Now, by using the first equation in (1.3), we estimate for any $1\leq i,k\leq N$ and $s\in[t-\tau(t),t]$,

[TABLE]

Then, Lemma 2.1 gives

[TABLE]

and due to the monotonicity property of the influence function $\psi,$ for $k\neq i,$ we deduce

[TABLE]

On the other hand, we find

[TABLE]

Then, from (2.9), we obtain

[TABLE]

Using the last estimate in (2.8), we have

[TABLE]

that, used in (2.6), concludes the proof. ∎

Lemma 2.3**.**

Let $u$ be a nonnegative, continuous and piecewise $\mathcal{C}^{1}$-function satisfying, for some constant $0<a<1$, the differential inequality

[TABLE]

Then we have

[TABLE]

with $\gamma=1-\tau_{0}^{-1}H(a\tau_{0}\exp(\tau_{0}))$, where $H$ is the product logarithm function, i.e., $H$ satisfies $z=H(z)\exp(H(z))$ for any $z\geq 0$.

Proof.

Note that the differential inequality (2.10) implies

[TABLE]

Then, the result follows from Halanay inequality (see e.g. [24, p. 378]). ∎

We are now ready to proceed with the proof of Theorem 2.1.

Proof of Theorem 2.1.

For $t>0$, we introduce the following Lyapunov functional for the system (1.3)–(1.5):

[TABLE]

where $R_{v}$ is given by (2.2) and the diameters $d_{X}(t)$, $d_{V}(t)$ are defined in (2.1). We first estimate $\mathcal{L}_{1}$ as

[TABLE]

for almost all $t\geq 0$, due to (1.1) and the first inequality in Lemma 2.2. Here $\mbox{sgn\,}(\cdot)$ is the signum function defined by

[TABLE]

Analogously, we also get

[TABLE]

for almost all $t\geq 0$. Then, from Lemma 2.2, (2.11) and (2.13), we have, for almost all $t>0$,

[TABLE]

On the other hand, since $h^{\prime}(t)=\alpha(\tau(t))\tau^{\prime}(t)\leq 0$, we have $D^{+}\left(h(t)\mathcal{L}(t)\right)\leq 0,$ namely

[TABLE]

Moreover, it follows from the assumption (2.3) that there exists a positive constant $d_{*}>d_{X}(-s)+R_{v}\tau_{0}$ such that

[TABLE]

This, together with (2.14) and (1.2), implies

[TABLE]

Now, observe that, if $d_{X}(t-s)\geq d_{X}(-s),$ then

[TABLE]

Similarly, when $d_{X}(t-s)<d_{X}(-s),$ then

[TABLE]

Thus, from (2.15), (2.16) and (2.17) we deduce that

[TABLE]

Note that, for $s\in[0,\tau(t)]$,

[TABLE]

due to Lemma 2.1 and the first inequality of Lemma 2.2. Analogously, we also find for $s\in[0,\tau(t)]$ that

[TABLE]

This gives

[TABLE]

Thus we get

[TABLE]

Combining this and (2.15), we obtain

[TABLE]

Since the left hand side of the above inequality is positive, we have

[TABLE]

We then again use (2.18) to find

[TABLE]

for $s\in[0,\tau(t)]$ and $t\geq 0$. Hence, by Lemma 2.2 together with the monotonicity of $\psi$ we have

[TABLE]

for almost all $t>0$, where $\psi_{*}=\psi(d_{*}+4R_{v}\tau_{0})$. We finally apply Lemma 2.3 to complete the proof. Note that, in order to have an exponential decay rate of $d_{V}$ independent of $N,$ it is sufficient to observe that $\beta_{N}\geq 1/2$ for $N\geq 3.$ ∎

3. A delayed Vlasov alignment equation

In this section, we are interested in the behavior of solutions to the particle system (1.3) as the number of particles $N$ goes to infinity. At the formal level, we can derive the following delayed Vlasov alignment equation from (1.3) as $N\to\infty$:

[TABLE]

with the initial data:

[TABLE]

where $f_{t}=f_{t}(x,v)$ is the one-particle distribution function on the phase space $\mathbb{R}^{d}\times\mathbb{R}^{d}$. Here the interaction term $F[f_{s}]$ is given by

[TABLE]

For the equation (3.1), we provide the global-in-time existence and uniqueness of measure-valued solutions and mean-field limits from (1.3) based on the stability estimate. We also establish the large-time behavior of measure-valued solutions showing the velocity alignment.

3.1. Existence of measure-valued solutions

In this part, we discuss the global existence and uniqueness of measure-valued solutions to the equation (3.1). For this, we first define a notion of weak solutions in the definition below.

Definition 3.1**.**

For a given $T>0$, we call $f_{t}\in\mathcal{C}([0,T);\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$ a measure-valued solution of the equation (3.1) on the time interval $[0,T)$, subject to the initial datum $g_{s}\in\mathcal{C}([-\tau_{0},0];\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$, if for all compactly supported test functions $\xi\in\mathcal{C}_{c}^{\infty}(\mathbb{R}^{d}\times\mathbb{R}^{d}\times[0,T))$,

[TABLE]

where $\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d})$ denotes the set of probability measures on the phase space $\mathbb{R}^{d}\times\mathbb{R}^{d}$ with bounded first-order moment and we adopt the notation $f_{t-\tau_{0}}:\equiv g_{t-\tau_{0}}$ for $t\in[0,\tau_{0}]$.

We next introduce the $1$-Wasserstein distance.

Definition 3.2**.**

Let $\rho^{1},\rho^{2}\in\mathcal{P}_{1}(\mathbb{R}^{d})$ be two probability measures on $\mathbb{R}^{d}$. Then $1$-Waserstein distance between $\rho^{1}$ and $\rho^{2}$ is defined as

[TABLE]

where $\Pi(\rho^{1},\rho^{2})$ represents the collection of all probability measures on $\mathbb{R}^{d}\times\mathbb{R}^{d}$ with marginals $\rho^{1}$ and $\rho^{2}$ on the first and second factors, respectively.

Theorem 3.1**.**

Let the initial datum $g_{t}\in\mathcal{C}([-\tau_{0},0];\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$ and assume that there exists a constant $R>0$ such that

[TABLE]

where $B^{2d}(0,R)$ denotes the ball of radius $R$ in $\mathbb{R}^{d}\times\mathbb{R}^{d}$, centered at the origin.

*Then for any $T>0$, the delayed Vlasov alignment equation (3.1) admits a unique measure-valued solution $f_{t}\in\mathcal{C}([0,T);\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$ in the sense of Definition 3.1. *

Proof.

The proof can be done by using a similar argument as in [13, Theorem 3.1], thus we shall give it rather concisely. Let $f_{t}\in\mathcal{C}([0,T];\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$ be such that

[TABLE]

for some positive constant $R>0$. Then we use the similar estimate as in [13, Lemma 3.1] to have that there exists a constant $C>0$ such that

[TABLE]

for $(x,v)\in B^{2d}(0,R),\;t\in[0,T]$. Subsequently, this gives

[TABLE]

and

[TABLE]

for $(x,v)\in B^{2d}(0,R),\;t\in[0,T]$. This, together with an application of [4, Theorem 3.10], provides the local-in-time existence and uniqueness of measure-valued solutions to the equation (3.1) in the sense of Definition 3.1. On the other hand, once we obtain the growth estimates of support of $f_{t}$ in both position and velocity, these local-in-time solutions can be extended to the global-in-time ones. Thus in the rest of this proof, we provide the support estimate of $f_{t}$ in both position and velocity. We set

[TABLE]

for $t\in[0,T]$, where supp${}_{x}f_{t}$ and supp${}_{v}f_{t}$ represent $x$- and $v$-projections of supp$f_{t}$, respectively. We also set

[TABLE]

We first construct the system of characteristics $Z(t;x,v):=(X(t;x,v),V(t;x,v)):[0,\tau]\times\mathbb{R}^{d}\times\mathbb{R}^{d}\to\mathbb{R}^{d}\times\mathbb{R}^{d}$ associated with (3.1),

[TABLE]

where we again adopt the notation $f_{t-\tau_{0}}:\equiv g_{t-\tau_{0}}$ for $t\in[0,\tau_{0}]$. The system (3.3) is considered subject to the initial conditions

[TABLE]

for all $(x,v)\in\mathbb{R}^{d}\times\mathbb{R}^{d}$. Note that the characteristic system (3.3)–(3.4) is well-defined since the $F[f_{s}]$ is locally both bounded and Lipschitz. We can rewrite the equation for $V$ as

[TABLE]

Then, arguing as in the proof of Lemma 2.1, we get

[TABLE]

due to (3.2). Using again a similar argument as in the proof of Lemma 2.1 and the comparison lemma, we obtain

[TABLE]

which further implies $R^{t}_{X}\leq R^{0}_{X}+tR^{0}_{V}$ for $t\geq 0$. This completes the proof. ∎

3.2. Stability & mean-field limit

In this subsection, we discuss the rigorous derivation of the delayed Vlasov alignment equation (3.1) from the particle system (1.3) as $N\to\infty$. For this, we first provide the stability of measure-valued solutions of (3.1).

Theorem 3.2**.**

Let $f^{i}_{t}\in\mathcal{C}([0,T);\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$, $i=1,2$, be two measure-valued solutions of (3.1) on the time interval $[0,T]$, subject to the compactly supported initial data $g^{i}_{s}\in\mathcal{C}([-\tau_{0},0];\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$. Then there exists a constant $C=C(T)$ such that

[TABLE]

Proof.

Again, the proof is very similar to [13, Theorem 3.2], see also [7, 9]. Indeed, we can obtain

[TABLE]

Then we have

[TABLE]

for $t\in[0,T)$. ∎

Remark 3.1**.**

Since the empirical measure $f^{N}_{t}(x,v):=\frac{1}{N}\sum_{i=1}^{N}\delta_{(x^{N}_{i}(t),v^{N}_{i}(t))}(x,v)$ associated to the $N$-particle system (1.3) with the initial data $g^{N}_{s}(x,v):=\frac{1}{N}\sum_{i=1}^{N}\delta_{(x_{i}^{0}(s),v_{i}^{0}(s))}(x,v)$ for $s\in[-\tau_{0},0]$, is a solution to the equation (3.1), we can use Theorem 3.2 to have the following mean-field limit estimate:

[TABLE]

where $C$ is a positive constant independent of $N$.

3.3. Asymptotic behavior of the delayed Vlasov alignment equation (3.1)

In this part, we provide the asymptotic behavior of solutions to the equation (3.1) showing the velocity alignment under suitable assumptions on the initial data. For this, we first define the position- and velocity-diameters for a compactly supported measure $g\in\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d})$,

[TABLE]

where supp${}_{x}f$ denotes the $x$-projection of supp$f$ and similarly for supp${}_{v}f$.

Theorem 3.3**.**

Let $T>0$ and $f_{t}\in\mathcal{C}([0,T);\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$ be a weak solution of (3.1) on the time interval $[0,T)$ with compactly supported initial data $g_{s}\in\mathcal{C}([-\tau_{0},0];\mathcal{P}_{1}(\mathbb{R}^{d}\times\mathbb{R}^{d}))$. Furthermore, we assume that

[TABLE]

Then the weak solution $f_{t}$ satisfies

[TABLE]

where $C$ is a positive constant independent of $t$.

Let us point out that the flocking estimate at the particle level, see Section 2 and Remark 2.3, is independent of the number of particles, thus we can directly use the same strategy as in [11, 13, 26]. However, we provide the details of the proof for the completeness.

Proof of Theorem 3.3.

We consider an empirical measure $\{g^{N}_{s}\}_{N\in\mathbb{N}}$, which is a family of $N$-particle approximations of $g_{s}$, i.e.,

[TABLE]

where the $x_{i}^{0},v_{i}^{0}\in\mathcal{C}([-\tau_{0},0];\mathbb{R}^{d})$ are chosen such that

[TABLE]

Note that we can choose $(x_{i}^{0},v_{i}^{0})_{i=1,\dots,N}$ such that the condition (2.3) is satisfied uniformly in $N\in\mathbb{N}$ due to the assumption (3.5). Let us denote by $(x^{N}_{i},v^{N}_{i})_{i=1,\dots,N}$ the solution of the $N$-particle system (1.3)–(1.5) subject to the initial data $(x_{i}^{0},v_{i}^{0})_{i=1,\dots,N}$ constructed above. Then it follows from Theorem 2.1 that there exists a positive constant $C_{1}>0$ such that

[TABLE]

with the diameters $d_{V}$, $d_{X}$ defined in (2.1), where $C_{1}>0$ is independent of $t$ and $N$. Note that the empirical measure

[TABLE]

is a measure-valued solution of the delayed Vlasov alignment equation (3.1) in the sense of Definition 3.1. On the other hand, by Theorem 3.2, for any fixed $T>0$, we have the following stability estimate

[TABLE]

where the constant $C_{2}>0$ is independent of $N$. This yields that sending $N\to\infty$ implies $d_{V}[f_{t}]=d_{V}(t)$ on $[0,T)$ for any fixed $T>0$. Thus we have

[TABLE]

Since the uniform-in-$t$ boundedness of $d_{X}[f_{t}]$ just follows from the above exponential decay estimate, we conclude the desired result. ∎

Acknowledgments

The first author was supported by National Research Foundation of Korea(NRF) grant funded by the Korea government(MSIP) (No. 2017R1C1B2012918 and 2017R1A4A1014735) and POSCO Science Fellowship of POSCO TJ Park Foundation. The research of the second author was partially supported by the GNAMPA 2018 project Analisi e controllo di modelli differenziali non lineari (INdAM).