Ring Polymer Molecular Dynamics of the C(^(1)D)+H_(2) Reaction on the Most Recent Potential Energy Surfaces

Ring polymer molecular dynamics(RPMD)calculations for the C(^(1)D)+H_(2)reaction are performed on the Zhang-Ma-Bian ab initio potential energy surfaces(PESs)recently constructed by our group,which are unique in very good descriptions of the regions around conical intersections and of van der Waals(vdW)interactions.The calculated reaction thermal rate coefficients are in very good agreement with the latest experimental results.The rate coefficients obtained from the ground˜a^(1)A′ZMB-a PES are much larger than those from the previous RKHS PES,which can be attributed to that the vdW saddles on our PESs have very different dynamical effects from the vdW wells on the previous PESs,indicating that the RPMD approach is able to include dynamical effects of the topological structures caused by vdW interactions.The importance of the excited˜b^(1)A′′ZMB-b PES and quantum effects in the title reaction is also underscored.

Three-Dimensional Diabatic Potential Energy Surfaces of Thiophenol withNeural Networks

Three-dimensional(3D)diabatic potential energy surfaces(PESs)of thiophenol involving the S0,and coupled 1ππ^(*) and 1πσ^(*) states were constructed by a neural network approach.Specifically,the diabatization of the PESs for the 1ππ^(*) and 1πσ^(*) states was achieved by the fitting approach with neural networks,which was merely based on adiabatic energies but with the correct symmetry constraint on the off-diagonal term in the diabatic potential energy matrix.The root mean square errors(RMSEs)of the neural network fitting for all three states were found to be quite small(<4 meV),which suggests the high accuracy of the neural network method.The computed low-lying energy levels of the S_(0) state and lifetime of the 0^(0) state of S_(1) on the neural network PESs are found to be in good agreement with those from the earlier diabatic PESs,which validates the accuracy and reliability of the PESs fitted by the neural network approach.

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.

Adiabatic Potential Energy Surfaces and Photodissociation Mechanisms for Highly Excited States of H_(2)O

Full-dimensional adiabatic potential energy surfaces of the electronic ground state X and nine excited states A,I,B,C,D,D',D'',E' and F of H_(2)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 D-E',E'-F,A-I and I-C 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.

Potential energy surfaces of ozone in the ground state

Equilibrium parameters of ozone, such as equilibrium geometry structure parameters, force constants and dissociation energy are presented by CBS-Q ab initio calculations. The calculated equilibrium geometry structure parameters and energy are in agreement with the corresponding experimental values. The potential energy function of ozone with a C2v symmetry in the ground state is described by the simplified Sorbie-Murrell many-body expansion potential function according to the ozone molecule symmetry. The contour of bond stretching vibration potential of an O3 in the ground state, with a bond angle （θ） fixed, and the contour of O3 potential for O rotating around O1-O （R1）, with O1-O bond length taken as the one at equilibrium, are plotted. Moreover, the potentials are analysed.

Full-Dimensional Potential Energy Surfaces of Ground(X^(2)A′)and Excited(A^(2)A″)Electronic States of HCO and Absorption Spectrum

In this work,high-fidelity full-dimensional potential energy surfaces(PESs)of the ground(X^(2)A′)and first doublet excited(A^(2)A″)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.

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.

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.

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.

Ground state properties and potential energy surfaces of 270Hs from multidimensionally-constrained relativistic mean field model

We study the ground state properties, potential energy curves and potential energy surfaces of the superheavy nucleus270 Hs by using the multidimensionally-constrained relativistic mean-field model with the effective interaction PC-PK1. The binding energy,size and shape as well as single particle shell structure corresponding to the ground state of this nucleus are obtained.270 Hs is well deformed and exhibits deformed doubly magic feature in the single neutron and proton level schemes. One-dimensional potential energy curves and two-dimensional potential energy surfaces are calculated for270 Hs with various spatial symmetries imposed. We investigate in detail the effects of the reflection asymmetric and triaxial distortions on the fission barrier and fission path of270 Hs. When the axial symmetry is imposed, the reflection symmetric and reflection asymmetric fission barriers both show a double-hump structure and the former is higher. However, when triaxial shapes are allowed the reflection symmetric barrier is lowered very much and then the reflection symmetric fission path becomes favorable.

