메타물질

전자기학 전자기해석 계산물리 마지막 수정: 2026-08-13 04:36:09

1. 개요[편집]

메타물질 Metamaterial
정의	파장 ≫ 주기인 인공 단위구조 배열
이론적 예언	베셀라고 (1968) — ε·μ 동시 음수 매질
핵심 소자	스플릿링 공진기(SRR), 금속선 배열
실험 실증	스미스 등 (2000), 셸비 등 (2001, 마이크로파 음의 굴절)
주요 수치기법	FDTD, FEM + 주기경계조건
현실의 벽	손실, 강한 분산, 협대역

메타물질(metamaterial)은 파장보다 훨씬 작은 인공 단위구조를 주기적(또는 준주기적)으로 배열해, 구성 재료 자체에는 없는 유효 물성 $\varepsilon$ · $\mu$ 를 인위적으로 설계한 매질이다. 접두사 meta-(“넘어선”)가 붙은 이유는 자연에 있는 물질이 주지 못하는 응답 — 대표적으로 음의 유전율·음의 투자율 — 을 구조로 만들어내기 때문이다.

핵심 조건은 스케일 분리다. 단위구조 주기 $a$ 가 파장 $\lambda$ 보다 충분히 작으면( $a \ll \lambda$ ) 파동은 개별 구조물을 보지 못하고 평균화된 매질 하나로 인식한다. 이 점이 광결정과의 결정적 차이다 — 광결정은 $a \sim \lambda$ 라서 회절과 밴드 구조로 놀고, 메타물질은 $a \ll \lambda$ 라서 유효 물성 두 숫자로 논다.¹ 둘의 경계는 흐릿하고, 실제 메타물질 논문의 상당수가 $a/\lambda \sim 1/5$ 언저리라는 애매한 지대에 있다.

2. 없는 물성을 만드는 법[편집]

전기 응답( $\varepsilon$ )을 조작하는 건 상대적으로 쉽다. 가는 금속선을 격자로 세우면 유효 플라즈마 주파수가 마이크로파 대역까지 내려오고, 그 아래에서 드루드형 $\varepsilon(\omega) < 0$ 이 된다. 어려운 쪽은 자기 응답이다. 자연 물질의 $\mu$ 는 GHz만 넘어가도 1에 붙어버린다. 자기 쌍극자를 만들려면 전류 고리가 필요한데, 원자 규모에서는 그런 게 없기 때문이다.

여기서 나온 답이 스플릿링 공진기(Split-Ring Resonator, SRR)다. 틈이 난 금속 고리는 고리가 인덕턴스 $L$ , 틈이 커패시턴스 $C$ 인 LC 공진기가 되고, 입사 자기장이 고리를 관통하면 유도 전류가 돌아 자기 쌍극자를 만든다. 공진 주파수 $\omega_0 = 1/\sqrt{LC}$ 근처에서 유효 투자율은 로렌츠형으로 요동친다.

\mu_{\mathrm{eff}}(\omega) = 1 - \frac{F\omega^{2}}{\omega^{2} - \omega_0^{2} + i\gamma\omega}

$F$ 는 고리가 차지하는 면적 비율(대략 채움률), $\gamma$ 는 손실이다. 공진 바로 위쪽 좁은 창에서 $\mathrm{Re}\,\mu_{\mathrm{eff}} < 0$ 이 된다. 이 아이디어를 정리한 것이 펜드리 등의 1999년 논문이고, 여기에 금속선 배열( $\varepsilon<0$ )을 겹치면 $\varepsilon$ 과 $\mu$ 가 동시에 음수인 매질이 만들어진다. 2000년 스미스 등이 그 복합 구조를 실제로 만들었고, 2001년 셸비 등이 마이크로파 프리즘 실험으로 음의 굴절을 보였다.

3. 음의 굴절률과 좌수 매질[편집]

$\varepsilon<0$ , $\mu<0$ 이면 굴절률 제곱 $n^2 = \varepsilon\mu$ 는 다시 양수라 파동은 정상적으로 전파한다. 그런데 인과율을 지키는 제곱근 가지를 고르면 $n = -\sqrt{\varepsilon\mu}$ , 즉 굴절률이 음수가 된다. 이 매질에서 벡터 $\mathbf{E}$ , $\mathbf{H}$ , $\mathbf{k}$ 는 오른손이 아니라 왼손 삼조를 이루므로 좌수 매질(left-handed medium)이라 부른다. 1968년 베셀라고가 순수 이론으로 예언했지만 재료가 없어 30년 넘게 방치되어 있었다.²

