A new scaling procedure to correct semiempirical MEP and MEP-derived properties

Summary The MNDO, AM1, PM3, and ab initio 6–31G^* and 6–31+G^* MEPs for 21 neutral and 12 charged molecules were computed in layers ranging from 1.2 to 2.0 times the van der Waals radii of atoms. Semiempirical and ab initio MEPs for each layer and two groups of layers were compared to gain insight into the relationships between semiempirical and ab initio MEPs. A detailed statistical study allowed us to obtain a new set of scaling coefficients able to correct the semiempirical MEPs to provide better representations of the ab initio values. The corrected semiempirical MEPs were used to obtain electrostatic charges, whose quality was tested by the comparison between semiempirical Coulombic MEPs and ab initio quantum mechanical MEPs.     

Interpolating moving least-squares methods for fitting potential energy surfaces: an application to the H2CN unimolecular reaction

Classical trajectories have been used to compute rates for the unimolecular reaction H2CN-->H+HCN on a fitted ab initio potential energy surface (PES). The ab initio energies were obtained from CCSD(T)/aug-cc-pvtz electronic structure calculations. The ab initio energies were fitted by the interpolating moving least-squares (IMLS) method. This work continues the development of the IMLS method for producing ab initio PESs for use in molecular dynamics simulations of many-atom systems. A dual-level scheme was used in which the preliminary selection of data points was done using a low-level theory and the points used for fitting the final PES were obtained at the desired higher level of theory. Classical trajectories were used on various low-level IMLS fits to tune the fit to the unimolecular reaction under study. Procedures for efficiently picking data points, selecting basis functions, and defining cutoff limits to exclude distant points were investigated. The accuracy of the fitted PES was assessed by comparing interpolated values of quantities to the corresponding ab initio values. With as little as 330 ab initio points classical trajectory rate constants were converged to 5%-10% and the rms error over the six-dimensional region sampled by the trajectories was a few tenths of a kcal/mol.     

Achieving fast convergence of ab initio free energy perturbation calculations with the adaptive force-matching method

This paper studies the possibility of improving the convergence of ab initio free energy perturbation (FEP) calculations by developing customized force fields with the adaptive force-matching (AFM) method. The ab initio FEP method relies on a molecular mechanics (MM) potential to sample configuration space. If the Boltzmann weight of the MM sampling is close to that of the ab initio method, the efficiency of ab initio FEP will be optimal. The difference in the Boltzmann weights can be quantified by the relative energy difference distribution (REDD). The force field developed through AFM significantly improves the REDD when compared with standard MM models, thus improving the convergence of the ab initio FEP calculation. The static dielectric constant ε_s of ice-Ih was studied with PW-91 through ab initio FEP. With a customized force field developed through AFM, we were able to converge ε_s to 80 ± 4 with 3,600 configurations. A similar ab initio FEP calculation with the TIP4P model would require 220 times more configurations to achieve the same accuracy. Our study indicates that the PW-91 functional underestimates ice-Ih ε_s by about 20%.     

Sensitivity of the H + CH3 → CH4 recombination rate constant to the shape of the CH stretching potential

Using a potential-energy surface obtained in part from ab initio calculations, the H + CH₃ → CH₄ bimolecular rate constant at T = 300 K is determined from a Monte Carlo classical trajectory study. Representing the CH stretching potential with a standard Morse function instead ofthe ab initio curve increases the calculated rate constant by an order of magnitude. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio stretching potentials.Two properties of the H + CH₃ α CH₄ potential-energy surface which significantly affect the recombination rate constant are the shape of the CH stretching potential and the attenuation of the H₃CH bending frequencies. Ab initio calculations with a hierarchy of basis sets and treatment of electron correlation indicate the latter is properly described [13]. The exact shape of the CH stretching potential is not delineated by the ab initio calculations, since the ab initio calculations are not converged for bond lengths of 2.0–3.0 Å [12]. However, the form of this stretching potential deduced from the highest-level ab initio calculations, and fit analytically by eq. (2), is significantly different from a Morse function. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio CH stretching potentials. This indicates that the actual CH potential energy curve lies between the Morse and ab initio curves. This is consistent with the finding that potential energy curves for diatomics are not well described by a Morse function [12].     

