금속 내 종방향 전자기 진동: 평면 및 구형 계면 위 표면 여기

초록

공간 분산 유전율의 통상적 정의를 수정했습니다. 벡터 헬름홀츠 방정식을 만족하는 전기장에 대한 축퇴 전자 기체 응답을 볼츠만 방정식 해법으로 구했습니다. 계산된 종방향 유전율 함수는 1952년 Klimontovich와 Silin, 1954년 Lindhard 결과와 일치합니다. 그러나 평면 전자기파 파동 벡터가 아닌 벡터 헬름홀츠 방정식 파라미터인 파동수 제곱에 의존합니다. 이 새 개념은 비국소 효과 시뮬레이션을 단순화합니다(예: 일반화된 로렌츠-미 이론, 푸리에 변환 불필요). 프레넬 계수 일반화로 종방향 전자기파 여기 허용. 이론 검증 위해 은·금 나노미터 구형 입자 소광 스펙트럼 계산. 이 입자들에 대해 일반화된 로렌츠-미 이론이 실험 데이터와 훌륭히 일치하는 청색 이동과 플라즈몬 공명 브로드닝 예측. 추가로 금 나노구체 직경 ≤2 nm에서 플라즈몬 공명 소멸 현상 비국소 이론으로 설명. 은에 대한 Klimontovich-Silin-Lindhard와 유체역학 유전율 함수 계산이 3-4 eV 광자 에너지 영역서 근접 결과. 고체 내 종방향 파 파동수 절댓값이 횡파보다 훨씬 큼 입증.

배경

평면 금속 표면 펨토초 레이저 펄스 조사 시 레이저 유도 주기적 표면 구조(LIPSS) 형성 흔함[1]. 일반 LIPSS 외 고공간주파수 LIPSS(HSFL) 관찰[1,2]. HSFL 공간 주기는 조사 파장 λ₀보다 현저히 작음. 예: 알루미늄 λ₀=0.8 µm서 20-200 nm 범위 추정[2,3]. 일반 LIPSS 파문 방향 레이저 편광 수직이나 HSFL은 수직이거나 때로 평행. 투명 유전체, 반도체, 금속서 유사 HSFL 형성. HSFL 기원 이차 고조파 생성, 특정 플라즈몬 모드 관여, 자기조직화, 유전체 불균일 파괴 시 국소장 증강 등 다양한 메커니즘 설명[2,3].

본 연구 목표: Λ≪λ₀ 짧은 주기 패턴 생성 파동 과정 탐색. 금속 내 종방향(L) 전자기파(플라즈마파) 특성 검토. 연구 단계: 첫째, 유전율 공간 분산 정의부터 시작. 통상 정의는 연구 매질이 균일 무한하지 않으면 무용함. 따라서 새로운 공간 분산 유전율 함수 ε 개념 제안. 이 함수는 두 벡터장 E(r,ω), D(r,ω) 간 직접 비례 관계 수립, 평면파 진폭 E(k,ω), D(k,ω) 아님. 결과적으로 ε은 평면파 파동 벡터 k 아닌 벡터 헬름홀츠 방정식 파라미터인 파동수 제곱 k² 의존. 둘째, 볼츠만 수송 방정식(완화시간 근사) 풀어 전도 전자 전자기 모드 응답 구해 횡방향/종방향 Lindhard 유전율 함수 도출. 종방향 Lindhard와 더 간단한 유체역학 함수가 넓은 파라미터 범위서 근접함 확인. 은·금 나노구체 광 소광으로 이론 예시. 비국소 미 이론이 귀금속 나노구체 크기 감소 시 플라즈몬 공명 청색 이동, 브로드닝, 최종 소멸 설명 최초 입증. 최종 개발 이론 모델로 종방향 모드 레이저 유도 표면 구조 형성 관여 가능성 검토. 전달 종방향파 고려해 프레넬 이론 수정.

방법

조각wise 균질 매질 내 전자기장 결정을 위해 고전 전자기 이론 적용. 매질 각 균질 영역 내 전기장 E를 벡터 헬름홀츠 방정식(VHE) 해로 가정:

$$ \Delta\,\mathbf{E} + k^{2}\, \mathbf{E}=0, $$ (1)

