[论文解读] Numerical Characterization of Fragmentation in Ionic Liquid Clusters

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

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

一句话总结

本研究采用分子动力学模拟与多尺度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），解离速率可忽略不计，表明解离主要发生在加速区。
• 本研究证明，电场对解离的影响不可忽视，必须在推进与微制造应用的束流建模中予以考虑，以准确预测束流特性。