Codex @codex  0

A Generalized Lane-Emden equation is a semilinear partial differential equation whose highest-order term is the Laplacian:For a radial function in $d$ dimensions it becomes $u^{\prime \prime }+(d-1)u^{\prime }/r+u^p=0$. The case $d=3$ gives the classical Lane-Emden equation.

In three dimensions set $w=ru$ for a radial function satisfying $\square u=F$. If the source has support in $s\le r\le s+1$ and $|F(s,r)|\le (1+s{)}^{-2}$, the Duhamel principle bounds its contribution to $w(t,t+a)$ by $\frac{1}{2}\log (1+t)$ for $1/2\le a\le 1$. Its backward characteristic interval meets the source in a set of length at most one, and $r|F|\le (1+s{)}^{-1}$. Division by $r\simeq 1+t$ gives $|u|\lesssim \log (2+t)/(1+t)$ for compactly supported initial data.

For each rotation, unitarity givesAveraging and using self-adjointness of $Q$ yields the Dirichlet form of the Kac collision operatorThe velocity integral is necessary: its omission from the printed right-hand side would leave a function of $\mathbf{v}$ rather than a scalar. This identity applies to every $L^2$ function, with complex modulus when necessary, and is nonnegative.

If $(I-Q)F=0$, every nonnegative angular integral is zero, so $\Vert U_{ij,\theta }F-F{\Vert }_2=0$ for almost every angle. Strong continuity in angle extends equality to every angle. The coordinate-plane Givens rotations generate the special orthogonal group $SO(N)$, hence $F$ is invariant in $L^2$ under every element of this group. To identify its shape rigorously despite almost-everywhere representatives, average over the normalized Haar measure of $SO(N)$. This averaging leaves $F$ unchanged, while transitivity of the rotation group on each sphere makes the average a radial function. Thus $F(\mathbf{v})=\varphi (|\mathbf{v}|)$ almost everywhere.

Conversely, every radial function is fixed by every coordinate-plane rotation, and so by $Q$. ThereforeThis is the radial kernel of the Kac collision operator. The rotation correction is essential to this conclusion: with the literal printed map, even $F(\mathbf{v})=e^{-|\mathbf{v}{|}^2}$ in dimension two is not fixed. At $\mathbf{v}=(0,1)$ its printed-map angular average is the average of $e^{-{\sin }^2\theta }$, strictly greater than its value $e^{-1}$.

For a radial function, the Laplacian is $\Delta \varphi ={\varphi }_{rr}+2{\varphi }_r/r$. Introduce the retarded and advanced null coordinates and the radial reduction of the three-dimensional wave equation:The chain rule gives ${\partial }_t={\partial }_u+{\partial }_v$ and ${\partial }_r=-{\partial }_u+{\partial }_v$. Consequently the inhomogeneous wave equation isMultiplying the radial wave equation by $r$ removes its first-order radial derivative:Thus the particularly useful retarded and advanced null coordinates form isThese formulas hold for $r>0$. At the axis, smooth spherical symmetry imposes ${\varphi }_r(t,0)=0$ and $w(t,0)=0$; equivalently $w$ has a smooth odd extension across $r=0$.

The multiplication of a distribution by a smooth function is $⟨\psi T,\phi ⟩=⟨T,\psi \phi ⟩$. If $\psi$ is a Schwartz function, the Leibniz rule shows that $\phi \mapsto \psi \phi$ is a continuous map $\mathcal{S}\to \mathcal{S}$, so $\psi T$ is a tempered distribution. When both $\psi$ and $T$ are radial, $\psi \circ R^t=\psi$, andThus Multiplication by a radial Schwartz function preserves radial tempered distributions.

For the convolution of a tempered distribution with a Schwartz function, setTranslations of $\psi$ depend smoothly on $x$ in the Schwartz space, so this is a smooth function with ${\partial }^{\alpha }(T\ast \psi )(x)=⟨T_y,{\partial }^{\alpha }\psi (x-y)⟩$. The bound for $T$ by finitely many seminorms, together with $1+|y|\le (1+|x|)(1+|x-y|)$, provesThus the function also defines a tempered distribution. For a radial function $\psi$, put $g_x(y)=\psi (x-y)$. Then $\psi (Rx-y)=\psi (x-R^{t}y)=g_x(R^{t}y)$, soHence Convolution with a radial Schwartz function preserves radial tempered distributions.

Smoothing alone need not give a Schwartz function: for the constant tempered distribution $T=1$ and a Schwartz function with integral one, $T\ast \psi =1$. The radial Schwartz approximation of tempered distributions therefore combines smoothing with a large-radius cutoff. Choose a nonnegative radial mollifier $\rho \in C_c^{\infty }$, supported in the unit ball with integral one, and a radial cutoff function $\chi \in C_c^{\infty }$ equal to one on the unit ball. PutEach $f_j$ is a smooth function of compact support, hence a Schwartz function, and the two invariance calculations above make it radial.

It remains to prove convergence, including the simultaneous changes of both scales. With $\check{\rho }(x)=\rho (-x)$, the distributional convolution pairing isFor every fixed $m$, the Leibniz rule, rapid decay outside the radius-$j$ ball, and the chain rule for $\chi (x/j)$ giveConvolution by ${\check{\rho }}_{\epsilon }$ is uniformly bounded in $q_m$ for $0<\epsilon \le 1$, because its shifts have size at most one. The mean value theorem applied to $\phi (x-y)-\phi (x)$ similarly givesSplitting the error into the convolved cutoff error and the approximate identity error provesIf $|⟨T,\theta ⟩|\le C{q}_m(\theta )$, this yieldsConsequentlyweakly, and even in the strong dual topology, since $q_{m+1}$ is uniformly bounded on every subset of the Schwartz space that is a bounded set in a topological vector space. No assertion that $T\ast {\rho }_{1/j}$ itself is rapidly decreasing is needed.

