dispersive | James Colliander

Zaher Hani Seminar: A Continuum Large Box Limit for the Cubic Nonlinear Schrödinger Equation

Thu, 21 Mar 2013 19:02:33 +0000

Analysis & Applied Math Seminar 2013-03-21

Abstract: Inspired by the general paradigm of weak turbulence theory, we consider the 2D cubic nonlinear Schrödinger equation on a box of size L with periodic boundary conditions. In an appropriate “large box regime” (L very large), we derive a continuum equation on ℝ2 that governs the dynamics of the discrete frequency modes over nonlinear time scales. This equation turns out to satisfy many surprising symmetries and conservation laws, as well as several families of explicit solutions. (This is joint work with Erwan Faou (INRIA, France) and Pierre Germain (Courant Institute, NYU)).

#Introduction#

Setup

2d cubic defocsusing or focusing NLS on a box of size $L$. Energy and mass conservation. NLS is GWP for small data in $H^s$ for $ s \geq 1$. We are not concerned with the existence issue. We are working in the setting of global-in-time solutions.

Physical and mathematical setup: weak nonlinearity.

Aim: understand out-of-equilibrium dnamics of small solutions. e.g. CKSTT
Take small data….nonlinear time scale is $\epsilon^{-2}$ where we imagine $ u \thicksim \epsilon v$ and we study $v$ with an $\epsilon^2$ coefficient and the data is of the size 1.
Fourier ansatz, transfer dynamics onto the coefficients. Comments on the $L^2$ norm dependence upon the box size parameter $L$.
Express the dynamics in terms of the $a_k (t)$.
Define the 4-frequency convolution hypersurface. He calls that $S_K$.
Interaction Representation. Conjugate by the fast linear dynamic….remove the linear dynamics. The new Fourier variable is called $\tilde{a}_k (t).$

All that has happened was a change of variables enabling us to view the dynamics on the Fourier side.

Weak turbulence “paradigm”

Aim: statistical description of out-of-equilibrium dynamics of small solutions (Zakharov 60s, Kolmogorov 50s)
RPA rand phase and amplitude.
$n(K,t) = {\mathbb{E}} |a_k (t)|^2$ is the wave spectrum or mass density.
propagation of chaos assumption. True at $t=0$, but not propagated.
Roughly, we have three main steps:
1. Statistical and time averaging.
2. large-box limit $L \rightarrow \infty$.
3. weak nonlinearity limit $\epsilon \rightarrow 0$ to arrive at a continuum equation for $n(K), ~ K \in {\mathbb{R}}^2.$

The Kolmogorov-Zakharov kinetic equation. Long convolution equation localized on the convolution hypersurvface and further localized on the resonant set.

Admits explicit stationary solutions called KZ spectra. These solutions are thought to offer some explanation to some cascade phenomena.
Non-rigorous, KZ spectra are not integrable, negliects some finite-size effects, some numerical discrepancies, the appearance of some coherent structures called “quasi-solitons” even in defocusing problems.

A New Limiting equation

Statistical averaging was causing problems. Let’s dispense with that but still take the large box and weak nonlinearity limits.

Resonant cutoff/normal forms transformation. He goes to the board and describes the separation. On the non-resonant portion, he makes a stationary phase type integration by parts, and then makes a crude estimate using the equation. This shows the non-resonant portion contributes at size $\epsilon^4 L^2$. Therefore, we concentrate our attention on the resonant terms.

(slide 12/37 and we are 20 minutes into the talk….)

He analyzes the convolution + resonance condition and identifies orthogonality properties based on the pythagorean relationship among frequencies.

Parametrization of rectangles in $\mathbb{Z}^2 / L$.

A lattice point $ J \in \mathbb{Z}^2 / L$ is called visible if $z = (p,q)/L$ with $gcd (|p|, |q|) =1. $ These points can be connected by a straight to the origin without hitting another point in the lattice.

Some new coordinates involving an $\alpha$ and $\beta$.

Co-prime equidistribution:

You can, in certain circumstances, replace sums by corresponding integrals with bounds. A classical number theory result establishes the density of visible lattice points in $\mathbb{Z}^2 / M$ is $\frac{6}{\pi^2}$. This lets you translate equidistribution into a co-prime equidistribution statement enabling us to replace sums by corresponding integrals with bounds.

Q:…interesting, I wonder to what extent similar ideas can be used on the sums we have omitted earlier in the argument. Perhaps those sums can also be represented as integrals with appropriate bounds?

Continuum limit.

Following these formal arguments leads to an integral equation resembling the KZ equation. Q: What are the differences/similarities with the KZ equation? One difference is that it preserves the Hamiltonian structure and has a positive definite Hamiltonian. He calls this equation $*$.

Symmetries lead to conserved quantities.

He writes the trilinear term in the equation as $\mathcal{T}(f,g,h)$.

Hamiltonian
Mass
Momentum
Position
Second momentum
Kinetic energy
Angular momentum

A scaling property.

Invariance under Fourier transform.

If $g$ solves $$ then $\hat{g}$ also solves $$.

Properties of the continuum equation

Boudedness properties and well-posedness. He reports on LWP and GWP properties of the equation $*$.

Gaussian family is a family of explicit stationary solutions. Gaussians are the unique maxima of the Hamiltonian functional.

Heavy tailed solutions.

Are there more?

Invariance of Harmonic oscillator eigenspaces. Hermite polynomials. The associated linear spans are invariant under the nonlinear flow $*$. The Hamiltonian of the harmonic oscillator is an integral of motion so the two flows commute and you get this easily.

Question: Is this equation $*$ completely integrable?

Rigorous Approximate Results

In analogy to the CKSTT cascade result, there is a reduction to an equation related to NLS. Can we transfer information from $*$ back to learn something about NLS?

Three difficulties:

Pass to the resonant sum.
Obtain good discrete to continuum error estimates.
Trilinear estimates on resonant sums.

….discussion of these issues…. identifies the small nonlinearity regime characterized by the condition $$ \epsilon^4 L^2 \ll \frac{\epsilon^2 \log L}{L^2}.$$

Möbius inversion formula.

Convergence Theorem: ….long statement. He gets a convergence statement on an interval that is longer than the nonlinear time scale by a factor $ M \leq \log \log L$.

Further Questions

Numerical study comparing NLS and $*$.
Other explicit solutions of $*$? Cascading solutions? Videos.
Is $*$ completely integrable?
Similar continuum limit for other equations?

GWP of Gross-Pitaevskii Equation on R4

Mon, 28 May 2012 18:55:42 +0000

Last week, I had a chance to visit Edinburgh in part to serve as the external examiner on the PhD Thesis (papers) of Tim Candy. Tim is now Dr. Timothy Candy and has an exciting research program to develop as a postdoc at Imperial.

It turned out I had lucky timing since my visit overlapped with a visit by Oana Pocovnicu. I had a chance to hear her speak about her recent work on the Gross-Pitaevskii equation. I took some notes during Oana’s talk and they appear below.

(joint work with Rowan Killip, Tadahiro Oh, and Monica Visan)

Edinburgh talk. 2012-05-21

Dynamics becomes more interesting with a nonvanishing condition at infinity.
This is the so-called energy critical case.

$$ i \partial_t u + \Delta u = (|u|^2 - 1)u, u(0) = u_0 $$

The modulus will tend to 1 as $ |x| \rightarrow 1$.

Literature

$R$
- Zhidkov 1987: introduced Zhidkov spaces.
- Gall 2004. gGWP in $X^1 (R)$
$R^2, R^3$
- Bethuel-Saut 1999 in $1+ H^1$.
- Gourbet 2007
- Gallo 2008
- Gerard 2006 in the energy space.
$R^4$
- Gerard 2006, small energy data such that $\nabla u \in L^2_t L^4_x.$

Remark: energy critical in $R^4$.

Gerard 2006 considered the energy space:

$$ E_{GP} = [ u = \alpha + v: |\alpha | =1, v \in \dot{H}^1, |v|^2 + 2 \Re (\overline{\alpha}v) \in L^2 (R^d)]. $$

Finite energy data do not have winding at spatial infinity. Therefore, to treat the finite energy case, it suffices to reduce the study to the setting where $u = 1 + v$ and $v$ satisfies…. She reduces the study to finite energy data so the set up excludes vortices right away.

Theorem (K-O-P-V): GP is GWP in the energy space $E_{GP} (R^4)$.

Two ingredients:

GWP of energy-critical defocusing NLS on $R^4$.
Perturbation theory: We will treat the equation as a perturbation off the cubic NLS.

Scaling Invariance

Dilation invariance of solutions for cubic NLS is described.
Dependence of $\dot{H}^s$ in terms of the scaling parameter $\lambda$.
critical, subcritical, supercritical.
Cubic NLS on $R^4$ is critical in $\dot{H}^1$. Quintic NLS on $R^3$ is also critical in $\dot{H}^1$.

Strichartz Estimates

Dispersive decay estimate
Strichartz Norm; supremum over the admissible pairs.
$N(I \times R^d)$ is the dual space of the Strichartz space $S(I\times R^d)$.
Homogeneous Strichartz estimate
Inhomogeneous Strichartz estimate
Admissible pairs on $R^4: (\infty, 2), (2,4), (6, \frac{12}{5})$.
By Sobolev embedding, we have some nice Strichartz containments.

Energy Critical NLS

LWP. Cazenave-Weissler 1989
GWP for small data. She then describes this by passing through Strichartz and identifies:
- If $\| \nabla e^{it \Delta } w_0 \|{L^6_t L^{12/5}_x}$ is small, we can close the argument.
- The smallness of this expression can be insured by shrinking $T$, but this depends upon the profile properties not just upon the norm of the data.
- GWP for small data follows.
Explains the blowup critereon showing that the spacetime $L^6$ norm controls the GWP+Scattering theory.

Main Results on defocusing energy-critical NLS

Bourgain 1999: GWP + Scattering, quintic NLS on $R^3$ with radial data.
- induction on energy
- localized Morawetz estimate
Grillakis 2000: global regularity for quintic NLS on $R^3$ with radial data.
CKSTT 2003: removed the radial assumption on $R^3$.
Ryckman-Visan 2007: GWP and scattering for cubic NLS on $R^4$.
Visan 2010: Simpler method for GWP+Scattering for cubic NLS on $R^4$, building on work of Dodson.
Kenig-Merle 2006: focusing energy-critical NLS on $R^3, R^4, R^4, R^5$. GWP+ Scattering for radial data with energy and kinetic energy smaller than those of the stationary solution.

Goal: prove existence of a global solution with control on the spacetime $L^6$.

Contradiction strategy.
Minimal blowup solution must exist.
Minimal blowup solutions mut be almost periodic. They are localized in physical and Fourier space.
Frequency localized Morawetz inequality. (only true for the minimal blowup solution). This is obtained by localizing in frequency the interaction Morawetz estimate.
This show that we have a smallness property on the spacetime $L^3$ norm on the high frequencies.
With some interpolation, we can then prove that the spacetime $L^6$ is bounded, contradicting the hypothesis.

Cubic NLS on $R^4$ (Visan)

(Original proof due to Ryckman-Visan but Visan recently simplified that following some ideas of Dodson.)

By contradiction and using concentration-compactness we have a minimal blowup solution.
There are only two scenarios. Rapid frequency cascade scenario; quasi-soliton scenario.

These are excluded using the long-time Strichartz estimates in the spirit of Dodson. The quasisoliton case is excluded using Morawetz.

Perturbation theory

Recalls the perturbation lemma from CKSTT, adapted to this problem.

She nicely describes the reduction to proving a local result on a time interval controlled by the energy. Once we have this type of local theory, we essentially convert the critical problem into one that behaves like the subcritical problem so GWP will follow.

Remarks on Proof

Subcritical quadratic terms in the Duhamel-Strichartz analysis on local intervals have a time factor. If this time factor is small enough, these subcritical terms can be absorbed. Oh, now I understand! The point here is that GP can be viewed as the energy-critical NLS plus some quadratic terms which don’t destroy energy conservation. This perspective guides the KOPV analysis. They show that the GP equation can be treated as a perturbation off the dilation invariant energy critical case.

Cubic-Quintic NLS with non-vanishing BC on $R^3$

They write $u=1+v$ and observe that $v$ satisfies energy critical NLS with subcritical lower order terms. The Hamiltonian is not sign definite so does not provide coercive control over the kinetic energy term. This is compensated for by using a lower order term $M(v)$, the $L^2$ norm of the real part of $v$. This quantity is not conserved. They show that it satisfies a Gronwall type estimate and that turns out to suffice.

Scattering for the GP equation in the case of large data

GP equation has traveling wave solutions that do NOT scatter.
Formation of traveling waves require a minimal energy in $R^d, d \geq 3$. Bethuel-Gravejat-Saut 2009, de Laire 2009.
Solutions with sufficiently small energy scatter. (Gustafson-Nakanish-Tsai 2006)
Can one prove scattering up to the minimal energy of a traveling wave?

Our goal is to fill in the gap. But, this problem does not seem too easy to attack, so we tried to apply these ideas on a simpler problem.

Killip-Oh-Pocovnicu-Visan

For a Cubic-Quintic NLS with zero boundary conditions (which has conserved mass and energy and has soliton solutions) the are working to show that if $v_0 \in H^1 (R^3)$ then scattering holds true if the mass is smaller than the mass of any soliton OR if it has positive energy, smaller than the enrgy of any solution.

(Final statement is a work in progress.)

Princeton Math Colloquium: Soliton Resolution for the Radial Energy Critical Focusing Nonlinear Wave Equation

Wed, 14 Mar 2012 18:45:45 +0000

(Apologies for typos or misquotations…; please help me improve the post. –J. Colliander)

Princeton Math Colloquium 2012-03-14

Presenter: Carlos Kenig

University of Chicago

Title: A case study for critical non-linear dispersive equations: the energy critical wave equation

Abstract: We will discuss recent work on the energy critical wave equation. The issues studied are global existence, scattering, finite time blow-up, universal profiles at blow-up and soliton resolution. This is viewed not as an isolated series of results, but as a way of approaching many similar critical non-linear dispersive equations.

We hope this is a model case for how nonlinear dispersive wave equations should proceed. This is a model case. Techniques of harmonic analysis have been introduced into the subject. The focus of those studies concentrated on local well-posedness. In the last fifteen years, there has been a new concentration on global-in-time, blowup aspects.

Linear wave equation

We start with a review of the linear wave equation.

$$ (LW) \partial^2_t w - \Delta w = h, data (w_0, w_1).$$

We write the soluton as $$ w(t) = S(t)(w_0, w_1) (t) + D(h)(t). $$

Finite speed of propagation. Picture of a cone.

An important estimate in these studies is the Strichartz estimate.

Critical focusing NLW

$$ \partial^2_t u - \Delta u = u^5, data (w_0, w_1). $$

Defocusing has $ - u^5 $. Reviews scaling and explains criticality.

Small data theory for NLW

Small (in $\dot{H}^1 \times L^2$) initial data evolves uniquely into a global-in-time solution which asymptotoically approaches (in same topology) into a linear solution as $ t \rightarrow \pm \infty$. Moreover, for any data, we have short time existence and a maximal time interval of existence.

This problem has an energy which is constant throughout the solution lifetime. Focusing vs. defocusing. There is competition, in the focusing case, between these two terms.

In the defocusing case, work of Struwe, Grillakis, Shatah-Struwe, Bahouri-Shatah proves that for any finite energy data, the solution exists globally and scatters. Another factor, Grillakis showed persistence of higher regularity.

In the focusing case, this fails. Levine 1974 showed that if $E(u_0, u_1) \leq 0$ then the solution can not exist globally in time (in either direction). (This is done by obstruction.) Recently, Kreiger-Schlag-Tataru 2009 constructed solutions for which $T_+ < \infty$ but for which the energy norm remains bounded. Also, in the focusing case, the elliptic equation admits a nonnegative solution $W$. This is a solution which is independent of time. This solution satisfies $\Delta u + u^5 = 0$. This solution is called the ground state. The reason for this terminology is that it arises as the optimizer in a Sobolev inequality and is charaterized variationally. This elliptic equation has been studied in connection with the Yanage problem in differential geometry. The formula for $W$ is explicit $$ W(x) = \frac{1}{(1 + |x|^2/3)^{1/2}}.$$ $W$ is the unique nonnegative solution of the elliptic equation (Gidas-Ni-Nirenberg 1979) and the only $\dot{H}^1$ solution (Pohozaev 1965). It does not scatter to a linear solution. It is “non-dispersive.” Donninger-Krieger 2012 have constructed global-ini-time solutions which are bounded in the energy class, are radial, and don’t scatter to either a linear solution or to $W$. There are other objects.

Recall a collection of recent results

Theorem (Kenig-Merle 2008): If $E(u) < E(w)$ then: 1. If the gradient is smaller, we have global existence and scattering. 2. If the gradient is bigger, we have breakdown in (both directions of) finite time. The case of equality is impossible under the hypothesis.

The proof of this result is an application of a general method. We call this a concentrated compactness rigidity method. Provides a brief summary of these ideas: small data, variational aspects, critical element extraction, rigidity.

Theorem (Duyckaerts-Merle 2008): There exists $W_-, ~ W_+$ radial which have the same energy as $W$. Then, there are three statements which provide a characterization of the dynamics.

(original statement of DM had an extra hypothesis which is now removed with inputs from Duyckaerts-Kenig-Merle 11, Krieger-Nakanishi-Schlag 11)

Now, we need to go beyond the energy threshold of $W$.

Existence of Type II blowup solutions

These solutions exist for a finite time but their critical norm remains bounded.

Krieger-Schlag-Tataru 2009

Describes a singularity formation process along a rescaled $W$ profile plus a continuous remnant.

We next show this is a “universal” phenomena.

