Morrow Zachary, Liu Chang, Kelley C T, Jakubikova Elena
Department of Mathematics , North Carolina State University , Raleigh , North Carolina 27695 , United States.
Department of Chemistry , North Carolina State University , Raleigh , North Carolina 27695 , United States.
J Phys Chem B. 2019 Nov 14;123(45):9677-9684. doi: 10.1021/acs.jpcb.9b08210. Epub 2019 Nov 1.
The potential energy surface (PES) describes the energy of a chemical system as a function of its geometry and is a fundamental concept in computational chemistry. A PES provides much useful information about the system, including the structures and energies of various stationary points, such as local minima, maxima, and transition states. Construction of full-dimensional PESs for molecules with more than 10 atoms is computationally expensive and often not feasible. Previous work in our group used sparse interpolation with polynomial basis functions to construct a surrogate reduced-dimensional PESs along chemically significant reaction coordinates, such as bond lengths, bond angles, and torsion angles. However, polynomial interpolation does not preserve the periodicity of the PES gradient with respect to angular components of geometry, such as torsion angles, which can lead to nonphysical phenomena. In this work, we construct a surrogate PES using trigonometric basis functions, for a system where the selected reaction coordinates all correspond to the torsion angles, resulting in a periodically repeating PES. We find that a trigonometric interpolation basis not only guarantees periodicity of the gradient but also results in slightly lower approximation error than polynomial interpolation.
势能面(PES)将化学体系的能量描述为其几何结构的函数,是计算化学中的一个基本概念。势能面提供了有关该体系的许多有用信息,包括各种驻点的结构和能量,如局部极小值、极大值和过渡态。构建具有超过10个原子的分子的全维势能面在计算上成本高昂且通常不可行。我们小组之前的工作使用多项式基函数的稀疏插值,沿着化学上重要的反应坐标(如键长、键角和扭转角)构建代理降维势能面。然而,多项式插值不能保持势能面梯度相对于几何结构角分量(如扭转角)的周期性,这可能导致非物理现象。在这项工作中,对于所选反应坐标均对应于扭转角的体系,我们使用三角函数基函数构建代理势能面,从而得到一个周期性重复的势能面。我们发现,三角函数插值基不仅保证了梯度的周期性,而且导致的近似误差比多项式插值略低。