Sfærisk harmonisk funksjon

Illustrasjon av hvordan verdiene til den sfærisk harmoniske funksjonen Mal:Nowrap for ℓ = 3 varierer på kuleflaten.

Sfærisk harmoniske funksjoner er spesielle funksjoner definert på overflaten til en kule. De kalles derfor ofte for kuleflatefunksjoner. Når Laplace-operatoren uttrykkes i kulekoordinater, er de egenfunksjoner til den vinkelavhengige delen av operatoren. På den måten spiller de en viktig rolle i løsning av mange partielle differensialligninger. Dette gjelder også for løsninger av Schrödinger-ligningen slik at funksjonene dermed bidrar til forståelse av egenskapene til atomer og molekyler.

Trigonometriske funksjoner er periodiske i vinkelkoordinaten på en sirkel. De kalles derfor alternativt for «sirkelfunksjoner». På tilsvarende vi er de sfærisk harmoniske funksjonene periodiske i de to koordinatene på en kuleflate. Betegnelsen «harmonisk» skyldes at de opptrer i løsningene av Laplace-ligningen når denne uttrykkes i kartesiske koordinater x, y og z. Løsningene av denne ligningen er da homogene polynom i disse tre variable og sies å være harmoniske. De sfærisk harmoniske funksjonene fremkommer når disse polynomene uttrykkes i kulekoordinater r, θ og φ som også omtales som «sfæriske koordinater».

Rent matematisk danner funksjonene grunnlaget for representasjoner av rotasjonsgruppen SO(3) i tre dimensjoner. Mer fysisk fremkommer de på tilsvarende vis ved kvantisering av dreieimpuls i kvantemekanikken.

Definisjon

De kartesiske koordinatene x, y og z  kan uttrykkes ved kulekoordinater $(r,\theta ,\varphi )$ som $x=r\sin \theta \cos \varphi ,$ $y=r\sin \theta \sin \varphi$ og $z=r\cos \theta .$ Laplace-operatoren i tre dimensjoner tar dermed formen

$\begin{array}{cc}{\mathrm{\nabla }}^2& =\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}+\frac{{\partial }^2}{\partial z^2}\\ & =\frac{{\partial }^2}{\partial r^2}+\frac{2}{r}\frac{\partial }{\partial r}+\frac{1}{r^2}[\frac{1}{\sin \theta }\frac{\partial }{\partial \theta }(\sin \theta \frac{\partial }{\partial \theta })+\frac{1}{{\sin }^2\theta }\frac{{\partial }^2}{\partial {\varphi }^2}]\end{array}$

Løsninger av Laplace-ligningen ${\mathrm{\nabla }}^{2}F(x,y,z)=0$ kan da skrives som et produkt $F(r,\theta ,\varphi )=R(r)Y(\theta ,\varphi )$ av en radiell del og en angulær del. Hvis denne siste delen skal være kontinuerlig og periodisk i begge sine koordinater, må den oppfylle den partielle differensialligningen

$[\frac{1}{\sin \theta }\frac{\partial }{\partial \theta }(\sin \theta \frac{\partial }{\partial \theta })+\frac{1}{{\sin }^2\theta }\frac{{\partial }^2}{\partial {\varphi }^2}]Y_{\ell m}(\theta ,\varphi )=-\ell (\ell +1)Y_{\ell m}(\theta ,\varphi )$

hvor ℓ er et heltall. Samtidig må indeksen m ligge i intervallet Mal:Nowrap og tar derfor Mal:Nowrap forskjellige verdier. Løsning av ligningen er de sfærisk harmoniske funksjonene.

Beregning

Grafisk fremstilling av realdelen til de første kuleflate-funksjonene. Rødt og grønt indikerer positive eller negative verdier, mens avstand fra origo angir størrelsen.

Fra formen til differensialligningen som definerer kuleflatefunksjonene, følger at dens løsninger må kunne skrives som et produkt av to funksjoner,

$Y_{\ell m}(\theta ,\varphi )={\mathrm{\Theta }}_{\ell m}(\theta ){\mathrm{\Phi }}_m(\varphi )$

Disse må igjen tilfredsstille hver sin egen differensialligning. Den første er

$\frac{d^2}{d{\varphi }^2}{\mathrm{\Phi }}_m+m^2{\mathrm{\Phi }}_m=0$

som har løsninger av formen ${\mathrm{\Phi }}_m(\varphi )=e^{im\varphi }$ eller reelle kombinasjoner av to slike. For at ${\mathrm{\Phi }}_m(\varphi )={\mathrm{\Phi }}_m(\varphi +2\pi ),$ må m  være et heltall. Samtidig må den andre funksjonen oppfylle

$\frac{1}{\sin \theta }\frac{d}{d\theta }(\sin \theta \frac{d{\mathrm{\Theta }}_{\ell m}}{d\theta })+[\ell (\ell +1)-\frac{m^2}{{\sin }^2\theta }]{\mathrm{\Theta }}_{\ell m}=0$