DKM 2009, 2010:

(Radial) Any Type II blowup (critical norm stays bounded) solution with just a little bit bigger size than $W$. Then the soution looks like a rescaled $W$ plus errors. All type II blowups look like rescaled $W$.

Corresponding nonradial result: Slightly larger (by $\eta_0$) in the $\sup$ sense as above. After a rotation and a translation, we find that the solution looks like a rescaled and translated $W_l$ plus a small error. The rescaling rate is understood. And we have $$\lim_{t \uparrow 1} \frac{x(t)}{1-t} = l e^1 .$$

So, the possible profiles in the nonradial case are $W$ or its Lorentz transformations.

There also exist explicit solutions which blow up in finite time. These are built with ODE techniques and are independent of $x$. You can then chop them up with finite speed of propagation. These solutions have exploding critical norm. These are called Type I solutions. Conceivably, there are solutions which are Type I along one sequence of times and Type II along a different sequence of times.

DKM 2011: $W_+$ turns out to be Type I.

Soliton resolution for radial solutions of NLW

For a long time, there has been a widespread belief in the math physics community that large, global in tie, solutions of dispersive equations, asymptotically in time, they decouple into a sum of finitely many modulated solitons, a free radiation term and a term that goes to zero at infinity. This is a kind of philosophy that guides the research. So far, this has only been proved for the integrable KdV equations. Eckhaus-Schurr carried this out using the completely integrable structure. Corresponding results for mKdV can be obtained via the Miura transform. Heuristic arguements for this conjecture, in the case of cubic NLS in 1D, were given by Ablowitz-Segur 1976 and Zakharov-Shabat 1971. These are all subcrticial equations, for which one expects that these decompositions are stable, unlike in critical equations.

For more general equatiosn, so far, results have been found for data close to the soliton, in subcritical nonlinearities due to several authors. Buslaev-Perelman 1992 NLS, Soffer-Weisntein 1990, Marterl-Merle 2001 for gKdV) For corresponding results with blowup there are works by Merle-Raphael, Martel-Merle. There have also been large solution results for critical equivariant wave maps onto the sphere due to Christodoulou-Tahvildar-Zadeh, Shatah-T-Z, Struwe. They sow convergence along some sequence of times converging to the blwup time locally in space to a soliton (harmonic map).

In the finite time blowup case, for the 1d nonlinear wave equation, Merle-Zaag have obtained reusults of this kind using a Lyapunov functional tool.

In critical elliptic problems, such as the ones mentioned earlier, in domains excluding a small ball, considering radial solutions, there have been obtained results on decompositions inot “towering bubbles” (the analog of a finite sum of modulated solitons), as the size of the ball goes to zero. Musso-Pistoia 2006

The first general results for radial solutions of NLW were for type II solutions and held only for a sequence of times. We now have the full soliton resolution for radial solutions of NLW.

Theorem: Let $u$ be radial solution of NLW. Then, one of the following holds true:

Type I blowup.
Type II blowup with a full decomposition as a finite sum rescaled $W$, with ordered scaling speeds.
Radiation plus a sum of modulated solitons.

Consequence: Any solution which exists globally in time is bounded.

Remark 1: When the existence time is finite, the limit of the norm exists. It is either divergent or it is bounded. There are no mixed asymptotics. We do not yet have solutions which require more than one bump but we expect that they exist. For 1d NLW similar constructions were made by Cote-Zaag 2011. See also Musso-Pistoia 2006.

Remark 2: Quantifications. Each bump absorbs some energy.

Discussion….no claim of stability. This is not generic.

Ideas of the proof (global case)

The fundamental new ingredient of the proof is the following dispersive property that all radial solutions to NLW (other than 0 and $\pm W$ ) must have: $ \exists ~ R >0,, ~ \eta > 0$ such that for all $ t \geq 0$ or all $ t \leq 0$ $$ \int_{|x| > R + |t|} |\nabla_{x,t} u (x,t)|^2 dx \geq \eta. $$

Notes on Nonlinear Dispersive Wave Equations Workshop in Oberwolfach

Mon, 04 Oct 2010 22:39:37 +0000

This page contains notes by J. Colliander taken at the workshop:

Nonlinear Waves and Dispersive Equations
Organizers:
Date: September 12th - September 18th, 2010

I apologize for any mistakes! If any of the speakers would like me to post (or link to) their slides, please send me the file. –Jim Colliander

Jérémie Szeftel: Alas, I missed this talk….thanks Air Canada. ↩

arXiv ↩

slides ↩

Sebastian Herr: Small data theory for energy critical periodic NLS ↩

(joint work with Tataru and Tzvetkov)

Energy critical NLS focusing or defocusing on a manifold M. Specific examples with Laplace Beltrami operator. Mostly intersted in manifolds with periodic geodesics. For example $\mathbb{T}^3$ or tori crossed with $\mathbb{R}^d$.

Target is LWP.

Warm-up remarks ↩

Warm up: $M = {\mathbb{R}^d}$. Strichartz, dual Strichartz, dispersive decay $\implies$ (Cazenave-Weissler) LWP.
Non-Euclidean cases: Asymptotically Euclidean and nontrapping metrics have been studied.
Failure of sharp Strichartz estimates on torus and on sphere.
Trapping creates geometric obstructions to dispersion.
Trapping can create instabilities and failure of Strichartz estimates.
Known estimates: Strichartz with a loss of derivatives.

$NLS^{\pm}_5 (T^3)$ ↩

Available estimates have some loss. The loss obeys the scaling but it is insufficient to control the quintic nonlinearity. We end up needing an $L^4$ estimate, which is unavailable.
Our strategy is to use multilinear, scale invariant versions of Strichartz estimates to better share the derivatives.
Use almost orthogonality wrt spacetime to reduce estimates to smaller scales.
Replacements/refinements of the $X^{s,1/2}$? We use the critical function spaces $U^p, V^p$.
We will need refinements of these spaces which are sensitive to finer than dyadic frequency localizations.

New Strichartz Estimates ↩

We have the Strichartz estimates on functions supported on cubes in Fourier space
For all rectangular sets of arbitrary orientation and center, we get a better bound!
This boils down to a classical estimate (Landau 24) for counting the number of lattice points on 6d ellipsoid.

Perturbative Analysis ↩

$U^p$: Definition involving all partitions of the line using $U^p$-atoms.
These are Banach spaces which embed into $L^\infty$.
$V^p$: We need another type of space. These are functions of finite $L^p$ variation over the partitions of the line.
$U^p \rightarrow V^p_{rc} \rightarrow L^\infty$ (Embeddings)
$\| u \|{U^p{\Delta} H^s} = \| e^{-it \Delta} u \| {U^p (R; H^s})$. (Similarly wrt $V^p$, as in Ginibre’s Asterisque.)
We choose then $p=2$ and call the resulting spaces $X^s$ and $Y^s$.
Properties: $U^2_{\Delta} H^s \rightarrow X^s \rightarrow Y^s \rightarrow V^2_{\Delta} H^s$ (Embeddings)
We define restrictions to smaller time intervals….
$X^s$ and $Y^{-s}$ have a nice duality relationship.

Trilinear Strichartz ↩

Refinement which generalizes Bourgain’s $p=6$ Strichartz estimate.

Sketch of proof ↩

Decompose the largest frequency $N_1$ annulus in cubes of the second largest frequency $N_2$.
We can replace $Y^0$ by $V^2_{\Delta} L^2$.
We deduce control on the quintic nonlinearity using the trilinear estimate. Some gain is obtained by playing with the exponent $p$ in the $U^p$ spaces, which he attributed to elementary properties of these atomic spaces.
This gain and some other slack in the other trilinear estimate allows one to sum up over the dyadic scales.

I am confused at this point? Do we have some derivative slack or are thing really tight? Since we are considering an $H^1$ critical problem, there can be no slack….I discussed this with Sebastian a bit after the talk. I was confused; there is no derivative slack.

Next, there is a new localization (the rectangle decomposition). The cubes are decomposed into almost disjoint strips of a certain width. The almost orthogonality is gained from the temporal frequency! (This reminded me of the ideas from Koch-Tzvetkov and later developed by Ionescu-Kenig)

Contraction estimate ↩

It is not necessary to use the rectangles to get this estimate. For the qunitic case, we can avoid the rectangles. For the cubic NLS, by duality you have a 4-linear estimate and by Cauchy-Schwarz you are reduced to bilinear estimates. For the cubic case, it is necessary to use the rectangle decomposition.

Remarks ↩

With similar ideas, they can treat the cubic case on $R^2 \times T^2$ or $R^3 \times T$.
This involves bilinear refinements instead of cubic refinements.
small data GWP for energy critical NLS on certain manifolds where arguments of the Euclidean setting fail.
Large data is a very interesting problem.
This is the first critical result for NLS on a compact manifold.
Quintic NLS on the 3-sphere? Strichartz estimates fail but possible to control second Picard iteration.
Cubic NLS on $T^4$.

Questions ↩

Flat waveguides?
$L^2$ critical case?

Benjamin Schlein: Effective evolution equations from many body quantum dynamics ↩

Resources: Schlein’s talk at ICMP 2009, Schlein’s Zurich Lectures

Introduction ↩

Consider $N$ particles moving in 3d. These particles can be described in quantum mechanics by a wave function $\Psi_N \in L^2 (R^{3N})$. The probability density $| \Psi_N (x_1, x_2, …, x_N)|^2$ represents the probability of finding particle 1 at location $x_1$ and so forth. Bosons are symmetric wrt particle interchange. Fermions are antisymmetric. We will restrict in this talk to Bosonic symmetry: For all permutations $\pi$, $$ \Psi_n (x_{\pi_1}, \dots, x_{\pi_N}) = \Psi (x_1, \dots, x_n)$$

The dynamics of the wave equation is governed by the Schrödinger equation $$ i \partial_t \Psi_N = H_N \Psi_N $$ Typically, $$ H = \sum_{j=1}^N (-\Delta{x_j} + V_{ext} (x_j)) + \lambda \sum^N V(x_i - x_j). $$ We have well-defined local dynamics. The problem is that we have way too many particles in typical physical systems. We want to find effective descriptions of the dynamics. In certain regimes, we can approximate this complicated but linear evolution using effective equations

Mean Field Regime

The particles interact with many other particles. The strength of each of these many interactions is small so that the effect of all of them is of order 1: $N \gg 1, \lambda \ll 1$. We will assume that $N \lambda \sim 1$. The dynamcis generated by the mean field Hamiltonian: $$ H^{mf} = \sum (-\Delta{x_j} + V_{ext} (x_j)) + \frac{\kappa}{N} \sum^N V(x_i - x_j). $$ We study the dynamics emerging from a product wave function: $$ \Psi_N (x_1, \dots, x_N) = \prod_{j=1}^N \phi (x_j) $$ Because of the interactions, we can’t expect that the product wave function remains of product form. But, in the mean field case, we might expect that $\Psi_N (t) \sim \phi(t)^{N}$. If we assume this, we obtain a self-consistent Hartree equation. Here is the heuristic step: $$ \frac{\kappa}{N} \sum^N V(x_i - x_j) \sim \frac{\kappa}{N} \sum^{i \neq j} \int V(x_i - y) |\phi(y)|^2 dy \sim \kappa (V * |\phi(t)|^2) (xj). $$

Reduced Densities

$$ \gamma_N (t) = |\psi_N (t)\rangle \langle \psi_N (t)|$$
Partial traces
- When we take partial traces, we lose some information. It is integrated out. However, we are only interested in the data that can be extracted based on measurements of finitely many particles.

Theorem: (Under suitable assumptions on $V$). Let $\phi \in H^1 (R^3), \Psi$ a pure product wave function, $\Psi_N (t)$ the linear evolution of the many body system. The for all fixed $k \in {\mathbb{N}}, t \in R$, the reduced density matrices converge to the projectors build on the $\phi$ evolutions where $\phi$ solves the Hartree equation.

The more singular the potential, the more difficult it is to prove the theorem.
Spohn 1980: proved this for bounded $V$.
Erdös, Yau 2000: $V(x) = \pm \frac{1}{|x|}$.
Rodnianski, Schlein 2008: $V(x) = \pm \frac{1}{|x|}$, gives quantitative convergence with control by $\frac{C}{N}e^{kt}$.
- The RS work was based on an approach by K. Hepp.
- The approach is based on a representation of the problem on Fock space.
- Coherent states and quantum field theory ideas.
Knowles, Pickl 2009: Improved to more singular potentials.
Grillakis, Machedon, Margetis 2009 I, II: Second order corrections to the mean field dynamics, giving norm convergence.

Boson Stars ↩

$N$ particle Hamiltonian $$ H = \sum^N \sqrt{1 - \Delta {xj}} - G \sum \frac{1}{|xi - xj|}$$
$N \gg 1, G \ll 1, NG = \kappa$
$\forall N \exists \kappa(N)>0$ (kappa critical) such that:
- $\inf \frac{\langle \Psi , HN \Psi \rangle }{\| \Psi\|_2^2} = 0 ~if~ \kappa \leq \kappa(N)$
- $\inf \frac{\langle \Psi , HN \Psi }{\| \Psi\|^2} = - \infty ~if~ \kappa \geq \kappa(N)$
- Lieb, Yau proved that $\kappa(N) \rightarrow \kappa^H$ as $N \rightarrow \infty$
Look at the corresponding effective field equation.
- For $\kappa \leq \kappa^H$, we have global well-posedness.
- For $\kappa \geq \kappa^H$, there exists finite time blowup solutions (Fröhlich-Lenzmann 2006).

Theorem (Michelangeli-Schein 2010): Let $\phi \in H^2 (R^3)$ and form product wave function $\Psi_N$ and let $\Psi_N(t)$ evolves according to the regularized Hamiltonain (where the singularity is tamed by adding a small positive term to denominator which vanishes as $N \rightarrow \infty$). If we have $H^{1/2}$ control on the nonlinear level by constant $k$ over a time interval $[0,T]$ then we have convergence. Moreover, if the nonlinear problem explodes then the energy per particle in the linear problem also blows up. (Hypotheses were a bit strange to me….I asked about it after the talk and need to look at the paper.)

Dynamics of Bose-Einstein Condensates ↩

Drop the external potential. Effective dynamics in this case is described by the Gross-Pitaevskii equation: $NLS_3^+$ The derivation of effective dynamics in this setting has only been established for the defocusing case.

Adrian Constantin: Camassa-Holm ↩

Physical Background ↩

2d water waves over a flat bed. He draws a curve above a flat bottom at $y = - h_0$ and the free surface is given by the graph $y = \eta (x,t)$. He writes the Euler equations, mass conservation, imposes reasonable boundary conditions. These are generally accepted to be the right model. I will work with one other assumption: $u_y - v_x = 0$: irrotational. There are various scales you can plug into the problem and then you can non-dimensionalize. The problem can then be written in terms of just two parameters $\epsilon$ and $\delta^2$ where: $$ \epsilon = \frac{2}{h_0}, $$ $$ \delta = \frac{h_0}{\lambda} $$ Small amplitudes $\epsilon \ll 0$ and $\delta$ is the shallowness parameter so shallow water wave theory means that $\delta$ small. Shallow water small amplitude is when $\delta \ll 1$ and $\epsilon = O(\delta^2)$. If you study this, you get KdV and BBM equations. In this regime, these model equations enjoy global existence. The nondimensional form of the KdV equation is $$ \eta_t + \eta_x + \frac{3 \epsilon}{2} \eta \eta_x + \frac{\delta^2}{6} \eta_{xxx} = 0. $$ Here is the emerging BBM: $$ \eta_t + \eta_x + \frac{3 \epsilon}{2} \eta \eta_x + \delta^2 (\beta + \frac{1}{6})\eta_{xxx} - \beta \delta^2 \eta_{xxt} = 0, \beta \geq 0. $$ KdV is completely integrable and has solitons. BBM has some nice analytic features but only 5 conserved quantities. These derivations are $O(\delta^4)$.

Where does Camassa-Holm come into this business? Since all the waves that are physically reasonable, we have global existence. We would like to have a simple model that captures the phenomenon of wave breaking: $\eta$ is bounded, $|\eta_x|$ becomes unbounded in finite time. (This is described as a desirable extension in the book Linear and Nonlinear Waves, by Whitham.)

Emergence of Camassa-Holm Equation ↩

Moderate amplitude (shallow water): $\epsilon = O(\delta), \delta \ll 1.$

Johnson found a path like this to see Camassa-Holm emerge. (“Unfortunately, the original derivation of that equation was not correct.” “They assume that $\epsilon$ is small and later they assume that $1/\epsilon$ is small…”) You can derive an equation for the horizontal velocity at a particular depth. If you do this derivation at depth $\frac{1}{\sqrt{2}} |h_0|$. Another equation called Degasperi-Procesi emerges when you consider this at a different depth: $$ u_t - u_{txx} + 3k u_x + 4 u u_x = 3 u_x u_{xx} + u u_{xxx} $$ Both of these equations are integrable! (I did not know about the D-P equation before…) Both of these equations have solution which break down, in the fashion of wave breaking described above.

For CH, we have some conservation laws which implies the solution stays in $L^\infty$. Fokas and Fuchsteinner found the CH equation in a list of 12 equations that are the only completely integrable equations. It is difficult to use the integrable systems machinery to study the wave breaking. For DP, if you start with data in a nice enough space (say $H^{3/2}$) and you can then prove the solution stays in $L^\infty$.

Tzvetkov Q: Can the solution be extended after the wave breaks? A: You might be able to extend the solution like shocks. But the relevance of the wave breaking event in CH, it is not clear whether the CH is a good approximation of the Euler equations. Therefore, even if the PDE theory for CH can be extended, this does not mean you have a relevant extension modeling the water wave problem.

