Gegenbauer polynomials and positive deniteness

Christian Berg

University of Copenhagen, Denmark

Orsay, Paris November 27, 2015

Based on joint work with Emilio Porcu, University Federico Santa Maria Valparaíso, Chile

Christian Berg Gegenbauer polynomials Overview

1. Presentation of the problem and main results 2. Reminder about Gegenbauer polynomials 3. The connection to 4. Positive denite functions on locally compact groups G d d 5. The class P(S , G) of positive denite functions on S × G 6. The class P(S∞, G) and its representations 7. Some indications of proof 8. Applications to some homogeneous spaces

Christian Berg Gegenbauer polynomials Presentation of the problem

In Geostatistics one examines measurements depending on the location on the earth and on time. This leads to Random Fields of 2 stochastic variables Z(ξ, u) indexed by (ξ, u) belonging to S × R, 2 where S the 2-dimensional sphereis a model for the earth and R is a model for time. If the variables are real-valued, one considers a basic probability space (Ω, F, P), where all the random variables Z(ξ, u) are dened as measurable mappings from Ω to R. The covariance of two stochastic variables X , Y is by denition

cov(X , Y ) = E((X − E(X ))(Y − E(Y )).

For n variables X1,..., Xn the covariance matrix

cov X X n [ ( k , l )]k,l=1 is symmetric and positive semi-denite.

Christian Berg Gegenbauer polynomials Isotropic and stationary covariance kernels

One is interested in isotropic and stationary random elds 2 Z(ξ, u), (ξ, u) ∈ S × R, i.e., the situation where there exists a continuous function f :[−1, 1] × R → R such that the covariance kernel is given as

2 cov(Z(ξ, u), Z(η, v)) = f (ξ · η, v − u), ξ, η ∈ S , u, v ∈ R.

Here ξ · η = cos(θ(ξ, η)) is the scalar product equal to cosine of the length of the geodesic arc (=angle) between ξ and η. We require with other words that the covariance kernel only depends on the geodesic distance between the points on the sphere and on the time dierence.

Christian Berg Gegenbauer polynomials First main result

2 We shall characterize the class P(S , R) of continuous functions f :[−1, 1] × R → R which are positive denite in the following sense: 2 For any n ∈ N and any (ξ1, u1),..., (ξn, un) ∈ S × R the matrix f u u n [ (ξk · ξl , l − k )]k,l=1 is symmetric and positive semi-denite. Theorem (B-Porcu 2015) 2 The functions f ∈ P(S , R) are precisely the functions

∞ ∞ X X f (x, u) = ϕn(u)Pn(x), ϕn(0) < ∞, n=0 n=0

where (ϕn) is a sequence of real-valued continuous positive denite functions on R and Pn are the on [−1, 1] normalized as Pn(1) = 1. The series is uniformly convergent.

Christian Berg Gegenbauer polynomials Generalizations

This Theorem can be generalized in various ways: 2 d The sphere S can be replaced by S , d = 1, 2,.... The additive group R can be replaced by any locally compact group G. 2 ∞ The sphere S can be replaced by the Hilbert sphere S . d We shall characterize the set P(S , G) of continuous functions f :[−1, 1] × G → C such that the kernel

f (ξ · η, u−1v)

d 2 is positive denite on (S × G) . Here d = 1, 2,..., ∞. If G = {e} is the trivial group we get a classical Theorem of Schoenberg from 1942 about positive denite functions on spheres.

Christian Berg Gegenbauer polynomials Gegenbauer polynomials

To formulate these generalizations we need to recall: (λ) The Gegenbauer polynomials Cn for λ > 0 are given by the

∞ 2 −λ X (λ) n (1 − 2xr + r ) = Cn (x)r , |r| < 1, x ∈ C. (1) n=0 For λ > 0, we have the classical orthogonality relation:

Z 1 n 2 21−2λ 1 x2 λ−1/2C (λ) x C (λ) x dx πΓ( + λ) (2) ( − ) n ( ) m ( ) = 2 δm,n. −1 Γ (λ)(n + λ)n!

Christian Berg Gegenbauer polynomials

For λ = 0 we use the generating function

∞ 1 − xr X (0) n = Cn (x)r , |r| < 1, x ∈ . (3) 1 − 2xr + r 2 C n=0 It is well-known that

(0) Cn (x) = Tn(x) = cos(n arccos x), n = 0, 1,...

are the Chebyshev polynomials of the rst kind.

Z 1  π if n 0 1 x2 −1/2T x T x dx 2 δm,n > (4) ( − ) n( ) m( ) = if n 0 −1 πδm,n = ,

(λ) Warning: λ → Cn (x) is discontinuous at λ = 0.

Christian Berg Gegenbauer polynomials More on Gegenbauer polynomials

Putting x = 1 in the generating functions yields

(λ) Cn (1) = (2λ)n/n!, λ > 0, Tn(1) = 1.

Recall that for a ∈ C

(a)n = a(a + 1) ··· (a + n − 1), n ≥ 1, (a)0 = 1.

It is of fundamental importance that

(λ) (λ) |Cn (x)| ≤ Cn (1), x ∈ [−1, 1], λ ≥ 0.

The special value λ = (d − 1)/2 is relevant for the d-dimensional sphere d d+1 S = {x ∈ R | ||x|| = 1}, d ∈ N (d−1)/2 because of the relation of Cn to spherical harmonics.

Christian Berg Gegenbauer polynomials Spherical harmonics

d d A spherical harmonic of degree n for S is the restriction to S of a d+1 real-valued harmonic homogeneous polynomial in R of degree n. d Hn(d) = {spherical harmonics of degree n} ⊂ C(S )

d is a nite dimensional subspace of the continuous functions on S . We have

(d)n−1 Nn(d) := dim Hn(d) = (2n + d − 1), n ≥ 1, N0(d) = 1. n!

The surface measure of the sphere is denoted ωd , and it is of total mass 2π(d+1)/2 ||ωd || = . Γ((d + 1)/2)

The spaces Hn(d) are mutual orthogonal subspaces of the Hilbert 2 d space L (S , ωd ), which they generate.

Christian Berg Gegenbauer polynomials Ultraspherical polynomials

2 d This means that any F ∈ L (S , ωd ) has an orthogonal expansion ∞ ∞ X 2 X 2 F = Sn, Sn ∈ Hn(d), ||F ||2 = ||Sn||2, n=0 n=0 2 d where the rst series converges in L (S , ωd ), and the second series is Parseval's equation. Here Sn is the orthogonal projection of F onto Hn(d) given as

Nn(d) Z Sn(ξ) = cn(d, ξ · η)F (η) dωd (η), d ||ωd || S where ((d−1)/2) ((d−1)/2) cn(d, x) = Cn (x)/Cn (1). These polynomials are called ultraspherical polynomials or d-dimensional Legendre polynomials. The 2-dimensional Legendre polynomials are the classical Legendre polynomials previously denoted Pn. Christian Berg Gegenbauer polynomials Orthogonality relation for ultraspherical polynomials

Specializing the orthogonality relation for the Gegenbauer polynomials to λ = (d − 1)/2:

Z 1 2 d/2−1 ||ωd || (1 − x ) cn(d, x)cm(d, x) dx = δm,n. −1 ||ωd−1||Nn(d)

(Dene ||ω0|| = 2). Note that |cn(d, x)| ≤ 1, x ∈ [−1, 1].

Christian Berg Gegenbauer polynomials Schoenberg's Theorem from 1942

d Let P(S ) denote the class of continuous functions f :[−1, 1] → R d such that for any n ∈ N and for any ξ1, . . . , ξn ∈ S the n × n symmetric matrix f n [ (ξk · ξl )]k,l=1 is positive semi-denite. Theorem (Schoenberg 1942)

d A function f :[−1, 1] → R belongs to the class P(S ) if and only if

∞ X f (x) = bn,d cn(d, x), x ∈ [−1, 1], n=0 for a non-negative summable sequence b ∞ given as ( n,d )n=0 N d Z 1 ||ωd−1|| n( ) 2 d/2−1 bn,d = f (x)cn(d, x)(1 − x ) dx. ||ωd || −1

Christian Berg Gegenbauer polynomials Positive denite functions on groups

Consider an arbitrary locally compact group G, where we use the multiplicative notation, and in particular the neutral element of G is denoted e . In the representation theory of these groups the following functions play an crucial role. A continuous function f : G → C is called positive denite if for any n ∈ N and any u1,..., un ∈ G the n × n-matrix f u−1u n [ ( k l )]k,l=1 is hermitian and positive semi-denite. By P(G) we denote the set of continuous positive denite functions on G.

Christian Berg Gegenbauer polynomials A generalization

d We shall characterize the set P(S , G) of continuous functions f :[−1, 1] × G → C such that the kernel

−1 d f (ξ · η, u v), ξ, η ∈ S , u, v ∈ G (5)

is positive denite in the sense that for any n ∈ N and any d (ξ1, u1),... (ξn, un) ∈ S × G the n × n-matrix f u−1u n (6) (ξk · ξl ), k l ) k,l=1 is hermitian and positive semi-denite. d Note that for G = {e} we can identify P(S , G) with P(Sd ).

