Ab initio ground‐state potential energy function and vibration‐rotation energy levels of imidogen,NH

The accurate ground‐state potential energy function of imidogen, NH, has been determined from ab initio calculations using the multireference averaged coupled‐pair functional (MR‐ACPF) method in conjunction with the correlation‐consistent core‐valence basis sets up to octuple‐zeta quality. The importance of several effects, including electron correlation beyond the MR‐ACPF level of approximation, the scalar relativistic, adiabatic, and nonadiabatic corrections were discussed. Along with the large one‐particle basis set, all of these effects were found to be crucial to attain “spectroscopic” accuracy of the theoretical predictions of vibration‐rotation energy levels of NH. © 2015 Wiley Periodicals, Inc.     

Theoretical Study on the Potential Energy Curve and Vibration-rotation Spectra of BeF Radical

The potential energy curve (PEC) of BeF(X2Σ+) radical is investigated by using the complete active space self-consistent field (CASSCF) method followed by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach over the internuclear separation range from 0.0522 to 2.0472 nm. The PEC is fitted to the analytic Murrell-Sorbie function, which is employed to accurately determine the spectroscopic parameters. The present D0, De, Re, ωe, ωeχe, αe and Be are 6.14 eV, 6.22 eV, 0.1372 nm, 1236.12 cm-1, 9.11 cm-1, 0.0175 cm-1 and 1.4651 cm-1, respectively. These parameters have been compared with those of previous investigations reported in the literature. With PEC determined at the present level of theory, a total of 75 vibrational states have been predicted for the first time by numerically solving the radial Schrdinger equation of nuclear motion using the Numerov method. For each vibrational state, the complete vibrational levels, classical turning points, inertial rotation and centrifugal distortion constants are determined for the first time. Comparing with the available experiments and other theories, we find that the present spectroscopic parameter and molecular constant results are more accurate and complete than the previous theoretical investigations.     

An analytic potential energy function for the amide–amide and amide–water intermolecular hydrogen bonds in peptides

An analytic potential energy function is proposed and applied to evaluate the amide–amide and amide–water hydrogen‐bonding interaction energies in peptides. The parameters in the analytic function are derived from fitting to the potential energy curves of 10 hydrogen‐bonded training dimers. The analytic potential energy function is then employed to calculate the N? H…O?C, C? H…O?C, N? H…OH₂, and C?O…HOH hydrogen‐bonding interaction energies in amide–amide and amide–water dimers containing N‐methylacetamide, acetamide, glycine dipeptide, alanine dipeptide, N‐methylformamide, N‐methylpropanamide, N‐ethylacetamide and/or water molecules. The potential energy curves of these systems are therefore obtained, including the equilibrium hydrogen bond distances R(O…H) and the hydrogen‐bonding energies. The function is also applied to calculate the binding energies in models of β‐sheets. The calculation results show that the potential energy curves obtained from the analytic function are in good agreement with those obtained from MP2/6‐31+G** calculations by including the BSSE correction, which demonstrate that the analytic function proposed in this work can be used to predict the hydrogen‐bonding interaction energies in peptides quickly and accurately. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

CN~+分子势能函数和光谱常数的多参考组态相互作用方法研究

采用多参考组态相互作用方法和aug-cc-p V5Z基函数组计算了CN+分子11∑+,21∑+,13∑+和13Π电子态的势能曲线。利用MS势能函数拟合得到了相应的解析势能表达式。在此基础上求解CN+分子的核运动薛定谔方程,获得了全部振动和转动能级,并用Dunham系数展开式拟合出了光谱常数,与目前仅有的11∑+,21∑+态的文献报道结果进行了比较。结果可对航天尾气及工业过程光谱方法监控提供参考。     

OUH体系的结构和分析势能函数

采用密度泛涵B3LYP方法优化出了OUH分子的各种结构，确定了最稳定构型和离解能，以及它们的谐性力常数，并导出双原子分子UH,UO的Murrell-Sorbie势能函数及其光谱数据。采用多体项展式方法，导出OUH(X^4A')基态分子的分析势能函数，获得OUH(X^4A')体系的势能面，考察了这个势能函数的基本性质，正确地复现出OUH分子的平衡结构特征，结果表明：U+OH,O+UH,H+UO的反应均为无阈能的放热能反应。为进一步探讨OUH体系的反应动力学过程打下了基础。     

Theoretical investigations on analytical potential energy function and spectroscopic parameters for the state b3Πu of dimer 7Li2

