구면조화함수 계산기란?
구면조화함수 \(Y_{n}^{m}(\theta,\phi)\)는 구면좌표계에서 라플라스 방정식 해의 각도 부분에 해당합니다. 구면 위에서 완전한 정규직교 기저를 이루며, 물리학과 공학 전반에 걸쳐 등장합니다. 양자역학의 원자 오비탈, 전자기학과 중력 이론의 다극 전개, 우주론의 우주배경복사(CMB), 지구물리학의 지자기 모델, 컴퓨터 그래픽스의 조명 모델 등이 대표적인 예입니다. 이 계산기는 선택한 차수 \(n\), 차원 \(m\), 극각(천정각) \(\theta\), 방위각 \(\phi\)에 대해 복소수 값을 갖는 구면조화함수를 계산합니다.
사용 방법
먼저 규약을 선택하세요. type A(Wolfram/Mathematica)는 위상 인자 \(e^{i m \phi}\)를 사용하고, type B(Maple)는 \(e^{i m (\phi+\pi)}\)를 사용하는데, 이 둘은 \((-1)^{m}\)만큼 차이가 납니다. 정수인 차수 \(n \ge 0\)과 \(-n \le m \le n\)을 만족하는 정수 차원 \(m\)을 입력하세요. 천정각 \(\theta\)와 방위각 \(\phi\)는 선택한 각도 단위(도 또는 라디안 — 두 각도에 동일한 단위가 적용됩니다)로 입력합니다. 표시할 유효숫자 자릿수를 지정한 뒤, 실수부와 허수부, 그리고 크기와 위상을 확인하면 됩니다.
계산 공식
\(x = \cos\theta\)로 둡니다. 정규화 상수는 $$N = \sqrt{\frac{2n+1}{4\pi}\cdot\frac{(n-m)!}{(n+m)!}}$$입니다. 연관 르장드르 함수 \(P_{n}^{m}(x)\)는 $$P_{m}^{m}(x) = (-1)^{m}(2m-1)!!\,(1-x^{2})^{m/2}$$에서 출발하는 점화식으로 계산합니다. 그러면 type A의 경우 $$Y = N \cdot P_{n}^{m}(x) \cdot \left(\cos(m\phi) + i\cdot\sin(m\phi)\right)$$가 됩니다. 음수 차원에 대해서는 $$P_{n}^{-m}(x) = (-1)^{m}\frac{(n-m)!}{(n+m)!}\,P_{n}^{m}(x)$$를 사용합니다. 극점(\(x = \pm 1\))에서는 \(m = 0\)인 경우만 0이 아닌 값을 가집니다.
계산 예시
기본값, type A: \(n=2\), \(m=1\), \(\theta=45^\circ\), \(\phi=0\)인 경우를 살펴봅시다. \(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}\)이므로 \(N = \sqrt{0.0663146} = 0.257516\)입니다. \(P_{2}^{1}(x) = -3x\sqrt{1-x^{2}} = -3\cdot 0.5 = -1.5\)이고, \(\phi=0\)이므로 위상은 1이 되어 $$Y = 0.257516 \cdot (-1.5) = -0.386274 + 0\,i,$$ 크기는 \(0.386274\)가 됩니다.
자주 묻는 질문
결과가 왜 복소수인가요? \(e^{i m \phi}\) 인자 때문에 \(Y\)는 일반적으로 복소수입니다. 다만 \(m\phi\)가 \(\pi\)의 정수배일 때(예: \(\phi=0\))는 허수부가 0이 됩니다.
참고서에는 왜 +0.386274로 나와 있나요? 일부 자료는 콘던–쇼틀리(Condon–Shortley) 위상 인자 \((-1)^{m}\)을 생략합니다. Wolfram/type A 규약은 이 인자를 \(P_{n}^{m}\) 안에 포함하므로 음수 값이 나옵니다.
극점에서는 어떻게 되나요? \(\theta=0\) 또는 \(180^\circ\)에서는 \((1-x^{2})^{m/2}\) 인자가 0이 되므로, \(m \ne 0\)인 모든 경우에 \(Y = 0\)입니다. 오직 \(m=0\)인 항만 살아남습니다.
