Time- and frequency-domain study for electron beams penetrating dielectric nanospheres: fingerprints of Cherenkov and transition radiation

技術要約：誘電体ナノスフィアを貫通する電子ビームの時系列および周波数領域研究

問題提起
高速電子ビームは、エバネッセント近接場を介してフォトニックモードに結合する、高度に局在化した広帯域励起源としてナノフォトニクスにおいて機能する。金属ナノ粒子におけるプラズモン励起を調べるための確立された手法として、電子エネルギー損失分光法（EELS）やカソードルミネセンス（CL）があるが、誘電体ナノ構造（具体的にはシリコン（Si）ナノスフィア）における励起メカニズムは依然として複雑である。誘電体において、高速電子は分極電流を駆動し、ミー型共鳴、遷移放射（TR）、そして特定の条件下ではチェレンコフ放射（CR）を励起する。中心となる課題は、電子速度がチェレンコフ閾値を超える場合において、CRとTRの相互作用をどのように切り分けるかである。球体のような曲率を持つ幾何学的形状では、表面の曲率によってマルチポール係数内に放射チャネルの混合が生じ、平面界面と比較して遠視野CLスペクトルの解釈を困難にする。本研究では、高速電子ビームが誘電体ナノ粒子内のフォトニックモードをどのように励起し、これらの励起が光学応答にどのように現れるかを、シリコンを主要な例として用いて検討する。

手法
著者らは、Siナノスフィア（半径 $R = 79$ nm、インパクトパラメータ $b=0$）を貫通する高速電子を解析するために、時系列および周波数領域の両方を利用した統合的な理論フレームワークを採用している。

1. 数値シミュレーション（時間領域）： 本研究では、マクスウェル方程式を解くために、不連続ガラーキン時間領域（DGTD）有限要素法を用いている。この手法により、ナノメートルおよびフェムト秒の範囲での時空間分解が可能となる。電子源は、$z$軸に沿って一定速度 $v_e = \beta c$ で移動するガウス電荷分布（$\sigma_e = 5$ nm）としてモデル化されている。「反跳なし（no-recoil）」近似を採用し、直線的な軌跡を仮定している。
2. 解析理論（周波数領域）： 相補的な解析計算が、ミー理論に基づき行われている。これには、CL発光確率および角度分解スペクトル放射強度のマルチポール展開が含まれる。
3. 材料モデル： 物理的メカニズムを分離するために、2つの異なるシリコン材料モデルを調査している：
   • 定数誘電率モデル： 材料共鳴による複雑さを排除し、CRとTRを透明に観察するために、非分散・無損失モデル（$\varepsilon = 16$）を使用している。
   • Drude–Lorentzモデル： バンド間遷移、バルクプラズモン、および材料減衰を捉えるために、実験データ（最大30 eVまで）に適合させた、分散性と損失を持つ現実的なモデルである。

主な観測量には、時間依存の角度分解ポインティング・フラックス $\mathcal{P}(\theta, t)$、および周波数領域の角度分解CL発光確率 $\Gamma_{CL}(\theta, \omega)$ が含まれる。

主な結果

• 定数誘電率領域（無損失Si）：

  • CRとTRの共存： 電子速度がチェレンコフ閾値（$\beta = 0.7, n=4$）を超えている場合、CRとTRの両方が励起される。時間領域のスナップショットは、バルク内で形成され、境界へと伝播し、そこで部分的に透過される明確なチェレンコフ・コーンを明らかにしている。
  • 遠視野における相互作用： この無損失領域では、CRとTRが遠視野放射において同等の役割を果たす。強力なCR信号が、角度分解CLマップにおける特徴的な干渉パターンを効果的に覆い隠してしまう。
  • 干渉メカニズム： CRが存在しないより低い速度（$\beta = 0.1$）では、TRが支配的となる。遠視野CLは、球体の入り口と出口の地点で放出される電場の重ね合わせから生じる、「二重スリット」実験に類似した明確な干渉パターンを示す。干渉条件 $\omega(\theta) \propto \frac{m\pi}{R} \frac{\beta}{1+\beta \cos \theta}$ が観察される。
  • サイズ依存性： 遷移点における近接場振幅は異なる球のサイズ間で同様であるが、より大きな球は、有効な分極体積が大きいため、遠視野においてより強い双極子放射をサポートする。

• 現実的なDrude–Lorentz領域（分散性Si）：

  • CRの抑制： 現実的な材料モデルを使用すると、強い吸収と分散（Lorentz極によってモデル化）が、特に高エネルギー領域において、コヒーレントなチェレンコフ・フロントを著しく抑制する。その結果、$\beta = 0.7$ の場合でも、時間領域のスナップショットにおいて明確なCRコーンは観察されない。
  • TRの支配： 現実的なモデルでは、理論的なチェレンコフ閾値を十分に超えた速度であっても、TRが遠視野CLの主要なメカニズムとなる。
  • 近接場と遠視野の対応関係： 本研究は、近接場のダイナミクスと遠視野放射の間の直接的な関連性を確立している。入り口と出口の遷移点における近接場の干渉が、遠視野の双極子放射パターンを支配している。これらの近接場の建設的干渉が、遠視野CLスペクトルに見られる共鳴位置を再現する。
  • マルチポール挙動： 高次のマルチポール（$\ell > 1$）は現実的なモデルにおいて強く減衰し、定数誘電率の場合に見られる鋭いピークと比較して、共鳴の広幅化をもたらす。エネルギー・角度マップにおける曲がった干渉パターンはTRに起因しており、CR閾値はこれらのパターンが不鮮明になる境界として機能する。

意義および主張
本論文は、組み合わせた時間領域および周波数領域のフレームワークが、高速電子が誘電体ナノ粒子を貫通する際の放射励起チャネルの直接的な視界を提供し、従来の周波数領域のみの解析では不明瞭であった相互作用メカニズムを明らかにするものであると主張している。

具体的には、著者らは以下のことを実証している：

1. 材料のリアリズムの重要性： 理想化されたモデルは、閾値以上ではCRとTRが同等の強度を持つことを示唆しているが、現実的な材料特性（分散と損失）はCRを抑制し、TRを誘電体ナノスフィアにおける遠視野CLの主要な寄与因子にする。
2. メカニズムの分離： 時間領域のポインティング・フラックスによってTRとCRの寄与を分離することで、それらの干渉がどのように遠視野の角度分布を形成するかを明確にしている。
3. 診断的有用性： 観察されたTR干渉パターンは、遠視野の特性を推論するために利用できる。具体的には、既知の電子速度に対して、これらの干渉パターンを用いてSi球のサイズを理論的に決定したり、あるいは電子ビームの位置を走査することで、任意のナノ粒子の正確な形状を回収したりすることが可能である。

本研究は、遷移点における近接場のダイナミクスが、特に双極子成分に関して、遠視野で見られる主要な特徴の主因であることを結論付けており、これは誘電体ナノ構造における電子と光の相互作用に関する洗練された理解を提供するものである。