金属における縦電磁振動：平面および球面界面での表面励起

概要

空間分散誘電率の一般的定義を改訂。ベクトルヘルムホルツ方程式満たす電場に対する縮退電子ガス応答をボルツマン方程式解で導出。得られた縦誘電関数はKlimontovichとSilin（1952）およびLindhard（1954）結果と一致。しかし、平面電磁波波数ベクトルではなく、ベクトルヘルムホルツ方程式パラメータである波数二乗に依存。この新概念は非局所効果シミュレーション簡素化（一般化Lorentz-Mie理論等、フーリエ変換不要）。フレネル係数一般化し縦電磁波励起許容。理論検証に銀・金ナノメートル球消光スペクトル計算。これら粒子で一般化Lorentz-Mie理論はプラズモン共鳴青方シフトと拡大予測、実験データと優秀一致。さらに非局所理論で直径≤2 nm金球プラズモン共鳴消失説明。銀用Klimontovich-Silin-Lindhardと流体力学誘電関数計算、光子エネルギー3-4 eVで近接結果。固体中縦波波数絶対値が横波より桁違い大きいこと示す。

背景

フェムト秒レーザーパルス平面金属表面照射でレーザー誘起周期表面構造（LIPSS）形成知られる。LIPSS以外に超微細リップル、高空間周波数LIPSS（HSFL）観測。HSFL空間周期は照射波長λ₀より著しく小さい。例：アルミニウムでλ₀=0.8 µm時20-200 nm範囲推定。通常LIPSSリップル配向レーザー偏光垂直、HSFL配向しばしば垂直、時々平行。類似HSFL透明誘電体、半導体、金属表面形成。HSFL起源：第二高調波生成、特定プラズモンモード関与、自己組織化、誘電体不均一破壊時局所場増強等諸説。

本研究目標：周期Λ≪λ₀パターン生成波動プロセス探索。金属中縦（L）電磁波（プラズマ波）特性検討。新規ステップ：空間分散誘電率定義から着手。媒質均一無限大でない場合一般定義無用。新空間分散誘電関数ε提案：2ベクトル場E(r,ω)とD(r,ω)直接比例関係確立、平面波振幅E(k,ω)とD(k,ω)でない。従いεは波数二乗k²（電場E(r,ω)ベクトルヘルムホルツ方程式パラメータ）依存、平面波波数ベクトルk非依存。導出に緩和時間近似ボルツマン輸送方程式解き、電磁モード伝導電子応答決定。横・縦Lindhard誘電関数獲得。縦Lindhardと単純流体力学関数広パラメータ範囲で近接確認。銀・金ナノ球光消光理論説明用例。非局所ミー理論で貴金属ナノ球サイズ減少時プラズモン共鳴青方シフト、拡大、最終消失説明初実証。新理論モデルレーザー誘起表面構造形成へ縦モード関与可能性検証応用：透過縦波考慮しフレネル理論修正。

方法

区分均質媒質電磁場決定に古典電磁気理論適用。均質領域電場Eをベクトルヘルムホルツ方程式（VHE）解と仮定：ΔE + k²E = 0（式1）。通常通り電場Eと磁場H接線成分媒質境界連続。電子金属内閉込め考慮、金属表面S法線方向電流密度j追加境界条件（ABC）：(j·n)|_r∈S = 0。

金属伝導電流決定に緩和時間近似ボルツマン輸送方程式（BTE）解く：∂f/∂t + v·∂f/∂r + (e/m)(E+v×B)·∂f/∂v = (f₀-f)/τ（式2）。f位相空間(r,v)単粒子分布関数、v微視的電子速度、e,m電子電荷質量、B磁束密度、f₀平衡分布、τ緩和時間。

以下空間分散誘電関数導出、平面金属表面光反射、貴金属ナノ球光散乱研究に利用。

結果と考察

不均質媒質中ε空間分散

文献空間分散誘電関数ε定義：D(ω,r)=ε₀∫ε(ω,r-r')E(ω,r')dr'（式3）。フーリエ変換でD(ω,k)=ε₀ε(ω,k)E(ω,k)（式4）、空間分散ε(ω,k)平面波波数ベクトルk依存。我々見解：式3無限均質体積でのみ明確、境界考慮区分不均質系ではk媒質毎異なり不適。我々アプローチ：平面波展開不使用。空間分散誘電率ベクトルヘルムホルツ式(1)特定解とD(ω,r)関係規定：D(ω,r)=ε₀ε(ω,k)E(ω,r)（式5）。E(ω,r)電場分布、単なる点rベクトルE非。

縦横誘電関数