Geometric viewpoint as a geodeisc on the diffeomorphism group ↩

There is this famous paper of Arnold that shows that Euler may be viewed as a geodesic flow on the diffeomorphism group. CH and KdV can be similarly interpreted as a geodesic flow on the Bott-Virasoro algebra. This geometry thing is very nice, very appealing. However, this geometric point of view does not give a useful consequence from the viewpoint of analysis.

The best result for CH is that when the solution does not change sign it stays global. This is built from Nöther’s theorem, which provides a different view on the CH equation

Write ${\mathcal{D}} = [ \phi: C^\infty ~orientation~ preserving~ diffeos ]$. This is a Lie group and the tangent space at the identity is $C^\infty (S)$. We can then move this tangent plane around using Lie algebra properties by right-translating. The geodesic equation looks like $\phi_t = u(t, \phi)$ where $u \in \mathcal{D}$.

If I do this for $L^2$, I get $u_t + 3 u u_x = 0$. However, the Riemannian exponential map $exp_R$ is not a local chart.
If I do this for $H^1$ ( I believe this is referencing the Riemannian structure imposed on $C^\infty$ diffeos) we get CH.
Consider the Bott-Virasoro Algebra $Vir = C^\infty \times {\mathbb{R}}$ and you do some Bott cycle thing which looks like a diffeo flow with a twist, you get KdV. (This is a result of Olsheyenko(?) and Khesin.)
DP equation also has some interpretation this way but it is more complicated.

Integrable Sturcture ↩

CH Lax Pair. This is an isospectral problem. For CH, we have $\psi_{xx} = \frac{1}{4} \psi - \lambda m \psi, m = u - u_{xx} + k$. This is a weighted spectral problem. If $u$ solves CH, then the eigenvalues of this equation are time independent.
DP Lax Pair. $\psi_{xxx} - \psi_{x} - m z^3 \psi = 0, z \in {\mathbb{C}}, m = u - u_x +k$. When $m$ is strictly positive, we can perform certain Liouville substitutions which allow us to recast this as a regular Sturm-Liouville problem.

This talk did not properly survey the literature. Instead, the talk was intended to highlight a physically relevant derivation of the CH and to show that this equation is mathematically and physically interesting.

Tzvetkov asks: Is there a Miura transformation? Seems to be no….although AC appeared to me to answer a different question.

Merle asks: Can you track the blowup using the integrable machinery here? We need estimates on the eigenvalues hold if $m>0$ and you can give examples which show that sign-changing $m$ breaks down the needed estimates on the eigenvalues.

Ponce asks: What is the best LWP theory for CH? Answer: Kato’s theory needs $H^{3/2}$. You have existence, uniqueness and continuous dependence in $H^1$, but the continous dependence is weaker.

Ponce asks: Is the peakon stable? A: Yes, this is a result of Molinet and El Dika.

Claudio Muñoz: Dynamics of gKdV solitons under perturbations by potentials in front of nonlinear term ↩

My computer ran out of battery….

arXiv: Muñoz on KdV

arXiv: Muñoz on NLS

Mihalis Dafermos: Superradiance, trapping and decay for waves on Kerr spactimes in the general subextremal case $|a| < M$. ↩

(joint work with Igor Rodnianski)

Kerr family $(0 \leq |a| \leq M)$ of metrics (in Boyer-Lindquist coordinates) $$ g_{M,a} = - \frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2 + \frac{\sin^2 \theta}{\rho^2}(a dt - (r^2 + a^2) d\phi)^2. $$ Here $\rho^2 = r^2 + a^2 \cos^2 \theta, \Delta = r^2 - 2 M r + a^2 = (r - r_{-})(r-r_{+}), r_{+} \geq r_{-}.$ This is a vacuum solution $(R_{\mu \nu} =0)$ and has Killing fields $\partial_t, \partial_{\phi}$. The domain of outer communications is $r > r_{+}$. The case $a=0$ is Schwarzschild 1916. The Kerr case is $a \neq 0$ and was discovered in 1963.

Penrose diagram for Kerr $(0 < |a| N M)$. What is a black hole? A spacetime has a black hole if the past of this infinity is not the entire spacetime. Both Kerr and Schwarzschild are expected to be unstable and the structure of the singularity should be something in between.

Penrose Diagrams (images taken from Dafermos-Rodnianski)

Penrose diagram of Schwarzschild spacetime

Penrose diagram of Kerr spacetime

Boundedness and decay for $\square_g \psi =0$ on Schwarzschild and Kerr ↩

These are natural questions from several points of view. One important application of these ideas is to address the stability properties of these solutions of the Einstein equations.

Current state of the art for the quantitative study of $\square_g \psi = 0$ ↩

Boundedness in general class of $C^1$ stationay axisymmetric spacetimes [DR].
“Integrated local energy decay” for exactly Kerr:
1. Slowly Rotating Case $|a| \ll M$ [DR], Tataru-Tohaneanu, Andersson-Blue
2. $|a| < M$, [DR] this talk!.
Pointwise-in-time decay from 1. and 2. (energy based method [DR] based on resolvent method of Tataru)

We typically think of proving decay as a two step process. We first prove that some spacetime integral of energy to the future an arbitrary hyperboloidal space-like hypersurface controlled by the energy on the hypersurface. This type of result has been shown in the slowly rotating case. Once we have items 1. and 2. (boundedness and integrated local energy decay) we can use the vector field method to get pointwise-in-time decay. These methods appear to be robust and might be applicable to nonlinear problems.

Review of the main features of Kerr Spacetimes ↩

Red-shift (associated to the event horizon)
Superradiance
Trapping (trapped null geodesics)

Red-shift

Two observes move in spacetime. You think of observer A emitting constant frequency signals and you imagine these being received by observer B so the frequency is shifted to the red. First discussed in 1939 by Oppenheimer-Snyder. Extremal case $a = M$: The red-shift factor at the horizon vanishes. The positivity of the surface gravity is a geometrical object underlying the red shift.

Superradiance

In Schwarzshild, the killing v.f. $\partial_t$ is timelike n the exterior becoming null on the horizon. Thus there is a conserved (by Nöther) non-negative definite energy by the time-like condition. The only subtlety is that the energy degenerates at the horizon.

In stationary perturbations of Schwarzschild, $\partial_t$ in general becomes spacelike near the horizon. this happens already for Kerr with $0 \neq |a| \ll M$. The corresponding energy is conserved but does not have a sign. For particle motion, this leads to the so-called Penrose Process. For waves, this leads to the phenomenon of superradiance (Zel’dovich).

In particular, using the conservation law associated to $\partial_t$ one cannot prove a priori boundedness, even away from the horizon. The energy radiated to null infinity might be bigger than the initial energy and this is called Superradiance.

For Schwarzschild, the only trouble is near the horizon because we have a useful energy control for radiated energy to null infinity. For Kerr, we don’t have that because of the superradiance phenomenon and this creates new difficulties. We need to prove boundedness and decay everywhere, not just near the horizon.

Trapping

On Schwarzschild, the photon sphere $r = 3M$ has the property that it contains null geodesics. These null geodesics thus neither escape to null infinity nor to the horizon. In Kerr, the behaviour persists, but it is more complicated! It is not obviously located in physical space but can be thought of more easily inside phase space. One can concentrate energy for arbitrarily large times near trapped null geodesics. One has to capture this to prove dispersive results. In particular, pointwise-in-time decay estimates for energy must lose derivatives (Ralston).

Proof of integrated local energy decay. ↩

We will only discuss the first energy. Higher order estimates require commutation with the redshift vector field, the Hawking v.f. and $\partial_t$.

The method of proof will exploit energy currents.

In the large $a$ case, the construction of these currents will need to frequency localized for two reasons.

To distinguish between non-superradiant and superradiant frequencies.
To degenerate at the correct value of r.

A convenient way of doing both at the same time is frequency localizing via Carter’s celebrated separation of the wave equation. Kerr geometry only has two killing fields. This is not enough to separate the equation. However, there is some extra symmetry there that helps you. In view of Ricci flatness, this separability is equivalent to separability of Hamilton-Jacobi equations and the existence of a Killing tensor. These three objects are devices to extract this hidden structure.

Separation.

Big display….can’t keep up with that. We are studying $\square_g \psi = F$ (where $F$ arises from cutoffs) and we take $\widehat{\Psi}$ and rewrite it using some structure of an oblate spheroidal metric in the $\theta, \phi$ variables. The content of what Carter noticed is that when you do this, you can show that there is a hidden ODE lurking in this decomposition…..

More big display…working pretty hard here, lots of indices….new coordinate $r^{*}$ so that things look more like the Regge-Wheeler coordinates in Schwarzschild. With this decomposition, we can identify the superradiant frequencies. The superradiant modes should be thought of as the modes which send infinite negative energy through the horizon.

Completely separated energy current identies (analogues of $\nabla^\mu (T_{\mu \nu}(\psi) (y \partial_{r^{*}})^\nu)$, etc.)

Lots of notation with symbols I don’t know how to make….

General Idea

From the above currents, produce integral identities with positive definite underlined bulk terms and (upon summation) we get the integrated decay except in regions which can’t be handled this way, basically because these frequency ranges are associated with trapping.

Kerr for small $|a|$. The constructions for all the other frequency ranges can be easily perturbed to yield positive definite bulks. The boundar terms, however, are not a priori controlled, this is the problem of superradiance. Since $|a|$ is small, this can be remedied by adding on a small amount of the redshift identity. Basically, for the small $|a|$ case, we have very little superradiance and can control it with the red shift. For large $|a|$, we need another idea.

The key observation seems to be that superradiant frequencies are not trapped. You can accomodate the superradiance using this idea by using the red shift.

Remark 1: In the small rotation case, the relationship between superradiance and red shift was the key idea. Remark 2: There are no trapped null geodesics which are orthogonal to $\partial_t$. This is a phenomenon identified by Alexakis-Ionescu-Klainerman in their works on uniqueness properties.

Some other important results

Positive and negative cosmological case
Ohter equations, like Dirac, Maxwell instead of wave equation. (Blue, Hafner, Finster et. al)

Open Problems ↩

Extremal case $a=M$ (recent results of S. Aretakis)
Higher dimensions (Schlue, Laul-Metcalfe)
Other measures of decay, Strichartz, …
Robust additional decay
Maxwell equations on Kerr (Blue) (Earlier work by Blue on Schwarzschild)
Equations of gravitational perturbation
Nonlinear stability of Kerr?

Stephen Gustafson: Dynamics on near-harmonic Schrödinger and Landau-Lifschitz maps ↩

(w. Nakanishi, Tsai) The paper that precedes the new stuff here is posted.

Landau-Lifschitz (30s), magnetizations $u(t,x) \in R^3$ with a constraint $|u(t,x)| = constant. The Landau-Lifschitz equation is:

$$ u_t = a_2 u \times \Delta u - a_1 u \times (u \times \Delta u), a_1 \geq 0. $$

Broader context:

$u(\cdot, t): R^2 \rightarrow S^2.$
energy $E(u) = \frac{1}{2} \int_{R^2} |\nabla u|2 dx$
heat flow $u_t = \Delta u + |\nabla u|^2 u = - E’(u)$
Schrödinger Map: $u_t = u \times \Delta u = J E’(u)$, $J$ is a complex structure.
Landau-Lifschitz is a combination of these equations
Also related to wave maps

Regurlarity vs. Singularity: energy critical problems ↩

Energy is scale invariant in $R^2$. $E = \int_{R^2} |\nabla u |^2 \geq 4 \pi |degree (u)|$ iff Harmonic map. So, what can you say?

Heat flow:

$E< 4 \pi \implies$ global smooth solutions Struwe 1985.
$E > 4 \pi \implies$ singularities may form, follows from Chang-Ding-Ye 92 via subsolution construction.

Wave map:

$E< 4 \pi \implies$ global smooth solutions Sterbenz-Tataru 09
$E > 4 \pi \implies$ singularities may form, follows from Chang-Ding-Ye 92 via subsolution construction. Rodnianski-Sterbenz 06, Kreiger-Schalg-Tataru 08, Raphael-Rodnianski 09

Schrödinger Map:

$E$ small $\implies$ global Bejenaru-Ionescu-Kenig-Tataru 08
Open: larger energy, singularity, for LL also.

Equivariant Maps ↩

Simples setting: near harmonic, equivariant maps.

$m \in Z^+$ is the degree
$(r, \theta) $ are polar coordinates
$R = {\hat{k}}\times$ (rotation about ${\hat{k}}$)

2-parameter family of harmonic maps (at minimal $E = 4 \pi m$). The energy is constrained. We work in a small energy shell above the $4 \pi m$ level. This is a restrictive class but known blowups are in this class. For heat flow case, we only have blowups with $m=1$ but for wave maps we have examples with $m \geq 1$.

Theorem Gustafson-Nakanishi-Tsai 09: For $m \geq 3$, solutions are global and converge to a (nearby) harmonic map (asymptotic stability):