Denne differensialligningen definerer de assosierte Legendre-polynomene $P_{\ell }^m(\cos \theta ).$ På denne måten finner man at de sfæriske harmoniske funksjonene kan skrives som

$Y_{\ell m}(\theta ,\varphi )=\sqrt{\frac{(2\ell +1)(\ell -m)!}{4\pi (\ell +m)!}}{P}_{\ell }^m(\cos \theta )e^{im\varphi },-\ell \le m\le \ell$

når de normeres ved integralet

$\underset{0}{\overset{2\pi }{\int }}d\varphi \underset{0}{\overset{\pi }{\int }}d\theta \sin \theta Y_{lm}^{*}(\theta ,\varphi )Y_{l^{\prime }{m}^{\prime }}(\theta ,\varphi )={\delta }_{l{l}^{\prime }}{\delta }_{m{m}^{\prime }}$

over hele kuleflaten. Da funksjonene er komplekse, er her konvensjonen

$Y_{\ell m}^{*}(\theta ,\varphi )=(-1{)}^m{Y}_{\ell ,-m}(\theta ,\varphi )$

benyttet. I litteraturen forekommer også andre konvensjoner.

Eksempel

De første sfærisk harmonisk funksjoner fra ℓ = 0 til og med ℓ = 3:

Yℓm	ℓ = 0	ℓ = 1	ℓ = 2	ℓ = 3
m = −3				$\sqrt{\frac{35}{64\pi }}{\sin }^3\theta e^{-3\mathrm{i}\varphi }$
m = −2			$\sqrt{\frac{15}{32\pi }}{\sin }^2\theta e^{-2\mathrm{i}\varphi }$	$\sqrt{\frac{105}{32\pi }}{\sin }^2\theta \cos \theta e^{-2\mathrm{i}\varphi }$
m = −1		$\sqrt{\frac{3}{8\pi }}\sin \theta e^{-\mathrm{i}\varphi }$	$\sqrt{\frac{15}{8\pi }}\sin \theta \cos \theta e^{-\mathrm{i}\varphi }$	$\sqrt{\frac{21}{64\pi }}\sin \theta (5{\cos }^2\theta -1)e^{-\mathrm{i}\varphi }$
m = 0	$\sqrt{\frac{1}{4\pi }}$	$\sqrt{\frac{3}{4\pi }}\cos \theta$	$\sqrt{\frac{5}{16\pi }}(3{\cos }^2\theta -1)$	$\sqrt{\frac{7}{16\pi }}(5{\cos }^3\theta -3\cos \theta )$
m = 1		$-\sqrt{\frac{3}{8\pi }}\sin \theta e^{\mathrm{i}\varphi }$	$-\sqrt{\frac{15}{8\pi }}\sin \theta \cos \theta e^{\mathrm{i}\varphi }$	$-\sqrt{\frac{21}{64\pi }}\sin \theta (5{\cos }^2\theta -1)e^{\mathrm{i}\varphi }$
m = 2			$\sqrt{\frac{15}{32\pi }}{\sin }^2\theta e^{2\mathrm{i}\varphi }$	$\sqrt{\frac{105}{32\pi }}{\sin }^2\theta \cos \theta e^{2\mathrm{i}\varphi }$
m = 3				$-\sqrt{\frac{35}{64\pi }}{\sin }^3\theta e^{3\mathrm{i}\varphi }$

Noen egenskaper

På en kuleflate med koordinater $\theta ,\varphi$ danner de sfærisk harmoniske funksjoner et «komplett sett». Det betyr at en vilkårlig funksjon $f(\theta ,\varphi )$ kan utvikles i en uendelig rekke som

$f(\theta ,\varphi )=\sum _{\ell =0}^{\mathrm{\infty }}\sum _{m=-\ell }^{\ell }{C}_{\ell m}{Y}_{\ell m}(\theta ,\varphi )$

hvor koeffisientene C_ℓm kan finnes fra integralet

$C_{\ell m}=\underset{0}{\overset{2\pi }{\int }}d\varphi \underset{0}{\overset{\pi }{\int }}d\theta \sin \theta f(\theta ,\varphi )Y_{\ell m}^{*}(\theta ,\varphi )$

Dette tilsvarer å utvikle en funksjon på sirkelen som en uendelig Fourier-rekke med ledd som involverer de trigonometriske funksjonene sinus og cosinus.

Addisjonsteorem

En variant av den vanlige cosinussetningen i planet eksistrer også på en kuleflate med sfærisk geometri. Der kan man angi et punkt med enhetsvektoren n i retning punktet (θ, φ). For to forskjellige retninger n og n'  som danner vinkelen γ, gjelder da setningen

$𝐧\cdot {𝐧}^{\prime }=\cos \gamma =\cos \theta \cos {\theta }^{\prime }+\sin \theta \sin {\theta }^{\prime }\cos (\varphi -{\varphi }^{\prime })$

Ved å benytte at kuleflatefunksjonene utgjør et fullstendig sett, kan man da utlede den viktige forbindelsen

$P_{\ell }(𝐧\cdot {𝐧}^{\prime })=\frac{4\pi }{2\ell +1}\sum _{m=-\ell }^{\ell }{Y}_{\ell m}(𝐧)Y_{\ell m}^{*}({𝐧}^{\prime })$

