Excitation of confined bulk plasmons in metallic nanoparticles by penetrating electron beams within a nonlocal analytical approach

技術要約：非局所的解析手法による、透過電子ビームを用いた金属ナノ粒子内の閉じ込められたバルク・プラズモンの励起

問題提起
閉じ込められたバルク・プラズモン（CBP）――金属ナノ粒子の体積内における電子密度の縦波的な集団振動――の、透過電子ビームによる励起は、局所応答近似（LRA）モデルでは正確に記述できない現象である。LRAの枠組みは、$\sim$10 nmより大きなナノ粒子における局在表面プラズモン（LSP）のモデル化には成功しているが、有限の電子圧縮性が顕著となる10 nm以下の領域では破綻する。LRAにおいて、透過プローブによるバルク・プラズモンの励起は、誘導ポテンシャルに物理的に不自然な対数発散を引き起こし、これを解決するために通常、アドホックな運動量カットオフが必要となる。さらに、LRAモデルは、プラズマ周波数を超える領域でのCBPの出現や、電子ビームのインパクトパラメータおよび運動エネルギーに対する励起の特異な依存性を説明することができない。非局所効果を組み込んだ、透過軌道における球形ナノ粒子の電子エネルギー損失（EEL）確率に関する透明性の高い解析的な表現が欠如している。

手法
著者らは、伝導電子を半径 $a$ の球形ナノ粒子内に閉じ込められた均質な自由電子ガスとして記述するために、線形流体力学モデルを採用している。このアプローチは、流体力学的分散パラメータ $\beta$ によって特徴付けられる電子ガスの有限の圧縮性を介して、空間分散を組み込んでいる。本研究では、透過型（$b < a$）および外部型（$b > a$）の電子ビーム軌道に対するEEL確率（$\Gamma_{EELS}$）の解析的な表現を導出している（ここで $b$ はインパクトパラメータ、$v$ は電子速度である）。

導出プロセスは以下の通りである：

1. 電子ビームを、$z$ 軸に沿って一定速度で移動する点電荷としてモデル化する。
2. 流体力学的境界条件（表面における法線電流の消失）に従い、金属内部の誘導電荷密度に関する非斉次ヘルムホルツ方程式を解く。
3. 外部および誘導された電荷密度とポテンシャルを球面調和関数を用いて展開し、多極展開解を得る。
4. 全EEL確率を、3つの明確な物理的寄与――バルク損失（$\Gamma_{bulk}$）、ベグレッツング効果（表面補正、$\Gamma_{Begr}$）、および外部損失（$\Gamma_{ext}$）――に分解する。
5. 標準的な自由電子系であるナトリウム（Na）ナノ粒子にこの導出された形式式を適用し、スペクトルのインパクトパラメータ、電子速度、およびナノ粒子サイズへの依存性を分析する。

主要な貢献と結果
本論文は、アドホックな運動量カットオフを用いることなく非局所効果を捉える、統一された解析的枠組みを提供しており、急速な多極収束を示している。主な知見は以下の通りである：

• CBPの出現とスペクトルのブルーシフト： LRAの予測ではプラズマ周波数において単一のバルク・プラズモン・ピークが示されるが、流体力学モデルは $\omega_p$ を超えるCBPに関連する複数のピークの出現を予測する。これらのCBPのスペクトル包絡線は、ナノ粒子サイズが減少するにつれて、また電子のインパクトパラメータが中心から表面に向かって増加するにつれて、顕著なブルーシフトを示す。
• インパクトパラメータ依存性と閾値： 本研究では、「バルク・プラズモン活性化閾値」（$e_{bt}$）を特定している。これは、CBPピーク強度が最大値の $1/e$ に低下するインパクトパラメータとして定義される。バルク損失が粒子内の経路長に比例するというLRAの予測とは対照的に、流体力学モデルは、効率的なCBP励起には電子がナノ粒子の深部まで浸透する必要があることを明らかにしている。閾値 $e_{bt}$ は、ナノ粒子半径が減少するにつれて減少する。
• 選択則： 多極解析を通じて、著者らはCBP励起を支配する選択則を特定している。
  • パリティ則： 中心軌道（$b=0$）では、位相整合条件により偶数-$l$ モードが主に励起される。$b$ が増加するにつれ、奇数-$l$ モードが活性化され、これが観察されたブルーシフトに寄与する。
  • 径方向節ルール： $n$ に関する微細な径方向節選択則が特定されている。中心軌道では、電子ビームとCBPモード間の破壊的干渉により、偶数 $n$ のモードは抑制され、奇数 $n$ のモードが励起される。
  • 速度依存性： 特定の $(l, n)$ モードの励起は、減少波長数 $q_e = \omega a / (\pi v)$ の特定の値で最大となり、これは電子ビームの速度が異なるCBPモードの干渉条件をチューニングすることを示唆している。
• サイズ効果： ナノ粒子半径が減少するにつれ、分散によりCBPモード間のエネルギー間隔が増大する。これによりモードの重なりが減少し、バルク・プラズモン包絡線が個別のピークへと分解される。また、活性化閾値 $e_{bt}$ も減少する。非常に小さなナノ粒子（$a < 1.5$ nm）では、インパクトパラメータに伴う支配的なピークエネルギーの進化は、滑らかなシフトではなく、急激なジャンプを示す。
• 項の物理的解釈： 解析的な式は、表面、バルク、およびベグレッツングの寄与を明示的に分離している。ベグレッツング項は、特にプローブが励起時とエネルギー損失時の両方において粒子内部に存在する成分（$\Gamma_{Begr,i}$）において、CBPの励起に極めて重要であることが示されている。

意義
本論文は、透過電子ビームによるCBP励起のメカニズムを理解するための統一的な解析的基礎を確立したと主張している。流体力学モデル内での閉じた形式の解析的表現を提供することで、本研究は5 nm以下の金属ナノ粒子に対するSTEM-EELS測定の設計および解釈のための実用的なツールを提供する。結果は、バルク・プラズモンの包絡線のブルーシフトやインパクトパラメータへの依存性に関して、実験および ab initio 計算で観察される定性的な挙動を再現している。本研究は、流体力学モデルが、局所的なアプローチが見落としている重要な非局所的傾向（閾値インパクトパラメータの存在や、表面付近での高次CBPの優先的な励起など）を捉えていることを強調しており、それによって、小さな金属系におけるバルク・プラズモンエネルギー領域の電子エネルギー損失に関する、より正確な物理的記述を提供している。