多原子反应体系的高精度拟合势能面

本文介绍了近几年来我们组构建多原子反应体系的高精度拟合势能面的进展。我们基于神经网络(NN)方法,成功构建了多原子气相分子体系和气相分子在金属表面解离的一系列势能面。这些势能面的拟合精度相当高,并且经过了严格的量子动力学测试,能广泛应用到动力学研究中。我们还提出了一种新的置换不变势能面的拟合方法,即基本不变量神经网络方法(FI-NN)。基本不变量的使用极大地减少了神经网络输入层多项式的个数,有效提高了势能面的计算速度。     

CASSCF方法对H2O＋激发态势能面的研究

采用Gaussian 98程序中的CASSCF（完全活性空间自洽场）计算了H2O＋激发态的势能面，并对计算结果进行拟合，得到了势能曲面的解析表达式.作者对活性空间和基组的选取做了分析，并将得到的势能面与广泛使用的MRD-CI结果进行了比较，结果基本一致.     

Exploring Potential Energy Surfaces for Aggregation‐Induced Emission—From Solution to Crystal

Aggregation‐induced emission (AIE) is a phenomenon where non‐luminescent compounds in solution become strongly luminescent in aggregate and solid phase. It provides a fertile ground for luminescent applications that has rapidly developed in the last 15 years. In this review, we focus on the contributions of theory and computations to understanding the molecular mechanism behind it. Starting from initial models, such as restriction of intramolecular rotations (RIR), and the calculation of non‐radiative rates with Fermi's golden rule (FGR), we center on studies of the global excited‐state potential energy surfaces that have provided the basis for the restricted access to a conical intersection (RACI) model. In this model, which has been shown to apply for a diverse group of AIEgens, the lack of fluorescence in solution comes from radiationless decay at a CI in solution that is hindered in the aggregate state. We also highlight how intermolecular interactions modulate the photophysics in the aggregate phase, in terms of fluorescence quantum yield and emission color.     

基于新一代密度泛函和神经网络势能面的量子反应动力学计算

介绍了近年来发展起来的新一代密度泛函XYG3及利用神经网络构造分子体系势能面的最新进展。以H₃和CH₅等体系为实例,表明基于高效准确的密度泛函电子结构计算,与神经网络高精度势能面构造的理想结合,可以得到可靠的化学动力学结果,并有望用于更大更复杂的体系。     

Theoretical Description of Water from Single-Molecule to Condensed Phase: Recent Progress on Potential Energy Surfaces and Molecular Dynamics

In this work, we review recent progress on the view of potential energy surfaces and molecular dynamics study of water and its related reactions in the last decade or so. Some important gas-phase reactions of water with radicals, chemisorbed dissociative dynamics of water on solid surfaces, and statistical mechanics and vibrational spectrum simulations of water from clusters to the condensed phase have been introduced. The recently developed machine learning techniques, such as the neural networks in a combination of permutational invariant polynomials or fundamental invariants, the atomic neural networks framework, the gaussian approximation potentials with the smooth overlap of atomic position kernel, as well as the many-body expansion framework for the construction of highly accurate potential energy surfaces, have also been discussed. Finally, some suggestions have been provided for further improvement of the potential energy surfaces and dynamics methods of water-related systems.     

Ne-HCl势能面和振转光谱的理论研究

利用量子化学计算方法CCSD(T)和大基组aug-cc-pVTZ加键函数3_s3_p2_d对Ne-HCl体系的分子间势能面进行了理论研究.结果表明,势能面上有两个势阱,分别对应于线性Ne-ClH和Ne-HCl构型.通过精确求解核运动方程发现,该从头算势能面分别支持5个(对Ne-HCl)和7个(Ne-DCl)振动束缚态.计算得到的振转跃迁频率与实值值吻合.     

CO2势能面和振动激发态的理论研究

采用振动自洽场-组态相互作用(SCF-CI)方法通过实验振动光谱优化了CO₂分子的势能函数,由该势能函数计算得到的纯振动光谱数据与实验值相比,所有能级的误差均在4cm^-1以内,均方根偏差为1.50cm^-1,所预测的Π态振转光谱也与实验值很接近.     

SO2基态势能面和振转能级的理论研究

在键长-键角坐标下,精确求解SO₂的核振动方程,并通过与实验数据比较来优化势能参数,由所得势能面计算得到38个振动能级,与实验值相比,均方根误差为0.93cm^-1.计算了³⁴SO₂的部分振动能级以及³²SO₂的J=6以下的部分振动能级,所得结果均与实验值较为吻合     

