ICARION：电场与中性气体环境中离子轨迹模拟的模块化框架

《Journal of the American Society for Mass Spectrometry》：ICARION: A Modular Framework for Ion Trajectory Simulation in Electric Fields and Neutral Gas Environments

　　

离子在电场和中性气体中的输运是离子迁移谱仪（ion mobility spectrometers, IMS）、质谱仪（mass spectrometers, MS）以及IMS-MS联用仪器运行的基础。对此类系统的数值建模需要在不同压力和场强区间内一致地处理场驱动运动、随机离子-中性碰撞以及反应。许多IMS-MS仪器运行于过渡区间，需要对碰撞进行显式粒子基处理。研究人员提出了ICARION（Ion Collision And Reaction IntegratiON），一个开源、模块化、基于粒子的模拟框架，用于对离子输运、碰撞和反应进行可复现的模拟。ICARION旨在支持离子动力学的机制研究、探索IMS、MS及相关技术中的新概念，以及为仪器开发和优化快速、可复现地筛选运行参数和几何结构。其模块化架构将多种可选择的场、碰撞和反应模型与结构化配置和数据处理相结合。具有不同电场、几何结构和环境条件的多个仪器域可以无缝地组合在一次耦合IMS-MS仪器的模拟中。这项工作是对ICARION的验证，而非引入新的物理模型。验证包括碰撞热化、双分子反应动力学、使用已有工具验证的均匀电场中的离子输运，以及耦合多域模拟。在这些案例中，ICARION重现了所选物理模型的预期行为，并能够在压力和场时间尺度上相差多个数量级的域之间进行连续模拟。虽然可以根据需要集成外部生成的场图，但计划中的内部场求解器将能够直接根据用户定义的电极构型计算电场。模块化架构还有助于集成额外的物理模型。

研究背景与意义

离子在电场和中性气体中的输运是离子迁移谱（ion mobility spectrometry, IMS）、射频（radio-frequency, RF）离子阱以及IMS-MS联用仪器运行的基础。在极低压力下，离子轨迹近似无碰撞且可确定；在高压力下，频繁碰撞导致可用统计漂移-扩散描述；而许多现代IMS-MS平台运行在中间压力区间，离子平均自由程与电极尺寸相当，必须用随机粒子基方法显式处理离子-中性碰撞。现有工具如SIMION、IDSimF、ITSIM和IMoS 2各有侧重，但耦合IMS-MS系统跨越高碰撞频率漂移区与近无碰撞质量分析区的模拟常依赖定制代码、工具链拼接或简化假设，可复现性和系统验证不足。为此，研究人员开发了ICARION，一个模块化、开源、基于粒子的离子碰撞与反应整合（Ion Collision And Reaction IntegratiON, ICARION）模拟框架。本文不提出新的微观相互作用模型，而是系统验证该框架在碰撞热化、双分子反应动力学、均匀电场中离子输运及IMS-四极杆耦合等场景中的数值一致性和可复现性。该研究成果发表在《Journal of the American Society for Mass Spectrometry》。

关键技术方法

研究人员基于模块化架构实现了ICARION仿真引擎，采用JSON声明式配置和HDF5自描述输出；在电场力作用下传播离子轨迹，并以Poisson过程随机采样碰撞和反应事件。碰撞模型包括有效单硬球散射（hard-sphere scattering, HSS）和几何分辨单接触硬球散射（geometry-resolved hard-sphere scattering, GR-HSS）。使用密度泛函理论（density functional theory, DFT）在ORCA软件中优化H3O+、质子化戊醛和质子化2,6-二叔丁基吡啶（2,6-DTBP）的分子几何结构，并用MobCal-MPI 2.0在氦气、298 K条件下计算碰撞截面（collision cross section, CCS）作为碰撞事件截面。并行计算采用OpenMP共享内存多线程，运行环境为WSL/GCC。

主要研究结果

随机碰撞热化

在无场条件下，研究人员对三种离子在0.2至2000 Pa压力和200至1000 K温度范围内进行了热化模拟。结果表明，离子系综的平动温度在约10至60个平均碰撞时间内收敛到背景气体温度；平衡速度分布与解析麦克斯韦速度分布在多个数量级上一致；热化时间τ95随压力呈τ95∝1/p标度。HSS与GR-HSS模型的热化时间比约为1.14至1.32，主要差异来源于碰撞事件截面σevent的不同，而非碰撞物理本质。

反应动力学

在伪一级条件下，研究人员模拟了H3O+向中性分子M的不可逆双分子质子转移反应。对三个固定反应速率系数k1=1×10^-9 cm³/s、k2=2×10^-9 cm³/s和k3=3×10^-9 cm³/s，模拟得到的反应物离子衰减和产物离子增长均与解析一级速率方程高度一致，证明ICARION对离子轨迹上反应事件的概率采样是正确的。

漂移管中的离子输运

在均匀电场漂移管模拟中，ICARION-HSS与IDSimF和SIMION在200 Pa以上、约3至10 Td条件下给出的约化离子迁移率偏差通常低于1%，验证了三种独立框架在相同碰撞物理和匹配输入参数下的数值一致性。低压下出现的表观偏差反映了非平衡输运效应，即碰撞弛豫时间占漂移时间比例增大，初始速度分布对表观迁移率仍敏感。GR-HSS模型给出的约化迁移率系统性低于HSS，迁移率比约为1.24至1.27，与σevent之比一致，说明差异主要来自碰撞事件率。

多域模拟：IMS-四极杆耦合

研究人员构建了由10 mbar氦气漂移管和10^-6 mbar RF四极杆质量过滤器组成的多域IMS-四极杆模型。在单一模拟中，离子可连续跨越两个压力相差7个数量级的域，无需手动转移离子系综。三种离子在域边界获得约90、171和251 μs的平均到达时间，展示出离子迁移率分离；四极杆不同RF/DC设置可选择性地透传不同离子，且相同随机数种子条件下漂移域轨迹完全可复现。

性能基准

在1000至50000个离子的系综上，ICARION在最快配置下的运行时间均短于IDSimF和SIMION；八线程时相对单线程加速约为2.8至5.7倍，十六线程出现共享内存扩展饱和。ICARION的峰值内存占用高于IDSimF，反映了其较短运行时间伴随较高的内存使用。

当前实现与未来方向

IC