여기서 Δ는 라플라스 연산자.

통상적 전기장 E, 자기장 H 접선 성분 매질 경계면 연속. 추가로 전자 금속 내 구속 고려해 금속 표면 S에서 전류 밀도 j 법선 성분 추가 경계 조건(ABC) 사용: (j·n)|_{r∈S}=0.

금속 내 전도 전류 결정을 위해 볼츠만 수송 방정식(BTE) 완화시간 근사형으로 풀기:

$$ \frac{\partial f}{\partial t}+\mathbf{v}\,\frac{\partial f}{\partial \mathbf{r}}+ \frac{e}{m}\,\left(\mathbf{E}+\mathbf{v}\times\mathbf{B} \right)\,\frac{\partial f}{\partial \mathbf{v}} = \frac{f_{0}-f}{\tau}, $$ (2)

여기서 f 위상공간(r,v) 단입자 분포함수, v 미시적 전자 속도, e,m 전자 전하·질량, B 자기 유도, f₀ 평형 분포함수, τ 완화시간.

아래 공간 분산 유전율 함수 유도 후 평면 금속 표면 빛 반사와 귀금속 나노구체 빛 산란 연구에 활용.

결과 및 고찰

이질 매질 내 ε 공간 분산

문헌서 공간 분산 유전율 ε 정의는 다음 관계식[4-6]:

$$ \mathbf{D}(\omega,\, \mathbf{r}) = \epsilon_{0} \iiint\limits_{-\infty}^{\infty} \! \mathbf{d} \mathbf{r}^{\prime}\, \epsilon\left(\omega, \, \mathbf{r}-\mathbf{r}^{\prime}\right)\,\mathbf{E}\left(\omega,\, \mathbf{r}^{\prime}\right), $$ (3)

