swarming models with repulsive-attractive effectsjcarrill/ravello/lecturejac-ravello1.pdf ·...
TRANSCRIPT
![Page 1: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/1.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Swarming Models withRepulsive-Attractive Effects
J. A. Carrillo
Imperial College London
Lecture 1, Ravello 2013
![Page 2: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/2.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Outline
1 MotivationsCollective Behavior ModelsVariationsFixed Speed models1st order Models
2 Outline of the course
3 Transversal Tool: Wasserstein DistancesDefinitionProperties
![Page 3: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/3.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Outline
1 MotivationsCollective Behavior ModelsVariationsFixed Speed models1st order Models
2 Outline of the course
3 Transversal Tool: Wasserstein DistancesDefinitionProperties
![Page 4: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/4.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Swarming by Nature or by design?
Fish schools and Birds flocks.
![Page 5: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/5.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Individual Based Models (Particle models)
Swarming = Aggregation of agents of similar size and body type generally moving ina coordinated way.
Highly developed social organization: insects (locusts, ants, bees ...), fishes, birds,micro-organisms (myxo-bacteria, ...) and artificial robots for unmanned vehicleoperation.
Interaction regions between individualsa
aAoki, Helmerijk et al., Barbaro, Birnir et al.
Repulsion Region: Rk.
Attraction Region: Ak.
Orientation Region: Ok.
Metric versus Topological Interaction
![Page 6: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/6.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Individual Based Models (Particle models)
Swarming = Aggregation of agents of similar size and body type generally moving ina coordinated way.
Highly developed social organization: insects (locusts, ants, bees ...), fishes, birds,micro-organisms (myxo-bacteria, ...) and artificial robots for unmanned vehicleoperation.
Interaction regions between individualsa
aAoki, Helmerijk et al., Barbaro, Birnir et al.
Repulsion Region: Rk.
Attraction Region: Ak.
Orientation Region: Ok.
Metric versus Topological Interaction
![Page 7: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/7.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Individual Based Models (Particle models)
Swarming = Aggregation of agents of similar size and body type generally moving ina coordinated way.
Highly developed social organization: insects (locusts, ants, bees ...), fishes, birds,micro-organisms (myxo-bacteria, ...) and artificial robots for unmanned vehicleoperation.
Interaction regions between individualsa
aAoki, Helmerijk et al., Barbaro, Birnir et al.
Repulsion Region: Rk.
Attraction Region: Ak.
Orientation Region: Ok.
Metric versus Topological Interaction
![Page 8: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/8.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Individual Based Models (Particle models)
Swarming = Aggregation of agents of similar size and body type generally moving ina coordinated way.
Highly developed social organization: insects (locusts, ants, bees ...), fishes, birds,micro-organisms (myxo-bacteria, ...) and artificial robots for unmanned vehicleoperation.
Interaction regions between individualsa
aAoki, Helmerijk et al., Barbaro, Birnir et al.
Repulsion Region: Rk.
Attraction Region: Ak.
Orientation Region: Ok.
Metric versus Topological Interaction
![Page 9: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/9.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
2nd Order Model: Newton’s like equations
D’Orsogna, Bertozzi et al. model (PRL 2006):dxi
dt= vi,
mdvi
dt= (α− β |vi|2)vi −
∑j 6=i
∇U(|xi − xj|).
Model assumptions:
Self-propulsion and friction termsdetermines an asymptotic speed of√α/β.
Attraction/Repulsion modeled by aneffective pairwise potential U(x).
U(r) = −CAe−r/`A + CRe−r/`R .
One can also use Bessel functions in 2Dand 3D to produce such a potential.
C = CR/CA > 1, ` = `R/`A < 1 andC`2 < 1:
![Page 10: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/10.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
2nd Order Model: Newton’s like equations
D’Orsogna, Bertozzi et al. model (PRL 2006):dxi
dt= vi,
mdvi
dt= (α− β |vi|2)vi −
∑j 6=i
∇U(|xi − xj|).
Model assumptions:
Self-propulsion and friction termsdetermines an asymptotic speed of√α/β.
Attraction/Repulsion modeled by aneffective pairwise potential U(x).
U(r) = −CAe−r/`A + CRe−r/`R .
One can also use Bessel functions in 2Dand 3D to produce such a potential.
C = CR/CA > 1, ` = `R/`A < 1 andC`2 < 1:
![Page 11: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/11.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
2nd Order Model: Newton’s like equations
D’Orsogna, Bertozzi et al. model (PRL 2006):dxi
dt= vi,
mdvi
dt= (α− β |vi|2)vi −
∑j 6=i
∇U(|xi − xj|).
Model assumptions:
Self-propulsion and friction termsdetermines an asymptotic speed of√α/β.
Attraction/Repulsion modeled by aneffective pairwise potential U(x).
U(r) = −CAe−r/`A + CRe−r/`R .
One can also use Bessel functions in 2Dand 3D to produce such a potential.
C = CR/CA > 1, ` = `R/`A < 1 andC`2 < 1:
![Page 12: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/12.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
2nd Order Model: Newton’s like equations
D’Orsogna, Bertozzi et al. model (PRL 2006):dxi
dt= vi,
mdvi
dt= (α− β |vi|2)vi −
∑j 6=i
∇U(|xi − xj|).
Model assumptions:
Self-propulsion and friction termsdetermines an asymptotic speed of√α/β.
Attraction/Repulsion modeled by aneffective pairwise potential U(x).
U(r) = −CAe−r/`A + CRe−r/`R .
One can also use Bessel functions in 2Dand 3D to produce such a potential.
C = CR/CA > 1, ` = `R/`A < 1 andC`2 < 1:
![Page 13: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/13.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Model with an asymptotic velocity
Classification of possible patterns: Morse potential. D’Orsogna, Bertozzi et al.model (PRL 2006).
![Page 14: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/14.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Model with an asymptotic speed
Typical patterns: milling, double milling or flocking:
![Page 15: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/15.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Velocity consensus model
Cucker-Smale Model (IEEE Automatic Control 2007):dxi
dt= vi,
dvi
dt=
N∑j=1
aij (vj − vi) ,
with the communication rate, γ ≥ 0:
aij = a(|xi − xj|) =1
(1 + |xi − xj|2)γ.
Asymptotic flocking: γ < 1/2; Cucker-Smale.General Proof for 0 < γ ≤ 1/2; C.-Fornasier-Rosado-Toscani.
![Page 16: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/16.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Collective Behavior Models
Velocity consensus model
Cucker-Smale Model (IEEE Automatic Control 2007):dxi
dt= vi,
dvi
dt=
N∑j=1
aij (vj − vi) ,
with the communication rate, γ ≥ 0:
aij = a(|xi − xj|) =1
(1 + |xi − xj|2)γ.
Asymptotic flocking: γ < 1/2; Cucker-Smale.General Proof for 0 < γ ≤ 1/2; C.-Fornasier-Rosado-Toscani.
![Page 17: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/17.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Variations
Outline
1 MotivationsCollective Behavior ModelsVariationsFixed Speed models1st order Models
2 Outline of the course
3 Transversal Tool: Wasserstein DistancesDefinitionProperties
![Page 18: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/18.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Variations
Leadership, Geometrical Constraints, and Cone of Influence
Cucker-Smale with local influence regions:dxi
dt= vi ,
dvi
dt=∑
j∈Σi(t)
a(|xi − xj|)(vj − vi) ,
where Σi(t) ⊂ 1, . . . ,N is the set of dependence, given by
Σi(t) :=
1 ≤ ` ≤ N :
(x` − xi) · vi
|x` − xi||vi|≥ α
.
Cone of Vision:
![Page 19: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/19.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Variations
Leadership, Geometrical Constraints, and Cone of Influence
Cucker-Smale with local influence regions:dxi
dt= vi ,
dvi
dt=∑
j∈Σi(t)
a(|xi − xj|)(vj − vi) ,
where Σi(t) ⊂ 1, . . . ,N is the set of dependence, given by
Σi(t) :=
1 ≤ ` ≤ N :
(x` − xi) · vi
|x` − xi||vi|≥ α
.
Cone of Vision:
![Page 20: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/20.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Variations
Roosting Forces
Adding a roosting area to the model:dxi
dt= vi,
dvi
dt= (α− β |vi|2)vi −
∑j 6=i
∇U(|xi − xj|)− v⊥i ∇xi
[φ(xi) · v⊥i
].
with the roosting potential φ given by φ(x) :=b4
(|x|
RRoost
)4
.
Roosting effect: milling flocks N = 400,Rroost = 20.
![Page 21: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/21.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Variations
Adding Noise
Self-Propelling/Friction/Interaction with Noise Particle Model:xi = vi,
dvi =
(α− β |vi|2)vi −∇xi
∑j 6=i
U(|xi − xj|)
dt +√
2σ dΓi(t) ,
where Γi(t) are N independent copies of standard Wiener processes with values inRd and σ > 0 is the noise strength. The Cucker–Smale Particle Model with Noise:
dxi = vidt ,
dvi =
N∑j=1
a(|xj − xi|)(vj − vi) dt +
√√√√2σN∑
j=1
a(|xj − xi|) dΓi(t) .
![Page 22: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/22.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Variations
Adding Noise
Self-Propelling/Friction/Interaction with Noise Particle Model:xi = vi,
dvi =
(α− β |vi|2)vi −∇xi
∑j 6=i
U(|xi − xj|)
dt +√
2σ dΓi(t) ,
where Γi(t) are N independent copies of standard Wiener processes with values inRd and σ > 0 is the noise strength. The Cucker–Smale Particle Model with Noise:
dxi = vidt ,
dvi =N∑
j=1
a(|xj − xi|)(vj − vi) dt +
√√√√2σN∑
j=1
a(|xj − xi|) dΓi(t) .
![Page 23: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/23.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Fixed Speed models
Outline
1 MotivationsCollective Behavior ModelsVariationsFixed Speed models1st order Models
2 Outline of the course
3 Transversal Tool: Wasserstein DistancesDefinitionProperties
![Page 24: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/24.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Fixed Speed models
Vicsek’s model
Assume N particles moving at unit speed: reorientation & diffusion:dXi
t = V it dt,
dV it =√
2 P(V it ) dBi
t − P(V it )
1N
N∑j=1
K(Xit−Xj
t)(V it − V j
t )
dt.
Here P(v) is the projection operator on the tangent space at v/|v| to the unit spherein Rd, i.e.,
P(v) = I − v⊗ v|v|2 .
Noise in the Stratatonovich sense: imposed by the rigorous construction of theBrownian motion on a manifold. Rigorous derivation: Bolley-Cañizo-C.
Main issue: phase transition? Degond-Liu-Frouvelle.
![Page 25: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/25.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Fixed Speed models
Vicsek’s model
Assume N particles moving at unit speed: reorientation & diffusion:dXi
t = V it dt,
dV it =√
2 P(V it ) dBi
t − P(V it )
1N
N∑j=1
K(Xit−Xj
t)(V it − V j
t )
dt.
Here P(v) is the projection operator on the tangent space at v/|v| to the unit spherein Rd, i.e.,
P(v) = I − v⊗ v|v|2 .
Noise in the Stratatonovich sense: imposed by the rigorous construction of theBrownian motion on a manifold. Rigorous derivation: Bolley-Cañizo-C.
Main issue: phase transition? Degond-Liu-Frouvelle.
![Page 26: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/26.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Fixed Speed models
Vicsek’s model
Assume N particles moving at unit speed: reorientation & diffusion:dXi
t = V it dt,
dV it =√
2 P(V it ) dBi
t − P(V it )
1N
N∑j=1
K(Xit−Xj
t)(V it − V j
t )
dt.
Here P(v) is the projection operator on the tangent space at v/|v| to the unit spherein Rd, i.e.,
P(v) = I − v⊗ v|v|2 .
Noise in the Stratatonovich sense: imposed by the rigorous construction of theBrownian motion on a manifold. Rigorous derivation: Bolley-Cañizo-C.
Main issue: phase transition? Degond-Liu-Frouvelle.
![Page 27: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/27.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Fixed Speed models
Vicsek’s model
Assume N particles moving at unit speed: reorientation & diffusion:dXi
t = V it dt,
dV it =√
2 P(V it ) dBi
t − P(V it )
1N
N∑j=1
K(Xit−Xj
t)(V it − V j
t )
dt.
Here P(v) is the projection operator on the tangent space at v/|v| to the unit spherein Rd, i.e.,
P(v) = I − v⊗ v|v|2 .
Noise in the Stratatonovich sense: imposed by the rigorous construction of theBrownian motion on a manifold. Rigorous derivation: Bolley-Cañizo-C.
Main issue: phase transition? Degond-Liu-Frouvelle.
![Page 28: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/28.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Fixed Speed models
Vicsek’s model
Assume N particles moving at unit speed: reorientation & diffusion:dXi
t = V it dt,
dV it =√
2 P(V it ) dBi
t − P(V it )
1N
N∑j=1
K(Xit−Xj
t)(V it − V j
t )
dt.
Here P(v) is the projection operator on the tangent space at v/|v| to the unit spherein Rd, i.e.,
P(v) = I − v⊗ v|v|2 .
Noise in the Stratatonovich sense: imposed by the rigorous construction of theBrownian motion on a manifold. Rigorous derivation: Bolley-Cañizo-C.
Main issue: phase transition? Degond-Liu-Frouvelle.
![Page 29: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/29.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
1st order Models
Outline
1 MotivationsCollective Behavior ModelsVariationsFixed Speed models1st order Models
2 Outline of the course
3 Transversal Tool: Wasserstein DistancesDefinitionProperties
![Page 30: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/30.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
1st order Models
1st Order Friction Model:
Edelshtein-Keshet, Mogilner (JMB 2000): Assume the variations of the velocity andspeed are much smaller than spatial variations, then from Newton’s equation:
md2xi
d2t+ α
dxi
dt+∑j 6=i
∇U(|xi − xj|) = 0
so finally, we obtain
dxi
dt= −
∑j 6=i
∇U(|xi − xj|) in the continuum setting V
∂ρ∂t + div (ρv) = 0v = −∇U ∗ ρ
Flock Solutions: stationary states xsi of the 1st order model are connected to
particular solutions of the Bertozzi etal 2nd order model of the form
xi(t) = xsi + tv0
with v0 fixed with |v0|2 = αβ
.
For which potentials do we evolve towards some nontrivial steady states/patterns?Is there any implication of the stability from first to 2nd order models?
![Page 31: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/31.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
1st order Models
1st Order Friction Model:
Edelshtein-Keshet, Mogilner (JMB 2000): Assume the variations of the velocity andspeed are much smaller than spatial variations, then from Newton’s equation:
md2xi
d2t+ α
dxi
dt+∑j 6=i
∇U(|xi − xj|) = 0
so finally, we obtain
dxi
dt= −
∑j 6=i
∇U(|xi − xj|) in the continuum setting V
∂ρ∂t + div (ρv) = 0v = −∇U ∗ ρ
Flock Solutions: stationary states xsi of the 1st order model are connected to
particular solutions of the Bertozzi etal 2nd order model of the form
xi(t) = xsi + tv0
with v0 fixed with |v0|2 = αβ
.
For which potentials do we evolve towards some nontrivial steady states/patterns?Is there any implication of the stability from first to 2nd order models?
![Page 32: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/32.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
1st order Models
1st Order Friction Model:
Edelshtein-Keshet, Mogilner (JMB 2000): Assume the variations of the velocity andspeed are much smaller than spatial variations, then from Newton’s equation:
md2xi
d2t+ α
dxi
dt+∑j 6=i
∇U(|xi − xj|) = 0
so finally, we obtain
dxi
dt= −
∑j 6=i
∇U(|xi − xj|) in the continuum setting V
∂ρ∂t + div (ρv) = 0v = −∇U ∗ ρ
Flock Solutions: stationary states xsi of the 1st order model are connected to
particular solutions of the Bertozzi etal 2nd order model of the form
xi(t) = xsi + tv0
with v0 fixed with |v0|2 = αβ
.
For which potentials do we evolve towards some nontrivial steady states/patterns?Is there any implication of the stability from first to 2nd order models?
![Page 33: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/33.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
1st order Models
1st Order Friction Model:
Edelshtein-Keshet, Mogilner (JMB 2000): Assume the variations of the velocity andspeed are much smaller than spatial variations, then from Newton’s equation:
md2xi
d2t+ α
dxi
dt+∑j 6=i
∇U(|xi − xj|) = 0
so finally, we obtain
dxi
dt= −
∑j 6=i
∇U(|xi − xj|) in the continuum setting V
∂ρ∂t + div (ρv) = 0v = −∇U ∗ ρ
Flock Solutions: stationary states xsi of the 1st order model are connected to
particular solutions of the Bertozzi etal 2nd order model of the form
xi(t) = xsi + tv0
with v0 fixed with |v0|2 = αβ
.
For which potentials do we evolve towards some nontrivial steady states/patterns?Is there any implication of the stability from first to 2nd order models?
![Page 34: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/34.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
1st order Models
1st Order Friction Model:
Edelshtein-Keshet, Mogilner (JMB 2000): Assume the variations of the velocity andspeed are much smaller than spatial variations, then from Newton’s equation:
md2xi
d2t+ α
dxi
dt+∑j 6=i
∇U(|xi − xj|) = 0
so finally, we obtain
dxi
dt= −
∑j 6=i
∇U(|xi − xj|) in the continuum setting V
∂ρ∂t + div (ρv) = 0v = −∇U ∗ ρ
Flock Solutions: stationary states xsi of the 1st order model are connected to
particular solutions of the Bertozzi etal 2nd order model of the form
xi(t) = xsi + tv0
with v0 fixed with |v0|2 = αβ
.
For which potentials do we evolve towards some nontrivial steady states/patterns?Is there any implication of the stability from first to 2nd order models?
![Page 35: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/35.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Mathematical Questions:What are the continuum models associated to these systems as the number ofindividuals gets larger and larger? Mean-field limits.
What is the good analytical framework to deal with the possible concentrationof mass in finite/infinite time in space or in velocity?
What is the good analytical framework to deal with particles and continuumsolutions at the same time?
How to deal with the stability of patterns, which perturbations?
![Page 36: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/36.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Mathematical Questions:What are the continuum models associated to these systems as the number ofindividuals gets larger and larger? Mean-field limits.
What is the good analytical framework to deal with the possible concentrationof mass in finite/infinite time in space or in velocity?
What is the good analytical framework to deal with particles and continuumsolutions at the same time?
How to deal with the stability of patterns, which perturbations?
![Page 37: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/37.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Mathematical Questions:What are the continuum models associated to these systems as the number ofindividuals gets larger and larger? Mean-field limits.
What is the good analytical framework to deal with the possible concentrationof mass in finite/infinite time in space or in velocity?
What is the good analytical framework to deal with particles and continuumsolutions at the same time?
How to deal with the stability of patterns, which perturbations?
![Page 38: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/38.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Mathematical Questions:What are the continuum models associated to these systems as the number ofindividuals gets larger and larger? Mean-field limits.
What is the good analytical framework to deal with the possible concentrationof mass in finite/infinite time in space or in velocity?
What is the good analytical framework to deal with particles and continuumsolutions at the same time?
How to deal with the stability of patterns, which perturbations?
![Page 39: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/39.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Schedule:
Lecture 2: Second order Models - Kinetic Equations for Swarming: measuresolutions - mean field limit with/without noise.
Lectures 3-4: First order Models - Aggregation Equations: derivation andmean-field limit, stability/instability of steady states for repulsive/attractivepotentials. Qualitative properties of Steady States.
Lecture 5-6: Second order Models - Kinetic Equations for Swarming: Flocksolutions: Stability. Mill Solutions: Instability. Hydrodynamic models: millsand double mills. Asymptotic Speed Models as friction limits.
![Page 40: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/40.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Schedule:
Lecture 2: Second order Models - Kinetic Equations for Swarming: measuresolutions - mean field limit with/without noise.
Lectures 3-4: First order Models - Aggregation Equations: derivation andmean-field limit, stability/instability of steady states for repulsive/attractivepotentials. Qualitative properties of Steady States.
Lecture 5-6: Second order Models - Kinetic Equations for Swarming: Flocksolutions: Stability. Mill Solutions: Instability. Hydrodynamic models: millsand double mills. Asymptotic Speed Models as friction limits.
![Page 41: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/41.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Schedule:
Lecture 2: Second order Models - Kinetic Equations for Swarming: measuresolutions - mean field limit with/without noise.
Lectures 3-4: First order Models - Aggregation Equations: derivation andmean-field limit, stability/instability of steady states for repulsive/attractivepotentials. Qualitative properties of Steady States.
Lecture 5-6: Second order Models - Kinetic Equations for Swarming: Flocksolutions: Stability. Mill Solutions: Instability. Hydrodynamic models: millsand double mills. Asymptotic Speed Models as friction limits.
![Page 42: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/42.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Outline
1 MotivationsCollective Behavior ModelsVariationsFixed Speed models1st order Models
2 Outline of the course
3 Transversal Tool: Wasserstein DistancesDefinitionProperties
![Page 43: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/43.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Definition of the distance1
Transporting measures:
Given T : Rd −→ Rd mesurable, we say that ν = T#µ, if ν[K] := µ[T−1(K)] for allmesurable sets K ⊂ Rd, equivalently∫
Rdϕ dν =
∫Rd
(ϕ T) dµ
for all ϕ ∈ Co(Rd).
Random variables:
Say that X is a random variable with law given by µ, is to sayX : (Ω,A,P) −→ (Rd,Bd) is a mesurable map such that X#P = µ, i.e.,∫
Rdϕ(x) dµ =
∫Ω
(ϕ X) dP = E [ϕ(X)] .
1C. Villani, AMS Graduate Texts (2003).
![Page 44: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/44.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Definition of the distance1
Transporting measures:
Given T : Rd −→ Rd mesurable, we say that ν = T#µ, if ν[K] := µ[T−1(K)] for allmesurable sets K ⊂ Rd, equivalently∫
Rdϕ dν =
∫Rd
(ϕ T) dµ
for all ϕ ∈ Co(Rd).
Random variables:
Say that X is a random variable with law given by µ, is to sayX : (Ω,A,P) −→ (Rd,Bd) is a mesurable map such that X#P = µ, i.e.,∫
Rdϕ(x) dµ =
∫Ω
(ϕ X) dP = E [ϕ(X)] .
1C. Villani, AMS Graduate Texts (2003).
![Page 45: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/45.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Definition of the distance1
Transporting measures:
Given T : Rd −→ Rd mesurable, we say that ν = T#µ, if ν[K] := µ[T−1(K)] for allmesurable sets K ⊂ Rd, equivalently∫
Rdϕ dν =
∫Rd
(ϕ T) dµ
for all ϕ ∈ Co(Rd).
Random variables:
Say that X is a random variable with law given by µ, is to sayX : (Ω,A,P) −→ (Rd,Bd) is a mesurable map such that X#P = µ, i.e.,∫
Rdϕ(x) dµ =
∫Ω
(ϕ X) dP = E [ϕ(X)] .
1C. Villani, AMS Graduate Texts (2003).
![Page 46: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/46.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Two piles of sand!
Energy needed to transport m kg of sand from x = a to x = b:
energy = m |a− b|2
W22 (ρ1, ρ2) = Among all possible ways to transport the mass from ρ1 to ρ2, find the
one that minimizes the total energy
W22 (ρ1, ρ2) =
∫Rd|x− T(x)|2 dρ1(x)
![Page 47: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/47.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Two piles of sand!
Energy needed to transport m kg of sand from x = a to x = b:
energy = m |a− b|2
W22 (ρ1, ρ2) = Among all possible ways to transport the mass from ρ1 to ρ2, find the
one that minimizes the total energy
W22 (ρ1, ρ2) =
∫Rd|x− T(x)|2 dρ1(x)
![Page 48: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/48.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Definition of the distanceKantorovich-Rubinstein-Wasserstein Distance p = 1, 2,∞:
Wpp (µ, ν)= infπ
∫∫Rd×Rd
|x− y|p dπ(x, y)
= inf(X,Y) E [|X − Y|p]
where the transference plan π runs over the set of joint probability measures onRd × Rd with marginals µ and ν ∈ Pp(Rd) and (X, Y) are all possible couples ofrandom variables with µ and ν as respective laws.
W∞(µ, ν) = infπ∈Π(µ,ν)
sup(x,y)∈supp(π)
|x− y| ,
Monge’s optimal mass transport problem:
Find
I := infT
∫Rd|x− T(x)|p dµ(x); ν = T#µ
1/p
.
Take γT = (1Rd × T)#µ as transference plan π.
![Page 49: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/49.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Definition of the distanceKantorovich-Rubinstein-Wasserstein Distance p = 1, 2,∞:
Wpp (µ, ν)= infπ
∫∫Rd×Rd
|x− y|p dπ(x, y)
= inf(X,Y) E [|X − Y|p]
where the transference plan π runs over the set of joint probability measures onRd × Rd with marginals µ and ν ∈ Pp(Rd) and (X, Y) are all possible couples ofrandom variables with µ and ν as respective laws.
W∞(µ, ν) = infπ∈Π(µ,ν)
sup(x,y)∈supp(π)
|x− y| ,
Monge’s optimal mass transport problem:
Find
I := infT
∫Rd|x− T(x)|p dµ(x); ν = T#µ
1/p
.
Take γT = (1Rd × T)#µ as transference plan π.
![Page 50: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/50.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Definition of the distanceKantorovich-Rubinstein-Wasserstein Distance p = 1, 2,∞:
Wpp (µ, ν)= infπ
∫∫Rd×Rd
|x− y|p dπ(x, y)
= inf(X,Y) E [|X − Y|p]
where the transference plan π runs over the set of joint probability measures onRd × Rd with marginals µ and ν ∈ Pp(Rd) and (X, Y) are all possible couples ofrandom variables with µ and ν as respective laws.
W∞(µ, ν) = infπ∈Π(µ,ν)
sup(x,y)∈supp(π)
|x− y| ,
Monge’s optimal mass transport problem:
Find
I := infT
∫Rd|x− T(x)|p dµ(x); ν = T#µ
1/p
.
Take γT = (1Rd × T)#µ as transference plan π.
![Page 51: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/51.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Definition
Three examples
W22 (δa, δb) = |a− b|2
W22 (ρ, δX0 ) =
∫|X0 − y|2 dρ(y)
= Var (ρ)
![Page 52: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/52.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
Outline
1 MotivationsCollective Behavior ModelsVariationsFixed Speed models1st order Models
2 Outline of the course
3 Transversal Tool: Wasserstein DistancesDefinitionProperties
![Page 53: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/53.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
Euclidean Wasserstein Distance
Convergence Properties
1 Convergence of measures: W2(µn, µ)→ 0 is equivalent to µn µ weakly-*as measures and convergence of second moments.
2 Weak lower semicontinuity: Given µn µ and νn ν weakly-* asmeasures, then
W2(µ, ν) ≤ lim infn→∞
W2(µn, νn).
3 Completeness: The space P2(Rd) endowed with the distance W2 is a completemetric space.
![Page 54: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/54.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
Euclidean Wasserstein Distance
Convergence Properties
1 Convergence of measures: W2(µn, µ)→ 0 is equivalent to µn µ weakly-*as measures and convergence of second moments.
2 Weak lower semicontinuity: Given µn µ and νn ν weakly-* asmeasures, then
W2(µ, ν) ≤ lim infn→∞
W2(µn, νn).
3 Completeness: The space P2(Rd) endowed with the distance W2 is a completemetric space.
![Page 55: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/55.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
Euclidean Wasserstein Distance
Convergence Properties
1 Convergence of measures: W2(µn, µ)→ 0 is equivalent to µn µ weakly-*as measures and convergence of second moments.
2 Weak lower semicontinuity: Given µn µ and νn ν weakly-* asmeasures, then
W2(µ, ν) ≤ lim infn→∞
W2(µn, νn).
3 Completeness: The space P2(Rd) endowed with the distance W2 is a completemetric space.
![Page 56: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/56.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
One dimensional Case
Distribution functions:
In one dimension, denoting by F(x) the distribution function of µ,
F(x) =
∫ x
−∞dµ,
we can define its pseudo inverse:
F−1(η) = infx : F(x) > η for η ∈ (0, 1),
we have F−1 : ((0, 1),B1), dη) −→ (R,B1) is a random variable with law µ, i.e.,F−1#dη = µ ∫
Rϕ(x) dµ =
∫ 1
0ϕ(F−1(η)) dη = E [ϕ(X)] .
![Page 57: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/57.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
One dimensional Case
Distribution functions:
In one dimension, denoting by F(x) the distribution function of µ,
F(x) =
∫ x
−∞dµ,
we can define its pseudo inverse:
F−1(η) = infx : F(x) > η for η ∈ (0, 1),
we have F−1 : ((0, 1),B1), dη) −→ (R,B1) is a random variable with law µ, i.e.,F−1#dη = µ ∫
Rϕ(x) dµ =
∫ 1
0ϕ(F−1(η)) dη = E [ϕ(X)] .
![Page 58: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/58.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
One dimensional Case
Wasserstein distance:
In one dimension, it can be checkeda that given two measures µ and ν withdistribution functions F(x) and G(y) then, (F−1 × G−1)#dη has joint distributionfunction H(x, y) = min(F(x),G(y)). Therefore, in one dimension, the optimal planis given by πopt(x, y) = (F−1 × G−1)#dη, and thus
Wp(µ, ν) =
(∫ 1
0[F−1(η)− G−1(η)]p dη
)1/p
= ‖F−1 − G−1‖Lp(R)
1 ≤ p ≤ ∞.
aW. Hoeffding (1940); M. Fréchet (1951); A. Pulvirenti, G. Toscani, Annali Mat. Pura Appl. (1996).
![Page 59: Swarming Models with Repulsive-Attractive Effectsjcarrill/Ravello/LectureJAC-Ravello1.pdf · icreauab Motivations Outline of the course Transversal Tool: Wasserstein Distances Swarming](https://reader034.vdocument.in/reader034/viewer/2022042214/5eb9fa6b2015db2ef05e5c36/html5/thumbnails/59.jpg)
icreauab
Motivations Outline of the course Transversal Tool: Wasserstein Distances
Properties
One dimensional Case
Wasserstein distance:
In one dimension, it can be checkeda that given two measures µ and ν withdistribution functions F(x) and G(y) then, (F−1 × G−1)#dη has joint distributionfunction H(x, y) = min(F(x),G(y)). Therefore, in one dimension, the optimal planis given by πopt(x, y) = (F−1 × G−1)#dη, and thus
Wp(µ, ν) =
(∫ 1
0[F−1(η)− G−1(η)]p dη
)1/p
= ‖F−1 − G−1‖Lp(R)
1 ≤ p ≤ ∞.
aW. Hoeffding (1940); M. Fréchet (1951); A. Pulvirenti, G. Toscani, Annali Mat. Pura Appl. (1996).