Ab initio potential energy surface and anharmonic vibration spectrum of NF_(3)^(+)

Potential energy surfaces(PESs), vibrational frequencies, and infrared spectra are calculated for NF_(3)^(+) using ab initio calculations, based on UCCSD(T)/cc-p VTZ combined with vibrational configuration interaction(VCI). Based on an iterative algorithm, the surfaces(SURF) program adds automatic points to the lattice representation of the potential function, the one-dimensional and two-dimensional PESs are calculated after reaching a convergence threshold, finally the smooth image of the potential energy surface is fitted. The PESs accurately account for the interaction between the different modes, with the mode q_(6) symmetrical stretching vibrations having the greatest effect on the potential energy change of the whole system throughout the potential energy surface shift. The anharmonic frequencies are obtained when the VCI matrix is diagonalized. Fundamental frequencies, overtones, and combination bands of NF_(3)^(+) are calculated, which generate the degenerate phenomenon between their frequencies. Finally, the calculated anharmonic frequency is used to plot the infrared spectra.Modal antisymmetric stretching ν_(5) and symmetric stretching ν_(6) exhibit a phenomenon of large-intensity borrowing. This study can provide data to support the characterization in the laboratory.

The effect of the attractive well of the potential energy surface for Ne-HC1 on rotationally inelastic partial wave cross sections

An interaction potential of the Ne-HC1 van der Waals complex is obtained by utilizing the Huxley analytic potential function to fit the accurate interaction energy data, which have been computed at the coupled cluster singles and doubles including connected triple excitations level and with the augmented correlation consistent polarized valence quintuple zeta basis set extended with a set of 3s3p2dlflg mid-bond functions [CCSD （T）/aug-cc-pV5Z-33211]. The close coupling calculation of state-to-state partial cross sections for collision of Ne with HC1 is first performed by employing the fitted interaction potential. This calculation is performed at the incident energies： 40, 60, 75 and 100 meV, separately. The effects of the long-range attractive and the short-range anisotropic interactions on the inelastic state-to-state partial cross sections are discussed in detail. Two maxima are present in the rotationally inelastic partial cross sections and they originate from different mechanisms.

Global Potential Energy Surface for the H＋CH4←→H2＋CH3 Reaction using Neural Networks

A globla potential energy surface （PES） for the H＋CH4←→H2＋CH3 reaction has been constructed using the neural networks method based on 47783 high level ab initio geometry points. Extensive quasi-classical trajectories and quantum scattering calculations were carried out to check the convergence of the PES. This PES, fully converged with respect to the fitting procedure and the number of ab initio points, has a very small fitting error, and is much faster on evaluation than the modified Shepard interpolating PES, representing the best available PES for this benchmark polyatomic system.

Theoretical Study on Dissociation Potential Energy Surface of Peroxynitric Acid

The lowest energy structures of peroxynitric acid have been studied with B3LYP/6-311＋ G（2d,2p） method. The potential energy surfaces （PES） along the O-N and O-O bonds have been scanned at CCSD（T）/aug-cc-pVDZ level, respectively. The calculated results show that on the O-N PES, the O3-N4 bond length of the loose transition state is 2.82A^° and the corresponding energy barrier is 25.6 kcal/mol, while on the O-O PES, the loose transition state with of O2-O3 bond length of 2.35A^° has the energy barrier of 37.4 kcal/mol. Thus the primary reaction path for peroxynitric acid is the dissociation into HO2 and NO2.

A Global ab initio Potential Energy Surface for F+H_(2)→HF+H

A global three dimensional potential energy surface for the F＋H2→HF＋H reaction has been developed by spline interpolation of about 15,000 symmetry-unique ab initio points, obtained from the multi-reference configuration interaction level with Davidson correction using the aug-cc-pV5Z basis set. In the entrance channel the spin-orbit coupling energy is also included.

Three-dimensional Potential Energy Surface and Bound States of the Ar2-Ne Complex

The first three-dimensional interaction potential energy surface （PES） of the Ar2-Ne complex is developed using the single and double excitation coupled cluster theory with noniterative treatment of triple excitations CCSD（T）. The aug-cc-pVQZ basis sets are employed for all atoms, including an additional （3s3p2d2flg） set of midpoint bond functions. The calculated single point energies are fitted to an analytic two-dimensional potential model at each of seven fixed rAr~ values. The seven model potentials are then used to construct the three- dimensional PES by interpolating along （r-re） using a sixth-order polynomial. The PES is used in the following rovibrational energy levels calculations. The comparisons of theoretical transition frequencies and spectroscopic constants with the experimental results are given.

