표면 플라스몬

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

1. 개요[편집]

표면 플라스몬
Surface Plasmon
정체금속-유전체 계면에 갇힌 전자 밀도파 + 전자기파의 결합 모드
정식 명칭표면 플라스몬 폴라리톤(SPP)
존재 조건TM(p) 편광 전용, εmεd < 0 이고 εm+εd < 0
분산광선(light line) 오른쪽 — 빛으로 직접 못 때린다
결합법크레치만·오토 프리즘, 격자, 근접장 프로브
한계금속 옴 손실 — 가시광에서 전파거리 수십 μm

빛을 파장보다 작게 가두는 방법은 있다. 대신 금속이 그 대가를 열로 받아 간다.

표면 플라스몬(surface plasmon)은 금속과 유전체의 계면을 따라 전파하면서 양쪽으로는 지수적으로 감쇠하는, 자유전자의 집단 진동과 전자기장이 결합한 표면파다. 장(field)이 함께 실려 있다는 점을 강조할 때는 표면 플라스몬 폴라리톤(surface plasmon polariton, SPP)이라고 부르고, 금속 나노입자처럼 전파하지 못하고 한 곳에 갇힌 경우는 국재 표면 플라스몬(localized surface plasmon, LSP)으로 구분한다.1

핵심은 하나다. 이 모드는 회절 한계와 무관하게 파장보다 훨씬 짧은 스케일로 장을 압축한다. 그래서 슈퍼렌즈의 소멸파 증폭, 표면증강 라만, SPR 바이오센서, 근접장 리소그래피가 전부 이 하나의 모드에 기대고 있다. 동시에 이 모드는 금속 안을 통과하므로 옴 손실을 피할 방법이 원리적으로 없다. 플라스모닉스의 20년은 “가둠은 좋은데 손실은 어쩔 거냐”는 질문 하나를 두고 돌아왔다.

2. 계면에 갇힌 파 — 왜 TM만 되는가[편집]

계면을 z=0z=0 에 두고 z<0z<0 을 금속(εm\varepsilon_m), z>0z>0 을 유전체(εd\varepsilon_d)로 놓자. xx 방향으로 ei(kx−ωt)e^{i(kx-\omega t)} 로 진행하고 zz 방향으로는 e−κj∣z∣e^{-\kappa_j |z|} 로 죽는 해를 찾는다. 각 영역에서 헬름홀츠 방정식이 요구하는 것은

κj2=k2−εjω2c2,Re κj>0\kappa_j^{2} = k^{2} - \varepsilon_j \frac{\omega^{2}}{c^{2}}, \qquad \mathrm{Re}\,\kappa_j > 0

이다. 여기서 편광에 따라 결과가 극적으로 갈린다.

TM(p) 편광 — 자기장이 H=(0,Hy,0)\mathbf{H}=(0,H_y,0) 이고 전기장이 xx–zz 평면에 있다. 맥스웰 방정식에서 Ex∝(1/ε) ∂zHyE_x \propto (1/\varepsilon)\,\partial_z H_y 이므로, HyH_y 와 ExE_x 의 연속 조건을 합치면

κdεd+κmεm=0⟺κmκd=−εmεd\frac{\kappa_d}{\varepsilon_d} + \frac{\kappa_m}{\varepsilon_m} = 0 \qquad\Longleftrightarrow\qquad \frac{\kappa_m}{\kappa_d} = -\frac{\varepsilon_m}{\varepsilon_d}

가 나온다. κd,κm\kappa_d,\kappa_m 은 정의상 양수이므로 εm\varepsilon_m 과 εd\varepsilon_d 의 부호가 서로 달라야 한다. 유전체가 εd>0\varepsilon_d>0 이니 금속은 εm<0\varepsilon_m<0 이어야 하고, 이는 플라스마 주파수 아래에서 금속이 실제로 갖는 성질이다.

TE(s) 편광 — 전기장이 E=(0,Ey,0)\mathbf{E}=(0,E_y,0) 일 때 같은 절차를 밟으면 Ey(κd+κm)=0E_y(\kappa_d+\kappa_m)=0 이 나온다. 두 κ\kappa 가 모두 양수이므로 해는 Ey=0E_y=0, 즉 TE 표면파는 존재하지 않는다. 이것은 근사도 관행도 아니라 경계조건이 내리는 판결이다. 플라스모닉 소자 설계도가 하나같이 “TM only”라고 적혀 있는 이유가 여기 있다.

3. 분산 관계와 운동량 불일치[편집]

위 조건을 κj\kappa_j 정의에 대입하면 SPP의 분산 관계가 닫힌 꼴로 나온다.

