Venturi S, Jaffe R L, Panesi M
University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.
NASA Ames Research Center, Moffett Field, California 94035-1000, United States.
J Phys Chem A. 2020 Jun 25;124(25):5129-5146. doi: 10.1021/acs.jpca.0c02395. Epub 2020 Jun 11.
This work introduces a novel methodology for the quantification of uncertainties associated with potential energy surfaces (PESs) computed from first-principles quantum mechanical calculations. The methodology relies on Bayesian inference and machine learning techniques to construct a stochastic PES and to express the inadequacies associated with the ab initio data points and their fit. By combining high fidelity calculations and reduced-order modeling, the resulting stochastic surface is efficiently forward propagated via quasi-classical trajectory and master equation calculations. In this way, the PES contribution to the uncertainty on predefined quantities of interest (QoIs) is explicitly determined. This study is done at both microscopic (e.g., rovibrational-specific rate coefficients) and macroscopic (e.g., thermal and chemical relaxation properties) levels. A correlation analysis is finally applied to identify the PES regions that require further refinement, based on their effects on the QoI reliability. The methodology is applied to the study of singlet (1A') and quintet (2A') PESs describing the interaction between O molecules and O atoms in their ground electronic state. The investigation of the singlet surface reveals a negligible uncertainty on the kinetic properties and relaxation times, which are found to be in excellent agreement with the ones previously published in the literature. On the other hand, the methodology demonstrated significant uncertainty on the quintet surface, due to inaccuracies in the description of the exchange barrier and the repulsive wall. When forward propagated, this uncertainty is responsible for the variability of 1 order of magnitude in the vibrational relaxation time and of factor four in the exchange reaction rate coefficient, both at 2500 K.
这项工作介绍了一种新颖的方法,用于量化与从第一性原理量子力学计算得到的势能面(PESs)相关的不确定性。该方法依靠贝叶斯推理和机器学习技术来构建一个随机PES,并表达与从头算数据点及其拟合相关的不足之处。通过结合高保真计算和降阶建模,经由准经典轨迹和主方程计算有效地向前传播得到的随机表面。通过这种方式,明确确定了PES对预定义感兴趣量(QoIs)不确定性的贡献。这项研究在微观(例如,振转特定速率系数)和宏观(例如,热和化学弛豫性质)层面上都进行了。最后应用相关性分析,根据PES区域对QoI可靠性的影响来识别需要进一步细化的区域。该方法应用于描述处于基态的O分子与O原子之间相互作用的单重态(1A')和五重态(2A')PESs的研究。对单重态表面的研究表明,动力学性质和弛豫时间的不确定性可忽略不计,这些结果与先前文献中发表的结果非常吻合。另一方面,由于交换势垒和排斥壁描述中的不准确,该方法在五重态表面上表现出显著的不确定性。当向前传播时,这种不确定性导致在2500 K时,振动弛豫时间变化1个数量级,交换反应速率系数变化4倍。