$$ { {| u(t) - H^{\mu} |{L^\infty}} } + { { { {a1}} E (u(t) - H^\mu) \rightarrow 0 ~(t \rightarrow \infty). $$

Remarks:

includes the pure Schrödinger map case $a_1 =0$
also for $m=2$ heat flow ($a_2 = 0$) in a symmetry sub-class:
- solutions are global and converge to a harmonic map family
- the parameters can drift, eg to give infinite-time blowup: $s(t) \rightarrow 0$. In particular asymptotic stability fails. (Here $s$ is the length scale of the harmonic map.)

He draws a picture. A $\delta$ neighborhood of the harmonic maps in the energy space, viewed as an infinite graph over $s>0$. For $m \geq 3$, solutions in the $\delta$ neigborhood move dynamically back down to the harmonic maps. For $m=2$, you can drift all over the place in the case of heat flow.

New results: global solutions for degree 2 (LL) with $a_1 > 0$. ↩

Setting, as above with dissipations.

Theorem GNT: For $m=2$, solutions are global and converge to the harmonic map family, but not to one particular harmonic map.

The harmonic map family parameter $\mu(t)$ does not have to converge in general. But it will converge if the initial perturbation has a slightly faster spatial decay: $$ u_{x_1} - u \times u_{x_2} \in |x| L^1 \implies \mu(t) \rightarrow \mu_\infty. $$
For $m=1$ and only for the heat-flow, finite time blowup can occur but not for more localized perturbations.

(We expect this holds for (LL) but we have no proof.)

Remark: New results of Bejenearu-Tataru for $m=1$ Schrödinger map of degree 1. Harmonic mapps are unstable, stable with more localization. You should think of the $m=1$ Schrödinger map case as the most delicate.

Standard “modulation theory” approach ↩

Take your solution and split it into the harmonic map piece plus a remainder. Rewrite things for this remainder term. You look at the linear part of this equation driving the remainder dynamics. Because of invariances of the equation, the linearized operator has zero modes. What would you do to kill the kernel? Choose the parameter at time $t$ so that the perturbation is orthogonal to the kernel of the linearized operator. This leads to an ODE on the parameter dynamics. It remains to get dispersive/diffusive estimates to prove that the parameter converges as time goes to infinity. Do we have these estimates?

Dispersive/diffusive estimates

The remainder $z$ is controlled by a derived quantity $q$, so the analysis is a bit indirect and we now fight to control this related quantity. The $L^2$ norm of $q$ measures the energy gap above $4 \pi m$. The (recast) remainder $q$ satisfies a reasonable nonlinear Schrödinger-(heat) equation. This equation has a potential which depends upon $m$. For $m>1$, we have a lower bound estimate on the potential $V$. In this case, we can get “Strichartz” estimates on $q$.

For $m \geq 4$, this standard approach works G-Kang-Tsai 08, [Guan-G-Tsai 08].
For $m \leq 3$, the orthogonality condition is incompatible with the desired $L^2_t$-decay estimates and the standard approach fails.
For $m \leq 2$,, the orthogonality condition makes no sense.
For $m=1$, we don’t even have an $L^2$-eigenfuction but rather a resonance.

A remedy for $m \leq 3 $ and its cost. ↩

Change the orthogoanlity condition. Instead of demanding that the remainder be orthogonal to the kernel of the linearized operator, you require that the remainder be orthogonal to a localized function (unrelated to the kernel). Of course, there is a penalty for this change. The parameter dynamics ODE transforms then to involve another term and we have different parameter dynamics. This extra term is analyzed in some way.
Solution 1: “Normal form” for $m=3$. By integrating by parts a few times, we get good control on a modified quantity $[\mu (t) - (\psi^s /s | q)]$. We need control on the correction term. For $m \geq 3$, the correction basically does nothing so things work as before.
Parameter drift for $m=2$ heat flow. In this case, the “normal form” correction need not be bounded. For the heat-flow case, it is possible to simplify up to converging errors and the nonintegrability of the correction term can be exploited to drive the blowup, blowdown and oscillation properties of the scale parameter $s(t)$.
Solution 2: Take $a_1 >0$ and exploit the dissipation “We need to somehow stop pretending that the Schrödinger and heat equations are the same…” In the dissipative case, we can extract some dissipation on the correction term. (Probably, this decay is not available in the Schrödinger case.) There are some factorization tricks where the operator is recast, some Duhamel tricks…and an iteration on dyadic time intervals where the time dynamics of the parameter $\mu(t)$ are updated. What emerges is an upper bound by $\log t$ on the parameter $\mu(t)$.

Conclusions ↩

Near harmonic dynamics for (LL) for degree $m \geq 3$.
For $m =2$, more complex behavior.
For $m=1$, do finite time singularities form? This is only known for the heat-flow.

Ioan Bejenaru: Near soliton evolution in 2d Schrödinger Maps ↩

(joint work with Tataru)

Much of this will be a deja-vu since it overlaps with Gustafson’s talk.

Schrödinger map, Heisenberg model in ferromagnietism or the conservative part of the Landau-Lifschitz equation. Energy Conservation, scale invariance $s_c = n/2$ which is the threshold for the well-posedness theory. The $n=2$ case is energy critical.

Main Question: Global regularity of smooth solutons?

Sulem-Sulem 86 established existence of local solutions for $s> [n/2] + 2$, and this was improved to $s> [n/2] +1 $ by McGahagan. An “easier” problem is the global regularity for “small” initial data. Chang-Shatah-Uhlenbeck 00, ….Bejenaru-Ionescu-Kenig-Tataru 08 established GWP of the SM for small data in the critical Sobolev space in $n=2$.

Small data now resolved. What happens for large data?

Large Data Theory ↩

The dynamics depend upon the target manifold. For the sphere target, the problem is called “focusing” and for the hyperbolic target, the problem is called “defocusing”. This terminology makes good sense for wave maps but is not as explicitly understood in the case of SM.

A key feature in these problems is played by the existence of solitons: $ u \times \Delta u = 0$, which are known as harmonic maps. there are no nontrivial finite energy harmonic maps for the hyperbolic target. There are nontrivial harmonic maps with finite energy for the $S^2$ target.

A SM which fails to be regular at one time, bubbles like a HM.

Main Conjecture: In the hyperbolic case, the problem is globally wellposed independent of the size of the data. In the spherical case, solutions emerging from data below $4 \pi$ will be globally wellposed while the problem with higher energy may develop singularities.

The above conjectre is known for the harmonic map flow Eells-Sampson 64, Struwe 85, Chang-Ding-Ye 92.

Singularity formation for the WM problem. Recent works RS08, KST08, RR10.

There is some progress on more general targets.

Equivariant Harmonic Maps on $S^2$. ↩

These are maps from the plan into the sphere. Think that the origin is mapped to the south pole. The point at infinity is mapped to the north pole. Think of the image of the positive x axis as a curve connecting the south and north pole. When you move in the domain around wrt theta one time, the curve connecting n and s pole moves around the sphere some number of times. Once you have these maps, you can fatten them up into a two parameter family of maps.

Basic setup for stability/instability ↩

Define the two parameter family of $m$-equivariant harmonic maps. If you have slightly more energy, then you float around the harmonic map family. But do you stay nearby a particular harmonic map or can you float far away? If you move far, the higher derivatives do not stay under control. We want to describe the trajectory of these maps.

Modulation Theory ↩

Linearize near a soliton, study the zero eigenvalue, and these solutions do not disperse. You want to get rid of this eigenvalue. There is room to do that because we have some choice about which soliton you linearize around. This approach has been developed by Gustafson-Kang-Tsai 06 and Gustafson-Nakanishi-Tsai 09. This has been pushed further recently but involves higher degree hypotheses. We have decided to concentrate on the $m=1$ case.

Theorem (Bejenaru-Tataru):

Let $m=1$ and $\gamma \ll 1$. The for each 1-equivariant initial data $u_0$ satisfying $ \| u0 - Q(0,1)^1 \| X \leq \gamma$, there exists a unique global solution $u$ so that $u - Q(0,1)^1 \in C(R,X)$ and $ \| u_0 - Q(0,1)^1\|_{C(R,X)} \lesssim \gamma$.

There exists a solution $u$ with the additional property that $ \| u(0) - Q(\alpha, \lambda)^1 \|$ …..ack……slide switched I could not keep up. It is an instability result with a large upper bound. It was not clear to me if the assertion was that this big drift really occurs but certainly this is suggested.

Ionescu-Gustafson-Bejenaru conversation: localizations of the perturbations can restore the stability for the heat flow case…

Frank Merle: Isolatedness of characteristic points for blow-up solutions of semilinar wave equation ↩

Background References (Incomplete)

(joint work with Hatem Zaag)

I want to give a talk about a series of works I have done with H. Zaag on the semilinear wave equation.

Semilinear Wave Equation, Blowup Surface ↩

$$ u_{tt} = \Delta u + |u|^{p-1}u$$ Here $p>1$. Let’s collapse to dimension 1. We have initial data $(u0, u1) \in H^1 \times L^2$.

Summary of the results:

Local Existence: We have local existence until a blowup time $[0,T)$.
Existence of blowup via ODE method. There is a more refined condition due to Levine: If a (not the same as mine) energy is negative then $T< \infty$.
The blowup phenomenon can be spatially localized. Therefore, as in the book of Alinhac, you can produce a blowup surface. The solution is well defined on all backwards cones behind the blowup surface.
Question: We want to understand the blowup surface. We don’t know anything about it besides that it is 1-Lipschitz.
A point is called non-characteristic on the blowup surface if the surface has smaller than slope 1 at that point so it does not touch the boundary of the light cone. Let us denote the set of characteristic points on the curve by $S$. The other points on the curve are non-characteristic and the set of such points is called $R$. Let us denote the blowup curve by $x \rightarrow T(x)$ so it is given by a graph $(x, T(x))$.

Caffarelli-Friedman 85: For $u_0 \geq 0, u_1 \geq 0$ and use monotonicity of the wave flow in 1 dimension to prove that $\partial_t u \geq (1 + \delta )|\partial_x u|$ and you can prove then that no characteristic points don’t exist.

This result is a bit misleading. We tried to prove the nonexistence of characteristic points and could not do it. So, we turned our attention to proving the existence of characteristic points.

Summary of Results ↩

Existence of characteristic points. There exist initial data $(u_0, u_1)$ which has $S$ nonempty.
$S$ is isolated, $R$ is open.
$T(\cdot)$ is $C^1$ on $R$.
The only way that a characteristic point can arise is like a “hat”. $T’$ from the right and from the left are well defined. (Alinhac has examples for quasilinear equations which can blowup at all points along a line segment of slope 1.) At points $S$ we have $T’$ of slope 1 on the right and slope -1 on the left.
At points along $R$, the solution is of one sign and points in $S$ are points where the solution changes sign.

Characteristic points are cusps along the graph of $T(\cdot)$.

A Lyapunov functional (Antonini-Merle).

He shows that the solution extends outside the light cone behind noncharacteristic points. This gives you the $T’$ well-defined (with same value from left and right) at a noncharacteristic point.

The talk was hard for me to type up and explain well….Frank emphasized that the proofs are quite intricate and not presentable in a linear fashion.

Ben Dodson: Defocusing $L^2$-Critical NLS ↩

References

Mass-Critical NLS ↩

$$ i u_t + \Delta u = \mu |u|^{4/d}u, u(0,x)= u_0 (x), x \in R^d$$

We concentrate on the defocusing case where $\mu =1$. This equation conserves the quantities

$M(u(t)) = \int |u(t,x)|^2 dx$
$E(u(t)) = \frac{1}{2} \int |\nabla u(t,x)|^2 dx + \frac{\mu d}{2(d+2)} \int|u(t,x)|^{\frac{2(d+2)}{d}} dx$

Strichartz Pairs $(p,q): \frac{2}{p} = d( \frac{1}{2} - \frac{1}{q} ), d \geq 3, p \geq 2.$

$$ A(m) = \sup [ { {| u |}{L^{2(d+2)/(d)}} } (R \times R^d): { {| {u0} |}_{L^2}} =M ] $$

Minimal Mass Blowup Solution Strategy ↩

Theorem (Tao-Visan-Zhang 08): If $u(t,x)$ is a minimal mass blowup solution on $I$ then $\exists ~x(t), \xi(t): I \rightarrow R^d, ~ N(t): I \rightarrow (0, \infty)$ and $$ u(t,x) = \frac{1}{(N(t))^{d/2}} e^{i x \cdot \xi(t)} Q_t ( \frac{x - x(t)}{N(t)}) $$ where $Q$ changes with time but ranges only in a precompact set. For any $\eta > 0, ~ \exists C(\eta) < \infty$ such that

$$ \int_{|x- x(t)| > \frac{C(\eta)}{N(t)}} |u(t,x)|^2 dx < \eta, $$

$$ \int_{|\xi- \xi (t)| > {C(\eta)}{N(t)}} |{\widehat{u}} (t,\xi)|^2 d\xi < \eta. $$

Theorem (Killip-Tao-Visan): To prove GWP it suffices to exlude three scenarios:

$N(t) \sim t^{-1/2}, ~t \in (0, \infty)$,
$N(t) =1 , ~ t \in (-\infty, \infty)$,
$N(t) \leq 1, {\liminf }_{ {t \rightarrow \pm \infty}} N(t) =0, ~ t \in (-\infty, \infty).$

Then he writes and doesn’t really explain……

$\int_{1}^{\infty} N(t)^3 dt < \infty$
$\int_{-\infty}^{\infty} N(t)^3 dt < \infty$

Collapse to $d=3$ for now.

Theorem (CKSTT 04): Interaction Morawetz Estimate

$$ | u(t) |^4_{ {L^4 (J \times R^3)}} \lesssim { {| u |^3_{L^\infty (L^2)}}} {| u |_{ {L^\infty(H^{1})}}} $$

He quotes some estimates from KVZ linking time integrated (over slabs) powers of $N(t)$ with Strichartz size on same slabs.

On LWP time intervals $J_k$ (defined by diagonal Strichartz norm of size $\epsilon$), we have $N(t_1) \sim N(t_2)$ on $J_k$.

Galilean Invariance Observations ↩

Using Duhamel formula, he claims that the galilean invariance $\xi(t)$ does not move too rapidly. This allows him to localize things near the frequency center and in this way tames the galilean invariance.

Planchon-Vega paper on interaction Morawetz describes why the interaction Morawetz estimate is galilean invariant. All these expressions involve galilean invariant right sides and left sides. He then explains that the Morawetz Action leading to the interaction estimate is galilean invariant.

This allows him to claim that $ i \partial_t (Iu) + \Delta (Iu) = |Iu|^{4/3} Iu + [|Iu|^{4/3} (Iu) - I(|u|^{4/3}u)]$ enjoys some nice control (if we could ignore the error term in square brackets). So, we turn our attention to the error term.

For $N \leq CK$, we have

$$| P_{ {|\xi - \xi(t)| > N }} u(t) |{ {L^2(L^6)}} \lesssim (\frac{K}{N})^{1/2} \rho(N) $$ where $\rho(N) \leq 1$ with $\lim{N \rightarrow \infty} \rho(N) = 0.$

$L^2_t$ interval decomposition induction argument ↩

Bust up $[-T, T]$ into small intervals on which we have good Strichartz control and….not clear what he is doing to me right here…. Sort the intervals into good and bad intervals where a bad interval is where $N(J_k) \geq \frac{\eta_1 (d) N}{2}$. He makes a crude estimate on the bad intervals and pays for them by adding up their contributions. On the good intervals, he changes the organization of the decomposition.

Either $$ \eta_1 (d) N \geq \sum_{J_k \subset G_j} N(J_k) \geq \frac{\eta_1 (d) N}{2} $$ or $G_j$ lies to the left of a bad interval or $G_j$ is ont he end of $[-T,T]$. This allows him to claim that the number of $G_j$ is bounded by $C(d) \frac{K}{N}$.

….not clear to me….but hopefully it will be after I work some more.

Decomposition of nonlinearity ↩

He expands the nonlinearity wrt the decomposition around the moving frequency center $\xi(t)$ and the moving spatial center $x(t)$. He dismisses some parts of the nonlinearity based on the induction hypothesis, and the smallness in $L^2$ on the frequency regime far from the moving galilean center. The bad term that remains needs further study.

Questions ↩

Colliander: What are the main new ideas beyond the works of Killip, Tao, Visan and Zhang? ###
- Galilean invariance taming trick.
- Barely slipping under the wire.
- Induction argument using $L^2_t$…
Colliander: And for lower dimensions?
- The critical spaces of Koch-Tataru $U^p, V^p$.
- Harder work on the decomposed nonlinearity due to the absence of the endpoint Strichartz estimate in $d=2$.
- $d=1$ is easier than $d=2$, which is a nightmare.

Postlude ↩

I had a nice conversation with Fabrice Planchon who reported having a longer discussion with Dodson in June. Fabrice suggested that the new elements are the galilean invariance trick, the induction argument exploiting the $L^2_t$ control on left and right side of the Duhamel estimate (only available in $d \geq 3$) and the role played by the time integrals of powers of $N(t)$. Technical difficulties in 3d emerge because the nonlinearity is not multilinear and the analysis there would be simpler if it were. In 2d, we have a nicer nonlinearity but the absence of the endpoint Strichartz estimate in 2d obstructs the $L^2_t$ induction argument. Dodson in fact uses the double endpoint! This conversation made me think that it might be a nice exercise to try to revisit the 2d argument under the false assumption that the forbidden endpoint Strichartz estimate and following the 3d strategy. Alternatively, there might be a streamlined (but incomplete) proof which exposes the strategy more cleanly if we assume somehow that the nonlinearity in the 3d case were multilinear.

Killip: Energy Supercritical Wave Equation in 3d ↩

Background References

Introduction ↩

$$ u_{tt} - \Delta u + u^7 = 0, ~ u: R \times R^3 \rightarrow R$$

$ E(u^\lambda) = \lambda^{-1/3} E(u)$ so that the energy does not control the small scale behavior. This is very alarming.

Definition: $$\cal{E} (t) = { {| u (t) |}^2_{ {\dot{H}^{7/6}}}}+ { {| u_t |}_{\dot{H}}^{1/6}}.$$

Theorem (Killip-Visan 2010): $\cal{E} (0 )< \infty $ then

$\cal{E}(t)$ diverges.
$u(t) - u^{\pm} (t) \rightarrow 0$ as $t \rightarrow \pm 0$, where $u^{\pm}$ is a solution of the linear wave equation.

Radial case was done by Kenig-Merle. Minimal blowup solutions have good spatial decay properties. This is really the main point of their work and ours. Two essential points in the KM work:

Radial Sobolev embedding: ${\dot{H}}^{7/6} \ni u \implies |u| \leq r^{-1/3}$.
$r u(r)= u_{out} (t-r) + u_{in} (t-r)$.
If the solution is small intially ${\cal{E}}(0) < \eta$ then scattering holds.

Scattering is equivalent to the finiteness of some spacetime Strichartz $L^{12} (R \times R^3)$.

Step 1: Minimal Criminal ↩

Keraani first proved the existence of minimal blowup solutions and there were used by Kenig-Merle. At each moment of time, this object has certain localization properties. It is frequency localized on a characteristic frequency scale $N(t)$ and is spatially localized near $x(t)$ at the Heisenberg dual scale $\frac{1}{N(t)}$. Here $N(t)$ defines a multiplier which captures 99% of the norm. Riesz’ interpretation of the Arzela-Ascoli theorem shows this object is precompact. This is why we call these objects almost periodic.

Ionescu: What does minimal mean?

Answer: Samllest ${\| {\cal{E}} (t)\|}{L^\infty{t}}.$

I can apply symmetries and subsequential limits to these minimal objects.

Step 2: Minimal Criminal satisfies one of three scenarios: ↩

$N(t) =1$ soliton-like
$N(t) \geq 1, N(t) \rightarrow \infty$ as $t \rightarrow \infty.$
Finite time blowup.

Step 3. No finite time blowup solutions. ↩

How could blowup occur? The norm lives on smaller and smaller sets. By finite speed of propagation, we can deduce that there is a point where concentration occurs. Suppose we have a minimal blowup solution. We then look at the backwards light cone. Outside the light cone, $u=0$ by minimality. We know that $u$ has finite 7/6 norm and it lives on a small set. But this means that the energy must go to zero and this means the solution is actually the zero solution so is not a finite time blowup solution.

Soliton and Cascade Solutions have finite energy.

Step 4. Solutions move more slowly than light speed. ↩

He draws a forward light cone. There is no energy at the apex. We know that the energy inside the ball defined by the light cone at time $T$ is bounded by $T^{1/3}$. Energy can come into the cone but nothing can go out due to light speed. This tells us that $$ \int_0^T \int_{|x| =t} |u|^8 dS dt \lesssim T^{1/3}. $$

This argument works well if $N(t)$ is not changing too fast. For varying $N(t)$, this can be shown to violate speed of propagation.
There are some other variations to get this nailed down.

Step 5. $L^p$ decay. ↩

We are worried that our super smooth function does not decay fast enough.

$\dot{H}^{7/6} \rightarrow L^9$ (embedding), but we can actually prove that it is in $L^6$. At any time, we can represent $u$ using a Duhamel formula: $$ u(0) = - \int_0^\infty \frac{\sin(t|\nabla|)}{|\nabla|} u^7 (t) dt. $$ You can use the energy-flux identity to turn this into the $L^6$ control. How? You split low into high and low frequencies. We are only afraid of the very low frequencies. $(u_l + u_h)^7$ so $u_l$ is small and some interpolations give you the control.

Why? There can be no other term at null infinity since we would be wasting stuff and this would not be minimal.

Step 6. A more quantitative $L^p$ estimate. ↩

$$ \int_{|x - x(t)| \geq R} |u|^8 dx \leq R^{-\gamma}. $$

Split the time interval $[0, \infty]$ into $[0, R/3]$. We can then set up a geometric bootstrap. Everything is fine on the short time interval. If we look far into the future, we get smallness in $L^\infty$ and can then interpolate against the $L^6$ control to get the target $L^8$ estimate.

Step 7. Climax $E(u) < \infty.$ ↩

We gain regularity.

Write the $H^1$ norm as an inner product: $\langle \nabla u(0), \nabla u(0) \rangle + \langle u_t, u_t \rangle$. Now play with Double Duhamel. We have a kernel decay like $|t - s|^{-1}$ which will not converge when integrated over $dt ds$. We would actually need $|t - s|^{-2 - \epsilon}.$ We introduce a Whitney ball. We decompose $R^3$ in Whitney balls w.r.t. the origin. The negative powers of $R$ gained above from the quantitative decay estimate allows us to sum over the Whitney balls.

This shows the energy is finite, after a lot of bookkeeping.

How do we use this to wrap things up and prove the theorem.

Step 8. Completion of Theorem ↩

No Soliton: $\frac{x}{|x|} \cdot p$ leads to the Morawetz identity which implies the estimate:

$$ \int \int \frac{|u(t,x)|^8}{|x|} dx dt \lesssim E(u).$$ This kills the soliton.

No Cascade:

Using the Whitney balls slack, we can in fact get tightness:

$$ \int \langle x - x(t) \rangle^\epsilon [|\nabla u|^2 + |u_t|^2] dx \leq \infty.$$

Questions/Comments: ↩

Nakanishi: Do you have the same result if you have bounded critical Sobolev norm hypothesis is only true in one direction of time?

Killip: If this nemesis existed, then I can time translate it to create a nemesis that I have just shown can not > exist. So, I believe this relaxed hypothesis can be made with the same conclusion.

Colliander: Peter Pang (an undergraduate at U. Toronto) has recently numerically simulated this problem in the radial case and observed that the critical Sobolev norm remains bounded and is not monotone in time.

Colliander: Can you relax the bounded critical norm hypothesis to one with very slow, say logarithmic, growth and maintain the scattering conclusion?

Killip: This makes my head spin. The minimal object approach, a la Kenig-Merle, is not amenable to this relaxation. It might be possible to approach this with the (more quantitative) gopher strategy of CKSTT.

Wilhelm Schlag: Global dynamics above the ground state energy ↩

(joint work with Kenji Nakanishi NLW, NLS)

Klein-Gordon and Schrödinger Equations ↩

$$u_{tt} - \Delta u + u = u^3, R^{1+3}$$

$$ i \partial_t \psi + \Delta \psi + |\psi|^2 \psi = 0, R^{1+3}$$

LWP in $H^1$. $T_* (| u(0)|_{\cal{H}}) >0 $ where $\cal{H} = H^1 + L^2$.

$E(u) = \int \frac{1}{2}(|\nabla u|^2 + |u_t|^2 + u^2) - \frac{1}{4}|u|^4 dx.$

If $E<0$ then you have finite time blowup.

Scattering set: $S_+ = [(u_0, u_1) \in \cal{H}: T_* = \infty, | u |_{ST} < \infty]$

$S_+$ is open, path connected, $S_+$ contains a small ball $B_\delta (0)$.

Questions and Answers ↩

Questions

$S_+ $ bounded in $ \cal{H} $.
$\partial S_+$: Is this smooth or very rough?
What is the dynamics of solutions on the boundary?
Does $\partial S_+$ separate regions of global existence versus finite time blouwp?

Answers

Recall $\exists ~ Q >0$ satisfying $-\Delta Q + Q = Q^3$.

Statements

Theorem (Nakanishi-Schlag 2010): (Radial Case for now)

$S_+$ is unbounded.
$\partial S_+ \cap [(u_0, u_1) \in \cal{H}: E(u_0, u_1) < E(Q, 0) + \epsilon^2]

Trichotomy: If you are slightly above Q, you either

Scatter to $Q$.
Scatter to 0.
blowup.

Computer Simulations ↩

(done with R. Donninger)

These were beautiful and provoke lots of ideas and wonder.

Structures in Phase Space ↩

$S_+ \cap surface$. Take $(Q+ Af, Bg), (Af, Bg)$, here with $f,g$ are radial fixed functions. Here $A,B$ are parameters and we draw a rectangle in (A,B) space and we color based on (numerical) GWP vs. Blowup.

$K(u) = \int | \nabla u|^2 + u^2 - u^4$

$ PS_{\pm} = [ (u_0, u_1) \in {\cal{H}}: E(u_0, u_1) < E (Q,0), K(u) \geq 0 (for +)]$ and $K(u) < 0$ for -.

PS denotes the Payne-Sattinger (1978) sets. What is up with these sets?

$-\Delta \phi + \phi = \phi^3, ~ J’ (\phi) = 0$ where $J(\phi) = \int \frac{1}{2} (|\nabla \phi|^2 + \phi^2) - \frac{1}{4} \phi^4 dx.$

$ K(\phi) = \langle J’(\phi) | \phi \rangle = 0.$

$\partial_{\lambda}|_{\lambda = 0} J(e^\lambda \phi) = K(\phi).$

Find the minimal height of the potential well. You do some mountain pass work.

$\inf [J(\phi): \phi \in H^1, \phi \neq 0, K(\phi) = 0)] = J(Q) = \inf [ J(\phi) - \frac{1}{4}K(\phi): \phi \in H^1, \phi \neq 0, K(\phi) \leq 0]$

Cor: $PS_{\pm}$ are invariant under the flow.

$PS_{+} \implies$ global existence.
$PS_{-} \implies$ finite time blowup.

$K(\phi) \geq 0 \implies K(\phi) \gtrsim \min (1, \| \phi \|_{H^1}^2 )$

Cor: $Q$ is unstable.

….as usual, Wilhelm is fast….deductions are rapid fire.

Ibrahim-Nasmoudi-Nakanishi proved that you not only have global existence in $PS_{+}$, but using the Bahouri-Gerard-Kenig-Merle compensated compactness machinery, you actually have scattering.

Final State Descriptions near $Q$ ↩

Theorem (Nakanishi-Schlag): ${\cal{Hrad}}^\epsilon = [ (u0, u1) \in {\cal{Hrad}} : {\cal{E}} (u0, u1) < J(Q) + \epsilon^2 ].$ Then, this set is a disjoint union of 9 nonempty sets. $\| (u, u_t) - (\pm Q, 0) \|_{\cal{H}} < C\epsilon.$

-: Scatter, Trapped by $\pm Q$, Finite time blowup
+: Scatter, Trapped by $\pm Q$, Finite time blowup

(Choptuik and Bizon have explored similar pictures in studying the GR setting.)

ack….too fast for me to type….grazing solutions…penetrating solutions…..exit mechanism….and now he is speeding up…..mind like a ferrari….beautiful phase space portraits

Jeremy Marzuola: Scattering and soliton stability in ${\dot{H}}^{-1/6}$ for quartic KdV ↩

(joint work with H. Koch)

The goal is to outline the ideas in this work.

The problem ↩

$$ \partial_t \psi + \partial_x (\partial_x^2 \psi + \psi^4) = 0 $$ with initial data $\psi_0$.

Quartic KdV is the first integer power gKdV that is not completely integrable. Also, we use multilinear estimates.

small data case: $| \psi0 |_{ {\dot{H}}^{-1/6}} \ll 1$.

$\psi0 = Q_{c} (x-x_0) + v_0, | v0 |_{ {\dot{H}}^{-1/6}}.$

Questions:

Scattering and GWP for small data (Yes)
Scattering and Asymptotic stability (Yes)
Existence of inverse wave operators (Almost)

Previous Results ↩

Pego-Weinstein 1994, Asymptotic stability with exponential weights.
Martel-Merle 2001-…, Asymptotic stability in energy space $H^1$ in a moving reference frame.
- Virial Identities
- Monotonicity properties
Côte 2006, Constructs multiple soliton solutions for gKdV.
Grünrock 2005, Multilinear estimates.
Tao 2006, Asymptotic stability in $H^1 \cap {\dot{H}}^{-1/6}$.

Function Spaces ↩

I don’t want to construct spaces in as much detail as done in the paper here. The convergence in the wave operators takes place in a Besov refinement of ${\dot{H}}^{-1/6}$.

$(U^p, V^p)$

These spaces are nicely presented in a paper by Hadac-Herr-Koch 2009. Tataru, Koch-Tataru.

Steps of Proof ↩

Improved linear estimates, there are many linear equations meriting detailed study.
- Airy $(\partial_t + \partial_x^3) \psi = 0.$
- The $u$ problem: $(\partial_t u + \partial_x ({\cal{L}} u))=0.$
- The $v$ problem: $(\partial_t v + ({\cal{L}} \partial_x v))=0.$
- Refined Kato smoothing estimates for Airy
- ${\cal{L}} = (-\partial_x^2) + c - p Q_c^{p-1})$ Refined (weighted) elliptic estimates for $\cal{L}$
- Virial identities (Martel-Merle) for the $v$ problem $\implies$ energy spaces for the linear evolution.

