[論文レビュー] On the Darwin-Howie-Whelan equations for the scattering of fast electrons described by the Schr\"odinger equation

Thomas Koprucki, Anieza Maltsi+1 more|arXiv (Cornell University)|Dec 21, 2020

Electron and X-Ray Spectroscopy Techniques参考文献 27被引用数 4

ひとこと要約

本稿は、走査型電子顕微鏡（TEM）における高速電子散乱のDarwin-Howie-Whelan（DHW）方程式について、厳密な数学的解析を提供し、有限ビーム近似の誤差見積もりを導出する。系統的行近似（18ビーム）では7桁の精度に達することが示され、2ビームモデルの1桁に比べて顕著に優れている。また、電子ホログラフィーなどの高精度なシミュレーションにおいて、Ewald球面の上下のビームを含める必要があることが示された。

ABSTRACT

The Darwin-Howie-Whelan equations are commonly used to describe and simulate the scattering of fast electrons in transmission electron microscopy. They are a system of infinitely many envelope functions, derived from the Schr\"odinger equation. However, for the simulation of images only a finite set of envelope functions is used, leading to a system of ordinary differential equations in thickness direction of the specimen. We study the mathematical structure of this system and provide error estimates to evaluate the accuracy of special approximations, like the two-beam and the systematic-row approximation.

研究の動機と目的

• TEM シミュレーションで用いられるDarwin-Howie-Whelan（DHW）方程式の数学的基盤を提供すること。
• 無限ビーム系を前提とした場合の、2ビームモデルや系統的行モデルといった有限ビーム近似の精度を分析すること。
• DHW方程式における包絡関数の無限系の切断によって生じる誤差を定量化すること。
• 多ビームシミュレーションにおける一般的なビーム選択ヒューリスティクスの妥当性を検証または反証する理論的誤差境界を確立すること。
• 電子ホログラフィーなどの高精度を要する定量的TEMアプリケーションにおいて、ビーム集合の選定を支援すること。

提案手法

• コラム近似の下で、時間に依存しないシュレーディンガー方程式からDHW方程式を導出し、厚さ方向に変調される包絡関数 ψg(z) に波動関数を分解する。
• 相互作用演算子のエルミート性に起因するハミルトニアン構造を持つ、zに関する1階常微分方程式として系をモデル化する。
• 散乱長と系のパラメータに基づく誤差見積もりを導入し、|α*k0|⁻² および α*z*/ℓ²_scatt を小さな摂動項として用いる。
• G1からG4までの異なるビーム集合Gの間で数値解を比較し、収束性と精度を評価する。基準として最大の集合G4を用いる。
• Ewald球面基準を用いて、励起誤差が最小で振幅が最大となる球面付近のビームを選択する。
• Juliaを用いた数値シミュレーションと照合して理論的誤差境界を検証し、ビーム振幅と励起誤差を可視化・定量化する。

実験結果

リサーチクエスチョン

• RQ1 有限ビーム集合に切断されたDHW方程式の数学的構造は何か？
• RQ2 2ビーム近似は、無限ビーム系と比較して、電子ビーム振幅をどの程度正確に予測できるか？
• RQ3 系統的行近似（gz = 0）は正しいビーム振幅をどの程度正確に捉えているのか？また、どのような状況で不十分になるか？
• RQ4 Ewald球面の上下のビームを含めることで、解の精度にどのような影響があるか？
• RQ5 理論的誤差境界は数値的に検証可能か？また、実用的TEMシミュレーションにおけるビーム選択にどのように寄与するか？

主な発見

• 2ビーム近似（G1）では、(0,0)および(1,0)モードの振幅において、わずか1桁の精度しか得られず、非常に粗い近似であることが判明した。
• 4ビームを用いた系統的行近似（G2）では精度が4桁に向上し、6ビーム（G3）に増やすことでさらに精度向上が得られなかった。
• Ewald球面の上下のビームを含めた（G4、18ビーム）近似では、精度が7桁に向上し、このような拡張が不可欠であることが示された。
• 全ビーム集合Gを用いた基準解は7桁の精度を達成しており、G4も同程度の精度を示し、高精度シミュレーションに適していることが検証された。
• 理論的誤差境界 |α*k0|⁻² ≈ 8.46e-6 および α*z*/ℓ²_scatt ≈ 0.017 は小さく、選定された設定において近似の妥当性が確認された。
• pyTEMにおけるビーム選択戦略（系統的行近似、|Ug| ≥ umin かつ |sg| ≤ es*）は、G3と同等のビーム範囲をカバーしており、有効であることが検証された。