IAS Workshop on Symplectic Dynamics 2: Wednesday

James Colliander

2012-03-14 12 min read lessons

IAS School of Mathematics

$Simonyi Hall$

Workshop web page

Happy Einstein Birthday!

Wednesday: 2012-03-13

9:00 - 10:00 Wilfrid Gangbo, Georgia Institute of Technology, “Lifting absolutely continuous curves from P(Td) to P2(Rd)” abstract
10:15 - 11:15 Jonatan Lenells, Baylor University, “Geometry of diffeomorphism groupos, complete integrability and optimal transport” abstract
11:30 - 12:30 David Ebin, SUNY, “Groups of diffeomorphisms and geodesics on them” abstract
2:30 - 3:30 Susan Friedlander, University of Southern California, “Well / Ill-posedness results for the magneto-geostrophic equations: the importance of being even”. abstract

Wilfrid Gangbo: Lifting absolutely continuous curves from $P (T d)$ to $P 2 (R d)$

$Wilfrid Gangbo$

(chalk talk; joint work with A. Tudorascu)

This work extends earlier work on the space of probability measures on the torus $P (T)$ to analogous results on $P (T d)$ . The earlier work used the embedding $P (T) \to L 2 (0, 1)$ but we don’t have this embedding in the higher dimensional case.

Let $P 2 (R d)$ be the set of Borel measures on $R d$ with finite second moment $\int | x | 2 μ (d x) < \infty$ . We say $μ 0 \sim μ 1$ if and only if $\int F d μ 0 = \int F d μ 1 \forall F \in C (T d), \forall F \in C (R d), \forall F (x + z) = F (x), z \in Z d$ .

I define $P (T) = P (T d) / \sim$ .

Let $γ$ be a measure on $R d \times R d$ which satisfy $π 1$ # $γ = μ 0$ and $π 2$ # $γ = μ 1$ . More generally, we write $W 22 (μ 0, μ 1) = inf γ \int R d \times R d | x - y | 2 γ (d x, d y) .$

Problem: Data: $v : (0, T) \times T d \to R d$ and $t \to σ t \in P (T d)$ . Assume that $\partial t σ t + \nabla \cdot (σ v) = 0$ (in the sense of distributions). Can we find $t \to^σ t \in P 2 (R d)$ and $ˆ v : (0, T) \times T d \to R d$ such that $\partial t ˆ σ + \nabla \cdot (ˆ σ ˆ v) = 0$ in the sense of distributions. Here $^σ t \sim σ t .$

Such a lift becomes important if I want to associate the rotation number. I want to write $d d t \int T d x d σ t = \int T d v t d σ t $ .$

Weak KAM: (A small fraction of what is known) $M = T d$ . Let $h : T * M \to R$ . Let $w (z 0, z 1) = z 0 (J z 1)$ where $J$ is the usual matrix satisfying $J 2 = - I d$ . Let $X h$ denote the associated Hamiltonian vector field: $˙ ϕ = X h (ϕ), ϕ 0 = (x 0, p 0) .$

The associated flow is denoted

ϕ t = (x t, p t) .

Existence of weak Lagrangian Tori: $― h : R d \to R$ is the effective Hamiltonian, $c \in R d$ .

$\exists u \in C (t d)$ with $h (x, c + \nabla u) = ― h (c)$ (viscosity)
$u c (x) = c \cdot x + u (t) .$ ( $\partial u c$ is invariant under $ϕ$ .)
$\forall x 0 \in T d \exists v 0$ such that if $(x t, p t) = ϕ t$ then $\forall T$ $u (x T) - u (x 0) = \int T 0 [l (x, ˙ x) + c \cdot ˙ x + ― h (c)] d t$ $lim t \to \infty^x t t = - \nabla ― h (c) .$

General Fact:

M \to

compact.

Specific to finite $d$ :

Given $x \in W 1, 2 (0, T, T d)$ and take two lifts $ˆ x, ˆ y \in W 1, 2, R d$ . We then find that $^x t -^y t = n \in Z d$

because we have

lim t \to \infty^x t t = lim t \to \infty^y t t

Obstacle in infinite dimensions:

Let $M 0 = P (T d)$ and in the sense of distributions we have $\partial t σ + \nabla \cdot (σ v) = 0.$

If $\nabla \cdot (σ w) = 0$ then $v + w$ is another velocity.

(wash board…)

Let $μ \in P 2 (R d)$ and define the tangent space $T μ P 2 (R d)$ and also the space $T μ P (T d)$ . These are defined with $L 2$ closures.

Pseudo symplectic form: ….going faster and I’m not keeping up with the typing.

Theorem (Gangbo-Kun-Pacuni 2011):

$Ω$ is a closed skew symmetric nondegenerate 2-form.
$\exists X H$ such that $- d H = Ω (X H, \cdot)$ .
$˙ f = X H (f) ⟺ \partial t f + \nabla x (v f) = \nabla v (f (\nabla V + \nabla W * ρ)) .$

This is a nonlinear Vlasov equation.

I want to state the analog of the weak KAM theorem in our context.

Theorem Let $― H$ be the effective Hamiltonian of $H$ restricted to $R d$ and let $c \in R d$ .

$\exists U : P (T d) \to R$ such that (in viscosity sense) $H (μ, c + \nabla w H) = ― H (c) .$
Given $σ 0 \in P (T d) \exists v 0 : T d \to R d$ such that if $f 0 = σ 0 δ v 0$ then $f t = σ t δ v t,$ $U (σ t) - U (σ 0) = \int T 0 [L (σ t, v t) + \int R d v t \cdot c d σ t + ― H (c)] d t .$
We also have $| 1 T \int T 0 d t \int R d v t d σ t + \nabla ― H (c) | \leq c o n s t \sqrt T .$

Corollary: If

(ˆ σ,^v t)

is an appropriate lift then

lim T \to \infty 1 T \int T 0 (\int R d x d σ t) d t = - \nabla ― H (c) .

Questions:

Mather: This has connections with fluids?

Answer: Kinetic theory. Consider the system $\partial 2 t x = 1 N \sum N j = 1 W (x i - x j) - \nabla V (x i)$ . When we consider the $N \to \infty$ limit, we can move the weak KAM theory from this $N$ particle system to the infinite particle case by moving to the setting of measures. This framework lets us prove convergence of discrete models to the PDE case.

Jonatan Lenells: Geometry of diffeomorphism groups, complete integrability and optimal transport

(pdf slides; Happy $Π$ day!; Einstein’s birthday)

(joint work with B. Khesin, G. Misiolek, S. Preston)

Outline

A new equation
Geometry of $D i f f (M)$
A sphere
Optimal Transport
Geometric Statistics

A new equation

ρ t + u \cdot \nabla ρ + 12 ρ 2 = - \int M ρ 2 d μ 2 μ (M) .

$M$ is a compact Riemannian manifold.
$μ (M)$ is the volume of $M$ .
This is an exciting equation because it is completely integrable for any $M$ .
This is a geodesic equation on $D i f f (M) / D i f f μ (M)$ .
Describes $˙ H 1$ -optimal transport
reduces to the Hunter-Saxton equation for $M = S 1$ . (derived in the context of liquid crystals in the early 90s.)

Euler-Arnold Equations. Summary of those ideas.

abc-metric. You can add some other terms to the original $L 2$ inner product involving $L 2$ inner products involvling codifferentials and the musical isomorphism. A lot of different equations arise as you take different values of the parameters. Writing down the associated abc Euler-Arnold equation generates a big equation which can be specialized into various equations. To obtain the $˙ H 1$ metric, we cancel away the terms associated with factors a and c. We simplify by setting $a = 0, b = 14, c = 0$ .

The equation induced by these choices has some degeneracy issues. These can be resolved by quotienting out part of the phase space. The function $u$ is not uniquely determined but its coset is uniquely determined. (Similar issues arise in Hunter-Saxton.) The equation we are considering here is a geodesic equation on a (quotiented) Diffeomorphism group.

Jacobian determinant.

A sphere

Theorem (Khesin-Misiolek-Lennels-Preston): The map which takes the coset

[η]

to its associated Jacoobian

\sqrt J a c μ η

is an isometry onto a subset of the sphere.

