A fully ab initio potential energy surface for ClH2 reactive system

An ab initio analytical potential energy surface called BW3 for the ClH2 reactive system is presented. The fit of this surface is based on about 1 200 ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H2 and HCl, and the endothermicity of the Cl + H2 abstraction reaction. For the Cl + H2 abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H + ClH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

A fully ab initio potential energy surface for C1H_2 reactive system

An ab initio analytical potential energy surface called BW3 for the CIH2 reactive system is presented. The fit of this surface is based on about 1 200 ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H2 and HCI, and the endothermicity of the Cl H2 abstraction reaction. For the Cl H2 abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H CIH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

Ab initio intermolecular potential energy surface of He-LiH

The intermolecular potential energy surface of He-LiH complex was studied using the full-electronic complete forth-order Miller-Plesset perturbation (MPPT) method.In ab initio calculations,the bond length of LiH was fixed at 0 159 5 nm.The potential has two local minima of Vm=-179.93 cm for the linear He LiH geormetrv at Rm=0.227 nm and Vm=-10.44 cm-1 for the linear He-HL1 geometry at Rm=0.516 nm The potemal exhibits strong anisotropy The analytic potential function with 31 parameters was determined by fitting to the calculated ab,mtio potentials The influence of variation of LiH bond length on the potential energy surface was also studied     

Complete basis set ab initio and hybrid density functional theory exploration of the potential energy surface in the reaction between an amino radical and nitrogen oxide

The potential energy surface for the reaction involving NH₂ plus NO was explored with a quadratic complete basis set ab initio approach and three hybrid density functional theory methods, the target being to accurately estimate activation barriers and the relative stability of the nitrogen–oxygen isomers. The computational results were compared with previously performed ab initio calculations and new, more accredited values for the NH₂NO rearrangement to HNNOH and for the HNNOH decomposition reaction were suggested. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 409–414, 1998     

General methodology in two dimensions for classical simulation of reactive and nonreactive events on ab initio potential energy surfaces

A completely general two-dimensional (2D) methodology for the classical simulation of reactive and nonreactive events on ab initio potential energy surfaces is introduced and tested. The methodology requires the minimum amount of information given a priori—geometries and energies at these geometries. From a list of ab initio geometries and energies, simulations may be executed and a distribution of outcomes obtained. The method introduced attempts a local approach at simulating the dynamics of the system, rather than a global analytic fit to the potential energy surface. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1431–1444, 1998     

DNA base trimers: empirical and quantum chemical ab initio calculations versus experiment in vacuo

A complete scan of the potential‐energy surfaces for selected DNA base trimers has been performed by a molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The resulting most stable/populated structures were then reoptimized at a correlated ab initio level by employing resolution of the identity, Møller–Plesset second‐order perturbation theory (RI‐MP2). A systematic study of these trimers at such a complete level of electronic structure theory is presented for the first time. We show that prior experimental and theoretical interpretations were incorrect in assuming that the most stable structures of the methylated trimers corresponded to planar systems characterized by cyclic intermolecular hydrogen bonding. We found that stacked structures of two bases with the third base in a T‐shape arrangement are the global minima in all of the methylated systems: they are more stable than the cyclic planar structures by about 10 kcal mol^?1. The different behaviors of nonmethylated and methylated trimers is also discussed. The high‐level geometries and interaction energies computed for the trimers serve also as a reference for the testing of recently developed density functional theory (DFT) functionals with respect to their ability to correctly describe the balance between the electrostatic and dispersion contributions that bind these trimers together. The recently reported M052X functional with a polarized triple‐zeta basis set predicts 11 uracil trimer interaction energies with a root‐mean‐square error of 2.3 kcal mol^?1 relative to highly correlated ab initio theoretical calculations.     

An ab initio potential energy surface and vibrational energy levels of ZnH2

A three‐dimensional potential energy surface of the electronic ground state of ZnH₂ (${X}^1\sum _g^ +$ ) molecule is constructed from more than 7500 ab initio points calculated at the internally contracted multireference configuration interaction with the Davidson correction (icMRCI+Q) level employing large basis sets. The calculated relative energies of various dissociation reactions are in good agreement with the previous theoretical/experimental values. Low‐lying vibrational energy levels of ZnH₂, ZnD₂, and HZnD are calculated on the three‐dimensional potential energy surface using the Lanczos algorithm, and found to be in good agreement with the available experimental band origins and the previous theoretical values. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