Ab initio centroid path integral molecular dynamics: application to vibrational dynamics of diatomic molecular systems

An ab initio centroid molecular dynamics (CMD) method is developed by combining the CMD method with the ab initio molecular orbital method. The ab initio CMD method is applied to vibrational dynamics of diatomic molecules, H2 and HF. For the H2 molecule, the temperature dependence of the peak frequency of the vibrational spectral density is investigated. The results are compared with those obtained by the ab initio classical molecular dynamics method and exact quantum mechanical treatment. It is shown that the vibrational frequency obtained from the ab initio CMD approaches the exact first excitation frequency as the temperature lowers. For the HF molecule, the position autocorrelation function is also analyzed in detail. The present CMD method is shown to well reproduce the exact quantum result for the information on the vibrational properties of the system.     

Molecular dynamics and metadynamics simulations of [2 + 2] photocycloaddition

Triplet state mechanism of [2 + 2] photocycloaddition forming a cyclobutane ring from two ethylenes is investigated in the context of photocatalysis. High‐level ab initio calculations are combined with ab initio adiabatic molecular dynamics and ab initio metadynamics for rare events modeling. In a photocatalytic scheme, a reactant reaches the triplet state either via intersystem crossing (ISC) or triplet sensitization. The model system adopts a biradical structure, which represents energy intersection with the ground state. The system either completes cyclization or undergoes fragmentation into two olefinic units. The potential and free energy surfaces of the cyclobutane/ethylenes system are mapped with multireference approaches describing possible reaction pathways. To obtain a full picture of a double bond photoreactivity, ab initio adiabatic dynamical calculations were used to estimate reaction yields and to model the effects of excess energy. The potential use of density functional theory based approaches for [2 + 2] photocycloaddition was investigated for future simulations and design of realistic photocatalytic systems. Dynamical aspects of [2 + 2] photocycloaddition via a triplet state manifold are investigated by combining ab initio multireference methods and ab initio molecular dynamics and metadynamics. The reaction pathways are studied for a model system of two ethylenes forming a cyclobutane ring to provide a basis for further studies on design of photocatalytic systems.     

Ab initio evaluation of bond orders and valence numbers in some substituted benzenes

Mayer's definition of bond order and valence number of an atom in a molecule in the ab initio SCF theory have been applied to calculate bond orders and valence numbers using different ab initio methods for some fluorobenzenes. The results have been discussed in the light of Mulliken population analysis. The traditional view that the valency of an atom is directly related to the atomic charge is found to be invalid in the present ab initio calculations.     

A hierarchical construction scheme for accurate potential energy surface generation: an application to the F+H2 reaction

We present a hierarchical construction scheme for accurate ab initio potential energy surface generation. The scheme is based on the observation that when molecular configuration changes, the variation in the potential energy difference between different ab initio methods is much smaller than the variation for potential energy itself. This means that it is easier to numerically represent energy difference to achieve a desired accuracy. Because the computational cost for ab initio calculations increases very rapidly with the accuracy, one can gain substantial saving in computational time by constructing a high accurate potential energy surface as a sum of a low accurate surface based on extensive ab initio data points and an energy difference surface for high and low accuracy ab initio methods based on much fewer data points. The new scheme was applied to construct an accurate ground potential energy surface for the FH(2) system using the coupled-cluster method and a very large basis set. The constructed potential energy surface is found to be more accurate on describing the resonance states in the FH(2) and FHD systems than the existing surfaces.     

Free energy perturbation study of water dimer dissociation kinetics