First Result:

$P_{Q’}^{\perp} \psi = \psi - \frac{\langle \psi, Q’ \rangle}{\langle Q’, Q’ \rangle} Q’$
${\tilde{P}}_{Q’}^{\perp} \psi = \psi - \frac{\langle \psi, Q \rangle}{\langle Q, {\tilde{Q}} \rangle} {\tilde{Q}}$ where ${\tilde{Q}} = x \cdot Q’ + \frac{2}{3} Q.$
$\cal{L} (\partial_x Q) = 0$
$\partial_x (\cal{L} Q’) = 0$
$\partial(\cal{L} \tilde{Q}) = Q’$
Variable coefficient operators (small modulations)
$U, V$ spaces/Littlewood-Paley.
Multilinear Estimtes
- Rely heavily upon the $L^6$ estimate: $\| u \|{L^6{t,x}} \leq \| |D|^{-1/6} u \|_{L^2}.$
- Bilinear Estimate…long expression hard to read….
- Example: $$\| \partial (v1 v2 v3 v4) \| ({ {\dot{Y}}^{-1/6}{\infty, T}}) \leq c \prod{j=1}^4 \| vj \| ({ {\dot{X}}^{-1/6}_{\infty, T}}).$$
Full nonlinear problem requires delicate modulation. If you do so, you can’t close the multilinear estimates. Instead, we only require orthogonality asymptotically, rather than at all times.
More multilinear estimates involving $Q, {\tilde{Q}}, Q’$.
GWP for small data/scattering in scaling spaces
Inverse wave operators.

Energy spaces ↩

Virial identity for the $v$-problem: $\eta (x) = -\frac{5}{3} \frac{Q’}{Q}$. Claim: $$ - \frac{d}{dt} \int \eta(x) v^2 dx + c \| [sech]^2 (\frac{3}{2} x) v \|_{H^1}^2 \leq 0.$$

So, we have some monotone decrease in this weighted space.

$$\partial_t \langle v, Q’ \rangle = \langle {\cal{L}} (\partial_x v), Q’ \rangle $$.

Kato Smoothing:

$\gamma_0 (x) = 1 + \int_{-\infty}^x (1 + |y|^2)^{-(1+\epsilon)/2} dy.$
$\gamma_\mu = \gamma_0 (\mu^{-1} (x - \mu^{-2} t))$

$ \frac{d}{dt} \int \gamma_\mu u^2 dx + \int (\gamma_\mu)’ (u_x^2 + \frac{1}{3 \mu^2} u^2) dx \leq 0.$

$ \partial_t \langle {\cal{L}}^{-1}v, v\rangle = 0.$

Lemma:

$$ E(v) = \int \gamma (x) (v_x^2 + v^2) dx + \lambda_E\int \eta(x) v^2 dx + \Lambda_E \langle {\cal{L}}^{-1} v , v\rangle

We define then our “natural” Energy spaces.

$X^s = L^\infty H^s \cap L^2 H^{s+1}_{\sqrt{\gamma’}}$
$Y^1 = L^1 H^1 + L^2 \sqrt{\gamma’} L^2$

We then build spacetime function spaces using the $U^2, V^2$ spaces (defined in S. Herr’s talk) based on these structures and the cubic dispersion relation….and not the linearized equation for the $v$ equation…..chalk coming too fast for me to write down…..ack.

Nonlinear Modulation ↩

$\psi (x,t) = Q_{c(t)} (x - y(t)) + w(x,t)$

$ \partial_t w + \partial_x (\partial_x^2 w + 4 Q^3 w) = \frac{\dot{c}}{c} {\tilde{Q}}(x-y) + ({\dot{y}} - c) (Q_c)’ (x-y) - \partial_x ( O(w^2)). $

Usually, we choose w $\perp Q, Q’$ through choice of $c, y$.

$ \frac{ {\dot{c}}}{c} \langle (Q_c) , (\tilde{Q}c ) \rangle = \langle w, (Qc ) \rangle.$

$ (\dot{y} - c^2) \langle (Q_c)’, (Q_c)’ \rangle = - \kappa < w, (Q_c)’>$

We then calculate: $$ \frac{d}{dt} \langle w, Q \rangle + \langle w, Q \rangle = O (w^2), $$

$$ \frac{d}{dt} \langle w, Q’ \rangle + \kappa \langle w, Q’ \rangle + O (w^2) = 0. $$

With this structure and the formalism of Tao, and some careful work, we can put it all together.

Postlude ↩

I had a nice follow-up conversation with Raphaël Côte. I wondered whether there were similar small data and remainder-atop-soliton scattering results for low power KdV equations. He pointed out that “clean” scattering does not hold in the small data case for the low power gKdV equations. Instead, there are modified scattering statements for data satisfying certain weighted conditions proved by Hayashi and Naumkin It is perhaps reasonable to expect corresponding statements about the error term in the asymptotic stability results around (multi)solitons. However, this is open for study.

Sijue Wu: Global and almost global wellposedness of the two and three dimensional full water wave equations ↩

arXiv: GWP of full 3D water wave problem

Introduction ↩

We are looking at the middle of the ocean. Let’s imagine infinite depth and no boundary. We have gravity pointing odwn and the density of the air is 0 and the density of the water is 1. We assume the water is inviscid, incompressible, irroational, surface tension is zero. The interface is called $\Sigma (t)$.

The motion of the fluid is described by the Euler equation ${v_t} - v \cdot \nabla v = (0, -1) - \nabla P$ in the interior $\Omega(t)$. We also have $div v = 0, curl v =0$, ….ack slide changed.

G.I. Taylor (1949) linearized about the flat interface and found that air above water is stable but water above air is unstable.

LWP ↩

LWP for arbitrary data [S. Wu 1997 (2d) 1999 (3d)]: Local existence in Sobolev spaces under the right Taylor stability condition.

Earlier Results:

Beal, How, Lowegrub 1992 formulated the Taylor sign condition: $ -\frac{\partial P}{\partial n} \geq c_0 > 0#.
Nalimov 1974 infinte depth
Yoshihara 1982.

The work has been extended in many directions. Iguch 2001, Ogawa and Tani 2002, Ambrose and Masmoudi 2005, Lannees 2005, Christodoulu and LIndblad 2003, Lindblad 2005, Coutand and Skholler 2005, Zhang and Zhang, Shatah and Zhang.

Global-in-time behavior ↩

What is the global in time behavior of the solution of the water wave equation?

We will focus on small and smooth data. This is reasonable since it is known that 90% of the waves on the ocean are smaller than 2m? I’d like to know the reference for this 90% claim. Maybe this is done using satellite data? Perhaps this remark motivates a probabilistic Cauchy theory which explains the infrequency of rogue waves?

….slides are changing fast….I can’t keep up so I will listen and make remarks wehn I can.

Quadratic interaction is too strong so the key idea is to use a change of variable which recasts the problem with a cubic nonlinearity.

A natural setting for studying 3D water wave is the Clifford Algebra and use Clifford analysis. The difficulties in 3D are that there is no Riemann mapping, the Clifford Algebra is noncommutative, products of analytic functions in 3D are not analytic. We find that in the 3D problem there is also a special structure allowing us to recast the problem so that quadratic problems disappear and the nonlinearity is cubic and higher orders in nature. It is not purely cubic, there are some quadratic terms but we can handle those as though they are cubic.

Statements ↩

Theorem: (2D) Assume initial wave is of small height, initial velocity is also $\epsilon$ small. Assume we have finitely many derivatives of f and g are in $L^2$. Then, there is a unique solution on a time interval $[0,e^{c/\epsilon}]$. During this time, the solution remains smooth and small.

Theorem: (3D) We assume less here. Suppose initial condition given as a graph. For data with small steepness (no smallness condition on the height) and possibly with infinite energy but also with small velocity on the interface, then the solution is uniquely defined and global-in-time, remains smooth and small.

It seems like we have a better result in 3D. But, in my opinion, these two results are equivalent, they are of equal strength: equally good/equally bad. We can view the 2D case inside the 3D problem and in that view we have an infinite energy 3D case. Maybe we can prove the 2D result under the small steepness condition.

Famous picture of Rogue wave with a ship in foreground.

Rogue waves are vastly massive waves (30m). Often appear in perfectly clear weather, wtithout warning. It’s exact causes are still unknown. Possible causes? Diffractive focusing (effect from caostline)? Focusing of currents? Nonlinear effects? We are avoiding wind and boundaries so we want to understand whether nonlinear effects can be explained as the source of rogue waves.

I am confused. The 3D result says that initial waves given as a graph over the bottom with small steepness remain small and smooth forever. So, this result does not explain or speak to the rogue wave phenomenon. Of course, it suggests that large initial steepness is required for a rogue wave to form within this model of the ocean. Again, this situation seems ripe to me for a probabilistic study of the Cauchy problem?

“Once you get the algebra part right, the analysis part just goes through without complication.”

We only need to know the fluid motion on the fluid interface. We therefore try to reduce the Euler equation to an equation on the fluid interface. This removes the difficulty of the free boundary.

Normal Forms Discussion ↩

The technical discussion seems to revolve around making a bilinear change of dependent variable with the goal of killing off the cubic terms. It doesn’t work….but when working in the right coordinate system with the right quantities, the nonlinearity of the 2D water wave equation is cubic and higher orders.

Nickolay Tzvetkov: On random data nonlinear wave equations ↩

Background References

(joint work with Nicolas Burq)

Framework ↩

Let $(M,g)$ be a Riemannian manifold of dimension $d=3$ with $\partial M = \phi$. We consider the cubic wave equation $$ (\partial_t^2 - \Delta_g) u + u^3 = 0 $$ with initial data $(u0, u1) \in H^s \times H^{s-1}$.

