Qu Chen, Conte Riccardo, Houston Paul L, Bowman Joel M
Department of Chemistry & Biochemistry, University of Maryland, College Park, Maryland 20742, USA.
Phys Chem Chem Phys. 2021 Apr 7;23(13):7758-7767. doi: 10.1039/d0cp04221h. Epub 2020 Sep 24.
We present a full-dimensional potential energy surface for acetylacetone (AcAc) using full and fragmented permutationally invariant polynomial approaches. Previously reported MP2/aVTZ energies and gradients are augmented by additional calculations at this level of theory for the fits. Numerous stationary points are reported as are the usual metrics to assess the precision of the fit. The electronic barrier height for the H-atom transfer is roughly 2.2 kcal mol. Diffusion Monte Carlo (DMC) calculations are used to calculate the ground state wavefunction and zero-point energy of acetylacetone. These together with fixed-node DMC calculations for the first excited-state provide the predicted tunneling splitting due to the barrier to H-transfer separating two equivalent wells. Simpler 1d calculations of this splitting are also reported for varying barrier heights including the CCSD(T) barrier height of 3.2 kcal mol. Based on those results the DMC splitting of 160 cm with a statistical uncertainty of roughly 21 cm, calculated using the MP2-based PES, is estimated to decrease to 100 cm for a barrier of 3.2 kcal mol. The fragmented surface is shown to be fast to evaluate.
我们使用全排列不变多项式方法和片段化排列不变多项式方法,给出了乙酰丙酮(AcAc)的全维势能面。此前报道的MP2/aVTZ能量和梯度,通过在该理论水平上进行额外计算以用于拟合。报告了大量的驻点以及用于评估拟合精度的常用指标。氢原子转移的电子势垒高度约为2.2千卡/摩尔。扩散蒙特卡罗(DMC)计算用于计算乙酰丙酮的基态波函数和零点能。这些结果与第一激发态的固定节点DMC计算一起,给出了由于氢转移势垒分隔两个等效阱而产生的预测隧穿分裂。还报告了针对不同势垒高度(包括CCSD(T)势垒高度3.2千卡/摩尔)的这种分裂的更简单一维计算。基于这些结果,使用基于MP2的势能面计算得出的DMC分裂为160厘米,统计不确定性约为21厘米,对于3.2千卡/摩尔的势垒,估计将降至100厘米。结果表明,片段化表面的评估速度很快。