이산 쌍극자 근사

1. 개요[편집]

이산 쌍극자 근사 Discrete Dipole Approximation (DDA)
제안	퍼셀 · 페니패커 (1973)
별칭	결합 쌍극자 방법(coupled dipole method)
대상	임의 형상·비균질 산란체
미지수	격자점마다 쌍극자 모멘트 3성분 → 3N
핵심 가속	블록 토플리츠 구조 + FFT, 반복 해법
유효성	|m|kd ≲ 0.5
대표 코드	DDSCAT, ADDA

이산 쌍극자 근사(Discrete Dipole Approximation, DDA)는 산란체를 격자 위에 놓인 유한 개의 편극 가능한 점쌍극자 배열로 치환하고, 각 쌍극자가 입사장과 다른 모든 쌍극자의 복사장을 함께 느낀다는 자기무모순(self-consistent) 조건을 풀어 산란 문제를 해결하는 방법이다. 1973년 퍼셀과 페니패커가 성간 먼지의 소광을 계산하려고 만들었고¹, 드레인과 플라타우의 DDSCAT 이후 사실상 비구형 산란 계산의 표준 도구 중 하나가 되었다.

미 산란이 완벽한 구에만 통하고 T-행렬법이 회전대칭체에서 진가를 발휘하는 데 비해, DDA의 팔림 포인트는 딱 하나다. 형상이 아무래도 상관없다. 검댕 응집체, 육각 얼음 기둥, 구멍 뚫린 나노입자, 굴절률이 내부에서 변하는 비균질체 — 격자에 채워 넣을 수만 있으면 그냥 돌아간다. 대가는 미지수 개수와 그에 비례하는 계산량이다.

2. 이산화 — 물체를 쌍극자로 바꾸기[편집]

부피 $V$ 의 물체를 격자 간격 $d$ 인 정육면체 격자로 채운다. 물체 내부에 들어오는 격자점 $N$ 개가 각각 하나의 쌍극자가 되고, $N d^3 = V$ 가 성립한다. 형상은 “이 격자점이 물체 안인가 밖인가”라는 불리언 배열로만 표현되므로, 형상 입력이 사실상 3차원 비트맵이다.

각 쌍극자에는 편극률 $\alpha_j$ 를 준다. 출발점은 정전기학의 클라우지우스-모소티 관계다. 유전율 $\varepsilon$ 인 물질을 간격 $d$ 의 입방 격자로 쪼개면(가우스 단위계)

$$\alpha^{\mathrm{CM}} = \frac{3 d^{3}}{4\pi}\,\frac{\varepsilon - 1}{\varepsilon + 2}$$

가 나온다. 이 값을 그대로 쓰면 답이 틀린다. 이유는 두 가지다.

• 복사 반작용(radiative reaction). 진동하는 쌍극자는 에너지를 복사하므로 편극률에 허수 보정 $-\tfrac{2}{3} i k^{3}$ 항이 붙어야 광학정리와 에너지 보존이 맞는다(드레인 1988).
• 격자 자체의 분산. 유한 간격 격자에 놓인 쌍극자 배열이 지지하는 파동의 위상속도는 연속 매질의 그것과 다르다. 이 차이를 $(kd)^2$ 차수까지 상쇄하도록 편극률을 손보는 것이 격자 분산관계 보정(Lattice Dispersion Relation, LDR)이다.

LDR 편극률은 드레인·굿맨(1993)이 유도한

$$\alpha^{\mathrm{LDR}} = \frac{\alpha^{\mathrm{CM}}} {1 + \dfrac{\alpha^{\mathrm{CM}}}{d^{3}}\left[(b_1 + m^{2} b_2 + m^{2} b_3 S)(kd)^{2} - \tfrac{2}{3} i (kd)^{3}\right]}$$

꼴이며, $b_1, b_2, b_3$ 는 고정 상수, $S$ 는 입사 방향과 편광에 의존하는 기하 인자다. 겉보기엔 잡스러운 보정 같지만 실제 효과는 크다 — 같은 격자 수에서 단면적 오차가 몇 배 줄어든다. “쌍극자만 촘촘히 박으면 되겠지”는 DDA에서 가장 흔한 초심자 함정이고, 편극률 처방이 그 절반을 해결한다.²

3. 자기무모순 방정식[편집]