Christian Berg Gegenbauer polynomials Simple properties

Proposition

d (i) For f1, f2 ∈ P(S , G) and r ≥ 0 we have rf1, f1 + f2, and d f1 · f2 ∈ P(S , G). d (ii) For a net of functions (fi )i∈I from P(S , G) converging pointwise to a continuous function f :[−1, 1] × G → C, we d have f ∈ P(S , G). d (iii) For f ∈ P(S , G) we have f (·, e) ∈ P(Sd ) and f (1, ·) ∈ P(G). d (iv) For f ∈ P(Sd ) and g ∈ P(G) we have f ⊗ g ∈ P(S , G), where f ⊗ g(x, u) := f (x)g(u) for (x, u) ∈ [−1, 1] × G. In d particular we have f ⊗ 1G ∈ P(S , G) and f 7→ f ⊗ 1G is an d embedding of P(Sd ) into P(S , G).

Christian Berg Gegenbauer polynomials d Characterization of the class P(S , G)

Theorem (B-Porcu 2015)

Let d ∈ N and let f :[−1, 1] × G → C be a continuous function. d Then f belongs to P(S , G) if and only if there exists a sequence P ϕn,d ∈ P(G) with ϕn,d (e) < ∞ such that

∞ X f (x, u) = ϕn,d (u)cn(d, x), n=0 and the above expansion is uniformly convergent for (x, u) ∈ [−1, 1] × G. We have