O3势能面和振动激发态的理论研究

由于O3在大气化学中的重要作用，近年来针对它的实验和理论研究皆较活跃．尽管对O3的高振动激发态已进行过广泛的研究［“］，得到过多种由实验光谱定出的势能面［‘，’］或由从头算得到的势能函数【‘］，但由于近年来又增加了一些新的高精度的振转光谱实验数据k，’」，而以     

次氟酸分子的势能函数及其振动光谱的同位素效应的理论研究

用三原子振动激发态的变分计算程序(TRIATOM)精确计算次氟酸分子H¹⁶OF的振动激发态的能级以及次氯酸分子中的H和O分别被D和¹⁸O取代后的H¹⁸OF,D¹⁶)OF和D¹⁸OF的同位素效应,理论计算值与已有的实验结果吻合较好。预测了一些尚未观测到的谱线频率及同位素效应,并确立了一个同位素位移的加和规则。     

Three-Dimensional ab initio Potential Energy Surface and Predicted Spectra for the CH4-Ne Complex

We present a new three-dimensional potential energy surface (PES) for CH\begin{document}$ _4 $\end{document}-Ne complex. The electronic structure computations were carried out using the coupled-cluster method with singles, doubles, and perturbative triples [CCSD(T)], the augmented correlation-consistent aug-cc-pVXZ (X = T, Q) basis sets were employed with bond functions placed at the mid-point on the intermolecular axis, and the energies obtained were then extrapolated to the complete basis set limit. Analytic intermolecular PES is obtained by least-squares fitting to the Morse/Long-Range (MLR) potential function form. These fits to 664 points have root-mean-square deviations of 0.042 cm\begin{document}$ ^{-1} $\end{document}. The bound rovibrational levels are calculated for the first time, and the predicted infrared spectra are in good agreement with the experimental values. The microwave spectra for CH\begin{document}$ _4 $\end{document}-Ne dimer have also been predicted for the first time. The analytic PES can be used for modeling the dynamical behavior in CH\begin{document}$ _4 $\end{document}-(Ne)\begin{document}$ _N $\end{document} clusters, and it will be useful for future studies of the collision-induced-absorption for the CH\begin{document}$ _4 $\end{document}-Ne dimer.     

Investigation of the Benzene–Naphthalene and Naphthalene–Naphthalene Potential Energy Surfaces: DFT/CCSD(T) Correction Scheme

The potential energy surfaces of the naphthalene dimer and benzene–naphthalene complexes are investigated using the recently developed DFT/CCSD(T) correction scheme [J. Chem. Phys. 2008 , 128, 114 102]. One and three minima are located on the PES of the benzene–naphthalene and the naphthalene dimer complexes, respectively, all of which are of the parallel‐displaced type. The stabilities of benzene–naphthalene and the naphthalene dimer are ?4.2 and ?6.2 kcal mol^?1, respectively. Unlike the benzene dimer, where the T‐shaped complex is the global minimum, the lowest‐energy T‐shaped structure is about 0.2 and 1.6 kcal mol^?1 above the global minimum on the benzene–naphthalene and the naphthalene dimer potential energy surfaces, respectively.     

范德华复合物Kr-C2H2的势能面和振转光谱的理论研究

采用[CCSD（T）]-F12方法和aug-cc-pVTZ基组,同时引入中心键函数（3s3p2d1f1g）构建了Kr-C₂H₂体系的高精度四维势能面.在构建势能面时考虑了分子间的振动方式及C₂H₂单体内的ν₁对称伸缩和ν₃反对称伸缩振动.将计算得到的四维势能面在Q₁方向和Q₃方向分别做积分得到C₂H₂单体分别处于振动基态和（ν₁,ν₃）=（1,1）激发态的平均势能面.计算结果表明,这2个平均势能面均存在2个等价的T型全局极小值和2个等价线性极小值.全局极小值的几何构型位于R=0.41 nm,θ=65.6°/114.4°,势阱深度为151.88 cm^-1.对径向部分采用离散变量表象法（DVR）,角度部分采用有限基组表象法（FBR）,并结合Lanczos循环算法计算了Kr-C₂H₂的振转能级和束缚态.计算结果表明,复合物在（ν₁,ν₃）=（1,1）区域的带心位移为-1.48 cm^-1,表现为红移,与实验值-1.38 cm^-1很接近;计算得到的红外跃迁频率也与实验值相吻合,说明得到的从头算势能面具有高精度.     