$j$ 번째 쌍극자가 느끼는 장은 입사장 $\mathbf{E}_{\mathrm{inc},j} = \mathbf{E}_0 e^{i\mathbf{k}\cdot\mathbf{r}_j}$ 에 나머지 모든 쌍극자가 만드는 장을 더한 것이다. 쌍극자 모멘트를 $\mathbf{P}_j$ 라 하면

$$\mathbf{P}_j = \alpha_j\left(\mathbf{E}_{\mathrm{inc},j} - \sum_{k \neq j} \mathbf{A}_{jk}\mathbf{P}_k\right) \quad\Longleftrightarrow\quad \sum_{k=1}^{N}\mathbf{A}_{jk}\mathbf{P}_k = \mathbf{E}_{\mathrm{inc},j}$$

가 된다($\mathbf{A}_{jj} = \alpha_j^{-1}$). 비대각 블록은 자유공간 그린 함수의 쌍극자–쌍극자 형태, 즉 $r = |\mathbf{r}_j - \mathbf{r}_k|$, $\hat{\mathbf{r}}$ 을 단위벡터라 할 때

$$\mathbf{A}_{jk} = \frac{e^{ikr}}{r}\left[k^{2}\left(\hat{\mathbf{r}}\hat{\mathbf{r}} - \mathbf{I}\right) + \frac{ikr - 1}{r^{2}}\left(3\hat{\mathbf{r}}\hat{\mathbf{r}} - \mathbf{I}\right)\right]$$

인 $3\times3$ 텐서다. 근거리항 $r^{-3}$, 중간항 $r^{-2}$, 원거리 복사항 $r^{-1}$ 이 한 식에 다 들어 있다. 결국 풀 것은 $3N$ 개 복소 미지수의 꽉 찬 선형계다. 이것이 DDA의 전부이자, 동시에 DDA를 비싸게 만드는 원흉이다.

4. 왜 FFT 없이는 못 쓰나[편집]

$N$ 은 금방 커진다. 크기 파라미터 $x = 10$ 짜리 입자를 파장당 격자 10개로 채우면 $N \sim 10^5$, 미지수 $3N \sim 3\times10^5$. 이걸 정공법으로 LU 분해하면 $O(N^3)$ 연산에 $O(N^2)$ 메모리 — 행렬을 저장하는 것부터 불가능하다. 그래서 DDA는 두 가지 관측 위에 서 있다.

첫째, 반복 해법. 행렬을 만들지 않고 행렬–벡터 곱만 제공하면 크리로프 부분공간법으로 풀 수 있다. 다만 $\mathbf{A}$ 는 복소 대칭($\mathbf{A}^{T} = \mathbf{A}$)이되 에르미트가 아니므로 교과서의 켤레기울기법이 그대로 적용되지는 않는다. 실제 코드는 복소 대칭계 전용 변형(COCG), 안정화 쌍공액구배(BiCGSTAB), QMR 같은 CG 계열을 쓴다.³ 조건수는 $|m|$ 이 클수록 나빠지고, 금속처럼 $|m|$ 이 큰 입자에서 반복수가 폭증하는 것이 이 방법의 고질병이다.

둘째, 균일 격자의 공짜 구조. $\mathbf{A}_{jk}$ 가 오직 $\mathbf{r}_j - \mathbf{r}_k$ 에만 의존하므로, 물체를 감싸는 직육면체 격자 위에서 $\mathbf{A}$ 는 3중 블록 토플리츠 행렬이다. 토플리츠는 두 배 크기의 순환행렬에 끼워 넣을 수 있고, 순환행렬–벡터 곱은 고속 푸리에 변환으로 대각화된다. 결과적으로 행렬–벡터 곱이 $O(N)$ 저장·$O(N \log N)$ 연산으로 떨어진다. 물체 밖 격자점은 $\mathbf{P}=0$ 으로 채워 넣으면 그만이다.

이 조합(반복법 + 3차원 FFT)이 DDA를 실용 도구로 만든 결정적 한 방이다. 대신 격자를 균일 정육면체로 고정해야 한다는 제약이 따라온다 — 국소 세분화를 하는 순간 토플리츠 구조가 깨져 FFT를 못 쓴다.

5. 단면적 뽑기 — 광학정리[편집]