An accurate three-dimensional potential energy surface for the He-Na_2 complex

An accurate three-dimensional potential energy surface(PES) for the He-Na2 van der Waals comple was calculated at the coupled cluster singles-and-doubles with noniterative inclusion of connecte triple(CCSD(T)) level of theory.A mixed basis set,aug-cc-pVQZ for the He atom and cc-pCVQZ for th sodium atom,and an additional(3s3p2d1f) set of midbond functions were used.The computed inte action energies in 819 configurations were fitted to a 96-parameter analytic potential model by leas squares fitting.The PES has two shallow wells corresponding to the T-shaped structure and the linea configuration,which are located at 12.5a0 and 14 a0 with depths of 1.769 and 1.684 cm-1,respectivel The whole potential energy surface exhibits weak anisotropy.Based on the fitted PES,state-to-stat differential cross sections were calculated.     

Molecular dynamics simulations of fluid methane properties using ab initio intermolecular interaction potentials

Intermolecular interaction energy data for the methane dimer have been calculated at a spectroscopic accuracy and employed to construct an ab initio potential energy surface (PES) for molecular dynamics (MD) simulations of fluid methane properties. The full potential curves of the methane dimer at 12 symmetric conformations were calculated by the supermolecule counterpoise‐corrected second‐order Møller‐Plesset (MP2) perturbation theory. Single‐point coupled cluster with single and double and perturbative triple excitations [CCSD(T)] calculations were also carried out to calibrate the MP2 potentials. We employed Pople's medium size basis sets [up to 6‐311++G(3df, 3pd)] and Dunning's correlation consistent basis sets (cc‐pVXZ and aug‐cc‐pVXZ, X = D, T, Q). For each conformer, the intermolecular carbon–carbon separation was sampled in a step 0.1 Å for a range of 3–9 Å, resulting in a total of 732 configuration points calculated. The MP2 binding curves display significant anisotropy with respect to the relative orientations of the dimer. The potential curves at the complete basis set (CBS) limit were estimated using well‐established analytical extrapolation schemes. A 4‐site potential model with sites located at the hydrogen atoms was used to fit the ab initio potential data. This model stems from a hydrogen–hydrogen repulsion mechanism to explain the stability of the dimer structure. MD simulations using the ab initio PES show quantitative agreements on both the atom‐wise radial distribution functions and the self‐diffusion coefficients over a wide range of experimental conditions. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

An ab initio MO study on the direct and migratory reaction mechanism for the formation of HF in H+ClF system

ＰｏｌａｎｙｉａｎｄｃｏｗｏｒｋｅｒｓｈａｖｅｓｔｕｄｉｅｄａｓｅｒｉｅｓｏｆｒｅａｃｔｉｏｎｓｂｅｔｗｅｅｎＨａｔｏｍｓａｎｄｉｎｔｅｒｈａｌｏｇｅｎｓＸＹ(Ｙｉｓｔｈｅｍｏｒｅｅｌｅｃｔｒｏｎｅｇａｔｉｖｅｈａｌｏｇｅｎａｔｏｍ)ｕｓｉｎｇｔｈｅｉｎｆｒａｒｅｄｃｈｅｍｉｌｕｍｉｎｅｓｃｅｎｃｅｍｅｔｈｏｄ[1,2]ａｎｄｈａｖｅｏｂｓｅｒｖｅｄｂｉｍｏｄａｌｅｎｅｒｇｙｄｉｓｔｒｉｂｕｔｉｏｎｓｆｏｒｔｈｅＨＹｐｒｏｄｕｃｔ;ｔｈａｔｉｓｔｏｓａｙ,ｔｈｅｔｏｔａｌａｖａｉｌａｂｌｅｅｎｅｒｇｙｆ…     

Organolithium complex as a gas sensing material for oxides from ab initio calculations and molecular dynamics simulations