The SAC‐CI (symmetry‐adapted‐cluster configuration‐interaction) method presented in Gaussian 03 program package is applied to investigate the adiabatic potential energy curves (PECs) of ⁷Li₂(b³Π_u). These calculations are performed at numbers of basis sets, such as 6‐311++G(3df,3pd), 6‐311++G(2df,2pd), 6‐311++G(df,pd), D95V++, D95(3df,3pd), D95(d,p), cc‐PVTZ, 6‐311++G and 6‐311++G(d,p). All the ab initio calculated points are fitted to the analytic Murrell‐Sorbie functions and then used to compute the spectroscopic parameters. The analytic potential energy function (APEF) for this b³Π_u state is reported. By comparison, the spectroscopic parameters reproduced by the APEF attained at 6‐311++G(2df,2pd) are found to be very close to the latest experimental findings. With the APEF obtained at the SAC‐CI/6‐311++G(2df,2pd) level of theory, a total of 62 vibrational states is found when J = 0. The complete vibrational levels, classical turning points, inertial rotation and centrifugal distortion constants for these vibrational states are also reported. The reasonable dissociation limit for this state is deduced using the calculated results at present. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The flying ice cube: Velocity rescaling in molecular dynamics leads to violation of energy equipartition

This article describes an unexpected phenomenon encountered during MD simulations: velocity rescaling using standard protocols can systematically change the proportion of total kinetic energy (KE) found in motions associated with the various degrees of freedom. Under these conditions, the simulation violates the principle of equipartition of energy, which requires a mean kinetic energy of RT/2 in each degree of freedom. A particularly pathological form of this problem occurs if one does not periodically remove the net translation of (and rotation about) the center of mass. In this case, almost all of the kinetic energy is converted into these two kinds of motion, producing a system with almost no kinetic energy associated with the internal degrees of freedom. We call this phenomenon “the flying ice cube.” We present a mathematical analysis of a simple diatomic system with two degrees of freedom, to document the origin of the problem. We then present examples from three kinds of MD simulations, one being an in vacuo simulation on a diatomic system, one involving a low resolution model of DNA in vacuo, and the third using a traditional all-atom DNA model with full solvation, periodic boundary conditions, and the particle mesh Ewald method for treating long-range electrostatics. Finally, we discuss methods for avoiding the problem. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 726–740, 1998     

Potential energy calculations on the 2∶1 inclusion compound between deoxycholic acid and (+)-camphor

The crystal structure of the 21 inclusion compound between deoxycholic acid and (+)-camphor has been solved by other authors. We have studied its crystal packing by potential energy calculations. Four arrangements of camphor, satisfactory from an energy point of view, were found and used to refine the crystal structure by a least-squares procedure. The results seem to indicate that more than one orientation is possible for the guest molecule and that the X-ray data do not allow us to establish which guest-molecule arrangements are actually present.     

Potential energy surface of the photolysis of isocyanic acid HNCO

