Search

A new method called adaptive force matching (AFM) has been developed that is capable of producing high quality force fields for condensed phase simulations. This procedure involves the parametrization of force fields to reproduce ab initio forces obtained from condensed phase quantum-mechanics/molecular-mechanics (QM/MM) calculations. During the procedure, the MM part of the QM/MM is iteratively improved so as to approach ab initio quality. In this work, the AFM method has been tested to parametrize force fields for liquid water so that the resulting force fields reproduce forces calculated using the ab initio MP2 and the Kohn–Sham density functional theory with the Becke–Lee–Yang–Parr (BLYP) and Becke three-parameter LYP (B3LYP) exchange correlation functionals. The AFM force fields generated in this work are very simple to evaluate and are supported by most molecular dynamics (MD) codes. At the same time, the quality of the forces predicted by the AFM force fields rivals that of very expensive ab initio calculations and are found to successfully reproduce many experimental properties. The site-site radial distribution functions (RDFs) obtained from MD simulations using the force field generated from the BLYP functional through AFM compare favorably with the previously published RDFs from Car–Parrinello MD simulations with the same functional. Technical aspects of AFM such as the optimal QM cluster size, optimal basis set, and optimal QM method to be used with the AFM procedure are discussed in this paper. [ABSTRACT FROM AUTHOR]

In this paper the structures of 4-aminophenol(H2O)1+ and 3-aminophenol(H2O)1+ clusters are investigated in molecular beam experiments by different IR/UV-double resonance techniques as well as the mass analyzed threshold ionization spectroscopy yielding both inter- and intramolecular vibrations of the ionic and neutral species. Possible structures are extensively calculated at the level of density functional theory (DFT) or at the ab initio level of theory. From the experimental and theoretical investigations it can be concluded that in the case of 4-aminophenol(H2O)1 one O[Single_Bond]H...O hydrogen-bonded structure exists in the neutral cluster but two structures containing either an O[Single_Bond]H...O or a N[Single_Bond]H...O hydrogen-bonded arrangement are observed in the spectra of the ionic species. This observation is a result of an intramolecular rearrangement reaction within the ion which can only take place if high excess energies are used. A reaction path via the CH bonds is calculated and explains the experimental observations. In the case of 3-aminophenol(H2O)1+ only one O[Single_Bond]H...O bound structure is observed both in the neutral and ionic species. Ab initio and DFT calculations show that due to geometrical and energetical reasons a rearrangement cannot be observed in the 3-aminophenol(H2O)1+ cluster ion. [ABSTRACT FROM AUTHOR]

Plane wave basis sets offer many advantages in ab initio molecular dynamics due to their efficiency and simplicity. In combination with hybrid density functionals, they become computationally expensive due to the evaluation of the Hartree-Fock exchange energy. The computational cost can be significantly reduced by screening the Kohn-Sham orbital products after localizing the orbitals in real space. However, such a procedure introduces apparent errors in the wavefunctions and nuclear forces resulting in unstable dynamics. It is shown here that a noise-stabilized dynamics approach can overcome this problem and at the same time permits using insufficiently converged wavefunctions for evaluating atomic forces. In this way, we achieve significant speed up even for a small system containing about 100 atoms. After benchmarking the accuracy and efficiency of this approach, we use it in combination with well-sliced metadynamics to compute the free energy barrier of formamide hydrolysis in alkaline aqueous medium. These results provide insight into the error of the Perdew-Burke-Ernzerhof functional in predicting the free energy barrier for hydrolysis reactions in water. [ABSTRACT FROM AUTHOR]

Ab initio molecular dynamics simulations with hybrid density functionals have so far found little application due to their computational cost. In this work, an implementation of the Hartree–Fock exchange is presented that is specifically targeted at ab initio molecular dynamics simulations of medium sized systems. We demonstrate that our implementation, which is available as part of the CP2K/Quickstep program, is robust and efficient. Several prescreening techniques lead to a linear scaling cost for integral evaluation and storage. Integral compression techniques allow for in-core calculations on systems containing several thousand basis functions. The massively parallel implementation respects integral symmetry and scales up to hundreds of CPUs using a dynamic load balancing scheme. A time-reversible multiple time step scheme, exploiting the difference in computational efficiency between hybrid and local functionals, brings further time savings. With extensive simulations of liquid water, we demonstrate the ability to perform, for several tens of picoseconds, ab initio molecular dynamics based on hybrid functionals of systems in the condensed phase containing a few thousand Gaussian basis functions. [ABSTRACT FROM AUTHOR]