An efficient approach is described for using accurate ab initio calculations to determine the rates of elementary condensation and evaporation processes that lead to nucleation of aqueous aerosols. The feasibility of the method is demonstrated in an application to evaporation rates of water dimer at 230 K. The method, known as ABC-FEP (ab initio/classical free energy perturbation), begins with a calculation of the potential of mean force for the dissociation (evaporation) of small water clusters using a molecular dynamics (MD) simulation with a model potential. The free energy perturbation is used to calculate how changing from the model potential to a potential calculated from ab initio methods would alter the potential of mean force. The difference in free energy is the Boltzmann-weighted average of the difference between the ab initio and classical potential energies, with the average taken over a sample of configurations from the MD simulation. In principle, the method does not require a highly accurate model potential, though more accurate potentials require fewer configurations to achieve a small sampling error in the free energy perturbation step. To test the feasibility of obtaining accurate potentials of mean force from ab initio calculations at a modest number of configurations, the free energy perturbation method has been used to correct the errors when some standard models for bulk water (SPC, TIP4P, and TIP4PFQ) are applied to water dimer. To allow a thorough exploration of sampling issues, a highly accurate fit to results of accurate ab initio calculations, known as SAPT-5s, as been used a proxy for the ab initio calculations. It is shown that accurate values for a point on the potential of mean force can be obtained from any of the water models using ab initio calculations at only 50 configurations. Thus, this method allows accurate simulations of small clusters without the need to develop water models specifically for clusters.     

An ab initio molecular fragment investigation of the electronic and conformational structure of formaldoxime

An ab initio molecular fragment study of formaldoxime is reported and compared with previous semiempirical and conventional ab initio studies.     

Distribution of reaction products (theory). Investigation of an ab initio energy-surface for F + H2 ⇀ HF + H

Bender, O'Neil, Pearson and Schaefer (BOPS) have computed ab initio energies for 232 collinear configurations of FHH, using extensive configuration interaction. We have fitted these points using an LEPS-type function. Comparison with semi-empirical surfaces for FHH shows that the general form of these surfaces is in good accord with the ab initio findings. Evidence is presented which indicates that the BOPS ab initio surface exhibits too great a drop in energy along the favoured route into the exit valley.     

Comparing ab initio density-functional and wave function theories: the impact of correlation on the electronic density and the role of the correlation potential

The framework of ab initio density-functional theory (DFT) has been introduced as a way to provide a seamless connection between the Kohn-Sham (KS) formulation of DFT and wave-function based ab initio approaches [R. J. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005)]. Recently, an analysis of the impact of dynamical correlation effects on the density of the neon atom was presented [K. Jankowski, K. Nowakowski, I. Grabowski, and J. Wasilewski, J. Chem. Phys. 130, 164102 (2009)], contrasting the behaviour for a variety of standard density functionals with that of ab initio approaches based on second-order M?ller-Plesset (MP2) and coupled cluster theories at the singles-doubles (CCSD) and singles-doubles perturbative triples [CCSD(T)] levels. In the present work, we consider ab initio density functionals based on second-order many-body perturbation theory and coupled cluster perturbation theory in a similar manner, for a range of small atomic and molecular systems. For comparison, we also consider results obtained from MP2, CCSD, and CCSD(T) calculations. In addition to this density based analysis, we determine the KS correlation potentials corresponding to these densities and compare them with those obtained for a range of ab initio density functionals via the optimized effective potential method. The correlation energies, densities, and potentials calculated using ab initio DFT display a similar systematic behaviour to those derived from electronic densities calculated using ab initio wave function theories. In contrast, typical explicit density functionals for the correlation energy, such as VWN5 and LYP, do not show behaviour consistent with this picture of dynamical correlation, although they may provide some degree of correction for already erroneous explicitly density-dependent exchange-only functionals. The results presented here using orbital dependent ab initio density functionals show that they provide a treatment of exchange and correlation contributions within the KS framework that is more consistent with traditional ab initio wave function based methods.     

Modeling electrochemical interfaces from ab initio molecular dynamics: water adsorption on metal surfaces at potential of zero charge

The origin of the potential difference between the potential of zero charge of a metal/water interface and the work function of the metal is a recurring issue because it is related to how water interacts with metal surface in the absence of surface charge. Recently ab initio molecular dynamics method has been used to model electrochemical interfaces to study interfacial potential and the structure of interface water. Here, we will first introduce the computational standard hydrogen electrode method, which allows for ab initio determination of electrode potentials that can be directly compared with experiment. Then, we will review the recent progress from ab initio molecular dynamics simulation in understanding the interaction between water and metal and its impact on interfacial potential. Finally, we will give our perspective for future development of ab initio computational electrochemistry.     

