I’ve just uploaded to the arXiv my paper “The high exponent limit $p \to \infty$ for the one-dimensional nonlinear wave equation“, submitted to Analysis & PDE.  This paper concerns an under-explored limit for the Cauchy problem

\displaystyle -\phi_{tt} + \phi_{xx} = |\phi|^{p-1} \phi; \quad \phi(0,x) = \phi_0(x); \quad \phi_t(0,x) = \phi_1(x) (1)

to the one-dimensional defocusing nonlinear wave equation, where \phi: {\Bbb R} \times {\Bbb R} \to {\Bbb R} is the unknown scalar field, p > 1 is an exponent, and \phi_0, \phi_1: {\Bbb R} \to {\Bbb R} are the initial position and velocity respectively, and the t and x subscripts denote differentiation in time and space.  To avoid some (extremely minor) technical difficulties let us assume that p is an odd integer, so that the nonlinearity is smooth; then standard energy methods, relying in particular on the conserved energy

\displaystyle E(\phi)(t) = \int_{\Bbb R} \frac{1}{2} |\phi_t(t,x)|^2 + \frac{1}{2} |\phi_x(t,x)|^2 + \frac{1}{p+1} |\phi(t,x)|^{p+1}\ dx, (2)

on finite speed of propagation, and on the one-dimensional Sobolev embedding H^1({\Bbb R}) \subset L^\infty({\Bbb R}), show that from any smooth initial data \phi_0, \phi_1, there is a unique global smooth solution \phi to the Cauchy problem (1).

It is then natural to ask how the solution \phi behaves under various asymptotic limits.  Popular limits for these sorts of PDE include the asymptotic time limit t \to \pm \infty, the non-relativistic limit c \to \infty (where we insert suitable powers of c into various terms in (1)), the small dispersion limit (where we place a small factor in front of the dispersive term +\phi_{xx}), the high-frequency limit (where we send the frequency of the initial data \phi_0, \phi_1 to infinity), and so forth.

Tristan Roy recently posed to me a different type of limit, which to the best of my knowledge has not been explored much in the literature (although some of the literature on limits of the Ginzburg-Landau equation has a somewhat similar flavour): the high exponent limit p \to \infty (holding the initial data \phi_0, \phi_1 fixed).  From (1) it is intuitively plausible that as p increases, the nonlinearity gets “stronger” when |\phi| > 1 and “weaker” when |\phi| < 1; the “limiting equation”

\displaystyle -\phi_{tt} + \phi_{xx} = |\phi|^{\infty} \phi; \quad \phi(0,x) = \phi_0(x); \quad \phi_t(0,x) = \phi_1(x) (3)

would then be expected to be linear when |\phi| < 1 and infinitely repulsive when |\phi| > 1 (i.e. in the limit, the solution should be confined to range in the interval [-1,1], much as is the case with linear wave and Schrödinger equations with an infinite barrier potential; though with the key difference that the nonlinear barrier in (3) is confining the range of \phi rather than the domain.).

Of course, the equation (3) does not make rigorous sense as written; we need to formalise what an “infinite nonlinear barrier” is, and how the wave \phi will react to that barrier (e.g. will it reflect off of it, or be absorbed?).  So the questions are to find the correct description of the limiting equation, and to rigorously demonstrate that solutions to (1) converge in some sense to that equation.

It is natural to require that \phi_0 stays away from the barrier, in the sense that |\phi_0(x)| < 1 for all x; in particular this implies that the energy (2) stays (locally) bounded as p \to \infty; it also ensures that (1) converges in a satisfactory sense to the free wave equation for sufficiently short times.  For technical reasons we also have to make a mild assumption that either of the null energy densities \phi_1 \pm \partial_x \phi_0 vanish on a set with at most finitely many connected components.  The main result is then that as p \to \infty, the solution \phi = \phi^{(p)} to (1) converges locally uniformly to a Lipschitz, piecewise smooth limit \phi = \phi^{(\infty)}, which is restricted to take values in [-1,1], with -\phi_{tt}+\phi_{xx} (interpreted in a weak sense) being a negative measure supported on \{ \phi=+1\} plus a positive measure supported on \{\phi = -1\}.  Furthermore, we have the reflection conditions

\displaystyle (\partial_t \pm \partial_x) |\phi_t \mp \phi_x| = 0.