$\mathbf{P}_j$ 가 다 구해지면 관측량은 합으로 떨어진다. 소광 단면적은 광학정리(전방 산란 진폭의 허수부 = 소광)의 이산 버전이다.

$$C_{\mathrm{ext}} = \frac{4\pi k}{|\mathbf{E}_0|^{2}} \sum_{j=1}^{N} \mathrm{Im}\!\left(\mathbf{E}_{\mathrm{inc},j}^{*}\cdot\mathbf{P}_j\right)$$

흡수는 각 쌍극자의 소산에서 자기 복사분을 빼면 된다.

$$C_{\mathrm{abs}} = \frac{4\pi k}{|\mathbf{E}_0|^{2}} \sum_{j=1}^{N} \left\{\mathrm{Im}\!\left[\mathbf{P}_j\cdot(\alpha_j^{-1})^{*}\mathbf{P}_j^{*}\right] - \tfrac{2}{3}k^{3}\left|\mathbf{P}_j\right|^{2}\right\}$$

산란은 $C_{\mathrm{sca}} = C_{\mathrm{ext}} - C_{\mathrm{abs}}$ 로 얻거나, 원거리장을 방향별로 합산해 미분 산란 단면적과 산란 행렬(뮐러 행렬)을 직접 만든다. 후자는 편광 정보까지 다 나오므로 리모트센싱 응용에서 필수다. 뒷항 $-\tfrac{2}{3}k^3|\mathbf{P}_j|^2$ 을 빠뜨리면 흡수가 없는 유전체에서도 $C_{\mathrm{abs}} > 0$ 이 나오는, 눈에 잘 안 띄는 버그가 생긴다.

6. 유효성 기준과 한계[편집]

DDA는 근사이므로 “언제 믿어도 되나”가 늘 따라붙는다. 실무 기준은 사실상 하나로 요약된다.

$$|m|\,k d \lesssim 0.5$$

즉 매질 내 파장당 격자점이 대략 10개 이상이어야 한다는 뜻이다. 이 조건에서 단면적 오차가 대체로 수 % 수준으로 들어온다. 여기에 두 가지가 더 붙는다.

• $N$ 이 형상을 묘사할 만큼 충분해야 한다. 파장 조건을 만족해도 격자 30개로 얼음 결정의 모서리를 흉내 낼 수는 없다.
• 표면 근처의 계단 오차(staircasing). 정육면체 격자로 곡면을 근사하니 경계에서 오차가 크고, 이게 $|m|$ 이 클수록 심해진다. 공명 근처의 플라스몬 응답처럼 표면장에 민감한 양은 수렴이 유난히 느리다.

기법	형상	비용 지배 요인	강한 곳
미 산란	구만	급수 항수 ~ x	정확해, 검증 기준
T-행렬법	회전대칭 유리	전개 차수	배향 평균, 다중체
DDA	임의	격자점 수 N	비균질·복잡 형상
FDTD · 유한요소 전자기	임의	전 영역 격자	광대역, 비선형·시변

$|m|$ 이 큰 금속(가시광 영역의 은·금)에서는 정확도와 수렴 둘 다 나빠져 DDA의 아성이 흔들린다. 이 영역에서는 유한요소 전자기나 표면 적분방정식(모멘트법) 쪽이 유리한 경우가 많다.

7. 변형과 확장[편집]

기본 골격이 단순한 덕에 파생형이 많다. 대부분 “격자 위 쌍극자 + 그린 함수”라는 뼈대를 유지한 채 그린 함수나 편극률만 갈아 끼운다.

• 주기 DDA. 무한히 반복되는 배열(격자 표면, 메타물질 층)은 단위 셀 안의 쌍극자만 미지수로 두고, 상호작용 텐서를 격자합으로 대체한다. 격자합이 조건부 수렴이라 에발트 합 같은 가속이 필요하다.
• 기판 위 입자. 자유공간 그린 함수를 반무한 매질의 그린 함수(조머펠트 적분)로 바꾸면 표면에 얹힌 나노입자를 다룰 수 있다. 그린 함수 평가 비용이 확 오르고 토플리츠 구조도 면내 방향에만 남는다.
• 편극률 처방의 갈래. LDR 외에도 필터링된 결합 쌍극자(FCD), 적분형 그린 텐서(IGT), 정확한 부피 적분 기반 처방 등이 있고, $|m|$ 이 큰 영역에서 서로 성능이 갈린다. 어떤 처방을 썼는지 밝히지 않은 DDA 결과는 재현이 안 된다.
• 근접장·힘·토크. $\mathbf{P}_j$ 를 갖고 있으면 맥스웰 응력 텐서 적분 없이도 광학 힘과 토크를 쌍극자 합으로 뽑을 수 있다. 광학 집게 시뮬레이션에서 쓰인다.
• 구현 계보. 드레인·플라타우의 DDSCAT(포트란, MPI/OpenMP)이 표준 참조 구현이고, 유르킨·훅스트라의 ADDA가 병렬·GPU와 다양한 편극률 처방을 갖춘 현대적 대안이다.