kspp(ω)=ωcεm(ω) εdεm(ω)+εdk_{\mathrm{spp}}(\omega) = \frac{\omega}{c}\sqrt{\frac{\varepsilon_m(\omega)\,\varepsilon_d}{\varepsilon_m(\omega)+\varepsilon_d}}

ksppk_{\mathrm{spp}} 가 실수이면서 유전체 쪽으로 속박되려면 kspp>εd ω/ck_{\mathrm{spp}} > \sqrt{\varepsilon_d}\,\omega/c 여야 하고, 이는 εm<0\varepsilon_m<0 에 더해

εm+εd<0\varepsilon_m + \varepsilon_d < 0

를 요구한다. 즉 금속의 유전율 실수부 크기가 유전체보다 커야 한다. 이 두 줄이 SPP 존재 조건의 전부다.

드루드 금속(ε∞=1)과 유전체 계면의 SPP 분산관계를 복소수로 풀어, ω–k 가지가 광선 오른쪽으로 휘어 ω_sp=ω_p/√(1+ε_d) 에 점근하는 것과 그 ω 에서 계면 양쪽 |E_z| 의 지수 프로파일을 같은 nm/px 축척으로 그린다. ħω_p=9 eV 로 환산하면 ε_m=−18(λ₀≈600 nm)에서 유전체 쪽 침투깊이 δ_d=394 nm, 금속 쪽 δ_m=21.9 nm 로 18배 차이가 난다. γ/ω_p 를 올리면 분산곡선이 ω_sp 못 미쳐 되돌아오고 전파길이 L=1/(2 Im k) 도 68.6 µm(γ/ω_p=0.005)에서 17.2 µm(0.02)로 짧아진다.

여기서 이 분야 전체를 규정하는 사실이 하나 튀어나온다. ksppk_{\mathrm{spp}} 는 언제나 같은 주파수에서 유전체 속을 달리는 빛의 파수보다 크다. 분산 곡선이 광선(light line) ω=ck/εd\omega = ck/\sqrt{\varepsilon_d} 의 오른쪽 아래에 통째로 놓인다는 뜻이고, 따라서

어떤 각도로 어떻게 때려도, 자유공간에서 온 평면파는 SPP를 여기시킬 수 없다.

운동량이 모자란다. 이 운동량 불일치(momentum mismatch)가 “금속에 레이저를 쐈는데 왜 아무 일도 안 일어나는가”의 답이고, 다음 절의 온갖 결합 장치들이 존재하는 이유다.

분산 곡선의 반대쪽 끝도 봐 두자. k→∞k\to\infty 극한에서는 εm+εd→0\varepsilon_m+\varepsilon_d\to 0 이 되어야 하므로, 무손실 드루드 모델 εm=ε∞−ωp2/ω2\varepsilon_m = \varepsilon_\infty - \omega_p^2/\omega^2 를 넣으면 곡선이 수평 점근선

ωsp=ωpε∞+εd\omega_{\mathrm{sp}} = \frac{\omega_p}{\sqrt{\varepsilon_\infty + \varepsilon_d}}

에 붙는다. 이것이 표면 플라스몬 주파수다. 배경 유전율을 무시한 교과서 표기(ε∞=1\varepsilon_\infty=1)로는 ωsp=ωp/1+εd\omega_{\mathrm{sp}}=\omega_p/\sqrt{1+\varepsilon_d} 이고, 진공 계면이면 그 유명한 ωp/2\omega_p/\sqrt{2} 가 된다. 다만 실제 은·금은 속박전자 배경 때문에 ε∞\varepsilon_\infty 가 5 안팎이라, ωp≈9\omega_p\approx 9 eV를 그대로 2\sqrt{2} 로 나눈 6.4 eV가 아니라 3.5~3.7 eV 근처에서 점근선이 잡힌다. 공식만 외우고 ε∞\varepsilon_\infty 를 빼먹으면 예측이 두 배 틀린다.

4. 결합 방법 — 모자란 운동량을 어디서 빌리나[편집]

