QEq方法在ReaxFF中的应用与局限性
QEq方法简介
QEq(电荷平衡方法)是ReaxFF中最常用的电荷分布方法，其主要特点包括：
-
基本原理：在每个分子动力学步骤中，系统会寻找使静电能量最小化的电荷分布
-
约束条件：整个系统必须保持电中性（总电荷为零）
-
核心思想：电子可以在系统中流动，直到每个原子的"电化学势"（或电负性）相等
-
数学表达：
min_q E_ele(q) = 1/2 ∑_i^N H_ii q_i^2 + ∑_i,j^N H_ij q_i q_j + ∑_i^N χ_i q_i subject to ∑_i^N q_i = 0其中q_i是每个原子i的电荷，χ_i是原子i的Mulliken电负性，H_ii与原子硬度相关
-
与EEM的比较：QEq与电负性均衡方法(EEM)相似，主要区别在于QEq使用类Slater电子密度替代点电荷相互作用
QEq方法的主要局限性
-
解离物种的非整数电荷：
- 当分子分解时，QEq预测出非整数的原子电荷
- 实际上，完全分离的原子或离子应该具有整数电荷
-
分子极化描述不准确：
- 不能正确描述分子在电场中的电荷重分布（极化）现象
-
长距离不合理的电荷转移：
- 根本问题：QEq假设系统是"完美导体"
- 含义：不考虑距离因素带来的电荷转移能量消耗
- 后果：相距很远的分子之间也会发生不合物理的电荷转移
-
在外部电场中的特殊问题：
- 电荷分布不随系统大小收敛
- 当电场电势差超过阈值时，电场力会超过分子间键合力，导致液体水分裂
- 电势阈值取决于系统的键合力强度（水约为6-7V）
- 对于相距较远无分子间键合力的两个水分子，任何电势差都会导致分子积累净电荷
-
QEq在电场中的数学问题：
- 电荷密度与所施加的电势呈线性关系
- 带电分布与所施加电势的线性关系导致非物理的负压力
- 在足够大的电势差下，这种负压力会导致系统分裂
ACKS2解决方案
ACKS2方法通过以下方式解决QEq的主要问题：
-
非全局电荷转移：
- 在大距离上阻止电荷转移
- 预测中性电荷在相距较远的水分子上
- 预测解离的氢和羟基分子上的中性电荷
-
理论基础：
- 从DFT严格推导
- 使用Kohn-Sham势的二阶近似
-
收敛性：
- 随着体系大小增加，电荷分布和水密度收敛到远离界面的体相值
-
局限性：
- 非全局电荷转移与短程库仑相互作用的组合导致电子屏蔽性质不准确
结论
QEq在外部电场中应用时存在固有问题，主要是因为它允许非物理的长程电荷转移。而ACKS2通过其非全局电荷转移特性解决了这些主要问题，但在电场屏蔽方面仍有不足。未来研究方向可能是将ACKS2电荷分布模型与长程静电相互作用结合，从而实现改进的电荷分布和准确的电场屏蔽性质。
注意：使用外加电场进行研究时谨慎使用Qeq方法！