$H^{1/2}(M)$ is the critical space for this problem. He sometimes denotes the problem with (*).

Theorem (deterministic theory):

The problem (*) is locally well-posed in $H^s \times H^{s-1}, ~ s \geq 1/2$ and globally for $s \geq 1$.
The problem (*) is ill-posed in $H^s \times H^{s-1}, ~ s \in (0, 1/2).$
1. For example, $\exists ~ (u_n (t))$ sequence of smooth solutions of (*) such that the initial data goes to zero in $H^s \times H^{s-1}$. But, $ \| (u_n(t), \partial_t u_n (t)) \|{L^\inftyT_ ; H^s \times H^{s-1}} = + \infty, ~\forall T>0.$ (inspired by work of Christ-Colliander-Tao)
2. Moreover, $\exists$ a single data $(u0, u1) \in H^s \times H^{s-1}$ such that $\forall ~T>0$, (*) has no solution in $L^\infty ([0,T]; H^s \times H^{s-1})$ satisfying the finite propagation speed. (instantaneous blowup inspired by work of Lebeau)

On $R^3$, there are refined global results for $s \in [3/4, 1]$ are due to Kenig-Ponce-Vega, Gallagher-Planchon, Bahouri-Chemin, Roy, …. Probably this can be transported to the torus (using finite propagation speed) but this is not written. OPEN QUESTION

Question: Can one still prove some form of well-posedness for $s< \frac{1}{2}$?

Idea: Yes, by randomizing the data.