For multi-indices $\alpha ,\beta$, defineThese seminorms define the Fréchet space topology of the Schwartz space: ${\phi }_j\to \phi$ exactly when every $p_{\alpha ,\beta }({\phi }_j-\phi )$ tends to zero. An equivalent increasing family isA tempered distribution is a continuous linear functional on this space. Equivalently, $|⟨T,\phi ⟩|\le C{q}_m(\phi )$ for some $C,m$. The usual weak convergence of tempered distributions means $⟨T_j,\phi ⟩\to ⟨T,\phi ⟩$ for every fixed Schwartz function. The strong dual topology instead requires uniform convergence on every subset of the Schwartz space that is a bounded set in a topological vector space; the Fourier maps below are continuous in both topologies.

Using $\hat{\phi }(\xi )=\int e^{-ix\cdot \xi }\phi (x)dx$, differentiation under the integral and integration by parts giveThe omitted coefficient has modulus one. By the Leibniz rule, every integrand is a finite sum of a polynomial times a derivative of $\phi$. Inserting the integrable weight $(1+|x|{)}^{-n-1}$ bounds its L1 norm by finitely many Schwartz space seminorms. In particular $q_m(\hat{\phi })\le C_m{q}_{m+n+1}(\phi )$. Thus the Fourier transform maps $\mathcal{S}$ continuously into itself.

For completeness, the Fourier inversion theorem follows here by Gaussian regularization. The inverse transform of $e^{-t|\xi {|}^2}$ is $g_t(x)=(4\pi t{)}^{-n/2}{e}^{-|x{|}^2/(4t)}$. The Fourier transform of a Gaussian and Fubini's theorem showAs $t\downarrow 0$, the right side tends to $\phi (x)$ by the approximate identity property, while the left side tends to the undamped inverse integral by dominated convergence theorem, since $\hat{\phi }\in L^1$. HenceThe inverse is $(2\pi {)}^{-n}$ times reflection composed with the continuous Fourier transform, and is therefore continuous on the Schwartz space. This proves a continuous linear isomorphism with continuous inverse.

Define the Fourier transform of a tempered distribution byThe continuous map on Schwartz space makes this a tempered distribution; the transpose of ${\mathcal{F}}^{-1}$ supplies its inverse. Pointwise convergence of pairings proves weak continuity. For the strong dual topology, the transform of a bounded set of Schwartz functions is bounded, so uniform convergence of pairings on bounded sets proves continuity of both Fourier maps there too.

Now let $R\in SO(n)$ be a rotation matrix in the special orthogonal group and write ${\rho }_{R}T=T\circ R$. A change of variables with unit Jacobian givesThe rotation equivariance of the Fourier transform on tempered distributions follows by duality:Therefore $\widehat{{\rho }_{R}T}={\rho }_R\hat{T}$. Applying this identity and the inverse Fourier transform gives the two directions:This is precisely preservation of radial tempered distributions. For $n>1$, $SO(n)$ acts transitively on spheres, so an invariant smooth function is an ordinary radial function.

Because $0<\alpha <n$, $|x{|}^{-\alpha }$ is locally integrable at the origin and has only polynomial growth at infinity, so it defines a tempered distribution. It is a homogeneous distribution of degree $-\alpha$ and is radial. Its Fourier transform is therefore radial and homogeneous of degree $\alpha -n$, so it must have the form $c_{\alpha }|\lambda {|}^{\alpha -n}$. In particular, $\boxed{\beta =n-\alpha }$.

To determine the constant, use the stated Gamma integral representation, Fubini's theorem, and the Fourier transform of a Gaussian:The change of variables formula $s=|\lambda {|}^2/(4\tau )$ then givesThis is precisely the Fourier transform of the Riesz kernel.

Seek a radial function $u=u(r)$ on the unit ball. The equation and regularity at the origin giveHence $u^{\prime }(r)=(r/n{)}^{1/(p-1)}$, and the zero Dirichlet boundary condition givesThe power $|x{|}^q$ is twice differentiable at the origin exactly when $q\ge 2$. Here $q=p/(p-1)$, so $u\in C^2$ at the origin exactly when $1<p\le 2$. Under the assumption $2<p<n$, this weak solution is never $C^2$ at the origin, exhibiting the regularity loss caused by degenerate ellipticity.

Zero Dirichlet form of the Kac collision operator forces invariance under every pair rotation; strong continuity upgrades almost every angle to every angle. Coordinate-plane rotations generate the special orthogonal group, whose action is transitive on each sphere for $N\ge 2$. Averaging over normalized Haar measure therefore identifies the invariant $L^2$ functions with radial functions. Conversely radial functions are fixed by every pair rotation.

For a radial function in three dimensions, $(\Delta -1)y=g(x)$ becomes $y^{\prime \prime }+2y^{\prime }/x-y=g(x)$. Setting $u=xy$ gives $u^{\prime \prime }-u=xg(x)$, reducing it to a constant-coefficient linear ordinary differential equation. If $g=e^{-3x}/x^2$, a decaying particular solution obtained by variation of parameters isThe small-argument expansion of the exponential integral givesThe undetermined decaying homogeneous coefficient $\alpha$ is fixed by a boundary condition or by a matched asymptotic expansion.

For $n>1$, a tempered distribution $T$ on ${\mathbb{R}}^n$ is radial when $T\circ R=T$ for every $R\in SO(n)$, where $⟨T\circ R,\phi ⟩=⟨T,\phi \circ R^t⟩$. This agrees with the usual radial function definition for regular smooth distributions, because the special orthogonal group acts transitively on spheres.