변환광학

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

1. 개요[편집]

변환광학 Transformation Optics
출발점	맥스웰 방정식의 좌표 형태 불변성
핵심 공식	ε′ = Λ ε Λᵀ / det Λ (μ′도 동일)
정식화	펜드리·슈리그·스미스 (2006) / 레온하르트 (2006)
선구	그린리프·라사스·울만 (2003, 전기 임피던스 역문제)
대표 산물	원통 클로크, 카펫 클로크, 굴절률 분포형 소자
구조적 한계	특이 물성 · 강한 분산 → 협대역

“빛이 이렇게 휘었으면 좋겠다”를 물성표로 번역하는 환율 계산기.

변환광학(transformation optics)은 맥스웰 방정식이 좌표 변환에 대해 형태 불변이라는 사실을 이용해, 원하는 광선·파면 거동을 좌표계의 왜곡으로 먼저 그린 다음 그 왜곡을 유전율·투자율 텐서 분포로 번역하는 설계 방법론이다. 방정식의 성질이지 새로운 물리가 아니라는 점이 중요하다 — 변환광학은 발견이 아니라 번역 규칙이다.

이 문서는 설계 방법론과 그 수치적 실현에 집중한다. 그런 물성을 실제로 만드는 인공 구조와 유효 물성 추출 문제는 메타물질에, $n<0$ 이라는 상태 자체의 물리는 음의 굴절률에 있다.

2. 형태 불변성 — 규칙 전체가 한 줄이다[편집]