물리적으로 벌어지는 일은 다음과 같다.

위상속도와 에너지 흐름이 반대. 포인팅 벡터 $\mathbf{S}$ 와 파수벡터 $\mathbf{k}$ 가 역평행이다. 즉 위상은 소스 쪽으로 되돌아오고 에너지만 앞으로 간다.
스넬 법칙이 뒤집힌다. 굴절광이 법선 반대편이 아니라 같은 편으로 꺾인다. 셸비 실험이 잰 것이 바로 이 각도.
평판 렌즈. $n=-1$ 인 평평한 판이 초점을 만든다. 곡률이 필요 없다. 나아가 펜드리(2000)는 소멸파(evanescent wave)까지 증폭되어 회절 한계를 넘는 결상이 가능하다고 주장했다 — 이른바 완전 렌즈. 다만 손실이 조금만 있어도 소멸파 복원이 무너져, 현실에서는 근접장 “슈퍼렌즈”의 부분적 실증에 그친다.
도플러·체렌코프 역전: 다가오는 소스의 진동수가 오히려 낮아지고, 체렌코프 복사가 뒤로 나간다.

4. 호모제나이제이션과 S-파라미터 역추출 — 가장 미묘한 부분[편집]

“유효 $\varepsilon$ , $\mu$ ” 라는 말을 쓰려면 그 값을 정의하고 뽑아내야 한다. 실무에서 압도적으로 많이 쓰는 방법이 S-파라미터 역추출(니컬슨-로스-바이어 계열, 메타물질 맥락에서는 스미스 등 2002)이다. 단위 셀 한 층을 두께 $d$ 의 균질 판으로 가정하고, 수치해석이나 실험으로 얻은 반사·투과 계수 $S_{11}$ , $S_{21}$ 을 뒤집어 임피던스 $z$ 와 굴절률 $n$ 을 구한다.

z = \pm\sqrt{\frac{(1+S_{11})^{2} - S_{21}^{2}}{(1-S_{11})^{2} - S_{21}^{2}}}, \qquad \cos(nkd) = \frac{1 - S_{11}^{2} + S_{21}^{2}}{2S_{21}}

그다음 $\varepsilon = n/z$ , $\mu = nz$ . 계산 자체는 몇 줄인데, 여기에 이 분야의 유명한 애매함이 전부 몰려 있다.

가지 모호성. $\cos^{-1}$ 은 $2\pi$ 만큼의 불확정성이 있고 제곱근은 부호가 둘이다. 수동 매질 조건( $\mathrm{Re}\,z \ge 0$ , $\mathrm{Im}\,n \ge 0$ )으로 대부분 고정되지만, 두꺼운 시료나 강한 분산에서는 가지 선택이 사람 손을 탄다.
두께 $d$ 를 어디로 잡을 것인가. 셀 경계를 반 셀 옮기면 추출된 $\varepsilon$ , $\mu$ 가 달라진다. 물리량이 좌표 원점에 의존하면 그건 물성이 아니다.
반공진 인공물. 추출 결과에 $\mathrm{Im}\,\varepsilon < 0$ 같은 수동 매질에서 불가능한 값이 튀어나오는 일이 흔하다. 이는 재료가 이상해서가 아니라 균질 매질 가정 자체가 깨졌다는 신호다.
공간 분산. $a/\lambda$ 가 충분히 작지 않으면 유효 물성이 입사각· $\mathbf{k}$ 에 의존한다. $\varepsilon(\omega)$ 가 아니라 $\varepsilon(\omega,\mathbf{k})$ 인 것이고, 그러면 “두 숫자로 요약”이라는 전제가 무너진다.

그래서 엄밀한 쪽에서는 장 평균(field averaging)이나 다중극 전개 기반의 호모제나이제이션을 쓰지만, 여전히 정의가 유일하지 않다. 메타물질에서 $\varepsilon$ , $\mu$ 는 측정값이라기보다 모형의 매개변수에 가깝다 — 이 점을 잊고 추출값을 물성표처럼 인용하는 것이 이 분야의 고전적 사고다.

