Radial decay and spherical Fourier spectra
Reconstructed and self-checked by GPT-6 Astra (OpenAI), Ultra reasoning effort, October 2026. Earlier edition: GPT-6.1 Sol (OpenAI), Ultra reasoning effort. Original exposition: CC0. Programme foundations retain their stated licences.
The Fourier transform of a radial integrable function in three dimensions can be computed by averaging a plane wave on a sphere. The resulting sine integral also computes sphere measures. A derivative of a radial delta requires a different step: we construct it as a limit of smooth functions and retain the derivative of the polar volume density. This gives whole distributional identities and exact checks using moments and Gaussians.
We use and complex bilinear pairings, so . The Schwartz Fourier foundation, F1–F5, proves the seminorm estimates, Gaussian constant, Fourier operations and transpose. The integration foundation, §§15.0–15.1, 15.6 and 16, proves convergence, absolute Fubini, linear changes, planar polar integration and generating-class uniqueness. The scalar foundation, §§13.1–13.5 and 13.7–13.10, proves calculus, exponential and trigonometric series, and smooth bumps.
For surface geometry we use exactly the latitude-and-area proof in U035, Theorem 1.1, and the orthogonal invariance and cap estimate immediately preceding U023, Lemma 2.1. Their graph measure and flux inputs are proved in U011, the graph construction and Theorem 2.1. The full measurable polar integration formula is proved in the angular foundation, A4. These earlier programme proofs are supplied with the lesson; external sources are not substitutes for them.
Lemma 1.1. Let be measurable and suppose . Then , its Fourier transform is radial and continuous, and, for ,
At zero its value is .
Proof. The supplied latitude proof uses
with and . On a hemisphere the graph density is ; planar polar radius satisfies . Thus , first on coordinate rectangles and then on all Borel sets by generating-class uniqueness. The equator and seam are null in graph charts, and the poles are null by the cap bound. Orthogonal invariance allows the latitude axis to be . Consequently,
The integral equals when , and its absolute value is at most .
The full polar formula now gives
. Absolute Fubini inserts (1.2) into its Fourier integral and proves (1.1). The sine integral is absolutely convergent: on , ; on , . The original sphere integral is bounded by , so dominated convergence against supplies its value and continuity at zero. Dominated convergence in the original Euclidean integral gives continuity everywhere. Orthogonal substitution proves radiality.
This is also the whole tempered transform. Indeed a bounded function defines a tempered distribution by the Schwartz estimate F1. For a Schwartz test , the double integral defining the transpose has absolute integral at most . Fubini therefore identifies it with integration against (1.1) and its continuous extension.
Theorem 1.2. In three dimensions,
Proof. For , the exponential primitive and integration by parts give
The boundary terms vanish because for : the exponential series bounds below by any fixed power, leaving . The same bound proves absolute integrability. Set ; the imaginary part of is . Lemma 1.1 gives the formula at .
For later moment checks, repeated integration by parts gives, for every integer and ,
Here , and , with the just-verified endpoints. In particular the physical mass is , which agrees with the continuous value at zero.
Theorem 2.1. For Euclidean surface measure on , ,
with smooth value at zero.
Proof. Graph scaling, proved in A4, multiplies surface measure by . Formula (1.2) therefore computes the plane-wave integral on this sphere. Its value at zero is its mass . Differentiating a plane wave introduces only products of coordinates of , whose absolute values are bounded by powers of on the sphere. Finite surface mass and dominated difference quotients thus give every coordinate derivative of the transform, and give continuity of those derivatives. This proves smoothness at zero without a polar-coordinate differentiability assumption.
For completeness, the scalar expression agrees there with
.
After , this is a series in ; factorial bounds dominate every differentiated series on a compact set. Finally, absolute Fubini with bound proves
.
Hence the computed function represents the entire transpose.
A radial delta derivative differentiates the volume density
We first construct all the delta derivatives needed here. Choose a smooth function , supported in , with integral one, and put . Such a function is supplied by the scalar bump construction followed by division by its positive integral. For , an integer, and , consider the smooth compactly supported function
.
Radial construction. For a smooth test define, near ,
This is smooth: stays in a compact positive interval and each derivative of the integrand is bounded by finitely many derivatives of on a compact annulus. Differentiation under the finite surface integral follows by dominated convergence. Polar integration and the ordinary substitution give
The integrations by parts have no endpoint terms, since the support lies strictly inside the positive axis. Uniform continuity of near proves the limit and its independence of the chosen unit-integral bump. This defines the notation
The chain and product rules bound this functional by
.
It is consequently a distribution of order at most , is supported on the sphere, and is tempered because this bound is controlled by finitely many Schwartz seminorms. It also acts on any smooth function near the sphere. Two extensions agreeing there have identical derivatives on the sphere and hence the same action. This proves cutoff independence directly. No general pullback theorem is being assumed.
Theorem 3.1. The distribution has whole Fourier transform
Proof. Since , the radial construction at gives
The first summand comes from differentiating the volume density ; the second comes from differentiating the test function.
Apply this to a plane wave, which is permitted by the established action on smooth functions near the sphere. For , (1.2) yields
At , the bracket is , giving . The plane-wave pairing is smooth in by differentiation under the finite sphere integral in (3.2). Equivalently, the cosine series becomes a locally differentiable series in .
To verify the whole transform, insert into (3.2). The function and its first -derivatives are absolutely dominated by , whose integral is finite by F1. Absolute Fubini gives
Thus no distribution supported at the frequency origin is left undetermined.
Exercises
Exercise 1 (foundation). For compute in dimension three. Normalize it to a probability density , give , and verify the zero-frequency mass.
Exercise 2 (intermediate). Compute the whole transform of by differentiating the decay parameter. Justify differentiation and verify the result at zero by a direct physical mass integral.
Exercise 3 (intermediate). With from Exercise 1, compute the transform, mass, mean and covariance matrix of the convolution . Prove the product rule in this case by absolute Fubini, and obtain the covariance both by physical moments and by the second-order Fourier expansion.
Exercise 4 (foundation). For compare the distributions and sphere surface measure on . Compute the former's whole transform and its value at zero.
Exercise 5 (intermediate). Define . Write its complete action on a smooth test near the sphere, and compute its Fourier transform and its zero-frequency value. Check the sign by differentiating the ordinary radial delta with respect to .
Exercise 6 (advanced). Define . Derive its action on a smooth test, prove it is a compact distribution of order at most two, and compute its entire Fourier transform. Retain the derivatives of the radial volume density and check the value at zero.
Exercise 7 (intermediate). Let . Compute its whole transform, total mass and leading term at frequency zero. Confirm the leading term by its signed second moments.
Exercise 8 (advanced). For and , evaluate by the physical pullback and independently by a one-dimensional Gaussian cosine integral. Identify the value of at which this pairing vanishes.
Solutions
Solution 1. Use in the proof of Theorem 1.2. Its inverse square has imaginary part , so
At zero the mass is . Hence
The density is nonnegative and has integral one, as required.
Solution 2. On a compact parameter interval , the derivative of is bounded by . Its Euclidean integral is . Dominated difference quotients therefore permit differentiating the whole Fourier integral. The minus derivative of the previous rational expression is
Its zero value equals the direct mass . Lemma 1.1 and its transpose argument ensure equality as distributions everywhere.
Solution 3. Define . Tonelli gives and finiteness almost everywhere; null exceptional values do not affect the density. Translation in the inner integral and absolute Fubini, with total bound , give
The radial moments are
Reflection in coordinate makes its first moment zero and makes mixed second moments zero. Coordinate permutations make the three diagonal moments equal, so .
The change is valid also for moments: , and . These yield integrable bounds for the first and second moment double integrals against . Expanding each coordinate product and using Fubini, the cross terms vanish by the zero means. Thus has mean zero and covariance .
Independently, finite second moments allow two frequency derivatives by dominated difference quotients, since each derivative is bounded by . Therefore
.
For , and is bounded near zero. Two applications of the fundamental theorem give . With the Hessian is , reproducing the covariance.
Solution 4. The radial construction, together with surface scaling, gives
Theorem 2.1 now gives the whole smooth transform
whose value at zero is . The normalization follows from , rather than from an assumed delta change-of-variables rule.
Solution 5. Apply the product and chain rules once to :
The construction has already proved compact support, order at most one and temperedness. Plane waves in this formula yield
For the sign check, put . Pairing the ordinary delta with gives ; differentiation in is . On each compact positive -interval all derivative bounds are uniform on a fixed annulus, so this is a valid distributional parameter derivative. Taking the minus -derivative of gives the formula above, also at zero. The full transpose is justified by the same finite-sphere Fubini bound as in Theorem 3.1.
Solution 6. For fixed , set . Writing gives
where . Integrating on the sphere gives , with positive sign. Its finite bound and support are explicit, and agree with the construction.
Taking two -derivatives of gives . Equivalently differentiate the ordinary delta plane-wave formula twice with :
At zero the value is , also equal to from a test equal to one near the sphere. Smoothness follows by differentiating the finite-sphere expression in , or by the series in . The result grows at most linearly. Substitution of into the displayed pairing and absolute Fubini against verify the entire tempered identity.
Solution 7. The masses and show that has mass zero. Linearity gives
Indeed the factorially convergent series gives , with the tail divided by bounded on . Symmetry makes the means and mixed moments zero, while scaling and equal diagonal moments give
One can check the quadratic coefficient directly from the plane-wave integral. The measure is even, so the sine term integrates to zero. Repeated scalar integration by parts in the Taylor remainder gives . The compact measure dominates the absolute value of every integral against , so its finite fourth moment bounds the integrated remainder by a constant times . The quadratic term is consequently
.
Solution 8. The real Gaussian proof F3 gives
On the unit sphere the radial derivative is times the Gaussian value. Inserting it in (3.2) gives
For a separate frequency calculation, take the even real part of the one-dimensional F3 formula:
For in a compact positive interval with lower endpoint , the differentiated integrand is bounded by , which is integrable by F3's polynomial Gaussian bound. Thus minus the -derivative gives
Use and multiply by , the product of the spherical polar factor and the coefficient in . The answer is again
.
All other factors are nonzero for , so the pairing vanishes exactly at .
References
- Michael E. Taylor, Fourier Analysis, Distributions, and Constant-Coefficient Linear PDE, §6, especially (6.2)–(6.5), p. 62: the freely accessible author notes supply the radial sphere-integral method. The author's symmetric Fourier normalization is converted here to the convention stated above. The lesson supplies its convergence, origin, density-derivative and exercise proofs explicitly.
- The supplied programme proofs linked in the introduction give the exact scalar, integration, graph, polar and Gaussian inputs. The original exposition here is CC0; each supplied foundation retains its own stated licence.