무전하·무전류 영역에서 맥스웰 방정식을 좌표 $x^i$ 에서 $x^{i'}$ 로 옮기면, 방정식의 모양은 그대로 유지되고 좌표 변환의 효과가 전부 물성 텐서로 흡수된다. 야코비 행렬을 $\Lambda^{i'}{}_{i} = \partial x^{i'}/\partial x^{i}$ 라 할 때

\varepsilon'^{\,i'j'} = \frac{\Lambda^{i'}{}_{i}\,\Lambda^{j'}{}_{j}\,\varepsilon^{ij}}{\det \Lambda}, \qquad \mu'^{\,i'j'} = \frac{\Lambda^{i'}{}_{i}\,\Lambda^{j'}{}_{j}\,\mu^{ij}}{\det \Lambda}

행렬로 쓰면 $\varepsilon' = \Lambda\,\varepsilon\,\Lambda^{\mathsf T}/\det\Lambda$ 다. 진공에서 출발하면 $\varepsilon=\mu=I$ 이므로 변환 후에도 $\varepsilon'=\mu'$ 가 되고, 그러면 매질의 임피던스가 진공과 같아 계면 반사가 원리적으로 0이다. 이것이 변환광학 설계가 흡수체가 아니라 “안내 장치”가 되는 이유다.

여기서 나오는 두 가지 사고 전환이 이 분야의 전부다.

장은 좌표를 따라간다. 원래 좌표계에서 직선으로 가던 광선은 새 좌표계의 좌표선을 따라 휘어 보인다. 구부러진 것은 빛이 아니라 격자다.
물성은 자유롭게 정하는 것이 아니라 계산되는 것이다. 사상 $x'(x)$ 를 정하는 순간 $\varepsilon'$ , $\mu'$ 는 강제된다. 설계 자유도는 물성이 아니라 사상에 있다.

3. 원통 클로킹 — 계산 예제이자 이 분야의 명함[편집]

반지름 $b$ 의 원판 내부를 껍질 $a<r'<b$ 로 밀어내는 압축 사상

r' = a + r\,\frac{b-a}{b}, \qquad \theta'=\theta,\qquad z'=z

를 넣으면 위 규칙이 곧바로 다음을 뱉는다.

\varepsilon_{r} = \mu_{r} = \frac{r'-a}{r'},\qquad \varepsilon_{\theta} = \mu_{\theta} = \frac{r'}{r'-a},\qquad \varepsilon_{z} = \mu_{z} = \left(\frac{b}{b-a}\right)^{2}\frac{r'-a}{r'}

원래 $r=0$ 이던 한 점이 반지름 $a$ 의 원으로 부풀었으므로, 그 안에는 장이 애초에 존재하지 않는다. 무엇을 넣든 바깥에서는 보이지 않는다. 2006년 펜드리·슈리그·스미스와 레온하르트가 같은 호에 각각 다른 경로로 이 결론에 도달했고, 같은 해 슈리그 등이 스플릿링 공진기 배열로 마이크로파 단일 주파수 2차원 클로크를 실증했다.¹

그런데 식을 보면 이 설계가 왜 완전할 수 없는지도 같이 보인다.

안쪽 경계에서 특이하다. $r'\to a$ 에서 $\varepsilon_\theta \to \infty$ , $\varepsilon_r \to 0$ 이다. 발산하는 곳은 클로크의 중심이 아니라 껍질의 안쪽 면이며, 유한한 값으로 잘라내는 순간 완전 은폐는 깨진다.
위상속도가 $c$ 를 넘는 지점이 생긴다. $\varepsilon_r\mu_r<1$ 인 영역이 필연이므로, 인과율(크라메르스-크로니히 관계)이 강한 주파수 의존성을 강제한다. 완전 클로킹은 본질적으로 단일 주파수 현상이고, 광대역 클로크는 물리적으로 금지되어 있다.
유한 대역 펄스는 은폐되지 않는다. 클로크를 돌아 나온 파는 직진한 파보다 먼 기하 경로를 갔는데도 같은 위상으로 도착해야 한다. 그런데 정보를 나르는 신호 선단(front)은 어떤 매질에서도 $c$ 를 넘지 못한다. 즉 완전 은폐는 단색 정상상태에서만 정의되는 개념이고, 펄스로 두드리면 반드시 흔적이 남는다.

실증 실험이 쓴 축소 물성(reduced parameters)은 이 벽을 우회한 타협이다. $E_z$ 편광에 대한 분산 관계만 보존하도록 $\mu_r = \left((r-a)/r\right)^2$ , $\mu_\theta = 1$ , $\varepsilon_z = \left(b/(b-a)\right)^2$ 로 재배분하면 만들어야 할 비자명 성분이 하나로 줄어든다. 대신 $\varepsilon\neq\mu$ 가 되어 임피던스 정합이 깨지고 껍질 바깥면에서 반사가 생긴다 — 실측 클로크의 성능이 이론값에 못 미치는 주된 이유다.

압축 사상 r′ = a + r(b−a)/b 가 만든 매질에서 광선방정식 d/ds(n dr/ds) = ∇n 을 RK4 로 적분한다. 이방성 원본은 편차 max 2.9e−7 b 로 숨긴 코어를 돌아 원래 직선을 되찾지만, 같은 사상을 등방 n(r) = b(r−a)/((b−a)r) 로 낮추면 최대 편향각이 143.21° 까지 벌어진다. 적분기는 루네버그 렌즈의 해석 초점으로 4차 수렴(오차비 15.95~16.01)을 확인했다.

4. 이방성을 없애는 법 — 준등각 변환과 카펫 클로크[편집]

위 물성표를 실제 재료로 만들려는 순간 벽에 부딪힌다. 위치마다 값이 다른 이방성 텐서를 요구하기 때문이고, 그러면 공진 구조(=손실+협대역)를 피할 수 없다. 그래서 나온 것이 이방성을 애초에 만들지 않는 사상을 고르자는 발상이다.

2차원에서 등각사상(conformal map)은 국소적으로 회전과 등방 확대만 하므로 $\Lambda\Lambda^{\mathsf T} \propto I$ 가 되고, 물성이 등방 굴절률 분포 $n(x,y)$ 하나로 줄어든다. 레온하르트의 접근이 정확히 이것으로, 복소해석의 해석함수 $w(z)$ 를 써서 $n = |dz/dw|$ 를 얻는다. 다만 등각사상은 경계 형상을 마음대로 지정할 수 없고(리만 사상 정리가 주는 자유도가 제한적이다), 이 방법의 은폐는 기하광학 극한에서만 엄밀하다.

준등각 변환광학(quasi-conformal TO, Li·펜드리 2008)은 이 제약을 수치적으로 완화한다. 원하는 경계 조건을 준 격자 생성 문제를 놓고 변형 Liao 범함수 같은 지표를 최소화해 이방성 비 $\alpha = \lambda_{\max}/\lambda_{\min}$ 가 1에 최대한 가까운 격자를 만든 뒤, 그 잔여 이방성을 무시하고 등방 $n$ 만 취한다. 결과가 유전체만으로 만들 수 있는 광대역·저손실 소자라는 점에서 변환광학이 실험실을 벗어난 결정적 계기가 됐다.

대표작이 카펫 클로크(ground-plane cloak)다. 평평한 거울 위의 융기를 평면 거울처럼 보이게 하는 사상으로, 융기 아래 공간을 숨긴다. 2009년 이후 마이크로파·근적외선·가시광에서 잇따라 실증됐고 실리콘 포토닉스 공정으로도 만들어졌다. 다만 공짜는 아니다.

준등각 근사가 남긴 잔여 이방성 때문에 반사 빔에 작지만 0이 아닌 측면 변위가 생긴다. 완전히 “평면 거울과 구별 불가”는 아니다.
자유공간에 떠 있는 물체가 아니라 지면 위 융기만 다룰 수 있다. 원통 클로킹 문제의 축소판이지 대체재가 아니다.
배경보다 굴절률이 큰 매질 속에서만 성립하도록 설계되는 경우가 많다 — 배경이 공기면 $n<1$ 영역이 요구되어 다시 분산이 필요해진다.

5. 광선 언어와의 접점 — 굴절률 분포형 소자[편집]

기하광학 극한에서 변환광학은 굴절률 분포 설계로 축약된다. 이 극한에서 광선은 광선 방정식

\frac{d}{ds}\left(n\,\frac{d\mathbf r}{ds}\right) = \nabla n

을 따르고, 이는 페르마 원리 $\delta\int n\,ds = 0$ 의 오일러-라그랑주 방정식이다. 즉 좌표 사상은 “광경로를 어떻게 다시 재느냐”의 문제로 바뀐다.

이 언어 안에서 변환광학은 고전 GRIN 소자 계보와 자연스럽게 이어진다.

루네버그 렌즈: $n(r) = \sqrt{2-(r/R)^2}$ 인 구. 평행 광선을 반대편 표면의 한 점에 모은다. 초점이 표면 위에 있으므로 광축이 없고, 모든 방향에서 동일하게 작동한다 — 위성 통신 안테나와 레이더 반사기에 실제로 쓰인다.
이턴 렌즈: 입사 광선을 온 방향으로 정확히 되돌려 보내는 역반사체. 중심에서 굴절률이 발산해 그대로는 만들 수 없고, 변환광학으로 중심 특이점을 유한 영역으로 펼치는 것이 표준 처방이다.
맥스웰 어안 렌즈: 구면 표면의 한 점을 대척점에 결상하는 고전 사상. 레온하르트의 등각 클로킹 구성에서 두 번째 리만 시트에 얹히는 매질이 이 계열이다.

역방향으로도 쓸모가 있다. 특이한 물성이 나오면 좌표계 탓을 해 보라는 것이 실무 감각이고, 실제로 완전정합층(PML)은 좌표를 복소수 방향으로 늘이는 변환의 물성 표현으로 정확히 재해석된다. PML의 이방성 텐서 형태가 클로크 물성과 닮아 보이는 것은 우연이 아니라 같은 공식의 두 사례다.

6. 수치 검증에서 반복해서 밟는 지뢰[편집]

변환광학 설계는 “계산이 잘 돌면 완벽하게 보이는” 부류의 문제라, 검증을 대충 하면 물리가 아니라 격자를 검증하게 된다.

특이 물성의 절단이 결과를 지배한다. 안쪽 경계를 $r=a$ 대신 $r=a+\delta$ 에서 끊거나 $\varepsilon_\theta$ 에 상한을 두는 것은 필수인데, 산란 단면적이 $\delta$ 에 민감하게 의존한다. $\delta$ 를 바꿔 가며 수렴 추세를 보고하지 않은 클로킹 계산은 신뢰할 수 없다. 이상적 원통 클로크가 물성의 작은 섭동에 극도로 취약하다는 것은 정식화 직후에 지적됐다.
FDTD에 상수로 넣으면 터진다. $\varepsilon_r<1$ 영역이 있으므로 비분산 상수는 초광속 위상속도를 뜻하고 시간영역 문제가 잘 정의되지 않는다. 드루드·로렌츠 분산 모형과 보조 미분방정식(ADE) 구현이 사실상 의무이며, 그러면 자동으로 협대역 여기만 검증하게 된다.
격자 정렬. 주파수영역 유한요소 전자기에서 $\varepsilon_\theta/\varepsilon_r = (r/(r-a))^2$ 는 안쪽 경계 근처에서 수천 배까지 벌어진다. 좌표선(원통이면 $r$ – $\theta$ )에 정렬된 사상 격자를 쓰지 않고 자유 삼각형 격자를 뿌리면, 요소 하나가 서로 다른 주축을 평균해 가짜 산란을 만든다. 은폐 성능이 격자 세밀화에 따라 단조롭게 좋아진다면 그건 물리가 아니라 이산화가 만든 은폐일 가능성이 높다.
PML과 클로크가 만나면 안 된다. 이방·분산 매질이 흡수 경계에 닿으면 표준 PML의 정합이 깨진다. 클로크를 진공층으로 둘러싸고 그 바깥에서 닫는 것이 안전하다.
평가 지표를 산란으로 잡아라. 장 분포 그림은 예쁘게 나오기 쉽다. 정량 지표는 벌거벗은 물체 대비 총 산란 단면적 감소율이고, 전방·후방 산란을 따로 보고해야 한다. 절차 자체는 검증 및 확인의 표준을 따르되, 여기서는 격자 수렴과 정칙화 파라미터 수렴 두 축을 모두 봐야 한다.

7. 다른 물리로의 이식 — 어디까지 형태 불변인가[편집]

같은 규칙이 통하려면 그 방정식이 좌표 변환에 대해 형태 불변이어야 한다. 이 조건은 물리마다 성패가 갈린다.

음향(유체). 압력에 대한 파동 방정식은 형태 불변이며, 대가로 관성 이방성 유체(밀도가 텐서, 체적탄성률은 스칼라)를 요구한다. 반대 극단에는 전단을 지지하지 않는 펜타모드 재료가 있고, 둘 사이는 한 개의 자유도를 가진 연속 계열로 이어진다는 것이 정리되어 있다. 실현은 음향 메타물질 소관.
탄성파. 여기서는 성립하지 않는다. 선형 탄성의 나비에 방정식은 일반 좌표 변환에서 형태를 유지하지 못하고, 변환된 방정식은 응력-속도가 결합된 윌리스형 구성식이 된다. “탄성 클로킹”이 음향보다 훨씬 어려운 근본 이유이며, 대칭성을 포기하거나 미소극성(Cosserat) 매질로 가는 우회로가 연구된다.
열전도. 확산 방정식도 형태 불변이라 이방성 열전도도 분포로 열 클로크·집열기를 설계할 수 있다. 파동이 아니라 위상을 맞출 필요가 없고, 그래서 분산 제약이 훨씬 약하다. 실제로 층상 복합재만으로 실증된 몇 안 되는 사례.
물결파. 수심이 변하는 선형 수면파 방정식에 같은 규칙을 적용해 원기둥을 우회하는 파를 만들 수 있고, 파장이 cm~m 급이라 실험실 수조에서 눈으로 확인 가능한 드문 클로킹이다.
정전·정상 전류. 라플라스 방정식은 당연히 형태 불변이고, 사실 이 갈래가 시간적으로 가장 앞선다. 전기 임피던스 단층촬영의 유일성 반례로서 특이 전도도 분포가 2003년에 이미 제시되어 있었다.²

8. 유산 — 망토는 못 만들었지만[편집]

투명 망토는 나오지 않았고, 위의 인과율 논증에 따르면 앞으로도 나오지 않는다. 그럼에도 이 방법론이 살아남은 이유는 설계 문제를 기하 문제로 바꿔 주기 때문이다.

메타표면·평면 광학. 원하는 위상 분포를 좌표 사상으로 유도한 뒤 국소 위상 소자로 이산화하는 절차가 표준이 됐다.
안테나·도파로. 벤딩 도파로의 곡률 손실 억제, 개구면 축소, 편평 반사판 설계처럼 “협대역이어도 상관없는” 마이크로파 문제에서는 실제로 쓰인다.
역설계와의 결합. 사상이 주는 물성 목표를 초기해로 삼고 위상 최적화·수반 감도 해석으로 제작 가능한 형상으로 옮기는 하이브리드가 흔하다. 좌표 사상은 최적화가 헤맬 초기해를 공짜로 주는 셈.
개념적 정리. 광결정의 음의 굴절, 쌍곡 매질의 이방성, 유효매질이론의 적층 극한 등 흩어져 있던 현상들을 “좌표를 어떻게 구부렸나”라는 단일 언어로 다시 읽게 해 줬다.³

9. 관련 문서[편집]

두 논문이 같은 학술지 같은 호에 나란히 실린 것은 우연이 아니라 편집부의 연출에 가깝다. 접근이 정반대였다는 점이 재미있는데, 한쪽은 전 벡터장을 다루는 물성 텐서 처방이고 다른 쪽은 복소평면 등각사상에 기반한 기하광학 구성이었다. 결과가 같은 방향을 가리켰으니 분야가 단숨에 굳었다. ↩
수학 쪽에서는 이것이 “클로킹”이 아니라 칼데론 역문제의 유일성 반례라는 이름으로 먼저 알려져 있었다. 물리학자들이 3년 뒤 같은 특이 텐서를 재발견해 망토라고 부르자 수학자들이 “우리가 먼저 했다”고 조용히 손을 든 사건. 학제 간 소통 부재의 교과서적 사례로 자주 인용된다. ↩
2006~2012년의 클로킹 논문 홍수는 규모 면에서 기억할 만하다. 사상 하나 고르고 야코비안 돌리면 논문 한 편이 나오다 보니, “새로운 형상의 클로크” 시리즈가 몇 년간 이어졌다. 결국 남은 것은 사상 목록이 아니라 준등각 변환 같은 제작 가능성 제약을 설계에 넣는 방법이었다. ↩