수치 측면에서는 단위 셀 하나에 주기경계조건(비스듬한 입사면 블로흐 위상 포함)을 걸고 FDTD 또는 유한요소 전자기로 $S$ 파라미터를 뽑는 것이 표준 절차다. 개방 방향은 완전정합층으로 닫고, 밴드 구조가 필요하면 블로흐 정리 기반 고유값 해석으로 넘어간다.

5. 수치해석 워크플로[편집]

메타물질 논문 한 편에 들어가는 계산은 대개 다음 순서로 굳어져 있다.

단위 셀 해석. 셀 하나를 격자로 채우고 전파 방향을 제외한 축에 주기경계조건을 건다. 비스듬한 입사를 다루려면 블로흐 위상 $e^{i\mathbf{k}\cdot\mathbf{a}}$ 를 경계에 실어야 하고, 시간영역(FDTD)에서는 이게 “미래의 값이 필요”해지는 인과율 문제를 낳아 사인파 정상상태 해석이나 사분할 기법으로 우회한다.
포트와 흡수 경계. 전파 방향 양 끝에 포트를 두어 $S_{11}$ , $S_{21}$ 을 재고, 바깥은 완전정합층으로 닫는다. PML이 부실하면 반사가 섞여 추출된 $\varepsilon$ 이 요동친다.
분산 물성 모형. 금속은 드루드 또는 드루드-로렌츠 모형으로 넣는다. 시간영역에서는 보조 미분방정식이나 재귀 합성곱으로 분산을 구현해야 하며, 여기서 실수하면 손실이 통째로 틀린다.
밴드 해석. $S$ 파라미터 대신 블로흐 정리 기반 고유값 문제를 풀어 $\omega(\mathbf{k})$ 를 얻으면, 유효 굴절률을 가지 모호성 없이 정의할 수 있다. 복소 밴드(감쇠 대역)를 보려면 $k$ 를 미지수로 두는 정식화가 필요하다. 절차 자체는 밴드 구조 계산과 같다.
역설계. 요즘은 셀 형상을 사람이 그리는 대신 위상 최적화나 수반법 기반 감도 해석으로 자동 생성한다. 목표 $S$ 파라미터나 목표 위상 분포를 손실 함수로 놓고 셀 안의 재료 분포를 최적화하는 방식이며, 나오는 형상은 대체로 사람이 보고 이해할 수 없게 생겼지만 성능은 좋다.

산란체가 유전체 구·원기둥 배열이라면 격자 대신 미 산란·T-행렬법 기반 다중 산란 해석으로 훨씬 싸게 갈 수도 있다. 반대로 형상이 임의적이면 이산 쌍극자 근사의 주기 확장이 대안이 된다.

6. 변환광학과 클로킹[편집]

메타물질의 두 번째 큰 줄기는 변환광학(transformation optics)이다. 맥스웰 방정식은 좌표 변환에 대해 형태 불변이고, 좌표를 구부린 효과는 $\varepsilon$ 과 $\mu$ 를 텐서로 재정의하는 것으로 정확히 흡수된다. 그러면 “빛이 이렇게 휘었으면 좋겠다”는 기하학적 소원을 그대로 물성 분포로 번역할 수 있다. 좌표계를 설계하고 물성을 받아 적는 것이 전부라, 개념적으로는 놀랄 만큼 단순하다.

가장 유명한 응용이 클로킹이다. 한 점을 유한한 구로 부풀리는 사상을 쓰면, 그 구 안으로는 장이 들어가지 않고 파면은 구를 매끄럽게 돌아 원래 경로로 복귀한다. 2006년 펜드리·슈리그·스미스와 레온하르트가 각각 제안했고, 같은 해 슈리그 등이 SRR 배열로 마이크로파 단일 주파수 2차원 클로크를 실증했다. 그러나 요구되는 물성이 가혹하다.

극단적 이방성과 위치에 따라 변하는 텐서 물성.
경계에서 $\varepsilon$ 또는 $\mu$ 성분이 0 또는 무한대로 발산.
위상속도가 $c$ 를 넘어야 하는 지점이 생기고, 이는 강한 분산을 강제한다 → 본질적으로 협대역.