金属誘電率導電率σ通じ表現：ε=ε_g + iσ/(ωε₀)（式6）。ε_g固体分極寄与、単純金属ε_g=1。σ決定に電流密度j=e∫v f dv=σE（式7）計算。先行研究と異なり波数ベクトルk導入せず、BTE解をv∇演算子含む無限級数形式で求めた：f=f₀ + e/(-iω+Γ) ∂f₀/∂ε [1 + v∇/(-iω+Γ)]⁻¹ v·E（式8）。Γ=1/τ。f₀零温度フェルミ-ディラック分布近似、式(7)ε積分後：j=ω_p²ε₀/(-iω+Γ) 3/4π ∬ u (1+ l u∇)⁻¹ (u·E) dΩ（式9）。ω_p²=e²n_e/(mε₀)プラズマ周波数、u=v/|v|、l=v_F/(-iω+Γ)、v_Fフェルミ速度。積分計算（式10-12）任意電場E依存j最終式（式13）導出。

式(1)解2タイプ：発散ゼロ（∇·E=0）と回転ゼロ（∇×E=0、式14）。平面波E∝exp[i(k·r-ωt)]で式14→k×E=0、縦波(L)示す。球体プロセスシミュレーション便利ベクトル球面調和関数L、M、N直交完全系使用。式14 L調和特定。L波波数分散則：ε^L(ω,k^L)=0（式15）。式(6)(13)から∇×E=0解縦誘電率：ε^L=ε_g - ω_p²/ω(ω+iΓ) 3/2 Φ(a²,1,3/2)（式16）。Φレルヒ超越関数、3/2 Φ(a²,1,3/2)=Σ 3/(2n+3) a^{2n}（式17）、a=k v_F/(ω+iΓ)。Kliewer and Fuchs定義と表記差のみ（式18）。恒等式（式19）用い書き換え可能（式20）。Γ=0でKlimontovichとSilin導出方程式（ランダウ減衰研究）と一致。式(16)(18)(20)等価誘電率通例Lindhard誘電関数と呼ばれる（[12]参照）が最初KlimontovichとSilin[8]取得。

横Lindhard誘電率[7]式(13)∇·E=0で導出。実ケースv_F k≪ωでDrude誘電率ε^T=ε_g - ω_p²/(ω²+iΓω)（式21）還元、多く金属実験データ整合。|a|<1で縦誘電率(16)流体力学誘電関数簡略化：ε^L(ω,k^L)=ε_g - ω_p²/(ω²+iΓω - 3/5(v_F k^L)²)（式22）。

平面金属表面平面電磁波反射

境界条件

平面電磁波誘電体-金属界面z=0入射、波数ベクトルxz平面内。誘電体媒質1電場入射E_iと反射E_r波、金属媒質2透過横波E_tと場合により縦波E^L成分。マクスウェル境界条件：電場磁場横成分z=0連続。電子金属外部放出せず、法線方向電流密度ゼロ：ẑ·j|_z=0=0（式23）。

マクスウェル境界条件全項x,y同一依存必要。結果：L波励起p偏光（入射波電場ベクトルE^(i)入射面平行）時のみ可能。横磁界(TM)波でプラズモンポラリトン生成可能、金属球類似。スネル則類似式導出：k_{1x}=k_{2x}=k_{2x}^L=k₁ sinθ₁（式24）。

反射透過係数

p偏光平面電磁波平面金属表面入射電場決定。入射波E_x^(i) x成分用い表現便利：反射波E_x(r)=-r E_x(i)、透過横波E_x(t)=t E_x(i)、透過縦波E_x^(a)=δ E_x(t)=t_L E_x^(i)（式25）。r反射係数、t、t_L透過係数。マクスウェル境界条件と式(23)ABC用い（式26形式）：r=-[(1+δ)ε₁k_{2z}-ε₂k_{1z}]/[(1+δ)ε₁k_{2z}+ε₂k_{1z}]=1-(1+δ)t（式27）、t=2ε₁k_{2z}/[ε₂k_{1z}+(1+δ)ε₁k_{2z}]（式28）、δ=(ε_g-ε)/ε_g · k_{2x}²/(k_{2z}k_{2z}^L)（式29）。δ=0でrはp偏光フレネル反射係数（例[4]式2.49）。同条件下tフレネル透過係数でない（t,r定義フレネルと異なる）。

金属ナノ球光消光

先行論文で著者一人がABC式(23)考慮しLorentz-Mie理論一般化。フレネル係数r類似、反射TMモードl次ミーコエフィシエントb_l：b_l=-[(1+δ_l)ε₁ k₂ψ'_l(k₂R)/ψ_l(k₂R) - ε₂ k₁ψ'_l(k₁R)/ψ_l(k₁R)]/[(1+δ_l)ε₁ k₂ψ'_l(k₂R)/ψ_l(k₂R) - ε₂ k₁ζ'_l(k₁R)/ζ_l(k₁R)]（式30）。δ_l=(ε^T-ε_g)/ε_g · l(l+1)j_l(k₂R)j_l(k₂^L R)/[ψ'_l(k₂R)k₂^L R j'_l(k₂^L R)]（式31）。ψ_l、ζ_lリッカティ-ベッセル/ハンケル関数、j_l球ベッセル関数、プライム引数微分。

