Fast Diffusion to Self-Similarity: Complete Spectrum, Long-Time Asymptotics, and Numerology

Abstract.The complete spectrum is determined for the operator on the Sobolev space W1,2ρ(Rn) formed by closing the smooth functions of compact support with respect to the norm Here the Barenblatt profile ρ is the stationary attractor of the rescaled diffusion equation in the fast, supercritical regime m the same diffusion dynamics represent the steepest descent down an entropy E(u) on probability measures with respect to the Wasserstein distance d2. Formally, the operator H=HessρE is the Hessian of this entropy at its minimum ρ, so the spectral gap H≧α:=2−n(1−m) found below suggests the sharp rate of asymptotic convergence: from any centered initial data 0≦u(0,x) ∈ L1(Rn) with second moments. This bound improves various results in the literature, and suggests the conjecture that the self-similar solution u(t,x)=R(t)−nρ(x/R(t)) is always slowest to converge. The higher eigenfunctions – which are polynomials with hypergeometric radial parts – and the presence of continuous spectrum yield additional insight into the relations between symmetries of Rn and the flow. Thus the rate of convergence can be improved if we are willing to replace the distance to ρ with the distance to its nearest mass-preserving dilation (or still better, affine image). The strange numerology of the spectrum is explained in terms of the number of moments of ρ.