Küresel Harmonik Fonksiyon Hesaplayıcı nedir?
Küresel harmonikler \(Y_{n}^{m}(\theta,\phi)\), Laplace denkleminin küresel koordinatlardaki çözümlerinin açısal kısmını oluşturur. Bir kürenin yüzeyinde tam ve dik (ortonormal) bir taban meydana getirir ve fizik ile mühendisliğin pek çok alanında karşımıza çıkar: kuantum mekaniğinde atomik orbitaller, elektromanyetik ve kütleçekim alanında multipol açılımları, kozmolojide kozmik mikrodalga arka plan ışıması, jeofizikte yer manyetik alanı modelleri ve bilgisayar grafiklerinde aydınlatma modelleri. Bu hesaplayıcı, seçtiğiniz \(n\) derecesi, \(m\) mertebesi, kutupsal (zenit) açı \(\theta\) ve azimut açısı \(\phi\) için karmaşık değerli harmoniği hesaplar.
Nasıl kullanılır?
Önce bir konvansiyon seçin: tip A (Wolfram/Mathematica) \(e^{i m \phi}\) fazını kullanır; tip B (Maple) ise \((-1)^{m}\) kadar farklılık gösteren \(e^{i m (\phi+\pi)}\) fazını kullanır. Tamsayı bir \(n \ge 0\) derecesi ve \(-n \le m \le n\) koşulunu sağlayan tamsayı bir \(m\) mertebesi girin. Zenit açısı \(\theta\) ile azimut \(\phi\) değerlerini seçtiğiniz açı biriminde (derece veya radyan — aynı birim her ikisi için de geçerlidir) yazın. Kaç anlamlı basamak gösterilmesini istediğinizi belirleyin; ardından reel ve sanal kısımları, büyüklüğü ve fazı okuyun.
Formül
\(x = \cos\theta\) olsun. Normalizasyon $$N = \sqrt{\frac{2n+1}{4\pi}\cdot\frac{(n-m)!}{(n+m)!}}$$ şeklindedir. İlişkili Legendre fonksiyonu \(P_{n}^{m}(x)\), $$P_{m}^{m}(x) = (-1)^{m}(2m-1)!!\,(1-x^{2})^{m/2}$$ başlangıç değerinden yola çıkan bir özyineleme (recurrence) ile elde edilir. Buna göre tip A: $$Y = N \cdot P_{n}^{m}(x) \cdot (\cos(m\phi) + i\cdot\sin(m\phi)).$$ Negatif mertebeler için $$P_{n}^{-m}(x) = (-1)^{m}\frac{(n-m)!}{(n+m)!}\, P_{n}^{m}(x)$$ kullanılır. Kutuplarda (\(x = \pm 1\)) yalnızca \(m = 0\) sıfırdan farklıdır.
Çözümlü örnek
Varsayılan değerler, tip A: \(n=2\), \(m=1\), \(\theta=45^\circ\), \(\phi=0\). Bu durumda \(x = \cos 45^\circ = 0.7071068\), \(\frac{2n+1}{4\pi} = \frac{5}{4\pi} = 0.3978874\), \(\frac{(n-m)!}{(n+m)!} = \frac{1}{6}\) olur; dolayısıyla \(N = \sqrt{0.0663146} = 0.257516\). \(P_{2}^{1}(x) = -3x\sqrt{1-x^{2}} = -3\cdot 0.5 = -1.5\). \(\phi=0\) olduğunda faz \(1\) olur ve $$Y = 0.257516 \cdot (-1.5) = -0.386274 + 0\,i$$ çıkar; büyüklük \(0.386274\)'tür.
Sıkça sorulan sorular
Sonuç neden karmaşık? \(e^{i m \phi}\) çarpanı, \(m\phi\) çarpımı \(\pi\)'nin tam katı olmadıkça (örneğin \(\phi=0\)) Y'yi karmaşık yapar; bu durumda sanal kısım sıfırlanır.
Başvuru kitabımda neden +0.386274 yazıyor? Bazı kaynaklar Condon–Shortley \((-1)^{m}\) fazını dahil etmez. Wolfram/tip A konvansiyonu bu fazı \(P_{n}^{m}\) içine alır ve negatif değeri verir.
Kutuplarda ne olur? \(\theta=0\) veya \(180^\circ\) iken \((1-x^{2})^{m/2}\) çarpanı sıfırdır; bu yüzden \(m \ne 0\) olan her durumda \(Y = 0\) olur, yalnızca \(m=0\) ayakta kalır.