Δk=kspp−(ω/c)sin⁡θ\Delta k = k_{\mathrm{spp}} - (\omega/c)\sin\theta 를 메워 주는 장치가 필요하다. 표준 해법은 셋이다.

  • 크레치만(Kretschmann) 배치. 굴절률 npn_p 인 프리즘 바닥에 금속 박막(보통 40~50 nm)을 직접 증착하고, 전반사 조건으로 프리즘 안에서 입사시킨다. 프리즘 쪽 소멸파의 면내 파수가 np(ω/c)sin⁡θn_p (\omega/c)\sin\theta 로 커져 있으므로, 박막 반대편(금속-시료 계면)의 SPP와 정합될 수 있다. 박막이 얇아야 소멸파가 반대편까지 새어 나간다. SPR 센서는 거의 전부 이 배치다.
  • 오토(Otto) 배치. 프리즘과 벌크 금속 사이에 수백 nm 두께의 공기·용액 틈을 둔다. 박막 증착이 필요 없고 금속 표면을 건드리지 않아 표면과학 연구에 좋지만, 틈 간격을 파장 스케일로 정밀 유지해야 해서 상용 기기로는 거의 안 쓴다.
  • 격자 결합. 주기 aa 인 요철을 새기면 역격자 벡터 G=2π/aG=2\pi/a 를 정수배로 빌릴 수 있다. 정합 조건은
kspp=ωcsin⁡θ±m2πak_{\mathrm{spp}} = \frac{\omega}{c}\sin\theta \pm m\frac{2\pi}{a}

이다. 파장·각도 스윕에서 반사율이 급락하는 지점으로 나타나며, 엄밀 결합파 해석이나 FDTD로 그대로 재현된다. 반대로 표면 거칠기·결함은 연속적인 GG 스펙트럼을 공급해 원치 않는 결합과 산란 손실을 만든다.

국소적으로 큰 kk 를 공급하는 근접장 프로브(NSOM 팁)나 나노 슬릿·안테나로 때리는 방법도 같은 원리다. 자유공간 파수 ω/c\omega/c 로는 도달할 수 없는 영역에 접근하려면 결국 파장보다 작은 무언가가 개입해야 한다.

5. 침투 깊이와 전파 거리 — 손실이라는 원죄[편집]

SPP의 실용성은 결국 세 길이가 결정한다. 은 계면 λ=633\lambda=633 nm, εAg≈−18+0.5i\varepsilon_{\mathrm{Ag}}\approx -18+0.5i, 공기(εd=1\varepsilon_d=1) 기준으로 계산하면 감이 온다.

길이식은/공기 @633 nm
유전체 침투깊이1/κd≈ 400 nm
금속 침투깊이1/κm≈ 23 nm
전파거리1/(2 Im kspp)≈ 60 μm

두 침투깊이의 비대칭이 SPP의 성격을 다 말해 준다. 장의 대부분은 유전체 쪽에 있고 금속 안으로는 겨우 20 nm 남짓 들어가는데, 그 20 nm에서 모든 손실이 발생한다. 전파거리 감쇠는 근사적으로

