[論文レビュー] Numerical Characterization of Fragmentation in Ionic Liquid Clusters

Madeleine Schroeder|arXiv (Cornell University)|Nov 9, 2021

Electrohydrodynamics and Fluid Dynamics参考文献 60被引用数 5

ひとこと要約

本研究では、分子動力学（MD）シミュレーションおよびマルチスケールN体モデルを用いて、電場下におけるイオン液体エレクトロスプライ射出源のイオンクラスター破壊を特徴付け、シミュレートされた実験的遅延電位分析（RPA）データとの強い一致を示した。主な貢献は、電場効果を考慮した物理的ベースの破壊モデルであり、近似的ベイズ計算を用いて実験的RPA曲線からビーム組成および温度を推定することで検証された。

ABSTRACT

Ionic liquid ion sources are a promising technology that can be used for many applications from space propulsion to focused ion beam microetching. Ionic liquid ion sources produce ion beams by extracting single ions and metastable solvated ion clusters from the surface of the ionic liquid and accelerating them using an electric field generated by applying a voltage between a sharp tip and a plate with an aperture. The solvated ion clusters often fragment in the electric field region, reducing the specific impulse and efficiency for propulsion applications and increasing the beam spot size for focused ion beam applications. Fragmentation behavior has previously been characterized in the region with no electric field. However, fragmentation in an electric field has not been investigated as experimental results are difficult to interpret. The goal of this work is to use various numerical methods to characterize fragmentation under the effect of an electric field. Molecular dynamics simulations are performed of various ionic liquid clusters under different conditions to determine the rate of fragmentation. Physics-based models are compared to the molecular dynamics results with the goal of deriving a new model that accounts for the effect of the electric field on fragmentation. Approximate Bayesian computational methods are employed to infer the temperature of different ionic liquid cluster types and the percentage of the beam composed of each species by comparing simulated retarding potential analysis curves to experimental ones. Finally, the results of multi-scale N-body simulations are postprocessed and compared to experimental data. Results show agreement between experimental data and N-body simulations using the fragmentation rates determined by molecular dynamics.

研究の動機と目的

• 電場がイオン液体クラスターの破壊速度に与える影響を理解すること。
• 電場効果を考慮した物理的ベースの破壊モデルを構築し、既存のアレニウスおよびショットキー・モデルを改善すること。
• 近似的ベイズ計算を用いて、実験的RPAデータとシミュレーション結果を一致させ、ビーム種の組成および温度を推定すること。
• MDから得た破壊率を用いたマルチスケールN体シミュレーションを、実験的RPA曲線と照合して検証すること。
• スリットサイズの変化が、加速領域における電場分布および破壊行動に与える影響を評価すること。

提案手法

• 異なる電場強度下でのイオン液体クラスターに対する分子動力学（MD）シミュレーションを実施し、破壊速度および経路を計算する。
• クラスターの解離を破壊と定義し、MD軌道からの平均寿命を計算する。
• 電場強度およびクラスター幾何学的形状を組み込んだ修正アレニウス型モデルを導出することで、破壊速度を予測する。
• MDから得た破壊率を用いてN体シミュレーションを実装し、ビームダイナミクスおよびRPA応答をシミュレートする。
• 近似的ベイズ計算（ABC）を適用し、シミュレートされたRPA曲線と実験データを一致させ、温度およびビーム組成を推定する。
• N体シミュレーション出力の後処理を行い、実験測定値と比較可能な合成RPA曲線を生成する。

実験結果

リサーチクエスチョン

• RQ1 電場強度は、エレクトロスプライ・アイオングラム源におけるイオン液体クラスターの破壊速度にどのように影響するか？
• RQ2 電場下での異なるイオン液体クラスター種において、支配的な破壊経路は何か？
• RQ3 クラスターの幾何学的形状およびエネルギー含量は、電場存在下での破壊結果にどのように影響するか？
• RQ4 スリットサイズの変化が、先端から抽出器までの領域における電場分布および破壊行動に及ぼす影響はどの程度か？
• RQ5 電場効果を含む物理的ベースの破壊モデルをN体シミュレーションと組み合わせることで、実験的RPAデータを正確に予測できるか？

主な発見

• MDから得た破壊率を用いたN体シミュレーションは、実験的RPAデータと優れた一致を示し、シミュレーションフレームワークの妥当性が裏付けられた。
• スリットサイズ（50–150 µm）の違いに対しても、先端先端における電場強度の変動は5%未満であり、発射点での破壊に及ぼす影響は最小限であった。
• 先端から抽出器までの領域、特にスリット付近では、電場分布に顕著な差が観察され、先端から500 µm離れた位置で電場強度に最大50%の変動が生じた。
• 近似的ベイズ計算により、実験的RPA曲線と一致させるシミュレーション曲線から、ビームの温度および質量組成が成功裏に推定された。
• 先端から500 µmを超えた低電場領域における破壊速度は、高電場領域に比べて無視できるほど小さく、破壊は主に加速領域で発生することが示された。
• 本研究では、破壊に及ぼす電場効果が非無視可能であり、宇宙推進およびマイクロファブリケーション応用におけるビーム特性の正確な予測のためには、ビームモデリングに電場効果を組み込む必要があることが明らかになった。