It turns out that the above conditions uniquely determine \phi, and one can even solve for \phi explicitly for any given data; such solutions start off smooth but pick up an increasing number of (Lipschitz continuous) singularities over time as they reflect back and forth across the nonlinear barriers \{\phi=+1\} and \{\phi=-1\}.  (An explicit example of such a reflection is given in the paper.)

[The above conditions vaguely resemble entropy conditions, as appear for instance in kinetic formulations of conservation laws, though I do not know of a precise connection in this regard.]

In the remainder of this post I would like to describe the strategy of proof and one of the key a priori bounds needed.  I also want to point out the connection to Liouville’s equation, which was discussed in the previous post.

— Strategy —

The top level strategy is based on compactness methods: first show that the family \phi^{(p)} of solutions is compact in a suitable topology, show that any limit point \phi of these solutions obeys the properties stated above, and then show that these properties uniquely determine the limit \phi.

The compactness (in the local uniform topology) is easy, coming from the Arzelà-Ascoli theorem and energy conservation (and relies crucially on the fact that the one-dimensional NLW (1) remains subcritical even as p \to \infty; I do not know what to predict in the critical case of two dimensions, or the supercritical case of higher dimensions).  The uniqueness is also not too difficult, as one can solve for \phi obeying these conditions by quite classical means based on the method of characteristics (i.e. following the solution along null rays) and the fundamental theorem of calculus.  The main difficulty is to show that a limit point \phi of the solutions \phi^{(p)} actually do satisfy all of the properties listed above.  For this, it turns out to be necessary to establish a number of a priori estimates on \phi^{(p)} and its first derivatives that will survive the passage to the limit.  Some of these estimates are based on the pointwise conservation laws for Liouville’s equation that I discussed in the previous post; but now I want to turn to another key estimate, namely the pointwise estimate

\displaystyle |\phi^{(p)}(t,x)| \leq 1 + \frac{\log p}{p} + O( \frac{1}{p} ) (4)

as p \to \infty, which of course will keep the limiting solution \phi^{(\infty)} confined to range in the interval [-1,1].  This bound turns out to be best possible, and is important in bounding the nonlinear effects of (1) (for instance, it implies that the nonlinear term in (1) is pointwise bounded by O(p) ).

The connection to the Liouville equation arises from the ansatz

\displaystyle \phi^{(p)} = p^{1/(p-1)} (1 + \frac{1}{p} \psi( p t, p x ) )

then (1) becomes (for positive \phi)

\displaystyle -\psi_{tt} + \psi_{xx} = (1 + \frac{1}{p} \psi)^p

which formally converges to the Liouville equation - \psi_{tt} + \psi_{xx} = e^\psi in the limit p \to -\infty.  For comparison, the estimate (4) corresponds in this ansatz to an upper bound \psi \leq O(1) on \psi.

— Confinement —

Now we prove (4).  To simplify the discussion let us assume that \phi = \phi^{(p)} is always non-negative (one can restrict to this case by sign reversal symmetry, finite speed of propagation, and energy conservation).  As in the previous post, it is convenient to work in null coordinates

\displaystyle u = t+x, v = t-x

in which case the equation (1) becomes

\displaystyle \phi_{uv} = - \frac{1}{4} \phi^p.  (5)

In particular, \phi_u is decreasing in the v direction, and \phi_v is decreasing in the u direction.  With fixed initial data, this gives us an upper bound

\displaystyle \phi_u, \phi_v \leq O(1) (6)

uniformly in p.

Now suppose that (t_0,x_0) is the first time and place in which (4) fails, in the sense that

\displaystyle \phi(t_0,x_0) = 1 + \frac{\log p}{p} + \frac{K}{p}

for some large constant K.  Then from (6) we know that \phi is also large just a bit before (t_0,x_0); in particular, we will have

\displaystyle \phi(t,x) \geq 1 + \frac{\log p}{p}

on the spacetime diamond

\displaystyle \{ (t,x): t_0 + x_0-\frac{cK}{p} \leq t + x \leq t_0 + x_0; t_0 - x_0-\frac{cK}{p} \leq t - x \leq t_0 - x_0 \}

where c is a small cosntant (independent of K).  Inserting this into (5), we obtain

\displaystyle \phi_{uv} \leq - c' p

on this diamond, for some other small constant c’ > 0.  Integrating this on this diamond and using the various bounds on \phi at the corners, one arrives at

\displaystyle -c'' \frac{K^2}{p} = O(\frac{K}{p})

for some other small constant c” > 0, which leads to a contradiction for K large enough (independently of p).