We have a general method to do this locally in time Burq-Tzvetkov 2008.
A very particular method for globally in time [Burq-Tzvetkov 2008]((http://arxiv.org/abs/0707.1448 “Random data Cauchy theory for supercritical wave equations II : A global existence result”)), exploiting invariant measures a la Bourgain.

Goal for today: General method for globally in time. We can skip this invariant measure business. But, if we are only PDE people, there is a method which allows us to globalize without relying upon the invariant measure aspects.

Randomized data on $T^3$ ↩

Starting from $(u0, u1) \in H^s \times H^{s-1}$ we form their Fourier series $$ u0 = \sum_{n \in Z^3} c^0_n e^{i n \cdot x}$$ (same for u1) and we define $$ u_0^\omega = \sum_n g_n^0 (\omega) c_n^0 e^{i n \cdot x} $$ with natural hypotheses on the random variables to ensure the data stays real valued. He also decomposes the $g(\omega)$ in real and imaginary parts is a system of i.i.d. random variables with a joint distribution $\mu$ satisfying $\exists ~c>0, ~ \forall ~ \gamma >0, \int_{-\infty}^{\infty} e^{\gamma x} d\mu (x) \leq e^{c \gamma^2}$.

Examples:

Gaussians: $d \mu (x) = e^{-x^2/2} \frac{dx}{2\pi}$
Bernoulli: $d \mu (x) = \frac{1}{2}( \delta_{-1} + \delta_1)$

The gaussians generate a dense set in $H^s$. Bernoulli does not but leaves the data on the same sphere in $H^s$.

Theorem: Let $M = T^3, (u0, u1) \in H^s \times H^{s-1}, ~ s \in [0,1]$. Then (*), with data $(u_0^\omega, u_1^\omega)$ is globally well-posed almost surely in $\omega$.

Consider the probability measure $\rho$ on $H^s \times H^{s-1}$ defined by the map: $ \omega \rightarrow $(u0^\omega, u1^\omega)$. Every function gives a different measure, so I have many measures.

Theorem (again): There exists a set $\Sigma$ such that $\rho (\Sigma) =1$ and such that $\forall ~ (v_0, v_1) \in \Sigma$ there is a unique global solution of (*) with data $(v0, v1)$ such that $$ (u, u_t) \in [Free ~Evolution ~of~ (v_0, v_1)] + C(R; H^1 \times L^2). $$ In addition, the solution satisfies the finite propagation speed and, moreover, if we denote by $\Phi(t)$ the constructed flow we have the following properties:

$\Phi (t) (\Sigma) = \Sigma$
$\forall (v_0, v_1) \in \Sigma$, $\| \Phi(t)(v_0, v-1)\|_{H^s \times H^{s-1}} \lesssim \langle t \rangle^{1-s/s +}, s>0.$ (Remark: The implicit constant here is a random variable.)
Measure same thing in $L^2 \times H^{-1}$ and we get the bound $e^{c t^2}$.

Steps in the proof ↩

Global existence step. (inspired by Paley and Zygmund)
Construction of the set $\Sigma$. (inspired by the invariant measure consideration by Bourgain)
Control on the flow for $s>0$. (inspired by the high/low frequency decompositon a la Gallagher-Planchon and by recent work by Colliander-Oh)
Control on the flow for $s=0$. Here the analysis degenerates. (inspired by the work of Yudovich on the Euler equations) “We can say that we have developed a probabilistic version of the Yudovich argument.”

On the proof of the Global existence step for $s>0$ ↩

Large deviation estimates. Consider $\square_g u_{lin}^\omega$ with the randomized data $(u0^\omega, u1^\omega)$. For $s>0, ~\delta > 0, ~\exists c>0, ~ \forall \lambda \geq 0$ we have the large deviation estimate $$ p ( \omega: \| \langle t \rangle^{-\delta} u_{lin}^\omega \|_{L^\infty (R \times T^3)} > lambda ) \leq \frac{1}{c}e^{-c \lambda^2}. $$ Of course, this is much better than what we can get from Strichartz.

We look for solutions as $u = u_{lin}^\omega + v$ and we study $\square_g v + (v + u_{lin}^\omega)^3 = 0$ with zero initial data. We have the energy $E(v) = \frac{1}{2} \int |\nabla v|^2 + |v_t|^2 + \frac{1}{4}\int v^4 dx. We then calculate $\frac{d}{dt} E(v) = \int \partial_t v (v^3 - (v + u_{lin}^\omega)^3).$ We are lucky because the $v^3$ terms cancel and by Gronwall we have global existence for $\omega$’s of big probability. Then, we make some intersections and do some measure theory to finish.

This argument gives exponential control. We revisit the analysis using the high/low frequency truncation ideas to improve to polynomial control.

Remark: We can prove similar results for ANY manifold by using a randomization due to Lebeau.

Questions ↩

Schlein: How is the set $\Sigma$ invariant?

Tzvetkov: The set \Sigma is of the form random orbit of the data plus smooth functions. Since the smooth functions have zero measure, we can throw them into \Sigma.

Ionescu: How do you see in the analysis that you are studying the defocusing question?

Tzvetkov: In the Gronwall business, we used the sign.

Pierre Germain: Global existence for coupled Klein-Gordon equations with different speeds ↩

Background References:

General Problem: Understand global existence and scattering for nonlinear dispersive equations with very nice data. ↩

We will assume the Cauchy data are small, smooth and localized. We will further restrict the problem to semilinear wave and Klein-Gordon equations in dimension 3.

NLW, $d=3$ ↩

$\square u = |u|^{p-1}u$.
- Above the Strauss exponent $ p > 1 + \sqrt{2}$.
- At the Strauss exponenent, finite time blowup was shown by [John-Schaeffer]
$\square u = |u_t^2 - |\nabla u|^2$.
- Null form structure observed by Christodoulu and Klianerman gives global existence.
$ \square u = |u_t|^2.$
- finite time blowup [John]
$\partial_t^2 u^i - c_i \Delta u^i = \sum Q^i_{jk} (Du^j, Du^k)$
- Global existnce if $Q^i_{jk}$ is a null form. [Yokoyama, Ohta, Katayama, Sogge, Metcalfe, ….]

NLKG ↩

$\partial_t^2 u - Delta u + u = |u|^{p-1}u.
- For $p>2$ (the Strauss exponenet), you have global existence [Strauss].
$\partial_t^2 u - Delta u + u = Q(u,u)$ or $Q(Du, Du)$.
- global existence [Klainerman], [Shatah]
What about different propagation speeds? $\partial_t^2 u^i - c_i \Delta u^i + u^i = \sum Q^i_{jk} (u^j, u^k)$
- This case has some difficulties and my new result addresses this issue.

All the results I quoted have been provd using the vector field method. How does this work? You find a bunch of vector field $(\Gamma_i)$ which commute with the linear part of the equation. Then you estimate $\Gamma^\alpha u.$ The method does not apply to KG with different speeds. You don’t have sufficiently many commuting vector fields to treat the multiple speed KG case.

There were some other methods used for these problems. In particular, Shatah used a normal forms method. Christodoulu used a change of variables method but most of the theory has been built on the vector field method.

NLKG with different speeds is a toy model for Euler-Maxwell, provided you restrict to high frequencies and ignore certain things.

Statement ↩

Theorem: $$ \partial_t^2 u - \Delta u + u = Q(u,v), ~ \partial_t^2 v - c^2 \Delta v + u = P(u,v) $$ with some initial data for the two equations. (No derivatives in the quadratic nonlinearities.) Assume that the data has some $L^2 $ weighted (power 1 ) control and is small enough and we also have $H^N$ smallness with a big enough N. Then there eists a global solution which furthermore scatters in $H^N \times H^{N-1}$ which means the nonlinear evolution converges to a linear solution as time goes to infinity.

The vector field method does not apply. Instead, we use a spacetime resonances method which we have applied to the water wave problem and to the NLS equation. This is a new instance where we can apply this method. The method was developed in collaboration with Shatah and Masmoudi.

Spacetime resonance method ↩

For the sake of exposition, consider $i \partial_t u + P(D) u =u^2$ emerging from data $u0$. Let $f(t) = e^{-it P(D)} u(t)$ and consider this new unknown function instead of $u$. Write the Duhamel formula for $\hat{f}$. What you find is that $$ \hat{f} (t, \xi) = \hat{u_0} (\xi) + \int_0^t \int e^{i s [P(\xi + \eta) - P(\xi) - P(\eta)]} \hat{f} (\eta, s) \hat{f} (\xi - \eta, s) d\eta ds. $$ We have a problem if the phase is stationary either in s. What can save us is the oscillations. This is what we call time resonances. Or, if the phase is stationary in $\eta$ and this is what we call space resonances. Of course, the worst situation is when we have stationarity in both senses and this is what we call spacetime resonances.

Method

If the phase factor (redenoted as) $\phi \neq 0$ an integration by parts in $s$ and push the nonlinearity to cubic. This is just the normal forms method seen on the Fourier side.
If $\partial_\eta \phi \neq 0$ you can integrate by parts in $\eta$ and you gain an $s$ in the denominator which is “always pleasant when you are trying to prove global eistence.” This is the vector field method seen in Fourier space.

He draws two graphs where $\phi$ and where $\partial_\eta \phi$ vanish on the $\xi, \eta$ Cartesian product. We use pseudo-product operators Coifman-Meyer to decompose in the $(\xi, \eta)$ space. $$ \mathcal{F} ( B_{m(\eta, \xi)}) (f,g) (\xi)= \int m(\eta, \xi) \hat{f} (\eta) \hat{g}(\xi - \eta) d \eta. $$

Physical meaning

Time resonances are “standard resonances” in the dynamical systems sense.
Space resonances are when waves of different frequency move with the same group velocity (….not really explained)

Application to our problem ↩

You get a lot of different phase functions:

$\phi (\xi, \eta) = \langle \xi \rangle_l \pm \langle \eta \rangle_m \pm \langle \xi - \eta \rangle_n$ where $\langle x \rangle_\alpha = \sqrt{1 + \alpha^2 x^2}$ and $l,m,n$ are chosen among the two possibilities: 1 and $c$.
Look at the place where both $\phi$ and where $\partial_\eta \phi$ vanish.
- Sometimes this set is empty.
- Sometimes this set has the form $[ |\xi | = R, \eta = \lambda \xi]$ for real numbers $R, \lambda$.
- Actually, such a set is generic for interactions between waves with a dispersion relation $p(|\xi|)$ which depends only on the frequency size. Thus, the method can be applied to other settings.

He redraws the graph of the zero level sets for $\phi$ and $\partial_\eta \phi$. He then excises around the point where these sets intersect using a cutoff using a pseudo-product operator with a symbol $m$ which is increasingly singular along the set of simultaneous vanishing. This is a bit annoying because there are no general estimates for such pseudo-products. We would need to estimate the boundedness of $B_m: L^p \times L^p \rightarrow L^r$ where $m$ is singular along $[ |\xi | = R, \xi = \lambda \eta]$ for real parameters $R, \lambda$. The Coifmann-Meyer calculus requires nicer properties on $m$. In contrast, there is work by Lacey-Thiele on the bilinear hilbert transform which does have a singularity in the 2-multiplier but does not apply to our case.

We use that you are at the Strauss exponent so that rough estimates are enough to succeed.

In the theorem, we need to assume that resonances are separated. Look at the spacetime resonance set $\cal{R} = [\phi = 0] \cap [\partial_\eta \phi = 0]$ and project onto $\xi$, which I call “outcome frequencies”. If you project onto $\eta, \xi - \eta$ you get what I call “source frequencies”. We need to assume that $[outcome] \cap [source] = \phi$. This is generically true for different speeds $c$. In particular, we have this property true for all but a discrete set of speeds $c$.

There is alast point wihich is a bit problematic: Spacetime resonaces at $\infty$. $\phi, \partial_\eta \phi \rightarrow 0$ at $\infty$. To overcome this difficulty, we rely upon the high regularity $H^N$ hypothesis using Strichartz estimates. We then separate the analysis into low and high frequencies.

Questions ↩

Koch: Gain from modulation versus gain from bilinear estimate. Dualize the argument and you can recast as a condition on the nonvanishing of $\partial_\xi \phi$.

Postlude: ↩

After the talk, I learned from Pierre that he had written an expository article on the spacetime resonances method. T

Oana Ivanovici: Dispersive Estimates on convex domains ↩

(joint work with Fabrice Planchon)

Introduction ↩

Consider a domain $\Omega$ of dimension $d \geq 2$. We consider the wave equation $\partial_t^2 u - \Delta u = 0$ with initial data and vanishing on $\partial \Omega$.

Consider, for point of reference versus later statements, the situation where $\Omega = R^d$. Take $u_0 = \delta_a$ and $u_1 = 0$. Then the solution is given by the Green’s function $$ u_{a, R^d} (t,x) = \int \cos (t |\xi|) e^{i \xi \cdot (x-a)} d \xi. $$

Dispersive Estimates: $$ | \psi (h D_t) u_{a, R^d} |_{L^\infty} \leq C(d) h^{-d} \min (1, (\frac{h}{|t|})^{\frac{d-1}{2}}). $$

We are interested in the case where $\partial \Omega \neq 0$. We must confront reflected waves, glancing rays and waves which travel along the boundary.

Let $\Omega$ be a strictly convex domain. In particular, we will consider $\Omega$ to be the Friedlander domain. $\Omega = [(x,y) \in R^d: x>0, y \in R^{d-1}]$ with the associated Laplacian $\Delta = \partial_x^2 + (1+x) \Delta_y$. This is very close to the laplacian on the disk. Then, she draws the half space and describes the bicharacteristics as a bunch of circles bouncing along the floor.

Theorem: Take $a>0$ small so that $(a,0) \in \Omega$ (in the interior but close to $\partial \Omega$). $\exists ~T>0$ such that $\exists ~ C>0$ such that $\forall ~ h \in (0,1]$ we have $$ | \psi (h D_t) u (t, x, y) |_{L^\infty (\Omega)} \leq C(d) h^{-d} \min (1, (\frac{h}{|t|})^{\frac{d-2}{2} + \frac{1}{4}}). $$ The way to study this is to consider the set of points you can reach from the point $a$ upon traveling for time $T$. The method of proof involves a decomposition of the data in terms of wave packets which hit the boundary a certain number of times. The worst packets are localized in small cones that are almost parallel to the boundary.

…rapid discussion of some frequency localzations…lots of glancing rays pictures….subsequent reflections are denoted by $u_j$. Each reflection involves a loss of 1/6 derivative and there can be manyreflections ccumulating until a total loss of 1/4 derivative. After that there will be no more regularity loses.

Applications ↩

$\implies$ Spectral projector and Strichartz estimates. Smith and Sogge studied similar problems using a reflection across the boundary idea. For dimensions $d \geq 3, the reflection method does not have a chance to get optimal regularity losses. First, you don’t see the dispersion tangential to the boundary. Also, their study only captures the loss from one reflection but does not resolve the accumulated losses. The loss of 1/4 derivative happens at a special time after many reflections.

Works by Blair-Smith-Sogge are improved in this work. She draws some Strichartz diagrams and shows that her new dispersive estimate implies a wider range of valid Strichartz exponents.

We will soon see that the only possible losses are 1/6 or 1/4.

Cusp solutions hugging the boundary ↩

This result was announced at a conference in Evian by G. Lebeau. Lebeau explained the geometrical features of the argument but the analytical details were not written down. Fabrice and I are writing those down….

To demonstrate the loss, she writes the boundary and draws data that looks like a cusp. $$ u_0 (x,y) = \int e^{i \frac{\eta}{h} 9\frac{\xi^3}{3} + (x-a)\cdot \xi + y} \psi (\xi) \phi (\eta) d\xi d\eta. $$ The wave starts localized within $a$ of the boundary. AFter some time $t \sim 2 \sqrt{a}$ the cusp is upside down wrt boundary and then at the time $t = 4 \sqrt{a}$ the cusp reappears and the singularities only appear at these specific locations near the boundary. The situations is studied with $a \thicksim h^{1/2}$. For $a$ smaller than this power, we would not be able to repeat the construction for many reflections. It will degenerate. The caustic in this case is the line sliding along the boundary passing through the cusps. along the caustic, the intensity of light is much brighter. At points along the caustic, oscillatory integrals don’t enjoy good bounds.

Proof ↩

$$ u_h (z) = \frac{1}{h^{1/2}}\int e^{\frac{i}{h} \phi (z, \xi)} \sigma (z, \xi, h) d\xi, \xi \in R $$ Everyone knows that the number and degeneracy of the critical points of the phase function control the asymptotics of this guy as $h \rightarrow 0$.

Degenerate critical points: Let $k$ denote the *order of the caustic of u_h$ be defined by $\inf_{k’} [k’: | u_h | \sim O (h^{-k’})]$

Example 1: Let $\phi_F (z, \xi) = \frac{\xi^3}{3} + z_1 \xi + z_2.$ Here $z_1 = -\xi^2, ~z_2 = - \frac{2}{3} \xi^3$. So we have a fold. This type of phase function corresponds to $k = 1/6$. She draws a sideways parabola and projects it down onto a cusp.

Cusp type integral: $\phi_C (z,\xi) = \frac{\xi^4}{4} + z_1 \frac{\xi^2}{2} + z_2 \xi + z_3.$ (This has order 1/4) (Pearcy-type integral)

$\partial \phi:~ z_2 + 2 z_1 \xi + \xi^3 = 0
$\partial^2 \phi: ~ 2 z_1 + 3 \xi^2 = 0
$\partial_\eta (\eta \phi_c ): ~ z_3 _ \xi z_1 + z_2 \frac{\xi^2}{2} + \frac{\xi^4}{4}=0.

* Swallowtail:* $\phi_s (z, \xi) = \frac{\xi^5}{5} + z_1 \frac{\xi^3}{3} + z_2 \frac{\xi^2}{2} + z_3 \xi + z_4.$

We have a degenerate critical point of order 4…..ack….I am running out of battery and this is really nice stuff…

Axel Grünrock: Cauchy Problem for higher order KdV and mKdV equations ↩

I am interested in the question of optimal local well-posedness.

Background References

Grünrock: Paper described in this talk

Equations ↩

KdV hierarchy

Lax 1968 introduced the hierarchy of higher order KdV equations. $$ \partial_t u + \partial_x G_j (u) = 0 $$ We will refer to this as (hoKdV-j), the higher order KdV equation. $$ \langle Gj (u), v \rangle = \frac{d}{d\epsilon} Hj (u+\epsilon v)|_{\epsilon = 0} $$

where

$$ Hj (u) = \int P_j (u, \partial_x u, \dots, \partial_x^j u) dx. $$ These are the Hamiltonians of KdV.

$P_{-1} (u) = u$
$P_0 = - \frac{1}{2} u^2$
$P_1 (u) = - \frac{1}{2} u_x^2 - u^3$

The iteration procedure then defines the hierarchy:

$G_1 (u) = u_{xx} - 3 u^2 \implies u_t + \partial_x^3 u = 6 u u_x$
$ u_t + \partial_x^5 + 5 \partial_x ( \partial_x^2 u^2 - (\partial_x u)^2 - 3u^3)=0$
$ u_t + \partial_x^7 u - 7 \partial_x (\partial_x^4 u^2 - 2 \partial_x^2 (\partial_xu )^2 _ (\partial_x^2 u)^2 - 10 u \partial_x (u \partial_x u + 5 u^4) = 0.$
….

We can thus define some general structure of the higher order KdV equations based on rank properties where $rank_{KdV} = degree + \frac{1}{2}~ derivatives ~in~ x = j+2$. We find that $|\rho| = 2 (j-k) + 3.$

For all the equations in the hierarchy, we have the same scaling critical regularity of $s_c = - \frac{3}{2}$.

There is a second shared property for all the equations in the hierarchy. The Hamiltonians in the KdV hierarchy are all in involution with respect to the Poisson bracket: $$ { H_k, H_l } := \langle G_k (u), \partial_x G_l (u) \rangle, ~\forall k, l \geq -1.$$ We can therefore calculate that $$ \frac{d}{dt} H_k (u) = \langle G_k (u), \partial_t u \rangle = - \langle G_k (u), \partial_x G_l (u) \rangle = 0. $$

mKdV hierarchy

A similar tower or hierarchy of equations may be built around the mKdV equation using the Miura map: $v \rightarrow v_x + v^2$.

Sequence of ${\tilde{H}}j (v) = H(j-1) (v_x + v^2)$. This spawns ${\tilde{G_j}}(v)$ by writing $$ \partial_t v + \partial_x {\tilde{G}}_j (v) = 0 $$ which we denote by (homKdV-j). What can we say about the structure of the nonlinear terms in the mKdV hierarchy of equations.

The rank condition for KdV hierarchy is transferred via the Miura map into a rank condition for the mKdV hierarchy.

nonlinear terms in mKdV hierarchy are all odd in $v$, so no quadratic terms.
$|l| = 2 (j-k) + 1$
We thus find that the mKdV hierarchy enjoys a joint scaling invariance corresponding to $s_c = - \frac{1}{2}$.

Earlier Results ↩

(Incomplete)

1979 Saut: Existence of persistent solutions of hoKdV-j and homKdV-j in $H^j$ using the energy method which works equally well in the periodic or nonperiodic setting.
1993 Ponce: hoKdV-2, LWP in $H^s (R)$ provided that $s> \frac{7}{2}$ and, combining the LWP result with conservation laws, he obtained GWP for $s \geq 4$.
2008 Kwon: LWP for hoKdV-2 for $s> \frac{5}{2}$ and GWP for $s\geq 3$ using a refined Energy method developed by Koch-Tzvetkov for treating Benjamin-Ono.
1993/4 Kenig-Ponce-Vega: $\exists ~s_0 = s_0 (j)$ and $m - m(j)$ such that $\forall ~ s \geq s_0$, hoKdV-j is LWP in $H^s (R) \cap L^2 (|x|^m dx)$
- Corresponding results for homKdV-j. It was remarked there that the weights are not necessary for treating the cubic and higher power cases.
1995 Linares: homKdV-2 is GWP in $H^s (R)$ provided $s \geq 2$.
2008 Kwon: LWP improved down to $s\geq - 3/4$ and thus GWP in $H^1$.
2008 Pilod: Without the weights in the data spaces, one has ill-posedness in the hoKdV-j hierarchy, ~$\forall ~j \geq 2$. In particular, he showed that the flow map can not be $C^2, ~ \forall s\in R$. The argument involves an interaction between high and very low frequencies. Higher order Sobolev regularity is not beneficial at all.

Killip: Is there a contradiction here with the positive result of Kwong vs. Pilod?

Grünrock: Kwon uses energy methods so obtains continuous dependence, not $C^2$ dependence of the flow map.

New Results ↩

Data spaces: $\| f \|{\hat{H}s^r} =\| \langle \xi \rangle^s \hat{f} \|{L^{r’}\xi}, ~ \frac{1}{r} + \frac{1}{r’} = 1.$ Here $1 < r \leq 2. We have $H^{s,r} \subset {\hat{H}}_s^r.

Spacetime spaces: $ | u |{X{s,b}^{r,p}} = | \langle \xi \rangle^s \langle \tau - \phi (\xi) \rangle^b \hat u|{L^{r’}\xi (L^{p’}\tau)}.$ Here we have $\phi (\xi) \sim \xi^{2j + 1}$.

What are the crucial estimate we need that will lead to local well-posedness?

Ingredients (tools)

Smoothing estimates
- linear: $\| D_x^{\frac{2j-1}{3r} u \|{L^r{tx}} \lesssim \| u \|{X^r{0b}}$ if $b > \frac{1}{r}, ~ \frac{4}{3}< r \leq 2$ ~(fails for $r \leq \frac{4}{3}$.)
- triliner estimates with the same gain order (up to $\epsilon$).
- bilinear refinement: For $b > \frac{1}{p}, ~ 1 < r \leq r_{1,2} \leq p \leq 2, ~ \frac{1}{r} + \frac{1}{p} = \frac{1}{r_1} + \frac{1}{r_2},$ $$ \| M_{j,p} (u,v) \|{ {\hat{L}^rx {\hat{L}}^p_x}} \lesssim \| u \|{X{0b}^{r_1, p}} \| u \|{X{0b}^{r_2, p}$$ We have an increasing gain of regularity with these estimates of gain $D_x^{\frac{2j}{p}}$ in the parameter $\frac{1}{r}$ or $\frac{1}{p’}$, respectively.
Resonance relation $(k=2)$. $$ \sum_{i=0}^2 \langle \tau_i - \xi_k^{2j+1} \rangle \gtrsim |\xi \xi_1 \xi_2| \times ( \xi_1^{2(j-1) } + \xi_2^{2(j-1)})$$. We have a gain: $D_x^{\frac{2j+1}{p’} -}$ since $\langle \tau_0 - \xi_0^{2j+1} \rangle^{b - 1 - \epsilon = \frac{1}{p’}}$. This gain is decreasing in $\frac{1}{r}$ or $\frac{1}{p’}$, respectively.

Statements

homKdV-j: He expresses the LWP results in the $(\frac{1}{r}, s)$ as lines leaving the vertical $s$ axis and all passing through the point $(1,0)$.

The results on $H^s$-scale for $j \geq 3$ are new. We obtain GWP in $H^s$ for $s \geq [\frac{j+1}{2}]$ (integer part). Thus, the use of $\hat{H}^{s,r}$ spaces lead to new insights.

Moreover, the results converge toward a nice statement which identifies a common joint space $\hat{L}^1$ which contains finite measures and which contains $L^1$. Unfortunately, the result at that endpoint is not yet established.

For KdV, he draws a similar picture. The lines do not appear to converge. we are far away from finding a joint space.

Questions ↩

Tataru: C^2 vs. mereley continuous dependence properties?

Staffilani: Periodic case?

Grünrock: No, I don’t have results there.

Colliander: $NLS_3$ in $\hat{L}^1$?

Postlude ↩

For me, fantastically interesting conversations with Koch, Grünrock, Tataru and Vega.

OPEN: Is there a space of functions wherein each equation in the mKdV hierarchy is GWP?
OPEN: The space ${\hat{L}}^1$ appears to be a natural candidate given the visual description Axel gave of his results.
Corresponding questions about cubic NLS in one space dimension? L. Vega points out that ${\hat{L}}^1$ can not do the job because of nonuniqueness results for NLS evolution emerging from the Dirac mass.
NLS has galilean invariance; mKdV does not so perhaps there is some hope for mKdV in ${\hat{L}}^1$?
I will ask Boris Khesin about whether the integrable hierarchy of equations containing cubic NLS is exposed nicely somewhere. It might be interesting to try and carry out an analogous study of the NLS hierarcy.

Selberg: Global existence for the Maxwell-Dirac system in two space dimensions ↩

(joint work with Piero d’Ancona)

The Maxwell-Dirac system (MD): $$ (-i \alpha^\mu \partial_\mu + M \beta) \psi = A_\mu \alpha^\mu \psi$$ $$ \square A_\mu = -\alpha \langle \alpha_\mu \psi , \psi \rangle$$

$B = \nabla \times A, ~ E = \nabla A_0 - \partial_t A.$ We are interested in evolution starting from data $\psi_0, E_0, B_0$ satisfying the constraints $\nabla \cdot E_0 = |\psi_0|^2, ~ \nabla \cdot B_0 = 0.$

We are using the Lorenz gauge condition: $\partial_\mu A^\mu = 0$.

2d: $\alpha^0 = I, ~ \alpha^1 = \sigma^1, ~ \alpha^2 = \sigma^2, ~\beta = \sigma^3$ where the $\sigma$’s are the Pauli matrices and the $\alpha$’s are called the Dirac matrices.

He decomposes the electric field into divergence free and curl free parts. We can then write $E_0 = E_0^{df} + \Delta^{-1} \nabla (|\psi_0|^2)$. We are restricting the motion to take place in the $x^1, x^2$ plane so the magnetic field must be perpindicular to that plane. All fields are independent of $x^3$. $A = (A_1, A_2, 0)$ and $B = (0, 0, \partial_1 A_2 - \partial_2 A_1)$. Given the initial consitions on the $E, B$ fields and the Lorenz gauge condition, we can specify the initial data for the potential $A$.

Maxwell-Dirac and Dirac-Klein-Gordon ↩

DKG: $(-i \alpha^\mu \partial_\mu + M \beta) \psi = \phi \beta \psi$….ack too fast. Best reference for this is the paper of Glassey-Strauss 1979.

Energy - no sign
Charge: $\int |\psi (t,x)|^2 dx = const.$
Scale invariant regularity: $\psi_0 \in {\dot{H}}^{d-3/2}, ~ E_0, ~ B_0 \in {\dot{H}}^{d-2/2}$.
- MD is critical is charge critical in 3d
- charge subcritical in 2d and 1d.

Of course, we would like to prove global regularity. A natural strategy is to prove low regularity LWP and then exploit conservation laws. However, this is not so clear yet…..

Results ↩

Global:

1d MD GWP: 1973 Chadam
3d MD global regularity for small data: 1993 Georgiev
3d MD stationary solutions: Esteban, Georgiev, Séré 1996 EGS
2d DKG GWP: Grunrock and Pecher
2d MD GWP: [d’Ancona and Selberg 2010]((http://arxiv.org/abs/1004.1715 “Global well-posedness of the Maxwell-Dirac system in two space dimensions”)))

Are there stationary solutions for 2d MD? Are there other obstructions to decay/scattering? Are there size thresholds for the 3d MD setting. I should study the [EGS] works….

Local theory in 3d:

Gross 1966
Bournaveas 1996
Masmoudi and Nakanishi 2004
d’Ancona, Foschi, Selberg: Complete null structure of DKG 2007 and MD 2010 and almost optimal LWP.

2d DKG ↩

Charge class data and $(\phi, \phi_t) \in H^{1/2} \times H^{-1/2} (R^2)$.

LWP known for such data.
To get the global result, we need to control $D(t)$, which is his notation for the $H^{1/2} \times H^{-1/2}$ size of the evolving solution $(u(t), \partial_t u(t))$.

Theorem (Grunrock-Pecher 2010): 2d DKG is LWP up to time $T>0$ s.t $T^{1/2} [1 + D(0)] \leq \epsilon$. Moreover, $\sup_{|t| \leq T } D(t) \leq D(0) + C T^{1/2}$ with $C$ dependent on charge constant.

The globalizing procedure follows a general argument introduced by Colliander, Holmer and Tzirakis 2008.

How does it go?

$T^{1/2} [1 + D(0)] = \epsilon/2$
$T^{1/2} \sim \frac{1}{D(0)}$
You iterate $n$ steps and accumlate errors until you grow until $n C T^{1/2} \sim D(0)$. This develops the solution onto a time interval of size $nT \sim 1$ so you have advanced the solution to a local interval whose length only depends upon the charge. Therefore, you can iterate this process to make it go global.

We want to apply this procedure to do the same for MD.

Theorem (d’Ancona and Selberg 2010): 2d MD is LWP up to time $T>0$ s.t. $T^{1/2} [1 + D_T (0] \leq \epsilon$ where $\epsilon$ depends upon the chage constant. Moreover, $$\sup_{|t| \leq T} D_T (t) \leq D_T (0) + C T^{1/2} \log (\frac{1}{T}).$$

Corollary: 2d MD is GWP.

The iteration procedure is more involved than the CHTz scheme due to a logarithmic loss. There is an intermediate iteration which reduces matters to a harmonic series! This was exposed nicely so I watched it without typing…..

What lies behind the proof? ↩

LWP, Subcriticality
Growth estimate for EM field

Key Points

Null structore of nonlinear terms
Refined bilinear estimates needed to exploit the structure
Subcriticality is crucial

He went on to describe the null structure and ideas in extracting the required quantitative slack in the local theory to run the globalization scheme.

Jason Metcalfe: Long time existence for nonlinear wave equations in exterior domains ↩

(many years of collaboration with Sogge and Nakamura)

At the quadratic level, all that work only with derivative terms and not terms involving the solution itself.

Let $K$ be a compact obstacle with $C^\infty$ boundary. We want to solve, in dimensions 3 and 4,

Problem $S$: $$\square u = |u|^p$$ in the exterior of $K$ with vanishing neumann condition and small initial data. Let’s call this problem $S$. Let’s assume here that $K$ is nontrappling

Problem $Q$: $$ \square u = Q(u, u’, u”) $$ with vanishing Dirichlet boundary conditions and with $K$ starshaped.

There are issues that make it difficult to work with the Klainerman vector field, especially the boosts and the scaling vector fields.

In this work, we will only work with $Z = [\partial_i, \Omega_{jk}= x_j \partial_k - x_k\partial_j]$.

Localized Energy Estimate: $$ [\log(2+T)]^{-1/2} | \langle x \rangle^{-1/2} u’ |{L^2{tx}} \lesssim | u’(0)|{2} + \int0^T | \box u(s, \cdot)|_2 ds. $$

A weighted Sobolev inequality $$ R^{(n-1)/2} | h |{L^\infty (\frac{R}{2}< |x| < R) \lesssim | Z^{\leq \frac{n+2}{2}} h |{L^2 (\frac{R}{4}< |x| < 2R)}. $$

KSS: $$ [\log(2+T)]^{-1/2} | \langle x \rangle^{-1/2} Z^{\leq 10} u’ |{L^2{tx}} \lesssim \epsilon + \int_0^T | Z^{\leq 10} (\partial_t u)^2|{L^2} dx \lesssim \epslion + | \langle x \rangle^{-1/2} Z^{\leq 10}u |^2{L^2_{tx}}. $$

Problem $S$: ↩

$p> p_c$ where $p_c > 0$ solves $(n-1) p_c^2 - (n+1) p_c - 2 = 0$:

$n=3 \implies p_c = 1 + \sqrt{2}$
$n=4 \implies p_c =2$.

Theorem (Hidono-Metcalfe-Smith-Sogge-Zhou): $n=3,4; ~p_c < p < \frac{n+3}{n-1}, \gamma = \frac{n}{2}- \frac{2}{p-1}.$ $$\sum_{|\alpha | \leq 2} ( \| Z^\alpha u(0, \cdot)\|{\dot{H}}^\gamma + \| Z^\alpha \partial_t u(0, \cdot)\|{\dot{H}}^{\gamma -1} ) < \epsilon $$ $\implies global existence.

…

Problem $Q$: ↩

No boundary case
- $n=3: \frac{c}{\epsilon^2}$ is the life-span (Lindblad and Hörmander)
- $n=4: exp(\frac{C/\epsilon})$ is the life-span (Lindblad and Hörmander)
$(\partial^2_u Q)(0,0,0)$ (this kills $u^2$ terms, but we are considering here the startshaped boundary.)
- $n=3:$ (in progress with a student John Helms)
- $n=4: \infty$ is the lifetime (Metcalfe-Sogge)

Scipio Cuccagna: The Hamiltonian structure of the nonlinear Schrödinger equation and the asymptotic stability of its ground states ↩

Reference

arXiv:0910.3797v5 Abstract: In this paper we prove that ground states of the NLS which satisfy the sufficient conditions for orbital stability of M.Weinstein, are also asymptotically stable, for seemingly generic equations. Here we assume that the NLS has a smooth short range nonlinearity. We assume also the presence of a very short range and smooth linear potential, to avoid translation invariance. The basic idea is to perform a Birkhoff normal form argument on the hamiltonian, as in a paper by Bambusi and Cuccagna on the stability of the 0 solution for NLKG. But in our case, the natural coordinates arising from the linearization are not canonical. So we need also to apply the Darboux Theorem. With some care though, in order not to destroy some nice features of the initial hamiltonian.

(This talk relates to the talk of Schlag.)

It seems to me this talk is also closely related to the talks of Marzuola and Muñoz.

We study the nonlinear Schrödinger equation:

$$ i u_t = -\Delta u + V(x)u + \beta (|u|^2) u, ~ in R^3$$

Results for this work do not work for $ i u_t = -\Delta u -|u|^{p-1} u$ with $p < 1 + \frac{4}{n}$.

We assume existence of a family of ground states. When they are gound states the look like you expect but he also had a graph involving nodes and I didn’t understand…

Notions of stability:

linear stability (i.e. Weinstein’s sufficient hypotheses for orbital stability)
- Only for ground states?
orbital stability
asymptotic stability
- $\lim_{t \rightarrow +\infty} \| u(t,x) - e^{i \theta(t)} \phi_{\omega+} (x) - e^{it\Delta} (h+)(x)\|{H^1x} = 0$
CONJECTURE: 1. $\iff$ 2. $\iff$ 3.
Theorem: 1. $\implies$ 3. generically.

Specifically, we prove nonlinear Fermi golden rule (terminology introduced by Soffer and Weinstein, Buslaev and Perlman used different terminology):

a. Some key coefficients are $\geq 0$; b. Generically they are $>0$.

One wants to prove that the remainder scatters. We have discrete and continuous modes. One wants to find a way to describe a mechanism of transfer from the discrete modes into the continuous modes. We want some way of writing the coordinates of the dynamics to reveal a damping effect in the discrete modes due to the transfer of the energy from the discrete modes into the continuous modes. The description of this transfer mechanism is the goal of the Fermi golden rule.

Asymptotic stability is analogous to showing that $u(t)$ solving an NLS-type equation is not only of the same siaze in $H^1$ for all time but also showing that the solution scatters. This is the analysis we want to do on the remainder. Eigenvalues obstruct asymptotic stability. That explains the preoccupation of Schlag with his proof of the nonexistence of eigenvalues in the gap.

Near ground states, we write the solution in a canonical way as a sum of a modulated ground state plus a remainder term. The NLS can be recast as a dynamical system of the phase and scaling parameter coupled to the (presumably dispersive) behavior. He then changes variables so that the system is expressed as a matrix equation in which the “Hamiltonian structure is obscured”. This is the standard way in the literature that the system is expressed. But somehow this way of writing it is wrong. (???)

He makes some assumptions about the absence of embedded eignevalues. He suggests this hypothesis is not necessary but is not certain….some discussion with Tataru.

He writes on the board a horizontal line and draws points at 0, and sa few eignevalues parametrized by $\omega$. He then draws wavy stuff over the right half starting at some point to the right of the eigenvalues representing the continuoys spectrum.

…slides are coming fast and they are too dense for me to type in real time….

Alexandru Ionescu: Uniquness theorems in general relativity ↩

General relativity…

Spacetimes ↩

Spacetimes $(M^4, g)$ are solutions of the Einstein vacuum equations $$R_{\mu \nu} = g^{\alpha \beta} R_{\alpha \beta \mu \nu} = 0$$

The metric is in 4 dimenions, it has 10 components. The Riemann tensor has 20 components. These are 10 equations for the 20 components.

Minkowski space: $(R^3 \times R, -dt^2 + dx^2 + dy^2 + dz^2)$.

Besides being Ricci flat, in fact this solution also has zero Riemann tensor and this condition completely characterizes the Minkowski space.

Schwarzschild spaces:

$ds^2 = -(1 - \frac{2m}{r}) dt^2 + (1 - \frac{2m}{r})^{-1} dr^2 + r^2 (d\theta^2 + (\sin \theta)^2 d\phi^2)$ where $(r, t, \theta, \phi) \in (2m, \infty)\times R \times (0, \pi) \times S^1$. It took several decades to realize that $r=2m$ is merely a coordinate singularity. This was realized with the Kruskal coordinates in which the metric may be expressed $ds^2 = F^2 (-dt^2 + (dx’)^2) + r^2 (d\theta^2 + (\sin \theta)^2 d\phi^2)$. The Kruskal picture is the region between the lobes of hyperboloid of two sheets. The region below the lobes and inside the $|y| = |x|$ cone regions containing the lobes is called the black hole. The domain of outer communication is outside the cone.

Kerr spaces:

$m$ is the mass of the black hole and $J$ is the angular momentum of the black hole. We assume $m>0, ~a = \frac{J}{m} \in [0, m)$ and let $r_+ = m+ (m^2 - a^2)^{1/2}. In Boyer-Lindquist coordinates $(r, t, \theta, \phi) \in (r_+, \infty)\times R \times (0, \pi) \times S^1 $$ -\frac{\rho^2 \Delta}{\Sigma^2} dt^2 + \frac{\Sigma^2 (\sin \theta)^2}{\rho^2 }( d\phi - \frac{2amr}{\Sigma^2}dt)^2 + \frac{\rho^2}{\Delta}(dr)^2 + \rho^2 (d\theta)^2$$ where

$\Delta -= r^2 + a^2 - 2mr$
$\rho^2 = ….$ slide changed….
$\Sigma^2 = …

For Minkowski, 20 of 20 components of the Riemann tensor vanish. For Schwarzshild 19 of the 20 componenents of the Riemann tensor vanish in the right coordinates. For Kerr, 18 of the 20 components vanish in the right coordinates.

Key properties of Kerr spacetimes: ↩

Solutions of the einstein vacuum equations $R_{\alpha \mu} = 0$;
Killing vector field $T = \partial_t$ timelike at “infinity”;
Killing vector field $Z - \partial_\phi$ wiht closed orbits;
Geometric properties: asymptotic flatness, smooth bifurcate sphere, global hyperbolicity;
Rigididty: Kerr spaces are real-analytic.

“No hair” theorems: such properties charcaterize the Kerr spaces (Carter, Robinson, Hawking-Ellis, Mazur, Bunting, Weinstein, Chrusciel-Costa). “We are trying to understand final states.”

Main Conjecture: If $(M^4, g, T)$ is a regular stationary vacuum, then the domain of outer communication of $M^4$ is isometric to the domain of outer communication of some Kerr spacetime of mass $m$ and angular momentum $ma$, $a \in [0, m)$.

What is regular in the conjecture? It took a long time to characterize what that means

There is a lot of supporting evidence.

Carter 1971: axially symmetric black holes have only 2 degrees of freedom
- Mathematically, an imprecise statement. It said there are “no bifurcations”
Robinson 1975: the uniqueness conjecture holds in the case of axially symmetric black holes
- global argument involving the whole space.
Hawking-Ellis 1973: the conjecture holds in the case of real-analytic spacetimes.

Hawking’s strategy is to define an additional Killing vector-field in the spacetime and reduce to the Carter-Robinson theorem. The assumption of real analyticity is not what you really want.

Theorem 1 (Ionescu-Klainerman): The conjecture holds provided that the scalar identity is assumed to be satisfied on the bifurcation sphere.

Theorem 2 (Alexakis-Ionescu-Klainerman): The conjecture holds proved that the spacetime is assumed to be “close” to a Kerr spacetime.

* Theorem 3 (Aliexakis-Ionescu-Klainerman):* Assume $\cal{N}, \underline{\cal{N}}$ are smooth, null, nnexpanding hypersurfaces in an Einstein vacuum $(O, g)$ which intersect transversally in a 2-sphere $Z$. Then there is an opern neighborhood $O’$ of $Z$ and a nontrivial Killing vector-field $K$ in $O’$ which is tangent to the null generators of $\cal{N}, \underline{\cal{N}}$.

This is a local version of Hawking’s Rigidity Theorem, without assuming analyticity of the spacetime.

Construct the Hawking v. K in the domoan of dependence of $ \cal{N}\cup \underline{\cal{N}}$ (Friedrich-Racz-Wald)
Extend the v.f to a full neighborhood of Z by solving a transport equation $[L, K] = cL

Key steps in our strategy:

We deine some tensors: $\pi_{\alpha \beta}, W_{\alpha \beta \mu \nu}.
Prove a system of wave/transport equations of the form:
- $\square_g W = {\cal{M}} (W, Dw, \pi, D\pi)$
- D_L pi ={\cal{M}} (W, Dw, \pi, D\pi)$
Use Carleman estimates and a unique continuation argument to conclude that $W, \pi$ vanish in a neighborhood of $Z$.

Model Theorem (I-Klainerman): Assume $\phi \in C^2 (M)$ and $A, B^l \in C^0 (M)$ for $l = 0, \dots, d$. Assume that $$ \square \phi = A \phi + \sum D^l \cdot \partial \phi…..ack slide change.

Unique Continuation: assume $\phi$ is smooth in $(O,g)$ and solves a wave equation $D^\alpha D_\alpha \phi = A \phi + B^\alpha D_\alpha \phi.$ Assume $\phi$ vanishes in the set $[h<0]$, where $h \in C^\infty (O), ~ \nabla h \neq 0.$ Does $\phi$ vanish in a neighborhood of $[h \leq 0]$?

Suppose we have $T(u)=0 $ in $B$. Suppose $u_1, u_2$ solutions in $B$ and $u_1 \sim u_2$ inside small set $A \subset B$. Basically, there are three possibilities:

lack of uniquneess: $u_1 = u_2$ inside $A$ but $u_1$ is far from $u_2$ in the big set $B$.
Well-posedness: If $u_1$ is close to $u_2$ in $A$ then $u_1$ is close to $u_2$ in $B$.
Unique continuation:
- If $u_1 = u_2$ in $A \implies u_1 = u_2$ in $B$.
- If $u_1$ is close to $u_2$ in $A$ we are unable to conclude that $u_1$ is close to $u_2$ in $B$

Hormander’s pseudo-convexity condition: Unique continuation holds if $X^\alpha X^\beta D_\alpha D_\beta h < 0$…ack slide change.

The method is based on Carleman Estimates.

Model theroem in Kerr spaces (I-Klainerman): Assume $W, A, B, C$ are smooth tensors in the Kerr space $K^4$, and $$ \square_{g_0} W = A \cdot W + B \cdot D W $$ $$ {\cal{L}}_T W = C \cdot W. $$ If $W=0$ on the horizon then $W=0$ everywhere outside.

T-conditional pseudoconvexity property

We would really like a tensor $\cal{S}$ (an analog of the Riemann tensor $\cal{R}$) which has the following properties:

It describes locally the Kerr spaces
It satisfies a suitable geometric equation of the form $$ \square_g \cal{S} = A \cdot W + B \cdot DW$$ $$ {\cal{L}}_T W = C \cdot W $$ We want to then uniquely continue the vanishing.

Mars-Simon found such a tensor. This is a tensor for Kerr which is analogous to the Riemann tensor for Minkowski. The Riemann tensor is a local quantity which characterizes the Minkowski space in the sense that when it vanishes, we know that we are in Minkowski space. Similarly, the Mars-Simon tensor characterizes Kerr. To go from local vanishing to conclude global vanishing, we need an analytic continuation.

More precise statement of Theorem 1: The domain of outer communication $E$ of a regular stationary vacuum $(M, g, T)$ is locally isometric to the domain of outer communication of a Kerr spacetime, provided that the identity $$ - 4 m^2 {\cal {F}}^2 = (1 - \sigma)^4 $$ holds on the bifurcation sphere $S_0$.

More precise statement of Theorem 2: The domain of outer communication $E$ of a regular stationary vacuum $(M,g,T)$ is isometric to the domain of outer communication of a Kerr spacetime, provided that the smallness condition $$ | (1 - \sigma) {\cal{S}} (T, e_\alpha, e_\beta, e_\gamma)| \leq {\overline{\epsilon}} $$ holds along a Cauchy hypersurface in $E$ for some sufficiently small ${\overline{\epsilon}}.$

Postlude ↩

I discussed with Alex whether one could (or should….) formulate a statement similar to Theorem 2 about Minkowski space using the Riemann tensor like: Suppose that the Riemann tensor is small on some (small? geometric conditions?) set $A$ inside a spacetime manifold $(M,g)$. Can one conclude that the Riemann tensor must therefore vanish on $A$ or perhaps on a bigger set $B$? One can view the [AIK] theorem 2 as a Liouville-type theorem: a smallness condition on the Mars-Simon tensor on a subset of $(M,g)$ with certain conditions implies that the Mars-Simon tensor vanishes. Is there a corresponding Liouville-type theorem where smallness of the Riemann tensor on an appropriate subset implies that the Riemann tensor actually vanishes?

dispersive | James Colliander

Zaher Hani Seminar: A Continuum Large Box Limit for the Cubic Nonlinear Schrödinger Equation

Setup

Weak turbulence “paradigm”

A New Limiting equation

Properties of the continuum equation

Rigorous Approximate Results

Further Questions

GWP of Gross-Pitaevskii Equation on R4

Princeton Math Colloquium: Soliton Resolution for the Radial Energy Critical Focusing Nonlinear Wave Equation

Presenter: Carlos Kenig

Linear wave equation

Critical focusing NLW

Small data theory for NLW

Recall a collection of recent results

Existence of Type II blowup solutions

Soliton resolution for radial solutions of NLW

Ideas of the proof (global case)

Notes on Nonlinear Dispersive Wave Equations Workshop in Oberwolfach

Table of Contents

Jérémie Szeftel: Alas, I missed this talk….thanks Air Canada. ↩

arXiv ↩

slides ↩

Sebastian Herr: Small data theory for energy critical periodic NLS ↩

Warm-up remarks ↩

$NLS^{\pm}_5 (T^3)$ ↩

New Strichartz Estimates ↩

Perturbative Analysis ↩

Trilinear Strichartz ↩

Sketch of proof ↩

Contraction estimate ↩

Remarks ↩

Questions ↩

Benjamin Schlein: Effective evolution equations from many body quantum dynamics ↩

Introduction ↩

Mean Field Regime

Reduced Densities

Boson Stars ↩

Dynamics of Bose-Einstein Condensates ↩

Adrian Constantin: Camassa-Holm ↩

Physical Background ↩

Emergence of Camassa-Holm Equation ↩

Geometric viewpoint as a geodeisc on the diffeomorphism group ↩

Integrable Sturcture ↩

Claudio Muñoz: Dynamics of gKdV solitons under perturbations by potentials in front of nonlinear term ↩

Mihalis Dafermos: Superradiance, trapping and decay for waves on Kerr spactimes in the general subextremal case $|a| < M$. ↩

Penrose Diagrams (images taken from Dafermos-Rodnianski)

Boundedness and decay for $\square_g \psi =0$ on Schwarzschild and Kerr ↩

Current state of the art for the quantitative study of $\square_g \psi = 0$ ↩

Review of the main features of Kerr Spacetimes ↩

Red-shift

Superradiance

Trapping

Proof of integrated local energy decay. ↩

Separation.

General Idea

Some other important results

Open Problems ↩

Stephen Gustafson: Dynamics on near-harmonic Schrödinger and Landau-Lifschitz maps ↩

Regurlarity vs. Singularity: energy critical problems ↩

Equivariant Maps ↩

New results: global solutions for degree 2 (LL) with $a_1 > 0$. ↩

Standard “modulation theory” approach ↩

Dispersive/diffusive estimates

A remedy for $m \leq 3 $ and its cost. ↩

Conclusions ↩

Ioan Bejenaru: Near soliton evolution in 2d Schrödinger Maps ↩

Large Data Theory ↩

Equivariant Harmonic Maps on $S^2$. ↩

Basic setup for stability/instability ↩

Modulation Theory ↩

Frank Merle: Isolatedness of characteristic points for blow-up solutions of semilinar wave equation ↩

Background References (Incomplete)

Semilinear Wave Equation, Blowup Surface ↩

Summary of Results ↩

Ben Dodson: Defocusing $L^2$-Critical NLS ↩

References

Mass-Critical NLS ↩

Minimal Mass Blowup Solution Strategy ↩

Galilean Invariance Observations ↩