A UNIVERSAL ANALYTIC POTENTIAL-ENERGY FUNCTION BASED ON A PHASE FACTOR

Using a field equation with a phase factor, a universal analytic potential-energy function applied to the interactions between diatoms or molecules is derived, and five kinds of potential curves of common shapes are obtained adjusting the phase factors. The linear thermal expansion coefficients and Young's moduli of eleven kinds of face-centered cubic (fcc) metals - Al, Cu, Ag, etc. are calculated using the potential-energy function; the computational results are quite consistent with experimental values. Moreover, an analytic relation between the linear thermal expansion coefficients and Young's moduli of fcc metals is given using the potential-energy function. Finally, the force constants of fifty-five kinds of diatomic moleculars with low excitation state are computed using this theory, and they are quite consistent with RKR (Rydberg-Klein-Rees) experimental values.

Structure and Properties of Hydrotalcite Using Electrostatic Potential Energy Model

The electrostatic potential energy model of hydrotalcites was based on the theory of crystallography. The anionic potential energy of MgAl-hydrotalcites, with 20 layers and 2107 anions per layer, was calculated, and the anionic stability of the hydrotalcites was investigated. The charge density of the layer and the distance between the adjacent anions varied with the molar ratio of Al^3＋/（Mg^2＋ ＋ Al^3＋）. Anionic potential energy depended on the charge and size of the anions. Calculation results remained consistent with thermal stability and the ion exchange ability reported. This model is able to predict anionic stability of the hydrotalcites.

Assessment of grid inherent vulnerability considering open circuit fault under potential energy framework

A potential energy framework for assessment of grid vulnerability was presented.In the framework,the branch potential energy function model was constructed.Two indexes,current vulnerability and forecasting vulnerability,were calculated.The current vulnerability was used to identify the current vulnerable area through calculating the distance between the current transmitted power and initial transmitted power;and the forecast vulnerability under variation of power injection was used to predict the vulnerable area of next step and verify the current vulnerable area.Numerical simulation was performed under variant operating conditions with IEEE-30 bus system,which shows that almost area of 90% overlaps between current vulnerable area and forecasting vulnerable area,the overlapped area is termed as inherent vulnerable area of grid.When considering N-1 contingency,the assessment results of this method proposed agree with those of optimal power flow.When considering N-2 contingency,optimal power flow fails to obtain correct results,while the method based on energy framework gives reliable results.

Analysis of Potential Energy Surface for Butanone Isomerization

The potential energy surfaces for butanone isomerization have been investigated by density function theory calculation. Six main reaction pathways are confirmed using the intrinsic reaction coordinate method, and the corresponding isomerization products are 1-buten-2-ol, 2-buten-2-ol, butanal or 1-buten-l-ol, methyl 1-propenyl ether, methyl allyl ether, and ethyl vinyl ether, respectively. Among them, there are three pathways through butylene oxide, indicating butylene oxide is an important intermediate product during butanone isomer ization. The calculated vertical ionization energies of the reactant and its products are in a good agreement with the experimental values available. From the consideration for the relative energies Of transition states and the number of high-energy barriers we infer that the reaction pathway butanone-＊l-buten-2-ol---2-buten-2-oi is the most competitive. The obtained results are informative for future studies on isomerization of ketone molecules.

Structure and analytical potential energy function for the ground state of the BC^x （x=0, -1）

In this paper, the electronic states of the ground states and dissociation limits of BC and BC- are correctly determined based on group theory and atomic and molecular reaction statics. The equilibrium geometries, harmonic frequencies and dissociation energies of the ground state of BC and BC- are calculated by using density function theory and quadratic CI method including single and double substitutions. The analytical potential energy functions of these states have been fitted with Murrell-Sorbie potential energy function from our ab initio calculation results. The spectroscopic data （αe, ωe and ωeχe） of each state is calculated via the relation between analytical potential energy function and spectroscopic data. All the calculations are in good agreement with the experimental data.

Ab initio calculations of accurate dissociation energy and analytic potential energy function for the second excited state B^1∏ of ^7LiH