Im kspp≈ωc(εm′εdεm′+εd)3/2εm′′2(εm′)2\mathrm{Im}\,k_{\mathrm{spp}} \approx \frac{\omega}{c}\left(\frac{\varepsilon_m' \varepsilon_d}{\varepsilon_m' + \varepsilon_d}\right)^{3/2}\frac{\varepsilon_m''}{2(\varepsilon_m')^{2}}

로 쓰이고, 손실 항 εm′′\varepsilon_m'' 에 선형이다. 금은 같은 파장에서 ε′′\varepsilon'' 가 은의 두 배 이상이라 전파거리가 10 μm 수준으로 떨어진다. 대신 금은 산화되지 않아 센서용으로는 금이 표준이다 — 물리가 아니라 화학이 재료를 정한 사례.

여기에 근본적인 상충이 있다. 분산 곡선을 ωsp\omega_{\mathrm{sp}} 쪽으로 밀어 kk 를 키울수록 가둠은 좋아지지만, 같은 이유로 장이 금속 안으로 더 들어가 손실이 급증한다. 가둠 × 전파거리는 대략 상수이며, 이 트레이드오프를 우회하려고 나온 것이 장거리 SPP(대칭 박막 모드), 유전체 로딩 도파로, 하이브리드 플라스모닉 도파로 같은 변종들이다. 이득 매질로 손실을 보상하는 스패서(spaser) 계열도 시도되었지만 실용 소자로 자리 잡지는 못했다.

6. 국재 표면 플라스몬 공명[편집]

금속 나노입자는 계면이 닫혀 있어 SPP가 전파하지 못하고 공진한다. 입자가 파장보다 훨씬 작으면 준정전(quasi-static) 근사가 통해 구의 분극률이

α=4πa3 εm−εdεm+2εd\alpha = 4\pi a^{3}\,\frac{\varepsilon_m - \varepsilon_d}{\varepsilon_m + 2\varepsilon_d}

로 나오고, 분모가 0에 가까워지는 프뢸리히 조건 εm=−2εd\varepsilon_m = -2\varepsilon_d 에서 공명이 일어난다. 무손실 드루드 구를 진공에 두면 ω=ωp/3\omega = \omega_p/\sqrt{3} 이다. 평면 계면의 −εd-\varepsilon_d 조건과 구의 −2εd-2\varepsilon_d 조건이 다르다는 점에 주의 — 계수는 기하가 정한다.

입자가 커지면 준정전 근사가 깨지고 미 산란의 엄밀 급수해가 필요하다. 실제로 미가 1908년에 그 이론을 만든 동기가 콜로이드 금 용액의 붉은색이었고, 그 색이 바로 LSPR이다.2 구가 아닌 형상(막대·삼각형·쉘·다이머)은 이산 쌍극자 근사나 T-행렬법, FDTD, 경계요소법으로 푼다. 공명 파장이 종횡비에 민감하게 끌려가므로 형상 설계가 곧 색 설계다.

두 입자를 수 nm 간격으로 붙이면 그 틈(gap)에 장이 극단적으로 몰린다 — 핫스팟. 아래 SERS가 여기 기대고 있다.

7. 응용 — SPR 센서와 SERS[편집]

SPR 바이오센서. 크레치만 배치에서 반사율 최소각 θSPR\theta_{\mathrm{SPR}} 은 금속 표면 위 수백 nm 안의 굴절률에 의해 결정된다. 표면에 항체를 고정하고 시료를 흘리면 결합량에 비례해 국소 굴절률이 오르고, 각도(혹은 공명 파장)가 이동한다. 라벨(형광 표지)이 필요 없고 결합·해리 곡선을 실시간으로 얻어 결합상수 ka,kdk_a, k_d 를 직접 뽑을 수 있다는 것이 결정적 장점이다. 1990년 상용기가 나온 뒤 항체·약물 스크리닝의 표준 장비가 됐다. 감도는 굴절률 단위로 10−610^{-6}~10−710^{-7} RIU 수준. 다만 재는 것은 질량이 아니라 굴절률 변화이므로, 벌크 굴절률 드리프트·온도·비특이 흡착을 참조 채널로 빼 주는 것이 실험의 절반이다.

SERS. 거친 은 전극에서 피리딘의 라만 신호가 비정상적으로 강하다는 1974년 관측에서 시작해, 1977년 두 그룹이 독립적으로 이를 표면 증강으로 해석하면서 표면 증강 라만 산란이라는 분야가 열렸다. 지배적 기구는 전자기적 증강이고 스케일링이 유명하다 — 라만 신호는 입사장과 산란장 양쪽에서 증강을 받으므로

EF∼∣ElocE0∣4\mathrm{EF} \sim \left|\frac{E_{\mathrm{loc}}}{E_0}\right|^{4}

즉 장의 네제곱이다. 핫스팟에서 ∣Eloc/E0∣∼102|E_{\mathrm{loc}}/E_0|\sim 10^2 이면 증강인자가 10810^8 이 된다. 초기 문헌의 101410^{14} 급 단일분자 주장은 과대평가로 정리됐고, 현재 합의는 평균 10410^410610^6, 핫스팟에서 10810^8101010^{10} 정도다. 여기에 화학적 증강(전하 이동)이 한두 자릿수 더 얹힌다. 어느 쪽이든 계산으로 EF를 보고할 때는 어떤 부피에 대해 평균했는지를 밝히지 않으면 숫자가 아무 의미가 없다.

이 밖에 근접장 리소그래피, 열보조 자기기록(HAMR)의 근접장 트랜스듀서, 플라스모닉 광열 치료, 광촉매의 열전자(hot electron) 주입이 모두 같은 물리에서 나온다. 재미있게도 손실이 미덕이 되는 응용(광열·열전자)에서는 플라스모닉스가 오히려 편하게 산다.

8. 수치해석에서 금속을 다룰 때의 함정[편집]

플라스모닉 구조 해석은 전자기해석 중에서도 함정이 유난히 많다. 원인은 셋이다 — 재료가 분산성이고, 특성 길이가 파장보다 두 자릿수 작고, 답이 공명이라 조금만 틀려도 크게 틀린다.

  • 분산 모형은 필수다. 시간영역 코드에 εr\varepsilon_r 을 음의 상수로 넣으면 갱신식이 즉시 발산한다. 비분산 음수 유전율은 비인과적이기 때문이다. 드루드 항 −ωp2/(ω2+iγω)-\omega_p^2/(\omega^2+i\gamma\omega) 을 보조미분방정식(ADE)이나 재귀합성곱(RC/PLRC)으로 넣는 것이 표준.
  • 드루드 하나로는 부족하다. 금은 약 2.4 eV(≈520 nm), 은은 약 3.8 eV 위에서 d밴드 띠간 전이가 시작되므로, 가시광 전체를 커버하려면 드루드 + 로렌츠 극점 1~2개(혹은 임계점 모형)를 얹어야 한다. 한 극점 드루드만 쓰고 청색 쪽 결과를 보고하는 것은 흔한 사고다.
  • 재료 데이터 출처를 반드시 밝힌다. 존슨-크리스티(Johnson & Christy, 1972)와 팔릭(Palik) 핸드북은 같은 금속에 대해 눈에 띄게 다른 ε′′\varepsilon'' 를 준다. 존슨-크리스티 쪽이 손실이 작아 공명이 날카롭게 나오고, 대체로 증착 박막 실측과 더 잘 맞는다는 평가가 많다. 어느 데이터를 썼는지에 따라 Q값과 증강인자가 배수로 갈리므로, 데이터 출처 없는 플라스모닉 논문은 재현이 불가능하다. 게다가 실제 박막의 광학상수는 증착 조건·결정립·표면 거칠기에 따라 문헌값과 또 다르다.
  • 격자는 침투깊이 기준으로 잡는다. 파장 기준 λ/20\lambda/20 같은 관행값은 여기서 무의미하다. 금속 침투깊이가 20 nm대이고 갭 핫스팟은 1~2 nm이므로, 그 영역에서 셀 크기가 나노미터급이어야 한다. 균일 격자로 하면 셀 수가 폭발하니 부분격자(subgridding)나 비정렬 유한요소 전자기로 간다.
  • 곡면 계단화가 특히 아프다. 나노구·나노막대의 Yee 격자 계단 근사는 표면 전하 분포를 왜곡해 공명 파장을 통째로 옮긴다. 부분셀 유전율 평균(subpixel smoothing)을 쓰면 크게 나아진다.
  • 경계·수렴 검사. 표면파가 완전정합층에 비스듬히 진입하면 흡수되지 않고 터지는 사고가 잦다. 계면을 PML 안까지 곧게 연장하거나 CFS/CPML 계열을 쓴다. 공명 QQ 가 높아 정상상태 도달에 시간스텝이 오래 걸리므로, 시간영역이면 신호를 충분히 길게 받고 필요하면 필터 대각화로 극점을 뽑는다.
  • 국소 모형의 유효 범위. 갭이 1 nm 이하로 내려가면 전자 스필아웃과 터널링 때문에 국소 드루드 응답이 무너진다. 비국소(유체역학) 보정이나 양자 보정 유전율을 넣으면 예측 증강인자가 눈에 띄게 낮아지는데, 그쪽이 실측에 가깝다. 고전 계산만 돌리고 “증강 101210^{12}“라고 쓰면 과장이다.
  • 주파수영역 FEM의 별도 함정. 계면에서 ε\varepsilon 부호가 바뀌면 변분 형식의 강제성(coercivity)이 깨져, 특정 유전율 비에서 이산해가 수렴하지 않는 구간이 존재한다. 슈퍼렌즈에서와 같은 문제다.

검증 순서는 단순하게 잡는 것이 좋다. 평면 계면의 SPP 분산 관계를 수치로 뽑아 위 해석식과 맞추고, 그다음 미 이론이 답을 아는 구 하나를 풀고, 그러고 나서야 실제 구조로 넘어간다. 검증 및 확인의 원칙이 이 분야만큼 뼈아프게 적용되는 곳도 드물다.3

9. 관련 문서[편집]

10. Footnotes[편집]

  1. 용어가 은근히 지저분하다. 엄밀히 플라스몬은 전자 밀도 진동의 양자, 폴라리톤은 그것이 광자와 섞인 결합 모드다. 그래서 계면을 달리는 놈은 SPP라고 불러야 맞지만, 현장에서는 그냥 “surface plasmon”으로 뭉뚱그린다. 논문 심사에서 이걸 지적하는 심사위원과 안 하는 심사위원의 비율은 대략 반반이다. ↩

  2. 고대 로마의 리쿠르구스 컵이나 중세 성당의 스테인드글라스 붉은색도 같은 원리다. 금 나노입자 콜로이드의 LSPR인데, 만든 사람들은 유전율의 부호 같은 건 당연히 몰랐다. 1600년 넘게 “레시피”로만 전승된 플라스모닉스인 셈. ↩

  3. 이 바닥에서 “우리 시뮬레이션은 실험과 잘 맞습니다”는 문장을 볼 때 확인할 것: 광학상수 출처, 격자 수렴 검사, 그리고 증강인자를 어떤 부피로 평균했는지. 셋 중 하나라도 없으면 그 그림은 물리가 아니라 메시를 그린 것일 수 있다. ↩