Decay of waves IIIa: nonlinear tails in Minkowski space redux
by Willie Wong
Before we move on to the geometric case, I want to flesh out the nonlinear case mentioned in the end of the last post a bit more. Recall that it was shown for generic nonlinear (actually semilinear; for quasilinear and worse equations we cannot use Duhamel’s principle) wave equations, if we put in compact support for the initial data, we expect the first iterate to exhibit a tail. One may ask whether it is possible that, in fact, this is an artifact of the successive approximation scheme; that in fact somehow it always transpires that a conspiracy happens, and all the higher order iterates cancel out the tail coming from the first iterate. This is rather unlikely, owing to the fact that the convergence to is dominated by a geometric series. But to just make double sure, here we give a nonlinear system of wave equations such that the successive approximation scheme converges after finitely many steps (in fact, after the first iterate), and so we can also explicitly compute the rate of decay for the nonlinear tail. While the decay rate is not claimed to be generic (though it is), the existence of one such example with a fixed decay rate shows that for a statement quantifying over all nonlinear wave equations, it would be impossible to demonstrate better decay rate than the one exhibited.
In this particular case we will consider the following system of wave equations on 1+3 dimensional Minkowski space, with the assumption of spherical symmetry in force.
The system as a whole, as an equation on the vector is nonlinear, since the right hand side depends nonlinearly on the vector itself. However, we recognize that taken individually the components of the vector decouple, and each component solves a (possibly inhomogeneous) wave equation with no nonlinearities. (So we cheat a little bit. But this is how we can make the successive approximation scheme converge after finitely many steps.)
Recall that under the spherical symmetry change of coordinates, writing again and , we can rewrite the equations as
Why this particular form on nonlinearities? It is known that generic nonlinearities, quadratic in , lead to formation of singularities in finite time, with the number of spatial dimensions is less than or equal to three. In dimension 1, we can explicitly illustrate this: consider the wave equation . Then along the curves of constant , the quantity verifies an equation of the form , whose solution is . So if is negative at any point when , the solution will develop a singularity in finite time.
In low spatial dimensions, the rates of dispersion are slower (the energy density is expected to drop as the inverse of the area of the wave front, which in dimensions grows like ), and compared to the positive feedback coming from a generic quadratic derivative nonlinearity, may not be enough to radiate away the self-interaction sufficiently fast. This motivated my thesis advisor, Sergiu Klainerman, to study, when he was doing his own PhD, the conditions sufficient to guarantee the quadratic nonlinearities are not too strong. This lead to the recognition that the so-called null structures play an important role in the study of wave equations in low spatial dimensions. To cut the story short, the form chosen for the nonlinearity above renders it a null form, so it should be better behaved (decay faster) than the most general quadratic type nonlinearities. By choosing this example, my goal is to illustrate the fact that even in the good scenarios, one expect to have some fairly strong nonlinear tail effects. (A detailed discussion of this topic would take us too far from the focus of this post. Perhaps I’ll write about it in a part IV to this series.) (That we need a and a is also related. If we only have the function and choose the nonlinearity for to be , we see that the null form cancellation is too strong, and we won’t be able to see the effect at the level of the first iterate.)
Anyway, back to the example.
The key observation is captured already in Figure 4 (which we recall again here above): when the red dot is “sufficiently far in the future”, the crosshatched region only intersects the “outgoing legs” of the cyan region (using the fact the incoming wave, after bouncing off the center axis, becomes outgoing also). In other words, where the crosshatched region and the cyan region intersects, the wave is purely outgoing.
Now suppose the cyan region is a representation of the support of and , and we are trying to evaluate at the red dot (the part coming from initial data can be neglected due to linearity). To do so we integrate in the intersection of the crosshatched and cyan regions. Now in that region, Equation 21′ implies that and are outgoing waves, that is there derivatives vanish. Furthermore the wave equation also means their derivatives are constant along lines of constant .
Using the constancy of along constant lines, we have that , and similarly . So we get that
Now, the area integrated over is determined both by the cyan and crosshatched regions. The cyan region coming from determins the limits of the integration, while the crosshatched region determines the limits. We will just assume that the limits are some fixed values and , which is completely determined by the compact support of the initial data. The limits, however, depends on where we will evaluate . If we choose to evaluate at point , we see that by the requirement on the central axis, the lower bound of the region is precisely , and the upperbound . Therefore we have that, for the red dot sufficiently far in the future where the picture in Figure 4 is a good representation,
Now, sufficiently far into the future, is sufficiently large, and hence is very large. This means that we can throw away the second term in the parentheses as negligible. For generic and (note that if there is actually a cancellation so that the contribution from the nonlinearity vanishes), we can find a lowerbound once become large enough (the upper bound is trivial)
And in Equation 25 we capture the behaviour of the nonlinear tail. For any fixed , the integral . So we have that
giving us the “interior decay” rate. On the otherhand, if we were to take the limit toward with fixed, we have that the total integral is controlled by , and so we have the “radiative decay” rate coming purely from the weight:
for any sufficiently large fixed .
In both Equations 26a and 26b, the “eventually vanishing” phenomena observed in the linear case for data with compact support are destroyed.
Lastly, just note that Equations 26a,b can be glued together with boundedness on bounded regions to get the global estimate
which captures a more generic (in the sense of stability under nonlinear perturbations) behaviour of spherically symmetric waves on 1+3 dimensional Minkowski space.