The reasonable dissociation limit of the second excited singlet state B1∏ of ^7LiH molecule is obtained. The accurate dissociation energy and equilibrium geometry of the B^∏ state are calculated using a symmetry-adaptedcluster configuration interaction method in full active space. The whole potential energy curve for the B1H state is obtained over the internuclear distance ranging from about 0.10 nm to 0,54 nm, and has a least-square fit to the analytic Murrell-Sorbie function form. The vertical excitation energy is calculated from the ground state to the B^1∏ state and compared with previous theoretical results. The equilibrium internuclear distance obtained by geometry optimization is found to be quite different from that obtained by single-point energy scanning under the same calculation condition. Based on the analytic potential energy function, the harmonic frequency value of the B^1∏ state is estimated. A comparison of the theoretical calculations of dissociation energies, equilibrium interatomic distances and the analytic potential energy function with those obtained by previous theoretical results clearly shows that the present work is more comprehensive and in better agreement with experiments than previous theories, thus it is an improvement on previous theories.

The Analytical Potential Energy Function of NH Radical Molecule in External Electric Field

The geometric structures of an Nit radical in different external electric fields are optimized by using the density functional B3P86/cc-PVSZ method, and the bond lengths, dipole moments, vibration frequencies and IR spectrum are obtained. The potential energy curves are gained by the CCSD （T） method with the same basis set. These results indicate that the physical property parameters and potential energy curves may change with the external electric field, especially in the reverse direction electric field. The potential energy function of zero field is fitted by the Morse potential, and the fitting parameters are in good accordance with the experimental data. The potential energy functions of different external electric fields are fitted adopting the constructed potential model. The fitted critical dissociation electric parameters are shown to be consistent with the numerical calculation, and the relative errors are only 0.27% and 6.61%, hence the constructed model is reliable and accurate. The present results provide an important reference for further study of the molecular spectrum, dynamics and molecular cooling with Stark effect.

Accurate potential energy function and spectroscopic study of the X^2Σ^+,A^2Ⅱ and B^2Σ^+ states of the CP radical

This paper calculates the equilibrium internuclear separations, the harmonic frequencies and the potential energy curves of the X^2∑＋, A^2П and B^2∑＋ states of the CP radical by the highly accurate valence internally contracted multireference configuration interaction method with correlation-consistent basis sets （aug-cc-pV6Z for C atom and aug-cc-pVQZ for P atom）. The potential energy curves are all fitted with the analytic potential energy function by the least-square fitting. Employing the analytic potential energy function, we determine the spectroscopic constants （Be, αe and ωeχe） of these states. For the X2∑＋ state, the obtained values of De, Be, αe, ωeχe, Re and ωe are 5.4831 eV, 0.792119 cm-1, 0.005521 cm-1, 6.89653 cm-1, 0.15683 nm, 12535.11 cm-1, respectively. For the A2H state, the present values of De, Be,αe, ωeχe, Re and We are 4.586 eV, 0.703333 cm-1, 0.005458 cm-1, 6.03398 cm-1, 0.16613 nm, 1057.89 cm-1, respectively. For the B2E＋ state, the present values of De, Be, αe, ωeχe, Re and We are 3.506 eV, 0.677561 cm-1, 0.00603298 cm-1, 5.68809 cm-1, 0.1696 nm, 822.554 cm-1, respectively. For these states, the vibrational states with the rotational quantum number J equals zero （J = 0） are studied by solving the radial nuclear Schr6dinger equation using the Numerov method. For each vibrational state, the vibrational level, the classical turning points, the rotational inertial and the centrifugal distortion constants are calculated. Comparison is made with recent theoretical and experimental results.

Generalized Available Potential Energy

The kinetic energy generation in either the dry or moist atmosphere may be estimated by the same relationships if we introduce the new concept of generalized available potential energy. The largest magnitude of generalized available potential energy and corresponding reference state of either dry or moist atmosphere are calculated in terms of the mitial conditions and entropy variation of the atmosphere. The obtained relationships are applicable for the statically unstable atmosphere as well. The generalized available potential energy associated with reversible processes reaches the maximum with respect to same initial state. While the generation of kinetic energy in irreversible processes is characterized by sudden changes. When the reference state is assumed to be saturated, we may predict the final temperature and moisture fields corresponding to provided initial state and entropy variation.

Assessment of wind energy potential in China

China wind atlas was made by numerical simulation and the wind energy potential in China was calculated. The model system for wind energy resource assessment was set up based on Canadian Wind Energy Simulating Toolkit (WEST) and the simulating method was as follows. First, the weather classes were obtained depend on meteorological data of 30 years. Then, driven by the initial meteorological field produced by each weather class, the meso-scale model ran for the distribution of wind energy resources according each weather class condition one by one. Finally, averaging all the modeling output weighted by the occurrence frequency of each weather class, the annual mean distribution of wind energy resources was worked out. Compared the simulated wind energy potential with other results from several activities and studies for wind energy resource assessment, it is found that the simulated wind energy potential in mainland of China is 3 times that from the second and the third investigations for wind energy resources by CMA, and is similar to the wind energy potential obtained by NREL in Solar and Wind Energy Resource Assessment(SWERA) project. The simulated offshore wind energy potential of China seems smaller than the true value. According to the simulated results of CMA and considering lots of limited factors to wind energy development, the final conclusion can be obtained that the wind energy availability in China is 700～1 200 GW, in which 600～1 000 GW is in mainland and 100～200 GW is on offshore, and wind power will become the important part of energy composition in future.