古典/一般化Lorentz-Mie理論予測実験データ比較。[15] Hilger, Tenfelde, Kreibig誘電体表面堆積銀ナノ粒子消光スペクトル研究。第1段階：平均径2、3.5、4 nm銀粒子ビーム生成、1ビーム粒子径分布決定、消光スペクトル記録、真空中銀球現象論式Γ=Γ_b+A v_F/RパラメータA=0.25推定。平均径〈D〉=2 nm銀球ビーム実験径分布（D=1-4 nm領域）用い消光スペクトル計算。理論調整パラメータなし。誘電関数定義にLynch and 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。

Fig. 2理論スペクトルKlimontovich-Silin-Lindhardと単純流体力学誘電関数両用い計算。驚くべき両計算近接結果、プラズモン共鳴領域|a|>1にもかかわらず。

ナノメートル銀球消光スペクトル極大値（Fröhlich、プラズモン、表面プラズモンポラリトン(SPP)共鳴）3.5→3.65 eVシフト知られる。非局所モデル実験データ優秀整合、局所(ミエ)理論ω≃3.5 eV極大（Fig. 2、Table 1）。

青方シフト考察：静電近似でb_1のみ消光断面積Q_ext寄与、式(30)近似式(32)用い簡略化。Q_ext極大条件：Re[2(1+δ₁)ε₁+ε₂]=0（式33）。δ₁項でLモード励起考慮、Fröhlich共鳴条件Re(2ε₁+ε₂)=0（式34）と異なる。実験[15]でピーク周波数ω_mと共鳴幅Δω殆ど〈D〉非依存。局所理論A=0.25でD減少に伴いプラズモン共鳴拡大予測（Table 1）。非局所理論ほぼ等幅、異なるピーク位置予測。全粒子寄与重ね合わせΔω実験データ顕著整合。非局所理論A=0でもビームプラズモン共鳴拡大予測興味深い。

ω>4 eVでFig. 2滑らか理論曲線実験狭ピーク群上方。帯間吸収支配、Fig. 1確認可能。観測スペクトル特徴連続バンドから離散準位構造移行量子サイズ効果帰属可能。銀球サイズ〈D〉=3.5 nm増加時吸収先ず相対的最大値比増大し等間隔ディップ系列プラトー形成、〈D〉=4 nmでわずか減少。

プラズモン共鳴青翼形成研究で超微小銀粒子消光スペクトル計算Fig. 3。注目：D=1 nmでプラズモン共鳴完全消失。先行金ナノ球実験研究[19]で同効果観測。[19] Fig. 9で径1.7、1.9、2.0、2.1、2.3、2.5 nm粒子実験スペクトル局所ミエ理論計算比較。整合性悪くプラズモン共鳴拡大・位置記述不能。粒子サイズ変動と誘電関数修正試み失敗。著者ら提案局所ミエ理論補正で異常広/抑制コヒーレント振動バンドフィット不能。Fig. 4非局所ミエ理論適用で劇的変化。調整パラメータ不使用。Johnson and Christy[20]複素屈折率表金誘電関数決定。他パラメータA=1、トルエン屈折率1.37[19]採用。

縦波波数

縦モード横モードと異なり波数桁違い大きい。例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+3/5 a²+3/7 a⁴+3/9 a⁶+... = 1/(1-ε^T/ε_g)（式35）。最簡単ε_g=1、Γ=0でω>ω_p金属横波縦波透明、ω<ω_pでk^L、k^T複素数。

固体透明時、p偏光斜入射平面表面で縦波励起可能。特徴：平坦表面で縦波生成可能（表面プラズモンポラリトン特殊努力必要と対比）。干渉パターン電磁場強度界面垂直方向変調（平行非）、固体剥離で界面平行面ボイド形成可能。ω_p定義よりω>ω_p条件低キャリア密度固体（半導体等）で成立可能。縦誘電率式縮退電子ガス導出ゆえ本ケース不検討。

結論

空間依存誘電関数定義に先行研究全て平面電磁波物質相互作用考慮。ナノオプティクス領域分化場空洞局在化、境界条件何らか考慮必要時、このアプローチ非建設的非厳密。ベクトルヘルムホルツ方程式満たす電場媒質応答計算で解決。導出空間分散誘電関数波数二乗（ヘルムホルツ方程式パラメータ）依存、平面波波数ベクトル非依存。

金属中縦波励起考慮修正フレネル反射係数報告。類似一般化ミーコエフィシエント先行済。本稿理論ナノメートル銀・金クラスタ光消光シミュレーション検証。銀粒子ビーム表面プラズモン共鳴3.5→3.65 eVシフトと幅実験データ優秀整合。さらに非局所モデル直径約2 nm金球プラズモン共鳴消失説明。L波平面入射波平坦表面励起可能、表面プラズモンポラリトンと主要差異。

金属電磁振動特性検討。縦波波数絶対値横波より遥か大発見。例：銀3.5 eVで波数絶対値比130。縦波実部波数対応波長7 nm。波数大き差平坦表面L波励起阻止。しかしLモードナノメートル銀・金粒子励起実証。