음의 굴절률

편집 역사 토론
전자기학 전자기해석 계산물리 마지막 수정: 2026-08-14 04:41:26

1. 개요[편집]

음의 굴절률
Negative refractive index
성립 조건ε < 0 이면서 동시에 μ < 0
별칭좌수 매질(LHM), 이중음성 매질(DNG)
예언베셀라고 (1967 러시아어 / 1968 영역)
실증셸비 등 (2001, 마이크로파 프리즘)
대표 현상역 굴절 · 역 도플러 · 역 체렌코프 · 평판 렌즈
수치해석 필수 조건분산 모형(드루드·로렌츠) 없이는 계산 불가

음의 굴절률(negative refractive index)은 유전율 ε\varepsilon 과 투자율 μ\mu 가 동시에 음수인 매질에서 굴절률의 물리적으로 옳은 가지가 n=−εμn = -\sqrt{\varepsilon\mu} 로 결정되는 상태를 말한다. 이런 매질 안에서는 E\mathbf{E}, H\mathbf{H}, k\mathbf{k} 가 오른손이 아니라 왼손 삼조를 이루므로 좌수 매질(left-handed medium)이라 부르고, 에너지가 가는 방향과 위상이 가는 방향이 반대가 된다.

자연 물질에는 없다. ε<0\varepsilon<0 은 플라즈마나 금속에서 흔하지만 μ<0\mu<0 은 GHz 이상에서 사실상 존재하지 않아서, 이 상태는 메타물질이라는 인공 구조가 나오기 전까지 30년간 사고실험이었다. 이 문서는 n<0n<0 이라는 상태 자체의 물리와 그것을 수치적으로 다룰 때 생기는 문제를 다룬다. 그런 물성을 만드는 구조(스플릿링 공진기, 금속선 배열)와 S-파라미터 역추출 이야기는 메타물질에, 부파장 혼합물의 균질화 일반론은 유효매질이론에 있다.

2. 부호는 선택이 아니라 강제된다[편집]

n2=εμn^2 = \varepsilon\mu 이고 ε,μ\varepsilon,\mu 가 둘 다 음수면 n2>0n^2>0 이라 파동은 정상적으로 전파한다. 그러면 nn 을 ++ 로 잡아도 되지 않느냐는 질문이 나오는데, 답은 안 된다이다. 수동 매질에는 아주 작더라도 손실이 있으므로 ε=−1+iδ\varepsilon = -1+i\delta, μ=−1+iδ\mu = -1+i\delta (δ>0\delta>0)로 두자. 그러면

εμ=(−1+iδ)2=1−δ2−2iδ,εμ≈1−iδ\varepsilon\mu = (-1+i\delta)^2 = 1 - \delta^2 - 2i\delta, \qquad \sqrt{\varepsilon\mu} \approx 1 - i\delta

인데 이 가지는 Im n<0\mathrm{Im}\,n<0, 즉 전파하면서 세기가 커지는 비물리적 해다. 에너지 보존이 요구하는 Im n≥0\mathrm{Im}\,n\ge 0 을 만족하는 가지는 n≈−1+iδn \approx -1+i\delta 뿐이다. 즉 음의 굴절률은 부호 규약이 아니라 수동성이 강제하는 결과다.

손실이 큰 실제 매질에서 ε,μ\varepsilon,\mu 의 실수부가 둘 다 음수여야만 하는 것도 아니다. 위상속도가 에너지 흐름과 반대가 될 조건은 드파인-라흐타키아 조건

ε′ ∣μ∣+μ′ ∣ε∣<0\varepsilon'\,|\mu| + \mu'\,|\varepsilon| < 0