Gas sensing study of C₂H₄Li complex toward oxides viz. CO, CO₂, NO, NO₂, SO, and SO₂ gas molecules has been carried out using ab initio method. Different possible configurations of gas molecule adsorption on C₂H₄Li complex are considered. The structural parameters of most stable configuration of gas molecule adsorbed complexes are thoroughly analysed. Electronic properties are studied using total density of states (DOS) plot. Charge transferred between the gas molecule and the substrate is studied using NBO charge analysis. Gas sensing of all the six gas molecules is possible at ambient conditions. Atom centred density matrix propagation (ADMP) molecular dynamics simulations confirmed that all the gas molecules remain adsorbed on C₂H₄Li complex at room temperature during the simulation. This study suggests that the C₂H₄Li complex acts as a novel gas sensing material for CO, CO₂, NO, NO₂, SO, and SO₂ gas molecules at ambient conditions, below room temperature as well as at high pressure.     

Benchmarking PES-Learn's machine learning models predicting accurate potential energy surface for quantum scattering

Machine learning (ML) models, neural networks, and Gaussian processes have been used to predict the potential energy surface taking C₂-He (both static and dynamic scenario) and NCCN-He collision systems. The surface is restricted to ∽125 points where traditional spline becomes inefficacious. Quantum dynamics is performed by solving close-coupling equation to compute cross sections benchmarking the performance of the ML models. The current study forms a basis for any future investigation of larger molecules where conventional fitting fails due to sparser ab initio points and cuts down the computational time without compromising on the quality of the surface.     

Lithium hydroxide phase transition under high pressure: an ab initio molecular dynamics study.

The high-pressure phase transition in the deuterated lithium hydroxide crystalline state has been studied by Car-Parrinello molecular dynamics simulations, in the constant-pressure, constant-temperature ensemble. The recently developed metadynamics approach has been applied to encourage the system to transform into different phases in an affordable simulation time. A previously not completely characterized high-pressure phase has been obtained. The structural and spectroscopic properties have been studied and compared with the neutron scattering, infrared and Raman measurements. It has been found that the calculated structure differs slightly from the experimental hypothesis, and that the presence of strong hydrogen bonds is the source of the red shift and of the characteristic features of the OD-stretching bands in both IR and Raman spectra.     

Accurate ab initio potential energy surface and vibration‐rotation energy levels of hydrogen peroxide

The accurate ground‐state potential energy surface of hydrogen peroxide, H₂O₂, has been determined from ab initio calculations using the coupled‐cluster approach in conjunction with the correlation‐consistent basis sets up to septuple‐zeta quality. Results obtained with the conventional and explicitly correlated coupled‐cluster methods were compared. The core–electron correlation, scalar relativistic, and higher‐order valence–electron correlation effects were taken into account. The adiabatic effects were also discussed. The vibration–rotation energy levels of the H₂O₂, D₂O₂, and HOOD isotopologues were predicted, and the experimental vibrational fundamental wavenumbers were reproduced to 1 cm^?1 (“spectroscopic”) accuracy. © 2012 Wiley Periodicals, Inc.     

Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Exploring the hydration of Pb2+: ab initio studies and first-principles molecular dynamics

Even though lead is a well-known toxicant widely scattered throughout the world since antiquity, its chemistry is poorly documented at the molecular level. Here we investigate the hydration of the Pb(2+) ion by means of first-principles molecular dynamics (Car-Parrinello molecular dynamics, CPMD). We found that the hydrated cation is heptacoordinated in a dynamically holodirected arrangement roughly corresponding to a fluxional distorted pentagonal bipyramid. The time-averaged Pb-O bond length is especially large and amounts to 2.70 A with an associated root-mean-square deviation of 0.26 A. This results from a dynamic exchange between short (<2.6 A), intermediate (2.6-3.0 A) and long (>3.0 A) Pb-O bonds. The latter very long Pb-O distance implies that the determination of the coordination number n(c) from experimental work may not necessarily yield values directly comparable to the theoretical value of n(c)=7, since not all experimental techniques would recognize such a long distance as a bond to the metal cation. Pronounced disorders are evidenced in the second shell, characteristic of a chaotropic cation, and exchanges between the first and second shells cannot be excluded on a timescale of a few tens of picoseconds.