N d Z 1 n( )||ωd−1|| 2 d/2−1 ϕn,d (u) = f (x, u)cn(d, x)(1 − x ) dx. ||ωd || −1

Christian Berg Gegenbauer polynomials d Relation between the classes P(S , G), d = 1, 2,... Note that 1 2 P(S , G) ⊃ P(S , G) ⊃ · · · d d−1 d−1 The inclusion P(S , G) ⊆ P(S , G) is easy, since S can be d considered as the equator of S . That the inclusion is strict is more subtle. The intersection ∞ \ d P(S , G) d=1 ∞ can be identied with the set P(S , G) of continuous functions f :[−1, 1] × G → C such that for all n f u−1v n [ (ξk · ξl , k l )]k,l=1 is hermitean and positive semi-denite for (ξk , uk ), k = 1,..., n ∞ from S × G, where ∞ ∞ X 2 2 S = {(xk ) | xk = 1} ⊂ ` . k 1 Christian Berg = Gegenbauer polynomials Schoenberg's second Theorem

When G = {e} Schoenberg proved in 1942: Theorem A function f 1 1 belongs to ∞ ∞ d if and :[− , ] → R P(S ) = ∩d=1P(S ) only if ∞ X n f (x) = bnx n=0

for a non-negative summable sequence bn. The convergence is uniform on [−1, 1].

Christian Berg Gegenbauer polynomials A characterization of P(S∞, G)

Theorem (B-Porcu 2015)

Let G denote a locally compact group and let f :[−1, 1] × G → C ∞ be a continuous function. Then f belongs to P(S , G) if and only if there exists a sequence ϕn ∈ P(G) with P ϕn(e) < ∞ such that

∞ X n f (x, u) = ϕn(u)x , n=0 and the above expansion is uniformly convergent for (x, u) ∈ [−1, 1] × G.

Christian Berg Gegenbauer polynomials Schoenberg coecient functions

d k For f ∈ P(S , G) we know that also f ∈ P(S , G) for k = 1, 2,..., d and therefore we have d expansions

∞ X f (x, u) = ϕn,k (u)cn(k, x), (x, u) ∈ [−1, 1] × G, n=0 where k = 1, 2,..., d. We call ϕn,k (u) the k-Schoenberg coecient functions of d f ∈ P(S , G). In the case where G = {e} the k-Schoenberg coecient functions are non-negative constants and they are just called k-Schoenberg coecients.

Christian Berg Gegenbauer polynomials More about Schoenberg coecient functions

There is a simple relation between these k-Schoenberg coecients / coecient functions:

d+2 d Suppose f ∈ P(S , G) ⊂ P(S , G). Then for u ∈ G, n ≥ 0 (n + d − 1)(n + d) (n + 1)(n + 2) ϕn d 2(u) = ϕn d (u)− ϕn 2 d (u). , + d(2n + d − 1) , d(2n + d + 3) + ,

In the case G = {e} these relations were found by Gneiting (2013) and extended to general G in B-Porcu(2015) (available on the ArXive).

Christian Berg Gegenbauer polynomials A result of J. Ziegel, 2014, and its generalization

Theorem (Johanna Ziegel, 2014)

d+2 Let d ∈ N and suppose that f ∈ P(S ). Then f is continuously dierentiable in the open interval (−1, 1) and there exist d f1, f2 ∈ P(S ) such that

f1(x) − f2(x) f 0(x) = , −1 < x < 1. 1 − x2 Theorem (B-Porcu 2015)

d+2 Let d ∈ N and suppose that f ∈ P(S , G). Then f (x, u) is continuously dierentiable with respect to x in ] − 1, 1[ and

∂f (x, u) f1(x, u) − f2(x, u) = , (x, u) ∈] − 1, 1[×G ∂x 1 − x2

for functions f f d G . In particular ∂f (x,u) is continuous 1, 2 ∈ P(S , ) ∂x on ] − 1, 1[×G. Christian Berg Gegenbauer polynomials Some remarks

The result of Ziegel is the analogue of an old result of Schoenberg about radial positive denite functions (Ann. Math. 1938): If f :[0, ∞[→ R is a continuous functions such that f (||x||) is n positive denite on R , then f has a continuous derivative of order [(n − 1)/2] on (0, ∞). The above extension of Ziegel's result to general G is needed in our ∞ extension of Schoenberg's result to P(S , G). ∞ If f ∈ P(S , G) we have a sequence of expansions, d = 1, 2,... ∞ ∞ X X n f (x, u) = ϕn,d (u)cn(d, x) = ϕn(u)x n=0 n=0 valid for (x, u) ∈ [−1, 1] × G. Here ϕn,d (u), ϕn(u) ∈ P(G). As part of our proof we obtain

lim ϕn,d (u) = ϕn(u) d→∞ for each n ∈ N0, u ∈ G. Christian Berg Gegenbauer polynomials Some indications of proof

Lemma

d Any f ∈ P(S , G) satises

f (x, u−1) = f (x, u), |f (x, u)| ≤ f (1, e), (x, u) ∈ [−1, 1] × G.

Lemma

Let K ⊂ G be a non-empty compact set, and let δ > 0 and an open neighbourhood U of e ∈ G be given. Then there exists a d partition of S × K in nitely many non-empty disjoint Borel sets, say Mj , j = 1,..., r, such that each Mj has the property

−1 (ξ, u), (η, v) ∈ Mj =⇒ θ(ξ, η) < δ, u v ∈ U.

Christian Berg Gegenbauer polynomials The crucial Lemma

Lemma

For a continuous function f :[−1, 1] × G → C the following are equivalent: d (i) f ∈ P(S , G). (ii) f is bounded and for any complex Radon measure µ on d S × G of compact support we have Z Z f (cos θ(ξ, η), u−1v) dµ(ξ, u) dµ(η, v) ≥ 0. d d S ×G S ×G Having established the Lemma, the next idea is to apply (ii) to the measure µ = ωd ⊗ σ, where σ is an arbitrary complex Radon measure of compact support on G.

Christian Berg Gegenbauer polynomials Application to some homogeneous spaces

In a recent manuscript Guella, Menegatto and Peron prove characterization results for isotropic positive denite kernels on d d0 0 d d0 S × S for d, d ∈ N ∪ {∞}. They consider the set P(S , S ) of 2 continuous functions f :[−1, 1] → R such that the kernel

d d0 K((ξ, ζ), (η, χ) = f (ξ · η, ζ · χ), ξ, η ∈ S , ζ, χ ∈ S ,

is positive denite in the sense that for any n ∈ N and any d d 0 (ξ1, ζ1),..., (ξn, ζn) ∈ S × S the matrix K n [ ((ξk , ζk ), (ξl , ζl ))]k,l=1 is positive semi-denite. They prove the following:

Christian Berg Gegenbauer polynomials Theorem of Gyuella, Menegatto, Peron, 2015

Theorem 0 2 Let d, d ∈ N and let f :[−1, 1] → R be a continuous function. d d0 Then f ∈ P(S , S ) if and only if

∞ X 0 f (x, y) = bfn,mcn(d, x)cm(d , y), x, y ∈ [−1, 1], n,m=0

P where bfn,m ≥ 0 such that bfn,m < ∞. The above expansion is uniformly convergent, and we have

0 Nn(d)σd−1 Nm(d )σd0−1 bfn,m = × σd σd0 Z 1 Z 1 0 2 d/2−1 2 d0/2−1 f (x, y)cn(d, x)cm(d , y)(1 − x ) (1 − y ) dx dy. −1 −1

Christian Berg Gegenbauer polynomials Reduction to our results

d0 The idea of proof is to consider S as the homogeneous space O(d 0 + 1)/O(d 0), where O(d 0 + 1) is the compact group of d0+1 0 orthogonal transformations in R and O(d ) is identied with the subgroup of O(d 0 + 1) which xes the point d 0+1 ε1 = (1, 0,..., 0) ∈ S . It is elementary to see that the formula of the Theorem denes a d d0 function f ∈ P(S , S ). d d0 Let us next consider f ∈ P(S , S ) and dene 0 F :[−1, 1] × O(d + 1) → R by

0 F (x, A) = f (x, Aε1 · ε1), x ∈ [−1, 1], A ∈ O(d + 1).

d 0 Then F ∈ P(S , O(d + 1)) because

−1 0 F (x, B A) = f (x, Aε1 · Bε1), A, B ∈ O(d + 1).

Christian Berg Gegenbauer polynomials Continuation of the proof 1

By our main Theorem

∞ X 0 F (x, A) = ϕn,d (A)cn(d, x), x ∈ [−1, 1], A ∈ O(d + 1), n=0 and N d Z 1 n( )σd−1 2 d/2−1 ϕn,d (A) = f (x, Aε1 · ε1)cn(d, x)(1 − x ) dx σd −1 belongs to P(O(d 0 + 1)). 0 The function ϕn,d is bi-invariant under O(d ), i.e.,

0 0 ϕn,d (KAL) = ϕn,d (A), A ∈ O(d + 1), K, L ∈ O(d ).

This is simply because f (x, KALε1 · ε1) = f (x, Aε1 · ε1).

Christian Berg Gegenbauer polynomials Continuation of the proof 2

0 The mapping A 7→ Aε1 is a continuous surjection of O(d + 1) onto d0 S , and it induces a homeomorphism of the homogeneous space 0 0 d0 O(d + 1)/O(d ) onto S . It is easy to see that as a bi-invariant function, ϕn,d has the form

ϕn,d (A) = gn,d (Aε1 · ε1) for a uniquely determined continuous function g 1 1 . n,d :[− , ] → R We have in addition g d0 , because for d0 n,d ∈ P(S ) ξ1, . . . , ξn ∈ S 0 there exist A1,..., An ∈ O(d + 1) such that ξj = Aj ε1, j = 1,..., n, hence

−1 −1 gn,d (ξk · ξl ) = gn,d (Al Ak ε1 · ε1) = ϕn,d (Al Ak ). It is now easy to nish the proof.

Christian Berg Gegenbauer polynomials Some references

C. Berg, E. Porcu, From Schoenberg coecients to Schoenberg functions, Preprint submitted to ArXiv. J. C. Guella, V. A. Menegatto and A. P. Peron, An extension of a theorem of Schoenberg to products of spheres arXiv:1503.08174. Thank you for your attention

Christian Berg Gegenbauer polynomials