Variational transition state theory study of the reactions Li+HF and Li+HCl on the BO potential energy surfaces

Variational transition state theory is used to study the kinetics of the reactions Li+HF and Li+HCl on the BO potential energy surfaces, fitted to ab initio values by Lagana et al. The saddle point, the location of generalized transition state and other kinetic features on the bottleneck of the po- tential energy surface are obtained. We point out that the location of generalized transition state tends to the direction in which the two atoms with larger reduced mass separate off. It can be determined whether a potential energy surface is accurate or not by a fine analysis of its bottleneck region and a comparison of kinetic calculation with experimental result. This is helpful for building a more accurate potential energy surface.

Calculation of multidimensional potential energy surfaces for even-even transuranium nuclei: systematic investigation of the triaxiality effect on the fission barrier

Static fission barriers for 95 even-even transuranium nuclei with charge number Z = 94-118 have been systematically investigated by means of pairing self-consistent Woods-Saxon-Strutinsky calculations using the potential energy surface approach in multidimensional（β;, γ, β;） deformation space. Taking the heavier （252）;f nucleus（with the available fission barrier from experiment） as an example, the formation of the fission barrier and the influence of macroscopic, shell and pairing correction energies on it are analyzed. The results of the present calculated β;values and barrier heights are compared with previous calculations and available experiments. The role of triaxiality in the region of the first saddle is discussed. It is found that the second fission barrier is also considerably affected by the triaxial deformation degree of freedom in some nuclei（e.g., the Z =112-118 isotopes）. Based on the potential energy curves, general trends of the evolution of the fission barrier heights and widths as a function of the nucleon numbers are investigated. In addition, the effects of Woods-Saxon potential parameter modifications（e.g.,the strength of the spin-orbit coupling and the nuclear surface diffuseness） on the fission barrier are briefly discussed.

Deconstructing Vibrational Motions on the Potential Energy Surfaces of Hydrogen-Bonded Complexes

Internal vibrations underlie transient structure formation,spectroscopy,and dynamics.However,at least two challenges exist when aiming to elucidate the contributions of vibrational motions on the potential energy surfaces.One is the acquisition of well-resolved experimental infrared spectra,and the other is the development of efficient theoretical methodologies that reliably predict band positions,relative intensities,and substructures.Here,we report size-specific infrared spectra of ammonia clusters to address these two challenges.Unprecedented agreement between experiment and state-of-the-art quantum simulations reveals that the vibrational spectra are mainly contributed by proton-donor ammonia.A striking Fermi resonance observed at approximately 3210 and 3250 cm^(−1)originates from the coupling of NH symmetric stretch fundamentals with overtones of free and hydrogen-bonded NH bending,respectively.These novel,intriguing findings contribute to a better understanding of vibrational motions in a large variety of hydrogen-bonded complexes with orders of magnitude improvements in spectral resolution,efficiency,and specificity.

A procedure of determining potential energy surfaces of triatomic molecules from inversion of spectroscopic data

In this paper, we have suggested an iterative procedure of optimization of the linear parameters in an analytic potential energy function for a triatomic molecule, by combining both variational and second order perturbation methods. The most important feature of this procedure is that the objective function is an analytical expression which can be optimized easily. The application to the water molecule is presented.

Differential diffusion quantum Monte Carlo method: determination of potential energy surfaces of molecules

A differential approach for self-optimizing diffusion Monte Carlo calculation was proposed in this paper, which is a new algorithm combining three techniques such as optimizing, diffusion and correlation sampling. This method can be used to directly compute the energy differential between two system in the diffusion process, making the statistical error of calculation be reduced to Order of 10?-5 hartree, and recover about more than 80% of the correlation. We employed this approach to set up a potential energy surface of a molecule, used a “rigid move” model, and utilized Jacobi transformation to make energy calculation for two configurations of a molecule having good positive correlation. So, an accurate energy differential could be obtained, and the potential energy surface with good quality can be depicted. In calculation, a technique called “post-equilibrium remaining sample” was set up firstly, which can save about 50% of computation expense. This novel algorithm was used to study the potential as molecular spectroscopy and the energy variation in chemical reactions.

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.

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.