로 완화되며(′' 은 실수부), 이 조건은 ε′<0,μ′<0\varepsilon'<0,\mu'<0 보다 넓다. 손실이 큰 메타물질 실험에서 μ′\mu' 이 겨우 0 언저리인데도 음의 굴절이 관측되는 경우를 이걸로 설명한다.

3. 무엇이 뒤집히는가[편집]

맥스웰 방정식의 평면파 형태 k×E=ωμH\mathbf{k}\times\mathbf{E} = \omega\mu\mathbf{H}, k×H=−ωεE\mathbf{k}\times\mathbf{H} = -\omega\varepsilon\mathbf{E} 에 음수 ε,μ\varepsilon,\mu 를 넣으면 삼조의 손잡이가 바뀐다. 포인팅 벡터 S=E×H\mathbf{S}=\mathbf{E}\times\mathbf{H} 는 여전히 에너지가 흐르는 방향이지만 k\mathbf{k} 와 역평행이 된다. 여기서 따라 나오는 것들.

  • 역 스넬 굴절. 계면에서 보존되는 것은 접선 파수 kxk_x 이고, n2<0n_2<0 이면 n1sin⁡θ1=n2sin⁡θ2n_1\sin\theta_1 = n_2\sin\theta_2 에서 θ2<0\theta_2<0 이 된다. 굴절 광선이 법선의 반대편이 아니라 같은 편으로 꺾인다. 2001년 셸비 등의 마이크로파 프리즘 실험이 잰 것이 정확히 이 각도다.
  • 역 도플러. 다가오는 광원의 관측 주파수가 오히려 낮아진다. 비선형 전송선에서 2003년에 실험적으로 확인됐다.
  • 역 체렌코프. 매질 내 위상속도보다 빠른 하전입자가 내는 복사가 뒤쪽 원뿔로 나간다. 도파관형 메타물질에서 2000년대 후반에 확인됐다.
  • 역 구스-핸헨 변위, 볼록 렌즈가 발산 렌즈로 동작하는 등 기하광학 규칙이 통째로 부호를 바꾼다.
  • 인과율은 멀쩡하다. 위상속도가 음수이거나 cc 를 넘어도 정보가 실리는 선단 속도(front velocity)는 여전히 cc 다. 2002년 발란주 등이 “군속도까지 고려하면 음의 굴절은 없다”고 제기한 반론은 유한 대역 펄스의 정상상태 도달 과정을 오해한 것으로 정리됐고, 이후 실험과 시간영역 계산이 정상상태에서의 음의 굴절을 재현했다.

4. 평판 렌즈와 완전 렌즈 논쟁[편집]

베셀라고가 1967년 논문에서 지적한 가장 인상적인 결과는 곡률 없는 렌즈다.1 두께 dd, n=−1n=-1 인 평판에 광원을 거리 d1<dd_1<d 에 두면, 광선이 판 안에서 반대로 꺾였다가 다시 꺾여 판 안쪽에 한 번, 판 바깥에 한 번 상을 맺는다. 광축도 없고 배율도 1이며, 광원이 판보다 멀면 상 자체가 생기지 않는다. 통상적인 렌즈와는 다른 종류의 물건이다.

2000년 펜드리는 여기에 훨씬 강한 주장을 얹었다. 통상 결상의 분해능을 제한하는 것은 소멸파(evanescent wave)다. 물체의 부파장 정보는 kx>ω/ck_x > \omega/c 성분에 실리는데 이들은 자유공간에서 지수적으로 죽어 상면에 도달하지 못한다. 그런데 ε=μ=−1\varepsilon=\mu=-1 인 판에서는 계면 표면파(표면 플라스몬-폴라리톤)가 공진해 소멸파가 지수적으로 증폭되고, 판 반대편에서 진폭이 정확히 복원된다. 전달함수가 모든 kxk_x 에 대해 1이 되므로 회절 한계가 사라진다는 것 — 이것이 완전 렌즈(perfect lens) 주장이다.

논쟁이 격렬했던 이유는 “증폭”이라는 단어 때문이었다. 에너지를 만들어내는 것 아니냐는 반론(‘t Hooft 2001, 가르시아·니에토-베스페리나스 2002)이 나왔지만, 소멸파는 실수 파워를 나르지 않으므로 진폭 증가가 에너지 증가가 아니라는 점에서 원리적 모순은 없다. 진짜 문제는 다른 데 있었다.2

  • 손실이 조금만 있어도 무너진다. ε=−1+iδ\varepsilon=-1+i\delta 이면 전달함수는 kxk_x 가 클수록 δ\delta 에 민감해지고, 복원 가능한 최대 파수는 대략 kmax⁡∼ln⁡(1/δ)/dk_{\max} \sim \ln(1/\delta)/d 로 로그로만 자란다. 분해능을 두 배 올리려면 손실을 제곱만큼 줄여야 한다는 뜻이라, 실질적으로 두께가 파장의 몇 분의 일인 근접장 영역에서만 이득이 남는다.
  • 완전한 ε=μ=−1\varepsilon=\mu=-1 은 특이점이다. 손실이 정확히 0이면 모든 소멸파 차수가 동시에 공진해 응답이 발산하고, 문제 자체가 비정칙(ill-posed)이 된다. 즉 손실은 결함이 아니라 문제를 잘 정의되게 만드는 정칙화다.
  • 결국 남은 것은 근접장 슈퍼렌즈다. kx≫ω/ck_x\gg\omega/c 인 정전기 극한에서는 TM 편광이 μ\mu 와 무관해져 ε=−1\varepsilon=-1 만 맞추면 된다. 은 박막이 자외선 영역에서 이 조건을 근사적으로 만족하므로, 2005년 팡 등이 365 nm 광원에서 60 nm 급 반주기 패턴을 결상해 보였다. 이어서 원통형 이방성 구조로 소멸파를 전파파로 변환해 원거리로 내보내는 하이퍼렌즈가 나왔는데, 이쪽은 유효매질이론의 적층 매질 이방성이 직접 쓰이는 사례다.

5. 분산과 손실은 옵션이 아니다[편집]

”ε=−1\varepsilon=-1, μ=−1\mu=-1 인 손실 없는 물질”은 좁은 주파수 창에서조차 존재할 수 없다. 이유는 열역학과 인과율 두 갈래다.

분산 매질에서 시간평균 에너지 밀도는 란다우-리프시츠 표현으로

u=12[∂(ωε)∂ω∣E∣2+∂(ωμ)∂ω∣H∣2]>0u = \frac{1}{2}\left[\frac{\partial(\omega\varepsilon)}{\partial\omega}|\mathbf{E}|^{2} + \frac{\partial(\omega\mu)}{\partial\omega}|\mathbf{H}|^{2}\right] > 0

이어야 한다. ε\varepsilon 이 상수 −1-1 이면 ∂(ωε)/∂ω=−1<0\partial(\omega\varepsilon)/\partial\omega = -1 < 0 이라 에너지 밀도가 음수가 된다. 따라서 ε<0\varepsilon<0 인 구간에는 반드시 강한 주파수 의존성이 있어야 하고, ∂(ωε)/∂ω>0\partial(\omega\varepsilon)/\partial\omega>0 을 만족하려면 ε\varepsilon 이 ω\omega 에 대해 가파르게 증가해야 한다. 여기에 크라메르스-크로니히 관계를 얹으면 가파른 실수부 변화는 근처에 큰 허수부, 즉 흡수를 동반한다. 결론은 하나다 — 광대역·저손실 음굴절률은 물리적으로 금지되어 있다. 메타물질의 협대역성과 손실은 설계 미숙이 아니라 정리의 귀결이다.

6. 수치해석 — 상수로 넣으면 터진다[편집]

FDTD로 음굴절률 매질을 다루려는 사람이 가장 먼저 하는 실수가 격자에 εr=−1\varepsilon_r=-1, μr=−1\mu_r=-1 을 상수로 채워 넣는 것이다. 그러면 갱신식이 즉시 발산한다.3 위 절의 이유로 비분산 음수 상수 자체가 비인과적이라 시간영역 문제가 잘 정의되지 않기 때문이다. 제대로 된 절차는 이렇다.

  • 분산 모형을 반드시 넣는다. 표준은 드루드 모델 쌍이다. ε(ω)=1−ωpe2/[ω(ω+iγe)]\varepsilon(\omega) = 1 - \omega_{pe}^2/[\omega(\omega+i\gamma_e)], μ(ω)=1−ωpm2/[ω(ω+iγm)]\mu(\omega) = 1 - \omega_{pm}^2/[\omega(\omega+i\gamma_m)] 로 두고, 목표 주파수 ω0\omega_0 에서 εr=μr=−1\varepsilon_r=\mu_r=-1 이 되게 하려면 무손실 극한에서 ωpe=ωpm=2 ω0\omega_{pe}=\omega_{pm}=\sqrt{2}\,\omega_0 로 잡으면 된다. 시간영역 구현은 보조 미분방정식(ADE)이나 (조각별 선형) 재귀 합성곱으로 한다. 지올코프스키·헤이먼의 초기 FDTD 연구가 이 관행을 정착시켰다.
  • 정상상태 도달이 아주 느리다. 음굴절률 대역은 공진 근처라 QQ 가 높다. 평판 렌즈의 초점이 형성되려면 소멸파 공진이 충분히 자라야 하는데, 이 과도 과정이 수만 시간스텝을 잡아먹는다. 계산을 일찍 끊고 “초점이 안 생긴다”고 결론 내리는 것이 이 분야의 고전적 오독이다.
  • 수치분산이 공진 조건을 깬다. 이(Yee) 격자의 수치분산 때문에 격자가 지지하는 실효 파수가 이론값과 미세하게 어긋난다. ε=−1\varepsilon=-1 공진은 그 어긋남에 극도로 민감해서, 셀 크기를 반으로 줄이면 결과가 정성적으로 달라지는 일이 생긴다. 격자 수렴 검사를 생략할 수 없는 대표적 문제다.
  • 경계 처리. 개방 경계는 완전정합층으로 닫는데, 분산 매질이 PML에 접하면 표준 PML의 성능이 떨어지므로 매질을 PML 앞에서 끊거나 분산 대응 PML을 쓴다.
  • 주파수영역에서도 편하지 않다. 유한요소 전자기에서 ε\varepsilon 의 부호가 계면을 사이에 두고 바뀌면 변분 형식의 강제성(coercivity)이 사라진다. 이 부호 변화 계수 문제는 대비비 ε2/ε1\varepsilon_2/\varepsilon_1 가 특정 임계 구간에 들어가면 실제로 잘 정의되지 않고, 임계 구간의 폭은 계면의 모서리 각도에 의존한다. 매끄러운 계면에서는 대비비 −1-1 하나만 병적이지만, 뾰족한 코너가 있으면 병적 구간이 넓어진다. T-강제성 이론이 이 조건을 정리했고, 실무적으로는 약간의 손실을 넣어 정칙화하거나 코너를 둥글리는 것으로 대응한다. 손실을 넣더라도 계는 부정부호가 되어 조건수가 나쁘고, 표준 반복 해법·전처리기가 잘 듣지 않는다.

검증은 어떻게 하나. 평판 계면에서 프레넬 방정식의 해석해와 반사·투과 계수를 비교하는 것이 1차 관문이고, 표면 플라스몬 분산 관계를 수치로 뽑아 해석식과 맞춰 보는 것이 2차 관문이다. 시간영역이라면 정상상태 도달을 확인한 뒤 위상 전선의 진행 방향과 포인팅 벡터 방향이 반대인지를 직접 그려 본다. 검증 및 확인의 표준 절차를 그대로 따르되, 이 문제에서는 격자 수렴과 정상상태 수렴 두 축을 모두 봐야 한다는 점이 특이하다.

7. 음굴절 ≠ 음의 굴절률[편집]

용어가 자주 뒤섞이므로 구별해 둘 필요가 있다. 광선이 법선 같은 편으로 꺾이는 현상(음의 굴절)은 n<0n<0 없이도 일어난다.

  • 광결정의 음의 굴절. 주기 a∼λa \sim \lambda 인 구조에서는 등주파수 곡면의 곡률이 뒤집힐 수 있고, 군속도가 파수벡터와 예각을 이루지 않는 대역이 생긴다. 노토미(2000), 루오 등(2002)이 정리한 “유효 굴절률이 음수가 아니면서 모든 입사각에서 음의 굴절”이 이 부류다. 여기서는 블로흐 정리와 밴드 구조 계산이 도구이고, ε\varepsilon·μ\mu 두 숫자로 요약하는 접근은 애초에 성립하지 않는다.
  • 이방성 적층 매질. 금속-유전체 다층으로 만든 쌍곡 매질은 ε∥\varepsilon_\parallel 과 ε⊥\varepsilon_\perp 의 부호가 달라 등주파수면이 쌍곡면이 되고, 그 결과 굴절이 음으로 나타난다. 이것 역시 n<0n<0 이 아니라 이방성의 효과이며 계산은 유효매질이론의 위너 극한값에서 출발한다.
  • 전송선 회로. 직렬 커패시터·병렬 인덕터로 구성한 이중 전송선은 후진파를 지지해 마이크로파 대역에서 좌수 거동을 보인다. 손실이 작고 광대역이라 실용 소자(위상 보상기, 누설파 안테나)로 실제 쓰이는 유일한 갈래에 가깝다.

즉 “음의 굴절을 봤다”는 관측만으로 ε<0\varepsilon<0, μ<0\mu<0 을 주장할 수 없다. 이 구별을 흐린 논문이 2000년대 초에 상당히 많았고, 메타물질 문서에서 다루는 S-파라미터 역추출의 애매함도 상당 부분 여기서 왔다.4 좌표 변환으로 물성 텐서를 설계하는 변환광학은 이 갈래들을 하나의 기하학적 언어로 묶으려는 시도로 볼 수 있다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 베셀라고의 논문은 러시아어 원문이 1967년, 영역이 1968년 Soviet Physics Uspekhi 에 실렸다. 30년 가까이 인용이 거의 없다가 2000년대에 폭발했는데, 정작 본인은 “재료를 만들 수 있다고는 생각하지 못했다”고 회고했다. 이론이 재료를 30년 기다린 셈. ↩

  2. 완전 렌즈 논쟁은 물리학에서 드물게 댓글(Comment)과 답변(Reply)만으로 서너 라운드가 오간 사례다. 결론적으로 원리적 오류는 없었고 실현 가능성이 과대평가되었다는 쪽으로 정리됐는데, 이런 결말이 가장 흔한 결말이기도 하다. ↩

  3. ”ε=−1\varepsilon=-1 을 넣었더니 FDTD가 폭발했어요”는 이 바닥 게시판의 단골 질문이다. 답은 언제나 같다 — 물리적으로 존재할 수 없는 매질을 넣었으니 수치가 그걸 알려준 것이다. 발산이 버그가 아니라 진단인 흔치 않은 경우. ↩

  4. 그래서 이 주제의 논문을 읽을 때 첫 번째로 확인할 것은 저자가 재는 “굴절률”이 위상 기준인지 군속도 기준인지, 그리고 그것이 단일 주파수 정상상태인지 펄스 응답인지다. 2000년대 초 논쟁의 절반은 서로 다른 양을 같은 이름으로 부른 데서 나왔다. ↩