New ab initio Potential Energy Surfaces for Cl（^2P3/2,^2P1/2）＋H2 Reaction

New global three dimensional potential energy surfaces for the Cl＋H2 reactive system have been constructed using accurate multireference configuration interaction calculations with a large basis set. The three lowest adiabatic potential energy surfaces correlating asymptotically with Cl（^2p）＋H2 have been transformed to adiabatic representation, which leads to a fourth coupling potential for non-linear geometries. In addition, the spin-orbit coupling surfaces have also been computed using the Breit-Pauli Hamiltonian. Properties of the new potential are described. Reaction dynamics based on the new potential agrees with the recent experimental results quite well.     

Adiabatic Potential Energy Surfaces and Photodissociation Mechanisms for Highly Excited States of H2O

Full-dimensional adiabatic potential energy surfaces of the electronic ground state \begin{document}$ \tilde X $\end{document} and nine excited states \begin{document}$ \tilde A $\end{document}, \begin{document}$ \tilde I $\end{document}, \begin{document}$ \tilde B $\end{document}, \begin{document}$ \tilde C $\end{document}, \begin{document}$ \tilde D $\end{document}, \begin{document}$ \tilde D' $\end{document}, \begin{document}$ \tilde D'' $\end{document}, \begin{document}$ \tilde E' $\end{document} and \begin{document}$ \tilde F $\end{document} of H\begin{document}$ _2 $\end{document}O molecule are developed at the level of internally contracted multireference configuration interaction with the Davidson correction. The potential energy surfaces are fitted by using Gaussian process regression combining permutation invariant polynomials. With a large selected active space and extra diffuse basis set to describe these Rydberg states, the calculated vertical excited energies and equilibrium geometries are in good agreement with the previous theoretical and experimental values. Compared with the well-investigated photodissociation of the first three low-lying states, both theoretical and experimental studies on higher states are still limited. In this work, we focus on all the three channels of the highly excited state, which are directly involved in the vacuum ultraviolet photodissociation of water. In particular, some conical intersections of \begin{document}$ \tilde D $\end{document}-\begin{document}$ \tilde E' $\end{document}, \begin{document}$ \tilde E' $\end{document}-\begin{document}$ \tilde F $\end{document}, \begin{document}$ \tilde A $\end{document}-\begin{document}$ \tilde I $\end{document} and \begin{document}$ \tilde I $\end{document}-\begin{document}$ \tilde C $\end{document} states are clearly illustrated for the first time based on the newly developed potential energy surfaces (PESs). The nonadiabatic dissociation pathways for these excited states are discussed in detail, which may shed light on the photodissociation mechanisms for these highly excited states.     

LaH分子结构和基态势能曲线的量子化学计算

用常规的单参考态HF、B3LYP、MPn、QCISD(T)方法在能量一致相对论有效势近似下计算了LaH分子平衡结构和基态势能曲线,考察了这些方法在计算远离平衡的金属氢合键体系势能时存在的不趋于离解极限的缺陷,提出了从B3LYP的平衡位置附近势能曲线拟合得到适用于整个空间范围的Murrell-Sorbie解析势能函数的计算方法,由此计算的振转常数和已有的实验光谱数据完全吻合.     

Full-Dimensional Potential Energy Surfaces of Ground (\begin{document}${\tilde X^2}$\end{document}A') and Excited (\begin{document}${\tilde{A}^2}$\end{document}A") Electronic States of HCO and Absorption Spectrum

In this work, high-fidelity full-dimensional potential energy surfaces (PESs) of the ground (\begin{document}$\tilde X^2$\end{document}A\begin{document}$'$\end{document}) and first doublet excited (\begin{document}$\tilde A^2$\end{document}A\begin{document}$"$\end{document}) electronic states of HCO were constructed using neural network method. In total, 4624 high-level ab initio points have been used which were calculated at Davidson corrected internally contracted MRCI-F12 level of theory with a quite large basis set (ACV5Z) without any scaling scheme. Compared with the results obtained from the scaled PESs of Ndengué et al., the absorption spectrum based on our PESs has slightly larger intensity, and the peak positions are shifted to smaller energy for dozens of wavenumbers. It is indicated that the scaling of potential energy may make some unpredictable difference on the dynamical results. However, the resonance energies based on those scaled PESs are slightly closer to the current available experimental values than ours. Nevertheless, the unscaled high-level PESs developed in this work might provide a platform for further experimental and theoretical photodissociation and collisional dynamic studies for HCO system.