[论文解读] Semiclassical description of Intermolecular Coulombic Electron Capture in solutions

Nicolas Sisourat|arXiv (Cornell University)|Mar 9, 2026

Photochemistry and Electron Transfer Studies被引用 0

一句话总结

该论文开发了一个利用 OpenMM 的半经典分子动力学方法来模拟水溶液中的 ICEC，并展示 ICEC 的量子产率依赖于电子能量和阳离子浓度，在高浓度下接近单位值。

ABSTRACT

In this work, we present a semiclassical approach to model Intermolecular Coulombic Electron Capture (ICEC) in aqueous solutions using molecular dynamics simulations with OpenMM. We investigate the behavior of an excess electron in the presence of cations (Fe$^{3+}$) in water, focusing on the influence of electron energy and cation concentration on the ICEC quantum yield. Our simulations reveal that the ICEC quantum yield approaches unity at higher concentrations and initial electron energies, while it decreases at lower concentrations due to electron energy loss before reaching the cation.

研究动机与目标

• 使理解在与辐射化学与生物学相关的水环境中 ICEC 的机理成为可能。
• 开发一个半经典分子动力学框架，以在围绕阳离子的水中模拟 ICEC。
• 量化电子能量和阳离子浓度如何影响 ICEC 的量子产率。
• 研究 ICEC 的可观测特征，例如 Fe2+/Fe3+ 平衡和电子能量分布。

提出的方法

• 使用 TIP3P 水模型，采用 AMBER 力场作为溶剂。
• 将电子和阳离子视为经典粒子，使用包含库仑和排斥成分的自定义非键相互作用。
• 在电子上引入阻尼力以模拟水中能量损失，并在 300 K 下应用朗之万动力学。
• 通过 Senk 等人（2024）的模型估算 ICEC 截面，并根据 r_i(t) 与 vΔt 计算 P_ICEC(t)。
• 在每种能量/浓度条件下运行 100 条轨迹，并分析 ICEC 产率 Φ = N_ICEC/N_tot。
• 分析电子-阳离子距离、电子能量和 ICEC 能量分布以解释动力学。

Figure 1: Snapshot of the simulation box: O are in red, H in white, Fe in pink and e in blue.

实验结果

研究问题

• RQ1 在水溶液中 Fe3+ 浓度对 ICEC 量子产率有何影响？
• RQ2 初始电子动能如何影响 ICEC 的可能性和时序？
• RQ3 溶剂诱导的能量损失在不同浓度下对 ICEC 的限制起到何种作用？
• RQ4 在溶剂-阳离子体系中，ICEC 最有可能在距离/能量的哪些范围出现？

主要发现

• 在较高阳离子浓度下，ICEC 量子产率接近单位值。
• 在较低浓度下，溶剂中无弹性碰撞造成的能量损失降低了 ICEC 的概率。
• ICEC 最可能在电子-阳离子距离较近且电子能量较高时发生。
• 由于溶剂相互作用，电子动能迅速衰减，从而限制了 ICEC 的时序。
• 分布显示高能电子有助于早期的 ICEC 事件，而后续事件涉及速度较慢的电子。
• 在研究的浓度范围内，ICEC 产率与阳离子浓度之间可用二次拟合来描述关系。

Figure 2: Distance (in nm) between the electron and the cation (in black) and electron kinetic energy (in red) as function of time for a typical trajectory at C=0.2mol/L. The dashed line indicates the ICEC threshold (i.e. below this line ICEC is not energetically possible).