The dissociation curves of the photolysis of the isocyanic acid HNCO→HN+CO corresponding to the ground state (S0), the first triplet excited state (TO and the first singlet excited state (S1) have been studied respectively at the UHF/6-311G** and CIS/6-311G** levels using ab initio method. The energy surface crossing points, S1/T1 T1/S0 and S1/S0, have been found and the characteristics of the energy minimum crossing point were given, based on which, the changes of the crossing points' geometries along the lower electronic energy surface and its end-result have been located according to the steepest descent principle. The computational result indicates thatthe photolysis of the isocyanic acid HNCO→HN + CO has three competitive reaction channels ((A)-(C)), and from the kinetic piont of view, channel (A) is the most advantageous.     

Spectroscopic parameter and molecular constant investigations for low‐lying electronic states of P ion

The potential energy curves (PECs) of three low‐lying electronic states of P ion, X²Π_u, A²Σ, and B²Σ, have been studied using the full valence complete active space self‐consistent field method followed by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach and MRCI with Davidson correction (+Q). The correlation‐consistent basis sets, aug‐cc‐pV5Z and aug‐cc‐pV6Z, are used and the total energies are extrapolated to the complete basis set limit. Using these PECs obtained with the MRCI+Q/56‐extrapolation, the spectroscopic parameters for these electronic states are determined and compared in detail with experimental data and those of previous studies reported in the literature. The comparison shows that excellent agreement exists between the present results and the available experiments. The first 40 vibrational states for the three electronic states are also computed when the rotational quantum number J equals zero. For each vibrational state, the vibrational level G(υ), inertial rotation constant B_υ, and centrifugal distortion constant D_υ are determined when J = 0, which are in good accord with the available measurements. © 2012 Wiley Periodicals, Inc.     

Investigations on adiabatic potential energy curve of Li2()

The SAC-CI method is used to investigate the spectroscopic properties of ⁷Li₂(). The adiabatic potential energy curves are calculated and fitted to the analytic Murrell–Sorbie function. The spectroscopic parameters reproduced by the potential attained at cc-PVTZ are found to be very close to the experiments. With the potential obtained at the SAC-CI/cc-PVTZ level of theory, a total of 62 vibrational states is found when J = 0. For each vibrational state, the vibrational level, classical turning points, inertial rotation and centrifugal distortion constants are calculated. Good agreement is obtained when they are compared with the available RKR data.     

OCS电子基态势能面与振动光谱的理论研究

本文采用键长-键角内标系下的自洽场-组态相互作用方法精确计算了OCS分子的振动高激发态能级,并结合实验观测到的振动能级利用非线性最小二乘法优化电子基态势能函数中的势能参数。由优化所得的势能面计算出的振动激发态能级与50个实验观测到的振动能级比较,标准偏差为0.08cm^-^1。此外,还用该势能面计算了OCS同位素分子的振动能级,计算结果与实验值也十分吻合。     

Potential energy surfaces and vibrational spectra for isotopomers of N_2O

ＴｈｅＮ2Ｏｍｏｌｅｃｕｌｅｐｌａｙｓａｖｅｒｙｉｍｐｏｒｔａｎｔｒｏｌｅｉｎｏｚｏｎｅｃｈｅｍｉｓｔｒｙｉｎｔｈｅａｔｏｍｏｓｐｈｅｒｅ.Ｉｔｓｒｏｖｉｂｒａｔｉｏｎａｌｓｐｅｃｔｒａｃａｎｐｒｏｖｉｄｅｕｓｅｆｕｌｉｎｆｏｒｍａｔｉｏｎａｂｏｕｔｔｈｅｉｎｔｒａｍｏｌｅｃｕｌａｒｄｙｎａｍｉｃｐｒｏｐｅｒｔｉｅｓ.Ｎ2Ｏｉｓａｌｉｎｅａｒｍｏｌｅｃｕｌｅａｔｔｈｅｅｑｕｉｌｉｂｒｉｕｍｇｅｏｍｅｔｒｙｆｏｒｔｈｅｅｌｅｃｔｒｏｎｉｃｇｒｏｕｎｄｓｔａｔｅ.Ｎｕｍｅｒｏｕｓｉｎｖｅｓｔｉｇａｔｉｏ…     

Theoretical study on the [Si, C, N, O] potential energy surface

The structures, energetics, spectroscopies, and stabilities of the doublet [Si, C, N, O] radical are explored at the density functional theory and ab initio levels. Sixteen isomers are located, connected by 29 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest lying isomer is a linear SiNCO 1 (0.0 kcal/mol) mainly featuring a cumulene | . Si = N = C = O. The second and third low-lying isomers are bent OSiCN 2 (8.8) and bent OSiNC 3 (11.1), respectively. All the three low-lying isomers 1, 2, 3, and another high-lying species 5 (75.4) with a linear SiCNO structure are shown to have considerable kinetic stability and may be experimentally observable. The predicted results of isomers 1 and 2 are consistent with the previous mass spectrometry experiments. Moreover, the fourth low-lying species SiOCN 4 (23.9) with bent structure is expected to be observable in low-temperature environments. The bonding nature of the five isomers 1, 2, 3, 4, and 5 is analyzed. The calculated results are compared with those of the analogous molecules C(2)NO and Si(2)NO. Implications in interstellar space and N,O-doped SiC vaporization processes are also discussed.     

Au2、Au3小团族分子的结构和势能函数

原子团簇的结构和性质研究是当今物理学和材料学中的一个热门课题,过渡金属团簇特别是Au团簇,由于其独特的物理和化学性质而被广泛地应用于催化反应、材料吸附[1-2]和光的吸收中[3]。近年来,人们用不同的理论方法研究金原子团簇。H儯kkinen等人利用GGA方法研究了中性和阴离Au2-10团簇的性质[4];Bravo-P啨rez等人采用从头计算的HF和post-HF方法研究Au2-Au6小团簇[5],这些计算结果与实验值相比,相差较大。由于金团簇电子结构的复杂性,对Au体系考虑旋—轨耦合和电子相关效应是很重要的,这种计算的不确定性对Au的影响比IB簇的其它金属团…     

Effect of a complex formation on the calculated low-pressure rate constant of a bimolecular gas-phase reaction governed by tunneling

Many important bimolecular hydrogen-transfer processes that take place in the atmosphere proceed via a potential energy minimum (hydrogen-bonded complex) that precedes along the minimum energy path the unique saddle point of the reaction, the one corresponding to the hydrogen transfer. It is clear that the one-step low-pressure rate constant of such a reaction does not depend on the existence of any complex along the minimum energy path below the reactant if the reaction takes place by thermal activation over a transition state that lies quite above the reactants (for instance 10 kcal/mol). However, we have quantitatively shown in this article that the scenario notoriously changes if the reaction involves significant tunneling. In this work, we have theoretically calculated the rate constants and their temperature dependence for the reaction HO+HOH→HOH+OH by means of a canonical variational transition state theory and a canonical unified statistical theory (when necessary). Multidimensional tunneling effects have been included with a semiclassical transmission coefficient. Two kinds of modified potential energy surfaces (PESs), obtained from an original ab initio potential energy surface, previously calculated by us, have been used. The Eckart-modified PESs serve to model the hydrogen-abstraction profiles with no complexes along the path, while the Gaussian-modified PESs model the energy profiles with two complexes along the path symmetrically distributed at each side of the abstraction saddle point. Our results show that the existence of those complexes reduces the thickness of the classically forbidden region for energies below the adiabatic barrier, and then tunneling is promoted and the reaction is accelerated. The effect of the complex formation in several kinetic magnitudes, as the Arrhenius parameters and the kinetic isotope effect has also been analyzed. © 1999 John Wiley & Sons, Inc. J Comput Chem 20: 1685–1692, 1999     

Potential energy surface and highly excited vibrational lines of NO2 via algebraic approach

The algebraic Hamiltonian of NO₂ is optimized using U(4) algebra via fitting to 102 observed vibrational lines. The RMS error of the fitting is 2.39 cm^?1. We calculated highly excited vibrational energy levels using this optimized Hamiltonian, and then obtained the potential energy surface for the electronic ground state by using the classical limit of the U(4) algebraic Hamiltonian. We also calculated the dissociation energies, the force constants etc. Our results are in good agreement with the other theoretical results. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Potential energy surface for linear triatomic molecules: An algebraic method

Recently, we proposed a new transformation between the angle of canonical coordinates and the bond angle to describe the bending motion in Potential Energy Surfaces (PES) of bent triatomic molecules. In this work we extend the transformation to include linear triatomic molecules. Results for the linear triatomic molecule N₂O are reported.     

Theoretical studies on the potential energy surfaces and vibrational energy levels of HXeF and HXeCl

The potential energy surfaces for the electronic ground state of the HXeCl and HXeF molecules areconstructed by using the internally contracted multi-reference configuration interaction with theDavidson correction(icMRCI Q)method and large basis sets.The stabilities and dissociation barriersare identified from the potential energy surfaces.The three-body dissociation channel is found to bethe dominate dissociation channel for HXeCl,while two dissociation channels are possible and com-petitive for HXeF.Based on the obtained potentials,vibrational energy levels of HXeCl and HXeF arecalculated using the Lanczos algorithm.Our theoretical results are in good agreement with the avail-able observed values.Particularly,the calculated fundamental frequency of the H—Xe stretching vi-bration including the Xe matrix effect of HXeCl is found to be 1666.6 cm-1,which is only 17.6 cm-1higher than the recently observed value of 1649 cm-1.     

Ab Initio Ground-State Potential Energy Function and Vibration-Rotation Energy Levels of Aluminum Monohydride

The accurate ground-state potential energy function of aluminum monohydride (AlH) has been determined from ab initio calculations using the multireference averaged coupled-pair functional (MR-ACPF) method in conjunction with the correlation-consistent core-valence basis sets up to septuple-zeta quality. The vibration-rotation energy levels of the two isotopologues, AlH and AlD, were predicted to near the “spectroscopic” accuracy. The importance of electron correlation beyond the MR-ACPF level of approximation, the scalar relativistic, spin-orbit, adiabatic, and nonadiabatic effects was discussed. © 2019 Wiley Periodicals, Inc.