À quoi sert le calculateur d'harmoniques sphériques ?
Les harmoniques sphériques \(Y_{n}^{m}(\theta,\phi)\) constituent la partie angulaire des solutions de l'équation de Laplace en coordonnées sphériques. Elles forment une base orthonormée complète à la surface d'une sphère et interviennent dans toute la physique et l'ingénierie : orbitales atomiques en mécanique quantique, développements multipolaires en électromagnétisme et en gravitation, fond diffus cosmologique en cosmologie, modèles géomagnétiques en géophysique ou encore modèles d'éclairage en synthèse d'images. Ce calculateur évalue l'harmonique à valeurs complexes pour un degré n, un ordre m, un angle polaire (zénithal) \(\theta\) et un angle azimutal \(\phi\) choisis.
Mode d'emploi
Choisissez d'abord une convention : la convention A (Wolfram/Mathematica) emploie le facteur de phase \(e^{i m \phi}\) ; la convention B (Maple) utilise \(e^{i m (\phi+\pi)}\), qui en diffère d'un facteur \((-1)^{m}\). Saisissez un degré entier \(n \ge 0\) et un ordre entier m tel que \(-n \le m \le n\). Indiquez l'angle zénithal \(\theta\) et l'azimut \(\phi\) dans l'unité d'angle retenue (degrés ou radians — la même unité s'applique aux deux). Réglez le nombre de chiffres significatifs à afficher, puis lisez les parties réelle et imaginaire, ainsi que le module et la phase.
La formule
Posons \(x = \cos\theta\). La normalisation vaut $$N = \sqrt{\frac{2n+1}{4\pi}\cdot\frac{(n-m)!}{(n+m)!}}.$$ La fonction de Legendre associée \(P_{n}^{m}(x)\) se construit par récurrence à partir de $$P_{m}^{m}(x) = (-1)^{m}(2m-1)!!\,(1-x^{2})^{m/2}.$$ On a alors, en convention A : $$Y = N \cdot P_{n}^{m}(x) \cdot \bigl(\cos(m\phi) + i\cdot\sin(m\phi)\bigr).$$ Les ordres négatifs s'obtiennent via $$P_{n}^{-m}(x) = (-1)^{m}\frac{(n-m)!}{(n+m)!}\,P_{n}^{m}(x).$$ Aux pôles (\(x = \pm 1\)), seul \(m = 0\) donne une valeur non nulle.
Exemple détaillé
Valeurs par défaut, convention A : \(n=2\), \(m=1\), \(\theta=45°\), \(\phi=0\). On a alors \(x = \cos 45° = 0{,}7071068\), \(\frac{2n+1}{4\pi} = \frac{5}{4\pi} = 0{,}3978874\), \(\frac{(n-m)!}{(n+m)!} = \frac{1}{6}\), donc $$N = \sqrt{0{,}0663146} = 0{,}257516.$$ \(P_{2}^{1}(x) = -3x\sqrt{1-x^{2}} = -3\cdot 0{,}5 = -1{,}5\). Avec \(\phi=0\), la phase vaut 1, d'où $$Y = 0{,}257516 \cdot (-1{,}5) = -0{,}386274 + 0\,i,$$ de module \(0{,}386274\).
FAQ
Pourquoi le résultat est-il complexe ? Le facteur \(e^{i m \phi}\) rend Y complexe, sauf lorsque \(m\phi\) est un multiple de \(\pi\) (par exemple \(\phi=0\)), cas où la partie imaginaire s'annule.
Pourquoi mon manuel de référence affiche-t-il +0,386274 ? Certaines sources omettent la phase de Condon–Shortley \((-1)^{m}\). La convention Wolfram (type A) l'inclut à l'intérieur de \(P_{n}^{m}\), ce qui donne la valeur négative.
Que se passe-t-il aux pôles ? Pour \(\theta=0\) ou \(180°\), le facteur \((1-x^{2})^{m/2}\) s'annule : \(Y = 0\) pour tout \(m \ne 0\), et seul \(m=0\) subsiste.
