En mathématiques, les harmoniques sphériques sont des fonctions harmoniques particulières, c'est-à-dire des fonctions dont le laplacien est nul. Les harmoniques sphériques sont particulièrement utiles pour résoudre des problèmes invariants par rotation, car elles sont les vecteurs propres de certains opérateurs liés aux rotations.

Les polynômes harmoniques P(x,y,z) de degré l forment un espace vectoriel de dimension 2 l + 1, et peuvent s'exprimer en coordonnées sphériques (r, θ, φ) comme des combinaisons linéaires des (2 l + 1) fonctions :

rlYl,m(θ,φ),  avec −l≤m≤+l.

Les coordonnées sphériques (r,θ,φ) sont, respectivement, la distance au centre de la sphère, la colatitude et la longitude.

Tout polynôme homogène est entièrement déterminé par sa restriction à la sphère unité S2.

Définition — Les fonctions sur la sphère obtenues par restriction de polynômes homogènes harmoniques sont des harmoniques sphériques.

C'est pourquoi la partie radiale de l'équation de Laplace, différente selon le problème étudié n'apparaît pas ici.

Les harmoniques sphériques sont utilisées en physique mathématique, dès qu'intervient la notion d'orientation (anisotropie) et donc de rotation (groupe de symétrie orthogonal SO(3)) et que le laplacien entre en jeu :

Elles sont également utilisées pour des problématiques de rendu 3D, notamment en réalité augmentée afin d'encoder des données d'éclairage ambiant.

Résolution de l'équation de Laplace

modifier

On cherche les fonctions Yl,m(θ,φ) sous la forme d'un produit de deux fonctions d'une seule variable :

Yl,m(θ,φ)=kPl,m(cos⁡θ) e+imφ

où k est une constante, qui sera fixée ultérieurement par la normalisation. L'équation aux valeurs propres devient une équation différentielle linéaire d'ordre deux pour la fonction Pl,m(cos θ) :

−1sin⁡θd dθ(sin⁡θdPl,m(cos⁡θ)dθ)+m2sin2θPl,m(cos⁡θ)=El,mPl,m(cos⁡θ)
Représentations des premières harmoniques sphériques "réelles" (combinaisons linéaires des Yℓ,m de même ℓ). Les parties en bleu correspondent aux valeurs négatives, celles en jaune aux valeurs positives des harmoniques.

On fait le changement de variable : θ↦x=cos⁡θ qui conduit à l'équation différentielle généralisée de Legendre :

−d dx[(1−x2)dPl,m(x)dx]+m2(1−x2)Pl,m(x)=El,mPl,m(x)

Les valeurs propres de cette équation sont indépendantes de m :

El,m=l(l+1) 

Les fonctions propres Pl,m(x) sont les polynômes associés de Legendre. Ils se construisent à partir des polynômes de Legendre Pl(x) qui sont les fonctions propres de l'équation différentielle ordinaire de Legendre, correspondant au cas m = 0 :

−d dx[(1−x2)dPl(x)dx]=l(l+1)Pl(x)

On a la formule génératrice d'Olinde Rodrigues :

Pl(x)=12ll!dl dxl[x2−1]l

On construit alors les fonctions propres Pl,m(x) par la formule :

Pl,m(x)=(−1)m[1−x2]m/2dmPl(x)dxm

soit explicitement :

Pl,m(x)=(−1)m2ll![1−x2]m/2dl+m dxl+m[x2−1]l

Remarque : il suffit en pratique de calculer les fonctions Pl,m(x) pour m ≥ 0, car il existe une relation simple entre Pl,m(x) et Pl, –m(x) :

Pl,−m(x)=(−1)m(l−m)!(l+m)!Pl,m(x)

Expression des harmoniques sphériques

modifier

On obtient alors l'expression inscrite plus bas. Une manière simple de retenir cette expression est la suivante :

Yl,0=Pl(cos⁡θ)⋅2l+14π,

où Pl(x) est le polynôme de Legendre de degré l.

On obtient ensuite :

J+Yl,m=(l2−m2)+(l−m)⋅Yl,m+1