This isometry lets them transport all the questions about the geodesic equation on this complicated Diff phase space into corresponding questions about geodesics on the sphere. Since we understand the sphere well, we can conjugate results there using the mapping to obtain explicit solution formulae for the geodesic equation. Magical integrability!

Preceding works.

Khesin-Misiolek 2003: Showed Hunter-Saxton may be viewed within the Euler-Arnold framework.

Lennels 2006: Recognized the image of the map as a portion of the sphere.

Optimal Transport

Optimal Transport. Wasserstein distance between two probability measures.

Moser 1965, Ebin-Marsden 1970, Otto 2001, also Benamou-Brenier.

The $˙ H 1$ optimal distance induces what they call the spherical Hellinger distance since it resembles the Hellinger distance used in probability theory

Geometric Statistics

Statistical model.

Fisher-Rao information metric.

Theorem (KMLP): The $˙ H 1$ metric coincides with the Fisher-Rao metric when restricted to any k-dimensional submanifold of the (quotiented) $D i f f$ .

This is another reason why we think this metric is important. It arises from many different points of view.

Summary

We found a new integrable PDE.
The PDE is a geodesic equation on a quotieted diff with $˙ H 1$ metric.
A sphere
One can understand what is going on using the optimal transportation point of view using the $˙ H 1$ metric.
This metric coincides with a basic metric arising in geometric statistics.

Open problems

Global weak solutions of the PDE? All solutions break in finite time because you hi the boundary of the diffeomorphism coset. However, there is no problem when you view the dynamics on the sphere. The motion along the great circle may be continued. This type of development has taken place already in the context of the Hunter-Saxton equation. The process appears to be more complicated in this more general context since the Jacobian can vanish.
Transfer results from geometric statistics into this diffeomorphism quotient. Then reinterpret these objects in the setting of PDE. Amari-Nagaoka 2000 alpha-connections, dual connections.
Develop an optimal transport theory based on the $˙ H 1$ theory.
Find a Lax pair.
Find a bi-Hamiltonian structure.
Analyze the associated two-component equation (c.f. Lennels-Zhao 2011). There is a 2-component Hunter-Saxton so that object suggests we might find a corresponding generalization. This has been observed by LZ.

David Ebin: Groups of diffeomorphisms and geodesics on them

$David Ebin$

(joint work with Stephen Preston)

Maps from a manifold to itself. Discussion of various topologies of such maps.

Volume preserving maps. Diffeomorphisms (Volumorphisms)
Even dimensional manifolds with a symplectic form. We can consider the maps which preserve the symplectic form. (Symplectomorphisms)
For odd dimensional manifolds, we can consider maps which preserve the contact form. (Contactomorphisms)

In all these cases, we can discuss the geodesics…..ack low battery.

Boothby-Wang fibration.

Susan Friedlander: Well / Ill-posedness results for the magneto-geostrophic equations: the importance of being even

(joint work with Vlad Vicol, Walter Rusin, Francisco Gancedo, Weiran Sun)

Homage to Oscar Wilde…

Amain theme is that there is a difference in behavior of solutions in Active Scalar Equations when the associated Fourier multiplier is even versus odd.

Active Scalar Equations; Incompressible Fluids

θ t + u \cdot \nabla θ = 0.

\nabla \cdot u = 0.

$u = O [θ], P D O$

$R d$ or $T d$ . Even or odd Fourier multiplier symbol. The results I’ll describe are not influenced by the presence of a physical boundary. The emphasis will be on examining the influence of the operator O on the properties of the PDE.

Consider $u j = \partial i T i j θ, \nabla \cdot u =$ . Here $T i j$ is a $d \times d$ Calderon-Zygmund operators.

ODD Symbol: Locally well-posed in Sobolev spaces. Commutator in energy estimates. **Chae et. al, Friedlander-Vicol.
EVEN Symbol: Lipschitz ill-posed in Sobolev spaces. Friedlander-Vicol; Nonuniqueness for $L \infty$ -weak solutions. Techniques from convext integration. Shuydkoy.

Recent reviews of results for certain active scalar equations. “Regularity and blowup for active scalars” Kiselev 2010.

