다중극 전개

전자기학 계산물리 수치해석 마지막 수정: 2026-08-14 04:21:18

1. 개요[편집]

다중극 전개 Multipole expansion
대상 커널	1/\|r − r′\| (라플라스), e^ikR/R (헬름홀츠)
전개 기저	르장드르 다항식 → 구면조화 함수
차수별 이름	단극자 · 쌍극자 · 사중극자 · 팔중극자
수렴 조건	원천을 감싸는 구(브리유앵 구) 바깥
절단 오차	~ (r′/r)^p+1
가장 큰 후손	고속 다중극자법

다중극 전개(multipole expansion)는 국소적으로 뭉쳐 있는 원천 분포가 멀리서 만드는 장을, 원천의 모양 대신 몇 개의 숫자(다중극 모멘트)로 요약해 거리의 역멱급수로 펼치는 방법이다. 핵심 아이디어는 한 줄이다 — 충분히 멀리서 보면 어떤 복잡한 전하 덩어리도 점전하 하나로 보이고, 조금 더 정밀하게 보면 점전하 + 쌍극자로 보이고, 더 정밀하게 보면 여기에 사중극자가 붙는다.

수학적 실체는 커널 $1/|\mathbf{r}-\mathbf{r}'|$ 을 원천 좌표와 관측 좌표로 분리하는 것이다. 이 “변수 분리”가 다중극 전개의 진짜 값어치다. 분리가 되는 순간 원천 $N$ 개와 관측점 $M$ 개의 상호작용을 $N\times M$ 번 계산할 필요가 없어지고, 원천 쪽 합 한 번( $p$ 개 모멘트)과 관측 쪽 평가 한 번으로 끝난다. 고속 다중극자법이라는 20세기 대표 알고리즘 전체가 이 관측 하나에서 자란 것이다.¹

2. 커널 쪼개기[편집]

$r' = |\mathbf{r}'| < r = |\mathbf{r}|$ 이고 두 벡터의 사잇각을 $\gamma$ 라 하면, 생성함수 정의에 의해 르장드르 다항식이 그대로 나온다.

\frac{1}{|\mathbf{r}-\mathbf{r}'|} = \frac{1}{\sqrt{r^{2}-2rr'\cos\gamma+r'^{2}}} = \sum_{l=0}^{\infty}\frac{r'^{\,l}}{r^{\,l+1}}\,P_{l}(\cos\gamma)

문제는 $\gamma$ 가 두 점의 각도에 동시에 의존한다는 것 — 아직 변수 분리가 안 끝났다. 여기에 구면조화 함수 덧셈 정리

P_{l}(\cos\gamma) = \frac{4\pi}{2l+1}\sum_{m=-l}^{l} Y_{lm}^{*}(\theta',\varphi')\,Y_{lm}(\theta,\varphi)

