In the setting of Probability Measures on Euclidean Space and Random Vectors: Standing Notation, with Differential Calculus and Convexity on Euclidean Open Sets: Standing Notation in force, fix a dimension q. Let gs for 0<s≤1 be the Gaussian kernels of The Gaussian Kernels on Euclidean Space: Scaling, Derivatives up to Order Three, the Convolution Identity, Moments and Tails; let gs∗μ be the Gaussian smoothing of μ∈P(Rq) at scale s, with Lebesgue measure λq, integrable and bounded as there; and let series of real numbers, their convergence, absolute convergence and sums be those of Series of Real Numbers. The third index of a third-order partial derivative is written l here, the letter k being the series index. For k∈N put
sk=4−kandwk=2−k(q+8),
the multiplicative inverses of the natural powers 4k and 2k(q+8), where 2=1+1 and 4=2+2; thus 0<sk≤4−1 and 0<2sk≤2−1 for every k. Then the following hold.
1. (The kernel)¶ For every z∈Rq the series ∑k=1∞wkg2sk(z) converges absolutely. The function
K:Rq→R,K(z)=k=1∑∞wkg2sk(z),
called the multiscale Gaussian kernel, is of class C3 on Rq and even, K(−z)=K(z); for all z∈Rq and i,j,l∈[q] its partial derivatives are the absolutely convergent sums
∂iK(z)=k=1∑∞wk∂ig2sk(z),∂j∂iK(z)=k=1∑∞wk∂j∂ig2sk(z),∂l∂j∂iK(z)=k=1∑∞wk∂l∂j∂ig2sk(z);
and there are nonnegative real numbers M0,M1,M2,M3, depending only on q, such that for all z∈Rq and i,j,l∈[q],
0≤K(z)≤M0,∣∂iK(z)∣≤M1,∣∂j∂iK(z)∣≤M2,∣∂l∂j∂iK(z)∣≤M3.
2. (The potential of a measure)¶ Let ν∈P(Rq). For every x∈Rq the function y↦K(x−y) is bounded and Borel, so that
(K∗ν)(x)=∫RqK(x−y)ν(dy)
is a real number. The function K∗ν:Rq→R, called the potential of ν, is of class C3 on Rq, with
∂i(K∗ν)(x)=∫Rq∂iK(x−y)ν(dy),∂j∂i(K∗ν)(x)=∫Rq∂j∂iK(x−y)ν(dy),∂l∂j∂i(K∗ν)(x)=∫Rq∂l∂j∂iK(x−y)ν(dy)
for all x and i,j,l∈[q], the integrands being bounded Borel functions of y; it satisfies 0≤K∗ν≤M0, ∣∂i(K∗ν)∣≤M1, ∣∂j∂i(K∗ν)∣≤M2 and ∣∂l∂j∂i(K∗ν)∣≤M3 on Rq; and for every x∈Rq the series below converges absolutely and
(K∗ν)(x)=k=1∑∞wk(g2sk∗ν)(x).
3. (The kernel pairing)¶ Let μ,ν∈P(Rq). The function z↦K(pr1(z)−pr2(z)) on Rq+q is bounded and Borel, with the coordinate projections and product measure of that clause, and the kernel pairing
K(μ,ν)=∫Rq+qK(pr1(z)−pr2(z))(μ⊠ν)(dz)
satisfies 0≤K(μ,ν)≤M0, K(μ,ν)=K(ν,μ), and
K(μ,ν)=∫RqK∗νdμ=k=1∑∞wk∫Rq(gsk∗μ)(gsk∗ν)dλq,
where each product (gsk∗μ)(gsk∗ν) is integrable and the series has nonnegative terms and converges.