Meso-scale available gravitational potential energy in the world oceans

The pitfalls of applying the commonly used definition of available gravitational potential energy （AGPE） to the world oceans are re-examined. It is proposed that such definition should apply to the meso-scale problems in the oceans, not the global scale. Based on WOA98 climatological data, the meso-scale AGPE in the world oceans is estimated. Unlike previous results by Oort et al. , the meso-scale AGPE is large wherever there is a strong horizontal density gradient. The distribution of meso-scale AGPE reveals the close connection between the baroclinic instability and the release of gravitational potential energy stored within the scale of Rossby deformation radius.

Accelerating the Construction of Neural Network Potential Energy Surfaces： A Fast Hybrid Training Algorithm

Machine learning approaches have been promising in constructing high-dimensional potential energy surfaces （PESs） for molecules and materials. Neural networks （NNs） are one of the most popular such tools because of its simplicity and efficiency. The training algorithm for NNs becomes essential to achieve a fast and accurate fit with numerous data. The Levenberg-Marquardt （LM） algorithm has been recognized as one of the fastest and robust algorithms to train medium sized NNs and widely applied in recent NN based high quality PESs. However, when the number of ab initio data becomes large, the efficiency of LM is limited, making the training time consuming. Extreme learning machine （ELM） is a recently proposed algorithm which determines the weights and biases of a single hidden layer NN by a linear solution and is thus extremely fast. It, however, does not produce sufficiently small fitting error because of its random nature. Taking advantages of both algorithms, we report a generalized hybrid algorithm in training multilayer NNs. Tests on H＋H2 and CH4＋Ni（111） systems demonstrate the much higher efficiency of this hybrid algorithm （ELM-LM） over the original LM. We expect that ELM-LM will find its widespread applications in building up high-dimensional NN based PESs.

Ab initio calculation on accurate analytic potential energy functions and harmonic frequencies of c^3∑g^＋ and B^1-Пu states of dimer 7Li2

The comparison between single-point energy scanning （SPES） and geometry optimization （OPT） in determining the equilibrium geometries of c^3∑g^＋ and B^1-Пu states of dimer 7Li2 is made at numerous basis sets by using a symmetryadapted-cluster configuration-interaztion （SAC-CI） method in the Gaussian 03 program package. In this paper the difference of the equilibrium geometries obtained by SPES and by OPT is reported. The results obtained by SPES are found to be more reasonable than those obtained by OPT in full active space at the present SAC-CI level of theory. And the conclusion is attained that the cc-PVTZ is a most suitable basis set for these states. The calculated dissociation energies and equilibrium geometries are 0.8818 eV and 0.3090 nm for c^3∑g^＋ state, and 0.3668 eV and 0.2932 nm for B^1-Пu state respectively. The potential energy curves are calculated over a wide internuclear distance range from about 2.5α0 to 37α0 and have a least-squares fit into the Murrell-Sorbie function. According to the calculated analytic potential energy functions, the harmonic frequencies （We） and other spectroscopic data （ωeXe, Be and αe） are calculated. Comparison of the theoretical determinations at present work with the experiments and other theories clearly shows that the present work is the most complete effort and thus represents an improvement over previous theoretical results.

Time-dependent Wave Packet Quantum Scattering Study of Reaction S（3p）＋H2→HS＋H on a New ab initio Potential Energy Surface 3A’

A new potential energy surface is presented for the triplet state 3At of the chemical reaction S（3P）＋H2 from a set of accurate ab initio data. The single point energies are computed using highly correlated complete active space self-consistent-field and multi-reference configuration interaction wave functions with a basis set of aug-cc-pV5Z. We have fitted the full set of energy values using many-body expansion method with an Aguado-Paniagua function. Based on the new potential energy surface, we carry out the time-dependent wave packet scattering calculations over the collision energy range of 0.8-2.2 eV. Both the centrifugalsudden approximation and Coriolis Coupling cross sections are obtained. In addition, the total reaction probabilities are calculated for the reactant H2 initially in the vibrational states v=0-3 （j=0）. It is found that initial vibrational excitation enhances the title reaction.