We examined reactivity of H2 on Ru(0001) using molecular beam techniques and we compared our results to experimental results for similar systems. The dissociative adsorption of H2 on Ru(0001) is similar to that on Pt(111) and Ni(111), although on ruthenium nonactivated adsorption is strongly suggested. However, we find no clear signature of a steering- or precursor-based mechanism that favors nonactivated reaction paths at low kinetic energy. In comparison to Pd(111) and Rh(111) our results indicate that a universal mechanism enhancing reactivity at low energy does not have a mass dependence. In addition, we have compared our results to predictions of reactivity for H2 on Ru(0001) from six-dimensional dynamical calculations using two different generalized gradient approximation functionals. It leads us to conclude that the PW91 functional yields a more accurate value for the minimum energy path but does not impose enough corrugation in the potential. The revised-Perdew-Burke-Ernzerhof (RPBE) functional appears to behave slightly better at higher energies, but we find significant quantitative disagreement. We show that the difference is not due to different energy resolutions between experiment and theory. However, it may be due to a dependence of the reactivity on rotational state or on omission of relevant dimensions in the theoretical description. [ABSTRACT FROM AUTHOR]

We present a density functional theory study of interactions between spherical colloidal particles in amphiphile solutions. Theory is found to be in good agreement with previously published molecular dynamics simulations. It is used to analyze the effect of the amphiphile solution bulk density, the chain length, and the solvent mole fraction on the potential of mean force between the particles. The general features of the potential of mean force are rationalized in terms of formation of layers and bilayers of amphiphilic molecules in the intercolloidal gap. Theory yields the same general trends as observed in simulations and in experiments. In particular, the computed mean force changes its character from repulsive to attractive and back to repulsive as the solvent mole fraction is gradually increased. [ABSTRACT FROM AUTHOR]

The authors use ab initio molecular dynamics and the density functional theory+U (DFT+U) method to compute the hydration environment of the manganese ion in manganese (II) and manganese (III) porphines (MnP) dispersed in liquid water. These are intended as simple models for more complex water soluble porphyrins, which have important physiological and electrochemical applications. The manganese ion in Mn(II)P exhibits significant out-of-porphine plane displacement and binds strongly to a single H2O molecule in liquid water. The Mn in Mn(III)P is on average coplanar with the porphine plane and forms a stable complex with two H2O molecules. The residence times of these water molecules exceed 15 ps. The DFT+U method correctly predicts that water displaces NO from Mn(III)P–NO, but yields an ambiguous spin state for the MnP(II)–NO complex. [ABSTRACT FROM AUTHOR]

We report an ab initio-based global potential energy surface for H+CH4 that describes the abstraction and exchange reactions. The PES, which is invariant with respect to any permutation of five H atoms, is a fit to 20 728 electronic energies calculated using the partially spin-restricted coupled-cluster method (RCCSD(T)) with a moderately large basis (aug-cc-pVTZ). A first set of quasiclassical trajectory calculations using this PES are reported for the H+CD4→HD+CD3 reaction at collision energies of 0.65 and 1.52 eV and are compared to experiment and recent direct dynamics calculations done with density functional theory. [ABSTRACT FROM AUTHOR]

The physicochemical properties of two molten salts, namely, KCl and NaCl, have been studied with a molecular-dynamics approach using a density-functional-based tight-binding (DFTB) model. The obtained results have been compared with a number of previously reported simulations, carried out on smaller systems and using classical force-field techniques. A good agreement has been found for both structural parameters and macroscopic properties, such as self-diffusion coefficients. Furthermore, our DFTB results are very close to the available experimental data. From a more general point of view, our results demonstrate the applicability of DFTB as an efficient tool in the modeling of melts. At the same time, the quality of the obtained results supports the use of this as a reliable alternative to the more expensive ab initio dynamics approaches, if accurate parameters are provided. [ABSTRACT FROM AUTHOR]