SQG Equation: $u = R ⊥ θ$ , symbol $i (k 2, - k 1) | k |$ .

Constantin-Majda-Tabak 1994, Resnick 1995 (Chicago thesis; unpublished), Wu, Cordoba, Chae, Iyer, Ju, Fefferman

The SQG equation had been known in the geophysics community before its introduction to the mathematical community by Constantin et.al.

Local existence for smooth initial data BUT global existence of smooth solutions is OPEN (just as it is open for 3D Euler). Cordoba and Fefferman have ruled out the existence of certain solution scenarios.

“Modified” SQG Equation (Okhitani)

Insert a power of

(- Δ) 1 / 2 = Λ

in the map

θ \to u

so that

u = \nabla ⊥ Λ β - 2 θ

where $1 < \beta \leq 2. Chae, Constantin, Cordoba, Ganceda, Wu 2011. Local existence of smooth solutions in

H s

, global existence of weak solutions.

Note: result holds more generally when the symbol is ODD and order $\leq 1$ .

IPM equation: singular integral operator with EVEN symbol

Darcy’s law.

$u = R ⊥ R 1 θ .$

Cordoba-Gancedo-Orive 2007

Regular initial data, local existence, weak solutions, SQG and IPM present different behaviours. Global existence of smooth solutions is OPEN.

Even symbol: $(k 1 k 2 | k | 2, - k 21 | k | 2) .$

There is very different behaviors among these equations for rough data. “Patch type initial data”

SIPM equations (even, unbounded)

u = R ⊥ R 1 Λ β θ

Friedlander-Gancedo-Sun-Vicol (2012)

Locally Lipschitz ill-posed in $H s, s > 2$ . Proved for $0 < β \leq 2$ in $T d \times [0, \infty]$
Locally well-posed for some “patch-type” weak solutions. Proved for $0 < β < 1$ in $R 2$ .

(Discussion: the notions of “wellposedness” changes between the previous two bullet points.)

Symbol: $k 1 k ⊥ | k | β - 2$

Magnetogeostrophic (MG) equations

Friedlander-Vicol 2011

Long symbol $M$ , even, unbounded, 3D. $u = M θ$ . Here $M$ is a vector operator that defines a 3-vector $u$ .

Cauchy problem is ill-posed in Hadamard sense in Sobolev spaces. There is no Lipschitz solution map.

Ill-posedness: singular, even, symbol

Active scalar equation. Special direction with index

d

, often associated with gravity. A list of many conditions on the

d

th component

S d

of the Fourier multiplier operator….allowing them to build eigenfunctions to show Lipschitz failure.

Definition: Locally Lipschitz $(X, Y)$ well-posed.

$$| \theta_1 (\cdot, t) - \theta_1 (\cdot, t)|X \leq K | \theta1 (\cdot, 0) -\theta_2 (\cdot, 0) |_Y. $$

The spaces $X, Y$ are often chosen to be $(H r, H S)$ .

Theorem: Under the many assumptions on $S d$ , the active scalar equation is Lipschitz $(H r, H s)$ -illposed for any $r > R, s \geq r + 1$ .

Linear problem

Linearize around

θ 0 = sin m x d

. Write out a Fourier series. Crank out a recurrence relation.

Continued fractions, characteristic equation. These ideas where used by Michalkin and Sinai to show unstable eigenvalues for the shear flow for Navier-Stokes equations.

Ill-posedness of the nonlinear problem

Follows a proof by contradiction.

Effects of dissipation: MG

Dissipation:

ν (- D e l t a) 1 / 2

Using De Giorgi techniques, Caffarelli-Vasseur proved critical SQG. What can we say about the MG equation?

Case $1 / 2 M γ < 1$ : LWP in $H s$ , for $s > 52 + (1 - 2 γ)$ . Well-prepared initial data.
Case $0 < γ < 1 / 2$ : Diffusion is too weak to overcome the continued fraction construction.
Case $γ = 1 / 2$ : Unique global solution when the initial data and source are small in a suitable sense, then there exists a unique golbal solution. However, if the data are large in this respect then we can run the ill-posedness construction. This reveals a very precise dichotomy.

notes conference IAS symplectic