8. 실무 체크리스트[편집]

DDA 결과를 논문이나 보고서에 넣기 전에 확인하는 항목은 사실상 정해져 있다.

• $|m|kd \le 0.5$ 를 만족하는가. $m$ 이 파장에 따라 변하므로 스펙트럼 전체에서 확인해야 한다. 짧은 파장 끝에서만 조용히 깨져 있는 경우가 흔하다.
• 격자를 1.5~2배 촘촘히 해서 단면적이 몇 % 움직이는가(격자 수렴 지수 검사와 같은 정신).
• 편극률 처방과 반복 해법의 수렴 판정 기준(잔차 노름)을 명시했는가. 잔차를 $10^{-3}$ 에서 끊으면 단면적은 맞아도 근접장은 틀린다.
• 구를 넣어 미 산란 해와 대조했는가.

9. 어디에 쓰나[편집]

• 천체·성간 먼지: 원래 태어난 목적. 비구형 실리케이트·탄소질 입자의 소광 곡선과 편광을 맞추는 데 쓴다.
• 대기 에어로졸과 얼음 구름: 육각 기둥·판, 검댕 응집체(프랙탈 차원 ~1.8)의 산란 위상함수를 만들어 복사 전달 방정식과 광학 두께 계산에 넘긴다. 구로 근사하면 복사 강제력 추정이 통째로 틀어지는 대표 사례.
• 플라스모닉 나노입자: 막대·별·다이머의 국소표면플라스몬 공명 파장과 근접장 증강을 예측한다. 다만 위에서 말한 정확도 제약을 반드시 확인할 것.
• 생체 세포·조직: 적혈구처럼 구가 아닌 세포의 산란 지문을 계산해 유세포 분석·광간섭 진단에 쓴다.
• 검증: 구에 대해 돌려 미 산란 해와 비교하는 것이 DDA 코드 신고식의 국룰이다. 여기서 3% 이내로 안 맞으면 격자를 의심하기 전에 편극률 처방부터 확인한다.

10. 관련 문서[편집]

• 미 산란 · T-행렬법 · 레일리 산란
• 맥스웰 방정식 · 그린 함수 · 광학 정리
• 고속 푸리에 변환 · 크리로프 부분공간법 · 희소행렬
• FDTD · 유한요소 전자기 · 모멘트법 · 고속 다중극자법
• 복사 전달 방정식 · 광학 두께 · 편광
• 메타물질 · 헬름홀츠 방정식

11. Footnotes[편집]

1. Purcell, E. M. & Pennypacker, C. R. (1973), ApJ 186, 705. 제목이 “Scattering and absorption of light by nonspherical dielectric grains”으로, 애초에 천문학 논문이다. 나노포토닉스에서 이 방법을 쓰는 사람 중 상당수는 자기가 쓰는 코드의 조상이 성간 먼지 계산이라는 걸 모른다. ↩

2. 격자를 두 배 촘촘히 하면 $N$ 이 8배, FFT 비용은 그 이상으로 늘어난다. 반면 편극률 처방을 CM에서 LDR로 바꾸는 것은 코드 몇 줄에 비용 0이다. 수치해석에서 “공짜 점심”에 가장 가까운 순간. ↩

3. 이 지점이 자주 헷갈린다. $\mathbf{A}$ 는 대칭이지만 복소수라서 에르미트가 아니고, 따라서 표준 CG의 수렴 이론(양정부호 에르미트)이 성립하지 않는다. “대칭이니까 CG 쓰면 되지”라고 했다가 잔차가 춤추는 걸 보고 나서야 매뉴얼을 읽게 되는 것이 이 바닥 통과의례. ↩