où

J+=eiϕ(∂∂θ+itan⁡θ⋅∂∂ϕ)

est l'opérateur « d'échelle montante ».

Pour m négatif, Yl,m=(−1)m⋅Yl,−m∗

Souvent cette base se note |lm⟩ :

toute fonction sur la sphère S2 pourra donc s'écrire :

f(θ,ϕ)=fl,m⋅|lm⟩

(en convention de sommation d'Einstein), les coefficients complexes f l,m jouant le rôle de composantes de f dans la base des |lm⟩ (on dit parfois coefficients de Fourier généralisés).

En chimie ou en géophysique, il arrive qu'on préfère utiliser les harmoniques sphériques « réelles » et des coefficients de Fourier réels.

Expression mathématique

modifier

Les harmoniques sphériques formant une base[pourquoi ?][réf. souhaitée] orthogonale sur la sphère unité, toute fonction continue[réf. nécessaire] f(θ,φ) se décompose en[précision nécessaire] une série d'harmoniques sphériques :

f(θ,φ)=∑l=0+∞∑m=−l+lClm⋅Ylm(θ,φ)

où l et m sont des indices entiers, Cm
l
est un coefficient constant et souvent en mathématiques prend le nom de coefficient de Fourier généralisé relativement à cette base.

Le développement en harmoniques sphériques est l'équivalent, appliqué aux fonctions angulaires, du développement en séries de Fourier pour les fonctions périodiques.

Ym
l
est la partie réelle d'une fonction complexe Ym
l

Ylm(θ,φ)=Re⁡(Ylm_(θ,φ))

Ym
l
est appelée « fonction associée de Legendre » et est définie par

Ylm_(θ,φ)=2⋅(l−m)!(l+m)!⋅Plm(cos⁡θ)⋅eimφ

où i est l'imaginaire et Pm
l
est le polynôme associé de Legendre :

Plm(X)=(−1)m2l⋅l!⋅(1−X2)m/2⋅∂m+l∂Xm+l[(X2−1)l]

On a donc

Ylm(θ,φ)=2⋅(l−m)!(l+m)!⋅Plm(cos⁡θ)⋅cos⁡(mφ)

On a par exemple :

  • P00(cos⁡θ)=1 (Y0
    0
    est isotrope) ;
  • P10(cos⁡θ)=cos⁡θ ;
  • P11(cos⁡θ)=−sin⁡θ ;
  • P31(cos⁡θ)=32⋅sin⁡θ⋅(−5⋅cos2θ+1) ;

Les fonctions Ym
l
(θ, φ)
présentent de plus en plus de symétries au fur et à mesure que l croît (sauf lorsque l = 0, puisque Y0
0
est une fonction constante et décrit donc une sphère).

Polynômes de Legendre

modifier

Pour les harmoniques circulaires, on utilise des polynômes Pl de la fonction cosinus :

Yl(θ)=Pl(cos⁡θ)

Les polynômes Pl utilisés sont les polynômes de Legendre :

Pl(X)=12l⋅l!⋅dldXl[(X2−1)l] (formule de Rodrigues, mathématicien français)

On obtient :

  • P0(cos⁡θ)=1  (fonction isotrope) ;
  • P1(cos⁡θ)=cos⁡θ  ;
  • P2(cos⁡θ)=12(3cos2θ−1) ;
  • P3(cos⁡θ)=12(5cos3θ−3cos⁡θ) ;

Harmoniques sphériques normalisées

modifier

Base orthonormale des harmoniques sphériques

modifier

Parmi les 2l +1 fonctions, l'habitude a été prise de sélectionner une base orthonormale sur la sphère 𝕊2 munie de la mesure

dμ=14πsin⁡θ dθ dϕ,

soit le produit scalaire (hermitien en fait) :

⟨f1∣f2⟩=14π∬S2f1∗f2sin⁡θ dθ dϕ

Les harmoniques sphériques sont les solutions de l'équation aux valeurs propres[1] :

−ΔYl,m(θ,φ)=l(l+1)Yl,m(θ,φ)

où l'opérateur laplacien s'écrit en coordonnées sphériques sur la sphère de rayon unité J2 :

Δf(θ,φ)=defJ2f=1sin⁡θ∂ ∂θ(sin⁡θ∂f∂θ)+1sin2θ∂2f∂φ2

Elles sont fonctions propres de l'opérateur J3=−i∂∂ϕ :

J3Yl,m=m⋅Yl,m

Celles-ci, une fois normées sur la sphère sont alors notées usuellement Yl,m(θ , φ), où les angles (θ , φ) sont les coordonnées sphériques sur la sphère de rayon unité, et l et m sont deux nombres entiers tels que 0 ≤ l et –l ≤ m ≤ +l

Normalisation

modifier

Les harmoniques sphériques constituent une base orthonormale de fonctions propres de l'opérateur laplacien sur la sphère de rayon unité S2 au sens où :

Elles sont orthogonales pour le produit scalaire suivant :

∬S2dΩ(θ,φ)Y‾l′,m′(θ,φ)Yl,m(θ,φ)=δl,l′δm,m′

Dans cette formule, dΩ(θ, φ) représente l'angle solide élémentaire :

dΩ(θ,φ)=sin⁡θ dθ dφ

Toute fonction f(θ , φ) suffisamment régulière admet un développement en série :

f(θ,φ)=∑l=0+∞∑m=−l+lal,mYl,m(θ,φ)

où les coefficients complexes al,m se calculent par :

al,m=∬S2dΩ(θ,φ)Y‾l,m(θ,φ)f(θ,φ)

Expression des harmoniques sphériques normalisées

modifier

Les harmoniques sphériques généralisées sont définies sur la sphère S3. La normalisation des harmoniques sphériques conduit à l'expression finale :

Yl,m(θ,φ)=(−1)12(m+|m|)(2l+1)4π(l−|m|)!(l+|m|)!Pl,|m|(cos⁡θ)eimφ

Forme « réelle » des harmoniques sphériques

modifier

Si m ≠ 0 les harmoniques sphériques Yℓ,m ont des valeurs complexes. Il est cependant possible, pour une valeur donnée de ℓ de définir des combinaisons linéaires des Yℓ,m qui soient réelles, tout en constituant toujours une base normalisée sur la sphère unité.

Il suffit pour cela de prendre les combinaisons linéaires suivantes :

Y~ℓm={i2(Yℓm−(−1)mYℓ−m)si m<0,Yℓ0si m=0,12(Yℓ−m+(−1)mYℓm)si m>0.={i2(Yℓ−|m|−(−1)mYℓ|m|)si m<0,Yℓ0si m=0,12(Yℓ−|m|+(−1)mYℓ|m|)si m>0.

Il est facile de vérifier que ces expressions sont bien normalisées à l'unité. Ces relations s'inversent sans difficulté pour donner :

Yℓm={12(Y~ℓ|m|−iY~ℓ,−|m|)si m<0,Y~ℓ0si m=0,(−1)m2(Y~ℓ|m|+iY~ℓ,−|m|)si m>0.

En substituant les expressions précédentes des harmoniques sphériques, on obtient les expressions générales suivantes :

Y~ℓm={2(2ℓ+1)4π(ℓ−|m|)!(ℓ+|m|)!Pℓ|m|(cos⁡θ)sin⁡|m|φsi m<0,(2ℓ+1)4πPℓ0(cos⁡θ)si m=0,2(2ℓ+1)4π(ℓ−m)!(ℓ+m)!Pℓm(cos⁡θ)cos⁡mφsi m>0.

Ces fonctions sont utilisées fréquemment en chimie quantique pour représenter les parties angulaires des différentes orbitales atomiques associées aux différents électrons du cortège électronique des atomes.

Représentations graphiques

modifier

Représentation sphérique

modifier

Si l'on utilise la représentation sphérique

ρ=ρ0+ρ1⋅Ylm(θ,φ)

alors la surface représentatrice est une sphère bosselée ; les bosses correspondent aux parties où Ym
l
est positif, les creux aux parties où Ym
l
est négatif. Lorsque θ et φ décrivent l'intervalle [0 ; 2π[, Ym
l
(θ, φ)
s'annule selon l cercles :

  • m cercles suivant un méridien, une iso-longitude (intersection entre un plan contenant Oz et la sphère) ;
  • l – m cercles suivant un parallèle, une iso-latitude (intersection entre un plan parallèle à Oxy et la sphère).

Le paramètre l est appelé le « degré », m est appelé l'« ordre azimutal ». Entre les cercles d'annulation, la fonction est alternativement positive ou négative.

Ci-dessous sont représentées quatre coupes de l'harmonique sphérique Y2
3
 :

Comme précédemment, on peut représenter la fonction par la courbe en coordonnées sphériques :

Y32
ρ=ρ0+ρ1⋅Y32(θ,φ)
les parties en blanc sont positives, en bleu négatives
ρ=|Y32(θ,φ)|2

Représentation en coupe

modifier

Les harmoniques sphériques peuvent être représentées de façon plus simple sans les ventres de vibration, en ne gardant que les nœuds, comme le montre le tableau suivant[2]. Ce sont les sphères de la figure du haut, projetées sur un plan vertical. On retrouve sur la dernière ligne les quatre sphères de la première figure ci-dessus où l = 3. Les quatre valeurs de m y varient de 0 à 3 en valeur absolue. Sur la figure ci-après, on distingue les valeurs négatives pour tenir compte de ce que la rotation peut se faire dans un sens ou dans l'autre pour m > 0. Pour montrer la concordance avec les harmoniques, leur plus simple expression est donnée sous chaque sphère.

On reconnaît les nombres quantiques secondaire l, correspondant aux sous-couches s, p, d, f et m, magnétique, de l'atome d'hydrogène. Le nombre quantique principal n n'apparaît pas car les modes radiaux sont différents selon le problème étudié, résonance acoustique, atome d'hydrogène ou autre.

Pour montrer la concordance avec la littérature, l’expression des harmoniques sphériques est donnée sous chaque sphère. Le nombre et la valeur des zéros des polynômes de Legendre associés, non normalisés, donne le nombre de parallèles et leur position sur l’axe vertical. L’exponentielle imaginaire exp(imϕ), de module unité, utilisée habituellement au lieu des sinus et cosinus, donne le nombre de méridiens. Les valeurs de l ≥ 4 ne s’observent que dans les états excités ou les atomes de Rydberg où la valeur habituelle de l est 50 et dont l'orbitale est représentée non par une sphère mais par un anneau[3].

Représentation cartésienne et polaire

modifier

On peut représenter les harmoniques circulaires de trois manières :

Trois premières harmoniques circulaires
Représentation cartésienne Représentations polaires (tracé manuel) Représentations polaires (tracé exact)
Y1
Y2
Y3

Autres harmoniques

modifier

Harmoniques circulaires

modifier

Dans le plan, la décomposition s'écrit :

f(θ)=∑l=0+∞Cl⋅Yl(θ)

Y0 est une fonction constante, la courbe représentatrice en coordonnées polaires r = Y0(θ) est donc un cercle de rayon r0.

Yl est une fonction invariante par une rotation d'un angle de 1/l+1 tour, c'est-à-dire que

Yl(θ+2πl+1)=Yl(θ)

on dit que Yl admet une symétrie d'ordre l + 1.

Harmoniques sphériques généralisées

modifier

Lorsque l'on considère l'orientation d'un objet dans l'espace, il faut faire appel à trois angles ; on utilise en général les angles d'Euler (ψ, θ, φ).

Considérons une fonction continue de l'orientation f(ψ, θ, φ) ; comme précédemment, cette fonction peut être décomposée en harmoniques sphériques généralisées

f(ψ,θ,φ)=∑l=0+∞∑m=−l+l∑n=−l+lClmn⋅Ylmn(ψ,θ,φ)

où Cmn
l
est une constante. La fonction Ymn
l
s'écrit :

Ylmn(ψ,θ,φ)=eimφ⋅Plmn(cos⁡θ)⋅einψ

Le polynôme Pmn
l
est le polynôme de Legendre généralisé

Plmn(X)=(−1)l−m⋅in−m2l⋅(l−m)!⋅[(l−m)!(l+n)!(l+m)!(l−n)!]1/2⋅(1−X)−n−m2⋅(1+X)−n+m2⋅∂l−n∂Xl−n[(1−X)l−m(1+X)l+m]

Quand X décrit l'intervalle [–1 ; 1], cette fonction Pmn
l
est soit réelle, soit imaginaire pure. Y00
0
(ψ, θ, φ)
est la fonction isotrope (symétrie sphérique).

D'après la loi de composition des rotations, on a :

Ylmn(ψ1+ψ2,θ1+θ2,φ1+φ2)=∑s=−l+lYlms(ψ1,θ1,φ1)⋅Ylsn(ψ2,θ2,φ2)

et en particulier

Plmn(cos⁡(θ1+θ2))=∑s=−l+lPlms(cos⁡θ1)⋅Plsn(cos⁡θ2)

On a de manière générale :

Plmn=Plnm=Pl−m−n

Par exemple pour l = 1 :

P1mn(cos⁡θ)
m n
-1 0 +1
-1 12(1+cos⁡θ) −i2sin⁡θ 12(cos⁡θ−1)
0 −i2sin⁡θ cos⁡θ −i2sin⁡θ
1 12(cos⁡θ−1) −i2sin⁡θ i2(1+cos⁡θ)

Pour l = 2 :

P2mn(cos⁡θ)
m n
-2 -1 0 +1 +2
-2 14(cos⁡θ+1)2 −i2sin⁡θ(cos⁡θ+1) −1232(1−cos2θ) −i2sin⁡θ(cos⁡θ−1) 14(cos⁡θ−1)2
-1 −i2sin⁡θ(cos⁡θ+1) 12(2cos2θ+cos⁡θ−1) −32isin⁡θcos⁡θ 12(2cos2θ−cos⁡θ−1) −i2sin⁡θ(cos⁡θ−1)
0 −1232(1−cos2θ) −32isin⁡θcos⁡θ 12(3cos2θ−1) −32isin⁡θcos⁡θ −1232(1−cos2θ)
1 −i2sin⁡θ(cos⁡θ−1) 12(2cos2θ−cos⁡θ−1) −32isin⁡θcos⁡θ 12(2cos2θ+cos⁡θ−1) −i2sin⁡θ(cos⁡θ+1)
2 14(cos⁡θ−1)2 −i2sin⁡θ(cos⁡θ−1) −1232(1−cos2θ) −i2sin⁡θ(cos⁡θ+1) 14(cos⁡θ+1)2

Utilisation en physique

modifier

Les harmoniques sphériques sont d'une grande utilité pour la représentation mathématique des orbitales électroniques en physique quantique. En effet les Ylm(θ,φ) représentent les fonction d'ondes des électrons dans le potentiel électrostatique créé par le noyau de l'atome. Dans ce cadre, les nombres l et m représentent les nombres quantiques azimutaux et magnétiques[4].

Notes et références

modifier
  1. ↑ On a introduit un signe moins pour avoir des valeurs propres positives. En effet, l'opérateur laplacien est un opérateur négatif au sens où, pour toute fonction ϕ lisse à support compact, on a :
    ∫ϕΔϕ=−∫‖gradϕ‖2
    Cette égalité se démontre en utilisant la relation Δ = div grad et en intégrant par parties.
  2. ↑ Bernard Schaeffer, Relativités et quanta clarifiés, Publibook, 2007
  3. ↑ Atomes circulaires : propriétés et préparation
  4. ↑ Claude Cohen-Tannoudji, Bernard Diu et Franck Laloë, Mécanique quantique, vol. 1, Hermann, coll. « Enseignement des sciences », (ISBN 978-2-7056-6074-1)

Voir aussi

modifier

Articles connexes

modifier

Liens externes

modifier

Bibliographie

modifier