We extend our previous formulation of time-dependent four-component relativistic density-functional theory [J. Gao, W. Liu, B. Song, and C. Liu, J. Chem. Phys. 121, 6658 (2004)] by using a noncollinear form for the exchange-correlation kernel. The new formalism can deal with excited states involving moment (spin)-flipped configurations which are otherwise not accessible with ordinary exchange-correlation functionals. As a first application, the global potential-energy curves of 16 low-lying ωω-coupled electronic states of the AuH molecule have been investigated. The derived spectroscopic parameters, including the adiabatic and vertical excitation energies, equilibrium bond lengths, harmonic and anharmonic vibrational constants, fundamental frequencies, and dissociation energies, are grossly in good agreement with those of ab initio multireference second-order perturbation theory and the available experimental data. [ABSTRACT FROM AUTHOR]

A (quasirelativistic) two-component density functional theory (DFT) approach to the computation of parity-violating energy differences between enantiomers is presented which is based on the zeroth-order regular approximation (ZORA). This approach is employed herein to compute parity-violating energy differences between several P and M conformations of dihydrogen dichalcogenides (H2X2 with X=O, S, Se, Te, Po), of which some compounds have recently been suggested as potential molecular candidates for the first experimental measurement of parity-violating effects in chiral molecules. The DFT ZORA results obtained in this work with “pure” density functionals are anticipated to deviate by well less than 1% from data that would be computed within related (relativistic) four-component Dirac–Kohn–Sham–Coulomb schemes. In our implementation of the ZORA slightly larger relative deviations are expected for hybrid functionals, depending on the amount of “exact” exchange. For B3LYP (20% exact exchange) differences are estimated to amount to at most 3% in hydrogen peroxide, 2% in disulfane, and 1% or less for the heavier homologs. Thus, the present two-component approach is expected to perform excellently when compared to four-component density functional schemes while being at the same time computationally more efficient. The ZORA approach will therefore be of particular interest for the prediction of parity-violating vibrational frequency shifts, for instance, in isotopomers of H2Se2 and H2Te2. [ABSTRACT FROM AUTHOR]

We propose approximations which go beyond the local-density approximation for the short-range exchange and correlation density functionals appearing in a multideterminantal extension of the Kohn–Sham scheme. A first approximation consists of defining locally the range of the interaction in the correlation functional. Another approximation, more conventional, is based on a gradient expansion of the short-range exchange-correlation functional. Finally, we also test a short-range generalized-gradient approximation by extending the Perdew–Burke–Ernzerhof exchange-correlation functional to short-range interactions. © 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The performance of density functional theory methods for the modeling of condensed aqueous systems is hard to predict and validation by ab initio molecular simulation of liquid water is absolutely necessary. In order to assess the reliability of these tests, the effect of temperature on the structure and dynamics of liquid water has been characterized with 16 simulations of 20 ps in the temperature range of 280–380 K. We find a pronounced influence of temperature on the pair correlation functions and on the diffusion constant including nonergodic behavior on the time scale of the simulation in the lower temperature range (which includes ambient temperature). These observations were taken into account in a consistent comparison of a series of density functionals (BLYP, PBE, TPSS, OLYP, HCTH120, HCTH407). All simulations were carried out using an ab initio molecular dynamics approach in which wave functions are represented using Gaussians and the density is expanded in an auxiliary basis of plane waves. Whereas the first three functionals show similar behavior, it is found that the latter three functionals yield more diffusive dynamics and less structure.© 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A new Kohn–Sham formalism is developed for studying the lowest molecular electronic states of given space and spin symmetry whose densities are represented by weighted sums of several reference configurations. Unlike standard spin-density functional theory, the new formalism uses total spin conserving spin-density operators and spin-invariant density matrices so that the method is fully spin-adapted and solves the so-called spin-symmetry dilemma. The formalism permits the use of an arbitrary set of reference (noninteracting) configurations with any number of open shells. It is shown that the requirement of degeneracy of the total noninteracting energies of the reference configurations (or configuration state functions) is equivalent to the stationary condition of the exact energy relative to the weights of the configurations (or configuration state functions). Consequently, at any molecular geometry, the weights can be determined by minimization of the energy, and, for given reference weights, the Kohn–Sham orbitals can be determined. From this viewpoint, the developed theory can be interpreted as an analog of the multiconfiguration self-consistent field approach within density functional theory. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]