Khabibrakhmanov Almaz, Fedorov Dmitry V, Tkatchenko Alexandre
Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg.
J Chem Theory Comput. 2023 Nov 14;19(21):7895-7907. doi: 10.1021/acs.jctc.3c00797. Epub 2023 Oct 24.
Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α and the dipole-dipole dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.
排斥性短程和吸引性长程范德华(vdW)力在扩展分子系统的行为中起着重要作用。在使用经验力场(应用于此类系统的最流行计算方法)时,vdW力通常由类似 Lennard-Jones 的势来描述,遗憾的是其预测能力有限。在此,我们提出了一种量子力学vdW势的通用参数化方法,该方法仅需要两个自由原子属性——静态偶极极化率α和偶极 - 偶极色散系数。这是通过从量子德鲁德振子(QDO)模型推导势的函数形式、采用平衡距离和结合能的标度律以及应用微观对应态定律来实现的。与21种稀有气体二聚体的从头算结合曲线相比,结果表明vdW - QDO势对于vdW结合能曲线是准确的。vdW - QDO势的函数形式在零距离和无限距离处都具有正确的渐近行为。此外,结果表明阻尼vdW - QDO势可以准确描述由第II族元素组成的二聚体中的vdW相互作用。最后,我们证明了分子中原子vdW - QDO模型在预测分子系统准确色散能方面的适用性。目前的工作朝着构建通用vdW势迈出了重要一步,这可能会使(生物)分子计算研究受益。