Qu Chen, Houston Paul L, Allison Thomas, Bowman Joel M
Independent Researcher, Toronto, Ontario M9B0E3, Canada.
Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.
J Chem Theory Comput. 2025 Apr 8;21(7):3552-3562. doi: 10.1021/acs.jctc.4c01793. Epub 2025 Mar 27.
Given the great importance of linear alkanes in fundamental and applied research, an accurate machine-learned potential (MLP) would be a major advance in computational modeling of these hydrocarbons. Recently, we reported a novel, many-body permutationally invariant model that was trained specifically for the 44-atom hydrocarbon CH on roughly 250,000 B3LYP energies (Qu, C.; Houston, P. L.; Allison, T.; Schneider, B. I.; Bowman, J. M. , , 9339-9353). Here, we demonstrate the accuracy of the transferability of this potential for linear alkanes ranging from butane CH up to CH. Unlike other approaches for transferability that aim for universal applicability, the present approach is targeted for linear alkanes. The mean absolute error (MAE) for energy ranges from 0.26 kcal/mol for butane and rises to 0.73 kcal/mol for CH over the energy range up to 80 kcal/mol for butane and 600 kcal/mol for CH. These values are unprecedented for transferable potentials and indicate the high performance of a targeted transferable potential. The conformational barriers are shown to be in excellent agreement with high-level ab initio calculations for pentane, the largest alkane for which such calculations have been reported. Vibrational power spectra of CH from molecular dynamics calculations are presented and briefly discussed. Finally, the evaluation time for the potential is shown to vary linearly with the number of atoms.
鉴于线性烷烃在基础研究和应用研究中的重要性,精确的机器学习势(MLP)将是这些碳氢化合物计算建模的一项重大进展。最近,我们报道了一种新颖的多体置换不变模型,该模型专门针对44原子的碳氢化合物CH,依据大约250,000个B3LYP能量进行了训练(Qu,C.;Houston,P. L.;Allison,T.;Schneider,B. I.;Bowman,J. M. , ,9339 - 9353)。在此,我们展示了这种势对于从丁烷CH到CH的线性烷烃的可转移性的准确性。与其他旨在实现普遍适用性的可转移性方法不同,本方法针对的是线性烷烃。对于丁烷,能量的平均绝对误差(MAE)在能量范围高达80 kcal/mol时为0.26 kcal/mol,对于CH,在能量范围高达600 kcal/mol时上升至0.73 kcal/mol。这些值对于可转移势来说是前所未有的,表明了一种有针对性的可转移势的高性能。构象能垒与戊烷的高水平从头算计算结果显示出极佳的一致性,戊烷是已报道此类计算的最大烷烃。给出并简要讨论了来自分子动力学计算的CH的振动功率谱。最后,势的评估时间显示与原子数呈线性变化。