를 넣으면 비로소 원천 좌표 $(r',\theta',\varphi')$ 와 관측 좌표 $(r,\theta,\varphi)$ 가 완전히 갈라진다. 전하밀도 $\rho$ 에 대해 적분하면 퍼텐셜은

\Phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\sum_{l=0}^{\infty}\sum_{m=-l}^{l} \frac{4\pi}{2l+1}\,\frac{q_{lm}}{r^{\,l+1}}\,Y_{lm}(\theta,\varphi), \qquad q_{lm} = \int \rho(\mathbf{r}')\,r'^{\,l}\,Y_{lm}^{*}(\theta',\varphi')\,d^3r'

가 된다. $q_{lm}$ 이 구면 다중극 모멘트이고, 차수 $l$ 마다 $2l+1$ 개씩 있다. 원천이 어떻게 생겼든 이 숫자들만 알면 바깥 장은 완전히 결정된다. 라플라스 방정식의 외부 해가 $r^{-(l+1)}Y_{lm}$ 으로 남김없이 전개된다는 사실의 물리 버전이라 보면 된다.

카르테시안 표기로 낮은 차수만 쓰면 익숙한 얼굴이 나온다.

\Phi \approx \frac{1}{4\pi\varepsilon_0}\left[ \frac{q}{r} + \frac{\mathbf{p}\cdot\hat{\mathbf{r}}}{r^{2}} + \frac{1}{2}\sum_{i,j}\frac{Q_{ij}\,\hat{r}_i \hat{r}_j}{r^{3}} + \cdots \right]

$q=\int\rho\,dV$ , $\mathbf{p}=\int\rho\,\mathbf{r}'\,dV$ , $Q_{ij}=\int\rho\,(3x'_i x'_j - r'^2\delta_{ij})\,dV$ . 여기서 $\mathbf{p}$ 가 전기 쌍극자 모멘트다. 사중극자 텐서를 무자취(traceless)로 정의하는 이유는 자취 부분이 $r^{-3}$ 장에 기여하지 않기 때문이고, 그래서 성분 9개 중 독립인 것은 대칭성·무자취를 빼면 5개 — 구면 표기의 $2l+1=5$ 와 정확히 맞는다. 차수가 올라가면 카르테시안 텐서는 성분이 폭발하는데( $3^l$ ) 실제 자유도는 $2l+1$ 뿐이라, $l\ge 3$ 부터는 사실상 아무도 카르테시안을 쓰지 않는다.

3. 원점을 옮기면 값이 바뀐다[편집]

다중극 모멘트는 좌표 원점에 의존한다. 원점을 $\mathbf{d}$ 만큼 옮기면

q' = q, \qquad \mathbf{p}' = \mathbf{p} - q\,\mathbf{d}

이고 사중극자는 $q$ 와 $\mathbf{p}$ 가 섞여 들어온다. 여기서 나오는 규칙이 하나 있다 — 0이 아닌 가장 낮은 차수의 모멘트만 원점에 무관하다. 전체 전하가 있는 계에서 “이 분자의 쌍극자 모멘트는 얼마인가”라는 질문은 원점을 말하지 않으면 무의미하고, 이온의 쌍극자 모멘트를 인용하는 논문이 종종 사고를 치는 지점이 여기다. 중성 분자( $q=0$ )에서 쌍극자가 잘 정의되는 것은 운이 좋아서가 아니라 정리다.²

이 좌표 의존성은 결점이 아니라 알고리즘의 자원이기도 하다. “원점을 옮기면 모멘트가 이렇게 변한다”는 규칙이 곧 고속 다중극자법의 M2M·L2L 이동 연산자이며, 트리 위에서 모멘트를 부모 상자로 접어 올리는 작업 전체가 이 변환의 반복이다.

4. 수렴 반경 — 가장 자주 무시되는 조건[편집]

전개는 무조건 수렴하지 않는다. 원점을 중심으로 모든 원천을 감싸는 최소 반경의 구를 $R_{\max}$ 라 하면, 외부 전개는 $r > R_{\max}$ 에서만 수렴한다. 지구물리에서는 이 구를 브리유앵 구라 부른다.

이게 왜 실전 문제인가. 지구 중력장을 구면조화로 전개하면 브리유앵 구는 적도 반지름을 갖는 구이고, 지구는 편평하므로 극지방 지표면은 브리유앵 구 안쪽에 있다. 즉 극 근처 지표에서 EGM2008 같은 고차 모형을 그냥 평가하면 수학적 수렴 보장이 없다. 실무에서는 그럼에도 잘 작동하지만, 그 “잘 작동함”은 정리가 아니라 경험이다. 길쭉한 소행성이나 혜성 핵의 중력장을 구면조화로 다루려다 발산을 만나는 사고는 실제로 일어나며, 그래서 다면체 중력 모형이나 질량 집중(mascon) 모형으로 갈아탄다.

반대로 원천보다 안쪽을 보고 싶으면 역할을 뒤집는다.

\frac{1}{|\mathbf{r}-\mathbf{r}'|} = \sum_{l}\frac{r^{\,l}}{r'^{\,l+1}}P_l(\cos\gamma) \qquad (r < r'_{\min})

이쪽이 국소 전개(local expansion, 정칙 전개)이고, 멀리 있는 원천 무리가 내 상자 안에 만드는 부드러운 장을 표현한다. 외부 전개(다중극)와 국소 전개(로컬)라는 두 종류가 있고 둘 사이를 잇는 변환이 M2L이라는 구조는, 결국 $r_<$ 와 $r_>$ 중 무엇을 분자에 두느냐의 차이일 뿐이다.

5. 파동 문제로의 확장[편집]

정전기만 다루는 이야기가 아니다. 헬름홀츠 방정식 커널도 게겐바우어 덧셈 정리로 같은 형태로 쪼개진다.

\frac{e^{ik|\mathbf{r}-\mathbf{r}'|}}{|\mathbf{r}-\mathbf{r}'|} = ik\sum_{l=0}^{\infty}(2l+1)\,j_l(kr_<)\,h_l^{(1)}(kr_>)\,P_l(\cos\gamma)

$r^{l}, r^{-(l+1)}$ 자리에 구면 베셀 함수 $j_l$ 과 제1종 한켈 함수 $h_l^{(1)}$ 이 들어앉았을 뿐 구조는 동일하다. 여기서 갈라져 나오는 것들이 많다.

미 산란의 계수 $a_n$ , $b_n$ 은 각각 전기·자기 다중극 진폭 그 자체다. 구가 작으면( $x\ll1$ ) $a_1$ 만 살아남는 것이 곧 레일리 산란의 쌍극자 근사이고, 크기가 커질수록 고차 다중극 공명이 순서대로 켜진다.
T-행렬법 은 입사 다중극 계수를 산란 다중극 계수로 보내는 선형 사상을 행렬 하나로 저장하는 방법이다. 즉 “다중극 전개로 표현된 산란체의 완전한 이력서”.
다중극 복사. 크기 $a$ 인 원천의 $l$ 차 복사 세기는 $(ka)^{2l}$ 로 억제된다. 안테나가 대부분 쌍극자 근사로 설명되는 이유이며, 원자에서 전기쌍극자 금지 전이가 사중극자 전이로 겨우 일어나 수명이 극단적으로 길어지는 이유이기도 하다. 중력파에서 사중극자가 최저 차수인 것도 같은 계보의 사실이다(질량 단극자는 보존, 쌍극자는 운동량 보존으로 복사하지 못한다).

6. 어디에 쓰나[편집]

분자 정전기. 힘장에서 원자에 점전하만 두는 대신 원자 중심에 쌍극자·사중극자까지 얹는 것이 다중극 힘장이다(스톤의 분산 다중극 해석, AMOEBA 계열). 방향성 있는 수소결합과 $\pi$ 스택을 훨씬 잘 재현하지만 분자동역학 비용이 몇 배로 뛴다. 반대로 점전하 모형은 다중극 전개를 0차에서 잘라 버린 근사인 셈이다.
주기계 정전기. 에발트 합산의 역공간 항은 결국 상자 전체의 다중극(특히 총쌍극자) 항을 어떻게 처리하느냐의 문제로 귀결되고, 표면항·틴폴(tinfoil) 경계 조건 논쟁이 여기서 나온다.
중력·측지. 지구 중력장은 $J_2 \approx 1.083\times10^{-3}$ 인 대칭 사중극자 항이 압도적이며, 이것 하나로 인공위성 궤도면의 승교점 세차가 설명된다. 궤도 역학에서 태양동기 궤도를 설계할 때 쓰는 것이 바로 이 $J_2$ 세차다.
적분방정식 가속. 경계요소법·모멘트법의 조밀 행렬에서 멀리 떨어진 블록은 저계수(low-rank)이고, 그 저계수성의 물리적 이유가 정확히 다중극 전개다. 계층 행렬·ACA 같은 대수적 기법은 같은 사실을 커널 공식 없이 수치적으로 재발견하는 방법이라 볼 수 있다.

7. 수치적 성질[편집]

절단 차수를 $p$ 로 두면 오차는 대략

\left|\Phi - \Phi_p\right| \lesssim \frac{A}{r-r'}\left(\frac{r'}{r}\right)^{p+1}

로 거리비의 기하급수로 줄어든다. 이 성질이 실무적으로 매우 좋다.

분리비 $r'/r$ 가 $1/2$ 이면 항 하나당 정확도가 약 0.3 자리씩 좋아진다. 유효숫자 $d$ 자리를 원하면 $p \approx 3.3d$ 정도가 필요하다는 얘기.
3차원 전개의 항 수는 $(p+1)^2$ 이므로 저장·평가 비용이 $p$ 의 제곱으로 는다. 이동 연산자를 정직하게 구현하면 $p^4$ 까지 오르고, 회전 기반 분해나 지수(평면파) 표현으로 $p^3$ · $p^2$ 까지 낮춘다.
저주파 파탄. 헬름홀츠 커널 전개는 $kr$ 이 작을 때 $h_l^{(1)}(kr_>)$ 가 지수적으로 커지고 $j_l(kr_<)$ 가 지수적으로 작아져, 유한 정밀도에서 곱이 자릿수 소실로 무너진다. 고주파에서는 반대로 $p \sim kD$ 로 항 수가 늘어난다. 그래서 파동 문제의 다중극 구현은 주파수 대역에 따라 정식화를 갈아 끼운다.
모멘트 계산 자체는 테일러 급수 절단과 달리 원천 분포 전체를 적분하므로, 국소적으로 뾰족한 분포에도 안정적이다. 대신 고차 모멘트는 $r'^l$ 가중이라 바깥쪽 원천에 극도로 민감해서, 상자 경계에 걸친 입자 하나가 고차 항을 오염시키는 일이 생긴다.³

8. 관련 문서[편집]

그린가드와 로클린이 1987년 논문에서 한 일은 새로운 물리를 발견한 게 아니라, 100년 넘게 물리학 교과서 4장에 있던 전개를 자료구조와 결합한 것이다. 수치해석에서 큰 도약이 종종 이런 식으로 온다 — 아무도 새로운 수학을 안 만들었는데 갑자기 복잡도가 한 차수 내려간다. ↩
이 함정의 고전적 사례가 주기계다. 주기 경계 아래에서는 “단위 셀의 쌍극자 모멘트”가 셀을 어떻게 자르느냐에 따라 달라져서, 결정의 자발 분극을 소박하게 정의하려는 시도가 통째로 실패한다. 이 문제를 베리 위상으로 해결한 것이 현대 분극 이론이고, 논문 제목에 “polarization is not what you think it is” 같은 표현이 붙던 시절이 실제로 있었다. ↩
그래서 다중극 전개는 “먼 곳은 쉽다”는 직관의 수학적 형식화다. 반대로 말하면 가까운 곳은 절대 안 쉽다는 뜻이기도 해서, 어떤 가속 알고리즘을 써도 근거리 상호작용은 결국 직접 계산해야 한다. 공짜 점심은 원천에서 충분히 멀 때만 나온다. ↩