Ab initio study of L-orientational defect in the hydrogen-bonded pattern of liquid water

Current status of the theory of orientational defects in H-bonded pattern of liquid water is briefly reviewed. Ab initio calculated water clusters from dimer to heptamer are thoroughly analyzed in terms of H-bonded pattern. New water heptamer structure of norbornane type is found via ab initio HF 6-311G** calculation. Its normal vibrations are properly assigned. Two ab initio water hexamers that refer to orientational defects are revealed. This is the first evidence of ab initio orientational defect in H-bonded patterns of water clusters. Some properties including normal vibrations of these defects, are studied.     

Modeling the global potential energy surface of the N + N2 reaction from ab initio data

A new global potential energy surface for the N + N2 exchange reaction has been built from ab initio data. To overcome the difficulty of carrying out ab initio calculations for a large set of geometries the alternative strategy of fitting the minimum energy paths of the surface and their angular dependence using a modified LAGROBO functional form has been adopted. In this way we have been able to reproduce all the main features of the potential using a fairly small set of ab initio values.     

On the accuracy of higher-order force constants calculated at the self-consistent field level of theory

Cubic and quartic internal coordinate force constants are often obtained by a least-squares fit of a set of ab initio energies to an expansion of the potential energy. The reliability of such a procedure can now be determined, since analytic Cartesian cubic force constants are available which, when used in a finite difference procedure, also yield ab initio Cartesian quartic force constants. It is shown that unless the points of the energy grid are situated close to the equilibrium point unreliable results are obtained. This is why ab initio SCF anharmonic constants are more accurate than previously believed.     

Multiple time scale Hartree–Fock molecular dynamics

This is the first application of a rigorous, established multiple time-step method to ab initio molecular dynamics. The resulting algorithm is conceptually simple and easy to implement, but very effective. It translates the large mass differences present in ab initio molecular dynamics into substantial savings in computer time while retaining high accuracy. This transforms ab initio molecular dynamics from a desirable but prohibitively expensive possibility into a viable method, at least for short-time phenomena in small systems or for otherwise inaccessibly complicated potential energy surfaces. © 1993 John Wiley & Sons, Inc.     

Ab initio STUDIES ON NITROXYL RADICAL MODELS OF ORGANIC FERROMAGNETS

The magnetic coupling between the unpaired spin of nitroxyl group (>N--O) and itsneighbor sp~2 carbon atom has been studied in detail by UHF/6--31G and MP4/6-31G methods.and a number of conjugate nitroxyl diradicals have also been studied by ab initio UHFmethod. A simple topological rule about the ferromagnetic coupling between two nitroxylgroups has been revealed by our calculation. Based on this rule as well as ab initio cal-culation, some nitroxyl radical models of organic ferromagnets (OFM) have been proposed,and the possible way of synthesis has also been suggested. Moreover, the ab initio UHF-crystal orbital (UHF--CO) method has been used to characterize the electronic band struc-ture of one of the models. The ab initio UHF--CO results show that this model is a verypromising candidate of OFM,     

分子间相互作用的研究——-价轨道模型势从头算能量分解法

研究分子间相互作用是了解液体,固体性质和结构以及气体性质的关键,也是研究化学和生物化学催化机制及化学反应途径的重要方面.因此,近几年来这个领域的理论和实验研究引起了人们广泛的兴趣并取得了长足的进步.特别值得指出的是Morokuma等提出的基于单行列式从头算的能量分解法,较好地解决了很多体系中分子间相互作用的本质问题,受到了理论化学界的普遍重视.但是这种方法很费计算机时间,对于较大的体系,特别是含有重原子的体系,应用受到了限制.因此简化这种方法,使它能够比较容     

Ab initio evidence for the stepwise mechanism of the McLafferty rearrangement of the butanal radical cation

The mechanism of the McLafferty rearrangement of the butanal radical cation to ethylene and vinyl alcohol cation is found, by ab initio calculations, to be stepwise. The results of a previous ab initio study are inconclusive because of symmetry restriction in their geometry optimization.