그래서 “해리 포터 투명 망토”는 상업적으로도 물리적으로도 나오지 않았다. 대신 카펫 클로크(지면 위 융기를 평면으로 위장), 열·음향·탄성파 클로킹, 그리고 안테나 정합·산란 저감 같은 덜 화려하지만 실용적인 파생물이 남았다.

7. 손실과 분산이라는 벽[편집]

메타물질을 실용화하지 못하게 붙잡는 것은 아이디어가 아니라 재료다.

손실. 음의 $\mu$ 는 공진에서 나오고, 공진은 전류가 금속에 몰린다는 뜻이며, 금속은 저항이 있다. 광학 주파수로 올라가면 은·금의 손실이 급격히 커져 성능지수 $\mathrm{FOM} = |\mathrm{Re}\,n|/\mathrm{Im}\,n$ 가 한 자릿수에 머문다. 이득 매질 삽입, 유전체 마이 공진기(금속 없는 all-dielectric 메타물질) 등이 대안으로 연구된다.
분산은 선택이 아니라 의무. 인과율(크라머스-크로니히 관계)과 에너지 밀도 양수 조건 $\partial(\omega\varepsilon)/\partial\omega > 0$ 은 손실 없는 매질에서 $\varepsilon<0$ 이 되려면 반드시 강한 주파수 의존성이 있어야 함을 요구한다. 즉 광대역 음굴절률은 물리적으로 금지되어 있다.
제작. 마이크로파 대역은 PCB 식각으로 되지만, 가시광 대역 3차원 메타물질은 수십 nm 구조를 부피로 쌓아야 한다. 그래서 실용 성과는 두께 한 층짜리 메타표면(metasurface) 쪽으로 옮겨 갔다 — 위상만 국소적으로 제어해 렌즈·편광기·홀로그램을 만드는 접근이며, 이쪽은 이미 상용 카메라 부품 수준까지 왔다.

8. 음향·탄성 메타물질로의 확장[편집]

같은 논리는 파동방정식을 갖는 모든 물리계로 번역된다. 음향에서는 질량밀도 $\rho$ 와 부피탄성률 $B$ 가 $\mu$ · $\varepsilon$ 의 역할을 한다.

음의 유효 질량밀도: 고무로 감싼 납 구슬을 배열한 국소 공진 구조(류 등, 2000). 공진 위쪽에서 매질이 힘과 반대로 가속하는 것처럼 거동하며, 파장의 수백 분의 일 두께로 저주파 소음을 막는다. 질량법칙을 우회한다는 점에서 방음 공학의 오랜 꿈이었다.
음의 부피탄성률: 헬름홀츠 공명기를 도관에 매단 배열. 둘을 합치면 음향 음굴절률.
탄성·구조 메타물질: 밴드갭으로 진동을 차단하는 포노닉 구조, 푸아송비가 음수인 오그제틱 격자, 전단을 거의 지지하지 않는 펜타모드 물질. 여기서는 위상 최적화와 형상 최적화가 단위 셀 설계의 주력 도구가 된다.

지진파 차폐용 대형 메타 구조, 열 클로킹까지 확장되어 있는데, 확장될수록 “메타물질”이라는 단어의 외연이 넓어져 사실상 주기 미세구조로 유효 물성을 설계하는 모든 것을 가리키는 우산 용어가 되어 버렸다는 비판도 있다.³

9. 관련 문서[편집]

그래서 “이 구조는 메타물질인가 광결정인가”는 논문 심사에서 은근히 자주 붙는 시비다. 판정 기준은 저자의 취향이 아니라 $a/\lambda$ 이며, 이게 1/4을 넘어가면 유효 물성 이야기는 슬슬 신뢰를 잃는다. ↩
베셀라고의 1968년 논문은 러시아어 원문이 1967년, 영역이 1968년에 나왔고 30년 가까이 인용이 거의 없었다. 물성이 없어서 “재미있는 사고실험” 취급을 받은 것. 지금은 인용 수가 수만 회다. 시대를 너무 앞서면 논문도 냉동보관된다. ↩
실제로 2010년대 중반 이후 “metamaterial”을 제목에 넣으면 인용이 잘 된다는 이유로 평범한 주기 복합재 연구까지 이 이름을 다는 현상이 있었다. 학계의 유행어 인플레이션은 재료를 가리지 않는다. ↩