med Legendre-polynom. Den kalles for addisjonsteoremet for sfærisk harmoniske funksjoner og har mange anvendelser.

Kvantemekanikk

En partikkel i posisjon r og med impuls p har i klassisk mekanikk dreieimpulsen L = r × p. I en kvantemekanisk beskrivelse blir impulsen erstattet med en operator p → -iħ ∇ der ħ  er den reduserte Planck-konstanten. På denne måten blir også dreieimpulsen en operator

$\widehat{𝐋}=-i\hslash 𝐫\times \mathrm{\nabla }$

Den kan kvantiseres ved bruk av komponentene ${\hat{L}}_{\pm }={\hat{L}}_x\pm i{\hat{L}}_y$ sammen med ${\hat{L}}_z=-i\hslash (x{\partial }_y-y{\partial }_x).$ Ved bruk av kulekoordinater tar disse tre operatorene formen

$\begin{array}{cc}{\hat{L}}_{\pm }& =\hslash e^{\pm i\varphi }(\pm \frac{\partial }{\partial \theta }+i\cot \theta \frac{\partial }{\partial \varphi })\\ {\hat{L}}_z& =-i\hslash \frac{\partial }{\partial \varphi }\end{array}$

Kvadratet av den totale dreieimpuls finnes da som

$\begin{array}{cc}{\widehat{𝐋}}^2& ={\hat{L}}_{-}{\hat{L}}_{+}+{\hat{L}}_z^2+\hslash {\hat{L}}_z\\ & =-{\hslash }^2[\frac{1}{\sin \theta }\frac{\partial }{\partial \theta }(\sin \theta \frac{\partial }{\partial \theta })+\frac{1}{{\sin }^2\theta }\frac{{\partial }^2}{\partial {\varphi }^2}]\end{array}$

Fra definisjonen av kuleflatefunksjonene ser man herav at de er egenfunksjoner til de to operatorene ${\widehat{𝐋}}^2$ og ${\hat{L}}_z,$

$\begin{array}{cc}{\widehat{𝐋}}^2{Y}_{\ell m}& ={\hslash }^2\ell (\ell +1)Y_{\ell m}\\ {\hat{L}}_z{Y}_{\ell m}& =\hslash m{Y}_{\ell m}\end{array}$

Indeksen ℓ er derfor et kvantetall som angir størrelsen tilden totale, kvantiserte dreieimpuls. Den blir derfor ofte omtalt som det «orbitale kvantetallet». Da indeksen m samtidig gir egenverdien til dreieimpulsen langs z-aksen som ofte sammenfaller med retningen til et ytre magnetfelt, kalles denne vanligvis for det «magnetiske kvantetallet».

Rotasjonsmatriser

En egentilstand til dreieimpulsoperatoren kan skrives som en ketvektor $|\ell ,m⟩.$ De sfærisk harmoniske funksjonene kan da betraktes som bølgefunksjoner

$Y_{\ell m}(\theta ,\varphi )=⟨\theta ,\varphi |\ell ,m⟩$

i en posisjonsbasis dannet av punkter på kuleflaten. En rotasjon av dette koordinatsystemet er nå gitt ved den kvantemekaniske rotasjonsoperatoren

$\hat{R}(\mathit{\alpha })=e^{-i\mathit{\alpha }\cdot \widehat{𝐋}/\hslash }$

hvor α = (α,β,γ) er de tre eulervinklene som definerer rotasjonen. Bølgefunksjonene til systemet vil da forandres og finnes fra

$\begin{array}{cc}{Y}_{\ell m}(\theta ,\varphi )& =⟨{\theta }^{\prime },{\varphi }^{\prime }|\hat{R}(\mathit{\alpha })|\ell ,m⟩\\ & =\sum _{m^{\prime }=-\ell }^{\ell }{Y}_{\ell m^{\prime }}({\theta }^{\prime },{\varphi }^{\prime })D_{m^{\prime }m}^{(\ell )}(\mathit{\alpha })\end{array}$

etter å ha innført rotasjonsmatrisene

$D_{m^{\prime }m}^{(\ell )}(\mathit{\alpha })=⟨\ell ,m^{\prime }|\hat{R}(\mathit{\alpha })|\ell ,m⟩$

til Eugene Wigner. De kan eksplisitt beregnes for hver verdi av det orbitale kvantetallet ℓ og har dimensjon (2ℓ +1)×(2ℓ + 1). For ℓ = 1 kan de finnes uten bruk av kvantemekanikk da de tilsvarer rotasjon av en tredimensjonal vektor.

Se også

• Laplace-ligning
• Legendre-polynom
• Kvantisert dreieimpuls

Eksterne lenker

• H. Haber, The Spherical Harmonics, Phys. 116C, UC Santa Cruz (2012).
• Professor M, The Spherical Harmonics, YouTube video
• R.G. Littlejohn, Orbital Angular Momentum and Spherical Harmonics, Physics 221A , UC Berkeley (2019).