여기서 ε₀ 전기 상수, D(ω,r) 점 r 각주파수 ω 변위 벡터 진폭, E(ω,r') 점 r' 전기장 진폭. 식 (3) 푸리에 변환으로

$$ \mathbf{D}(\omega,\, \mathbf{k}) = \epsilon_{0} \, \epsilon(\omega,\, \mathbf{k})\,\mathbf{E}(\omega,\, \mathbf{k}) $$ (4)

평면 전자기파 파동 벡터 k 의존 공간 분산 ε(ω,k) 얻음. 우리 의견으론 식 (3) 무한 균일 부피서만 모호하지 않음, 경계 고려해야 하는 조각wise 이질 시스템서 k 각 매질서 다름.

우리 접근법 평면파 전개 안 함. 공간 분산 투자율은 D(ω,r)과 벡터 헬름홀츠 식 (1) 특정 해 간 관계 결정:

$$ \mathbf{D}(\omega,\, \mathbf{r}) = \epsilon_{0} \, \epsilon(\omega, \, k)\,\mathbf{E}(\omega,\, \mathbf{r}). $$ (5)

여기서 E(ω,r) 점 r 벡터 E 아닌 전기장 분포 의미.

종방향 및 횡방향 유전율 함수

금속 투자율 통상 전도도 σ로 표현[4]:

$$ \epsilon=\epsilon_{\mathrm{g}}+\frac{i\,\sigma}{\omega\,\epsilon_{0}}, $$ (6)

ε_g 고체 분극 허용 유전율 부분, 단순 금속 ε_g=1. σ 결정을 위해 전류 밀도 계산:

$$ \mathbf{j}= e \iiint\limits_{-\infty}^{\infty} \! \mathbf{v}\, f\, \mathrm{d}\/ \mathbf{v}=\sigma\,\mathbf{E}, $$ (7)

여기서 dv = v/m dε dΩ, dΩ=sinθ dθ dϕ, v,θ,ϕ 속도 구면좌표. 선행 연구와 달리 파동 벡터 k 도입 안 하고 v∇ 연산자 vE 작용 무한 급수 형태 BTE 해 구함:

$$ f=f_{0} + \frac{e}{-i \omega+\Gamma}\,\frac{\partial f_{0}}{\partial \epsilon} \left[ 1+\frac{\mathbf{v}\,{\mathbf{\nabla}}}{-i \omega+\Gamma} \right]^{-\,1}\! \mathbf{v}\,\mathbf{E}, $$ (8)

Γ=1/τ. f₀ 영온 페르미-디락 분포 근사, 식 (7) ε 적분 후

$$ \mathbf{j}=\frac{\omega_{\mathrm{p}}^{2}\,\varepsilon_{0}}{-i \omega+\Gamma}\,\frac{3}{4 \pi} \iint \mathbf{u}\,\left(1+ l\,\mathbf{u} \mathbf{\nabla} \right)^{-\,1} (\mathbf{u}\,\mathbf{E})\, \mathrm{d}\,\Omega, $$ (9)

ω_p² = e²n_e/mε₀, ω_p 플라즈마 주파수, u=v/v v 방향 단위 벡터, l=v_F/(-iω+Γ), v_F 페르미 속도. 추가 적분 계산:

$$\begin{array}{*{20}l} &\textstyle \iint \mathbf{u}\, (\mathbf{u}\,\mathbf{E})\,\mathbf{d}\,\Omega = \frac{4\/\pi}{3}\,\mathbf{E} \end{array} $$ (10)

$$\begin{array}{*{20}l} &\textstyle \iint \mathbf{u}\,(\mathbf{u}\,\mathbf{\nabla})^{2 n-1} (\mathbf{u}\,\mathbf{E})\,\mathbf{d}\,\Omega = 0 \end{array} $$ (11)

$$\begin{array}{*{20}l} &\textstyle \iint \mathbf{u}\,(\mathbf{u}\,\mathbf{\nabla})^{2 n} (\mathbf{u}\,\mathbf{E})\,\mathbf{d}\,\Omega = \frac{4 \pi}{2 n+3} \\ &\times \, \Delta^{n-1} \left[ \mathbf{\nabla}\, (\mathbf{\nabla}\cdot\mathbf{E})-\frac{1}{2 n+1}\, \mathbf{\nabla}\times \mathbf{\nabla}\times\mathbf{E}\right] \end{array} $$ (12)

n 자연수. 임의 전기장 E에 대한 j 의존성 최종 얻음:

$$\begin{array}{*{20}l} \mathbf{j}&=\frac{\omega_{\mathrm{p}}^{2}\,\varepsilon_{0}}{-i \omega+\Gamma} \left\{\mathbf{E} + 3\, \sum\limits_{n=1}^{\infty} l^{\,2 n} \right. \\ &\quad\left.\times \frac{\Delta^{n-1}}{2 n+3} \left[ \mathbf{\nabla}\, (\mathbf{\nabla}\cdot\mathbf{E})-\frac{\mathbf{\nabla}\times \mathbf{\nabla}\times\mathbf{E}}{2 n+1}\right]\right\}. \end{array} $$ (13)

식 (1) 해 두 유형: 발산 자유 ∇·E=0 만족, 회전 자유 ∇×E=0 만족.

평면파 E∝exp[i(k·r-ωt)]서 식 (14) k×E=0 변환, 종파(L)임. 구형 물체 과정 시뮬레이션 편의상 벡터 구면 조화 L, M, N 직교 완비계 사용. 이 경우 식 (14) L 조화 지정. L파와 L 모드 파동수 다음 분산 법칙 결정:

$$ \epsilon^{\mathrm{L}}\left(\omega, \, k^{\mathrm{L}}\right)=0. $$ (15)

식 (6), (13)서 식 (14) 제약 만족 식 (1) 해가 다음 종방향 투자율 줌:

$$ \epsilon^{\mathrm{L}}=\epsilon_{\mathrm{g}}- \frac{\omega_{\mathrm{p}}^{2}}{\omega\,(\omega+i \Gamma)}\,\frac{3}{2} \, \Phi \left(a^{2},\, 1,\,\frac{3}{2} \right) $$ (16)

Φ 레르히 피 함수,

$$ \frac{3}{2} \, \Phi \left(a^{2},\, 1,\,\frac{3}{2} \right) =\sum\limits_{n=0}^{\infty} \frac{3}{2 n+3}\, a^{2\,n}, $$ (17)

a = k v_F/(ω+iΓ).

얻은 투자율 Kliewer와 Fuchs[7] 것과 표기만 다름:

$$ \epsilon^{\mathrm{L}}=\epsilon_{\mathrm{g}}+ \frac{\omega_{\mathrm{p}}^{2}}{\omega\,(\omega+i \Gamma)}\,\frac{3}{a^{2}}\left[1-\frac{1}{i a} \tan^{-1}(i a) \right] $$ (18)

항등식

$$ \frac{1}{i a}\tan^{-1}(i a)=\frac{1}{2}\ln\frac{1+a}{1-a} $$ (19)

식 (18) 재작성:

$$ \epsilon^{\mathrm{L}}=\epsilon_{\mathrm{g}}- \frac{\omega_{\mathrm{p}}^{2}}{\omega\,(\omega+i \Gamma)}\,\frac{3}{a^{2}}\left[1-\frac{1}{2 a} \, \frac{\ln(1+a)}{\ln(1-a)} \right]. $$ (20)

Γ=0 경우 Klimontovich와 Silin[8] 유도 식 형태(참조 [9], [10, 식 (40.17)], [11]). 식 (16), (18), (20) 동등 투자율 통상 Lindhard 유전율 함수라 칭함([12] 인용)지만 최초 Klimontovich와 Silin[8] 유도.

횡방향 Lindhard 투자율[7] ∇·E=0 조건 식 (13)서 구함. 실제 v_F k≪ω 경우 드루드 투자율로 환원:

$$ \epsilon^{\mathrm{T}}=\epsilon_{\mathrm{g}}- \frac{\omega_{\mathrm{p}}^{2}}{\omega^{2}+i\,\Gamma\/\omega}. $$ (21)

이 함수 다수 금속 실험 데이터와 일치[13]. |a|<1 경우 종방향 투자율 (16) 유체역학 유전율 함수로 단순화:

$$ \epsilon^{\mathrm{L}}\left(\omega,\,k^{\mathrm{L}}\right)=\epsilon_{\mathrm{g}}- \frac{\omega_{\mathrm{p}}^{2}}{\omega^{2}+i\,\Gamma \omega-\frac{3}{5}\,\left(v_{\mathrm{F}}\, k^{\mathrm{L}}\right)^{2}}. $$ (22)

평면 금속 표면 평면 전자기파 반사

경계 조건

평면 전자기파 평면 금속 표면 반사 시 금속 내 여기되는 L파 파동 벡터 k^L 방향과 진폭 결정.

파동 벡터 xz평면 내 z=0 유전체-금속 계면 입사 평면파 고려. 유전체 매질 1 전기장 입사 E_i, 반사 E_r파, 금속 2 전달 횡파 E_t, 경우 종방향 E^L 성분. 맥스웰 경계 조건 따라 전기장, 자기장 횡성분 z=0 평면 연속. 추가 전자 금속서 방출 안 됨, 전류 밀도 법선 성분 z=0서 0:

$$ \hat{\mathbf{z}}\,\mathbf{j}|_{z=0}=0. $$ (23)

ẑ z축 방향 단위 벡터.

맥스웰 경계 조건 모든 항 x,y 의존성 동일해야 함. 이 요건 여러 결과. 첫째, L파 여기 가능 조건: p-편광 입사파 전기 벡터 E^(i) 입사면 평행 경우만. 즉 플라즈몬 폴라리톤 횡자기(TM)파로 생성 가능. 금속 구서도 유사 효과[14]. 둘째, 스넬 법칙 유사 공식 조건서 유도:

$$ k_{1x}=k_{2x}=k_{2x}^{\mathrm{L}}=k_{1}\,\sin\theta_{1} $$ (24)

1x,2x 매질 1,2 벡터 x-투영, θ₁ 입사각.

반사 및 투과 계수

평면 p-편광 전자기파 평면 금속 표면 입사 형성 장 결정. 입사파 E_x^(i) x성분 표현 편리: 반사파 E_x(r)=-r E_x(i), 전달 횡파 E_x(t)=t E_x(i), 전달 종파

$$ E^{\mathrm{(a)}}_{x}=\delta\,E^{\mathrm{(t)}}_{x}=t_{\mathrm{L}}\,E^{\mathrm{(i)}}_{x} $$ (25)

r 반사 계수, t, t_L 투과 계수.

맥스웰 경계 조건과 식 (23) ABC 다음 형태로부터:

$$ \hat{\mathbf{z}}\,(\mathbf{D}-\epsilon_{0}\epsilon_{\mathrm{g}}\,\mathbf{E})|_{z=0}=0, $$ (26)

다음을 얻음:

$$\begin{array}{*{20}l} r&=- \frac{(1+\delta)\,\epsilon_{1}\,k_{2z}-\epsilon_{2}\,k_{1z}}{(1+\delta)\,\epsilon_{1}\,k_{2z}+\epsilon_{2}\,k_{1z}}=1-(1+\delta)\,t \end{array} $$ (27)

$$\begin{array}{*{20}l} t&=\frac{2\,\epsilon_{1}\,k_{2z}}{\epsilon_{2}\,k_{1z}+(1+\delta)\,\epsilon_{1}\,k_{2z}}, \end{array} $$ (28)

$$\begin{array}{*{20}l} \delta&=\frac{\epsilon_{\mathrm{g}}-\epsilon}{\epsilon_{\mathrm{g}}}\,\frac{k_{2x}^{2}}{k_{2z}\,k_{2z}^{\mathrm{L}}} \end{array} $$ (29)

δ=0 시 r은 p-편광 파 프레넬 반사 계수(예: [4] 식 2.49). 동일 조건서 t 프레넬 투과 계수 아님(우리 t,r 정의 프레넬과 다름).

금속 나노구체 빛 소광

선행 논문서 저자 중 한 명 로렌츠-미 이론 일반화해 식 (23) ABC 허용. 반사 TM 모드 l차 미 계수 b_l 유도:

$$ b_{l}=- \frac{(1+\delta_{l})\,\epsilon_{1}\, \frac{k_{2}\,\psi_{l}^{\prime}(k_{2} R)}{\psi_{l}(k_{2} R)} - \epsilon_{2}\, \frac{k_{1}\,\psi_{l}^{\prime}(k_{1} R)}{\psi_{l}(k_{1} R)}}{(1+\delta_{l})\,\epsilon_{1}\, \frac{k_{2}\,\psi_{l}^{\prime}(k_{2} R)}{\psi_{l}(k_{2} R)} - \epsilon_{2}\, \frac{k_{1}\,\zeta_{l}^{\prime}(k_{1} R)}{\zeta_{l}(k_{1} R)}}, $$ (30)

여기서

$$ \delta_{l}=\frac{\epsilon^{\mathrm{T}}- \epsilon_{\mathrm{g}}}{\epsilon_{\mathrm{g}}}\, \frac{l\,(l+1)\,j_{l}(k_{2} R)\,j_{l}(k_{2}^{\mathrm{L}} R)}{\psi_{l}^{\prime}(k_{2} R)\,k_{2}^{\mathrm{L}} R \,j_{l}^{\prime}\left(k_{2}^{\mathrm{L}} R\right)}, $$ (31)

ψ_l, ζ_l 각 l차 리카티-베셀, 리카티-한켈 함수, j_l 구면 베셀 함수, 프라임 인수 미분.

고전·일반화 로렌츠-미 이론 예측 실험 데이터 비교. [15] Hilger, Tenfelde, Kreibig 유전체 표면 증착 은 나노입자 소광 스펙트럼 연구. 1단계서 평균 직경 2, 3.5, 4 nm 은 입자 빔 생성, 한 빔 입자 크기 분포 결정, 소광 스펙트럼 기록, 진공 은 구체 현상적 공식 Γ=Γ_b+A v_F/R 파라미터 A=0.25 추정. 먼저 평균 직경 〈D〉=2 nm, 실험 크기 분포(D=1-4 nm) 은 구체 빔 소광 스펙트럼 계산. 우리 이론 조정 파라미터 없음. 유전율 함수 정의 위해 Lynch와 Hunter[16] 벌크 은 굴절률 표 사용(Fig. 1 참조). 식 (16), (21), (22) ω_p=9.17 eV, Γ_b=0.021 eV, v_F=1.39×10⁶ m/s, A=0.25 적용. 계산 결과와 실험 스펙트럼 Fig. 2 제시.

그림 1. Johnson and Christy(◇)[20], Lynch and Hunter(+)[16], Weber(△)[21], Hao and Nordlander(파선)[22], Drachev et al.(실선)[23] 은 유전율 함수 실수부(a), 허수부(b)

그림 2. [15] 관찰 은 나노미터 입자 빛 소광 스펙트럼과 국소·비국소 모델 계산. 모든 이론 스펙트럼 공통 상대 단위

Fig. 2 이론 스펙트럼 Klimontovich-Silin-Lindhard와 더 간단한 유체역학 유전율 함수로 계산. 놀라운 점은 |a|>1 플라즈몬 공명 영역서도 두 계산 근접 결과.

나노미터 은 구체 소광 스펙트럼 최대(프뢸리히[17], 플라즈몬, 표면 플라즈몬 폴라리톤(SPP)[15] 공명) 3.5→3.65 eV 청색 이동 알려짐[18]. 비국소 모델 실험 데이터와 훌륭히 일치, 국소(미) 이론 ω≃3.5 eV 최대(Fig. 2, Table 1).

플라즈몬 공명 청색 이동 계산 다음 고려로 뒷받침. 정전기 근사서 b_1만 소광 단면 Q_ext 기여, 식 (30) 다음 근사로 단순화:

$$ \frac{k_{2} R\,\psi_{l}^{\prime}(k_{2} R)}{\psi_{l}(k_{2} R)}\simeq l+1;\,\,\, \frac{k_{1} R\,\zeta_{l}^{\prime}(k_{1} R)}{\zeta_{l}(k_{1} R)} \simeq -\,l. $$ (32)

Q_ext 최대 조건:

$$ \Re [2\,(1+\delta_{1})\,\epsilon_{1}+ \epsilon_{2}]=0. $$ (33)

조건 (33) L 모드 여기 고려(δ_1 항), 프뢸리히 공명 조건[17]과 다름:

$$ \Re (2\,\epsilon_{1}+ \epsilon_{2})=0. $$ (34)

실험[15]서 소광 스펙트럼 피크 주파수 ω_m, 공명 폭 Δω 거의 〈D〉 독립. 이 Δω 특징 고전 미 이론과 상충. 국소 이론 D 감소시 플라즈몬 공명 브로드닝 예측(A=0.25) Table 1. 반면 비국소 이론 대략 동일 공명 폭, 다른 피크 위치. 전 입자 기여 중첩 Δω 값 실험 데이터와 놀라운 일치. 흥미롭게 비국소 이론 A=0에서도 빔 플라즈몬 공명 브로드닝 예측.

ω>4 eV서 Fig. 2 매끄러운 이론 곡선 좁은 밀집 실험 피크 매시보다 높음. 이 영역 밴드간 흡수 우세, Fig. 1로 확인. 관측 스펙트럼 특이성 연속 밴드→이산 준위 구조 전이 결과 가능. 양자 크기 효과 금 나노구체 광학 특성 연구서 앞서 발견[19]. 은 구체 크기 〈D〉=3.5 nm 증가 시 흡수 최대 대비 최초 증가 후 일련 작은 등간격 딥 가진 고원 형성. 〈D〉=4 nm서 흡수 약간 감소.

플라즈몬 공명 청색 날개 형성 연구 위해 초미세 은 입자 소광 스펙트럼 계산 Fig. 3 제시. Fig. 3 주목할 특징: D=1 nm서 플라즈몬 공명 완전 소멸. 이 효과 앞서 금 나노구체 실험 연구서 관찰[19]. 특히 [19] Fig. 9서 1.7, 1.9, 2.0, 2.1, 2.3, 2.5 nm 직경 입자 실험 스펙트럼 국소 미 이론 계산과 비교. 일치도 낮아 플라즈몬 공명 브로드닝과 위치 기술 실패[19]. 입자 크기 변화와 유전율 함수 수정으로 피팅 개선 시도 실패. [19] 저자들 비정상적 넓거나 억제된 집단 진동 밴드 제안된 국소 미 이론 수정으로 피팅 저항. Fig. 4서 비국소 미 이론 적용 시 상황 극적 변화. 조정 파라미터 없음. Johnson and Christy[20] 복소 굴절률 표 금 유전율 함수용. A=1, 톨루엔 굴절률 1.37 등 기타 파라미터 [19]서 채택.

그림 3. 직경 2.2, 1.8, 1.4, 1.0 nm 은 입자 정규화 소광 단면 국소(파선), 비국소(실선) 미 이론 계산. 입자 작을수록 곡선 낮음. 모든 이론 단면 공통 상대 단위

그림 4. 국소(파선), 비국소(실선) 미 이론 계산 흡수 스펙트럼과 [19] Fig. 9 추출 실험 데이터(점) 톨루엔 중 금 구체 D=2.5, 2.1, 1.7 nm(곡선 1·원, 곡선 2·사각형, 곡선 3·삼각형). 모든 이론 스펙트럼 4.12 eV서 1로 정규화, 수직 변위

종방향파 파동수

종방향 모드 횡파 대비 훨씬 높은 파동수 값 특징. 예: Fig. 2 계산서 Re(k₂^L) 해당 공간 주기 Λ=2π/Re k₂^L ω 3→4 eV 증가 시 9→2 nm 감소. 이 ω 구간 |k₂^L/k₂| 절댓값 130→100 감소, 식 (27) δ 파라미터 θ₁=π/4서 0.01→0.005 감소. 평면 은 표면 L파 여기 무시 가능 결론. 단 L 모드 나노미터 은 클러스터서 중요함 발견.

식 (16) -ω_p²/(ω²+iΓω) 항 식 (21) ε^T-ε_g로 치환해 분산 식 (15) 재작성:

$$ 1+\frac{3}{5}\,a^{2}+\frac{3}{7}\,a^{4}+\frac{3}{9}\,a^{6}+\dots= \frac{1}{1-\epsilon^{\mathrm{T}}/\epsilon_{\mathrm{g}}}. $$ (35)

ε_g=1, Γ=0 최단순 경우 식 (35) 금속 ω>ω_p서 횡파·종파 모두 투명, ω<ω_p서 k^L, k^T 모두 복소수 예측.

고체 투명 시 p-편광파 경사 입사 평면 표면서 종파 여기 가능. 이 효과 몇 특징: 첫째, 종파 평면 표면서 생성 가능, 표면 플라즈몬 폴라리톤 여기 위해 특수 노력 필요[4,5]와 대조. 둘째, 간섭 무늬서 전자기장 강도 계면 평행 아닌 수직 방향 변조. 따라서 고체 박리 기인 계면 평행 평면서 보이드 나타날 수 있음. ω_p 정의 따라 ω>ω_p 조건 저전류 캐리어 밀도 고체(예: 반도체)서 충족 가능. 여기선 축퇴 전자 기체 유도 ε^L 공식이라 이 경우 미검토.

결론

공간 의존 유전율 정의 위해 선행 연구자들 모두 평면 전자기파와 물질 상호작용 고려. 이 접근 나노광학서 장 공진기 내 국소화, 경계 조건 고려 필요시 비건설적·비엄밀. 벡터 헬름홀츠 방정식 만족 전기장 매질 응답 계산로 해결. 유도 공간 분산 유전율 함수는 헬름홀츠 방정식 파라미터 파동수 제곱 의존, 평면파 파동 벡터 비의존.

종방향파 여기 고려 수정된 프레넬 반사 계수 보고. 미 계수 유사 일반화 선행. 본 이론 나노미터 은·금 클러스터 빛 소광 시뮬레이션 검증. 은 입자 빔 표면 플라즈몬 공명 3.5→3.65 eV 이동과 폭 실험 데이터 훌륭히 일치. 추가 금 구체 직경 ~2 nm 플라즈몬 공명 소멸 비국소 모델 설명. 중요: L파 평면 입사파로 평면 표면서 여기 가능. 표면 플라즈몬 폴라리톤과 주요 차이.

금속 내 전자기 진동 특성 검토. 종방향파 파동수 절댓값 횡파보다 훨씬 큼 발견. 예: 은 3.5 eV 광자 에너지서 파동수 절댓값 비 130. 종방향파 실수부 파동수 7 nm 파장 대응. 파동수 큰 차이 평면 표면 L파 여기 방지. 단 L 모드 은·금 나노미터 입자서 여기됨 보임.