Statement

In three dimensions, the satisfies, for a test function ff,

RiRjf(x)=p.v.R33zizjδijz24πz5f(xz)dzδij3f(x).R_iR_jf(x)=\operatorname{p.v.}\int_{\mathbb R^3} \frac{3z_iz_j-\delta_{ij}|z|^2}{4\pi|z|^5}f(x-z)\,dz -\frac{\delta_{ij}}3 f(x).

Here δij\delta_{ij} is the Kronecker delta. In particular its off-diagonal kernel obeys Kij(z)Cz3|K_{ij}(z)|\le C|z|^{-3}.

Distributional Hessian

Differentiate Γ(z)=1/(4πz)\Gamma(z)=1/(4\pi|z|) twice outside zero. The displayed kernel results. Its angular mean is zero because the spherical average of zizj/z2z_iz_j/|z|^2 is δij/3\delta_{ij}/3, making the principal value converge for smooth inputs. Integration by parts across a small sphere contributes δijδ0/3-\delta_{ij}\delta_0/3 to ijΓ\partial_i\partial_j\Gamma. Fourier transformation yields the symbol ξiξj/ξ2-\xi_i\xi_j/|\xi|^2, proving the formula. Summing i=ji=j gives iRi2=I\sum_iR_i^2=-I, which also checks the local term's sign.

Operator bounds

The off-diagonal kernel and its first derivatives have the Calderón–Zygmund bounds. The Fourier symbol has absolute value at most one, giving L2L^2 boundedness. The Calderón–Zygmund theorem therefore applies for 1<p<1<p<\infty.