Ab initio based polarizable force field parametrization

Experimental and simulation studies of anion-water systems have pointed out the importance of molecular polarization for many phenomena ranging from hydrogen-bond dynamics to water interfaces structure. The study of such systems at molecular level is usually made with classical molecular dynamics simulations. Structural and dynamical features are deeply influenced by molecular and ionic polarizability, which parametrization in classical force field has been an object of long-standing efforts. Although when classical models are compared to ab initio calculations at condensed phase, it is found that the water dipole moments are underestimated by approximately 30%, while the anion shows an overpolarization at short distances. A model for chloride-water polarizable interaction is parametrized here, making use of Car-Parrinello simulations at condensed phase. The results hint to an innovative approach in polarizable force fields development, based on ab initio simulations, which do not suffer for the mentioned drawbacks. The method is general and can be applied to the modeling of different systems ranging from biomolecular to solid state simulations.     

Ab initio molecular dynamics of liquid 1,3-dimethylimidazolium chloride

Density-functional-based Car-Parrinello molecular dynamics (CPMD) simulations have been performed for the ionic liquid 1,3-dimethylimidazolium chloride, [dmim]Cl, at 438 K. The local structure of the liquid is described in terms of various partial radial distribution functions and anisotropic spatial distributions, which reveal a significant extent of hydrogen bonding. The cation-anion distribution simulated with the BP86 functional is in qualitative agreement with the structural model derived from neutron diffraction data for the liquid, whereas the theoretical cation-cation distribution shows less satisfactory accord. Population analyses indicate noticeable charge transfer from anions to cations, and specific CH...Cl hydrogen bonds are characterized in terms of donor-acceptor interactions between lone pairs on Cl and antibonding sigma(CH) orbitals.     

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.     

Developing ab initio quality force fields from condensed phase quantum-mechanics/molecular-mechanics calculations through the adaptive force matching method

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.     

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.     

Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching

A novel least-squares fitting approach is presented to obtain classical force fields from trajectory and force databases produced by ab initio (e.g., Car-Parrinello) molecular dynamics (MD) simulations. The method was applied to derive effective nonpolarizable three-site force fields for liquid water at ambient conditions from Car-Parrinello MD simulations in the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure includes a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The various parameterizations of the water force field differ by an enforced smooth cut-off applied to the short-ranged interaction term. These were obtained by fitting to the trajectory and force data produced by Car-Parrinello MD simulations of systems of 32 and 64 H(2)O molecules. The new water force fields were developed assuming both flexible or rigid molecular geometry. The simulated structural and self-diffusion properties of liquid water using the fitted force fields are in close agreement with those observed in the underlying Car-Parrinello MD simulations. The resulting empirical models compare to experiment much better than many conventional simple point charge (SPC) models. The fitted potential is also shown to combine well with more sophisticated intramolecular potentials. Importantly, the computational cost of the new models is comparable to that for SPC-like potentials.     

Derivation of class II force fields. VIII. Derivation of a general quantum mechanical force field for organic compounds

A class II valence force field covering a broad range of organic molecules has been derived employing ab initio quantum mechanical "observables." The procedure includes selecting representative molecules and molecular structures, and systematically sampling their energy surfaces as described by energies and energy first and second derivatives with respect to molecular deformations. In this article the procedure for fitting the force field parameters to these energies and energy derivatives is briefly reviewed. The application of the methodology to the derivation of a class II quantum mechanical force field (QMFF) for 32 organic functional groups is then described. A training set of 400 molecules spanning the 32 functional groups was used to parameterize the force field. The molecular families comprising the functional groups and, within each family, the torsional angles used to sample different conformers, are described. The number of stationary points (equilibria and transition states) for these molecules is given for each functional group. This set contains 1324 stationary structures, with 718 minimum energy structures and 606 transition states. The quality of the fit to the quantum data is gauged based on the deviations between the ab initio and force field energies and energy derivatives. The accuracy with which the QMFF reproduces the ab initio molecular bond lengths, bond angles, torsional angles, vibrational frequencies, and conformational energies is then given for each functional group. Consistently good accuracy is found for these computed properties for the various types of molecules. This demonstrates that the methodology is broadly applicable for the derivation of force field parameters across widely differing types of molecular structures. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1782-1800, 2001     

Application of Static Charge Transfer within an Ionic‐Liquid Force Field and Its Effect on Structure and Dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3‐dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are ±0.5 e. While the three‐dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly—with the positive result that for ionic charges of ±0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial‐charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic‐liquid models in which the ionic charges are assumed to be ±1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled‐charge sets may be of benefit in the simulation of these systems—especially when looking at properties beyond liquid structure—thus providing an alternative to computationally expensive polarisable force fields.     

Molecular force field for ionic liquids III: imidazolium, pyridinium, and phosphonium cations; chloride, bromide, and dicyanamide anions

This is the third set of parameters of a force field for the molecular simulation of ionic liquids, developed within the spirit of the OPLS-AA model and thus oriented toward the calculation of equilibrium thermodynamic and structural properties. The parameter sets reported here concern the cations alkylimidazolium, tetra-alkylphosphonium, and N-alkylpyridinium, and the anions chloride, bromide, and dicyanamide. The force field is built in a stepwise manner that allows the construction of models for an entire family of cations, with alkyl side chains of different length, for example. Due to the transferability of the present force field, the ions studied here can be combined with those reported in our two previous publications to create a large variety of ionic liquids that can be studied by molecular simulation. The parameters reported were obtained through different series of ab initio calculations concerning the geometry, force constants, torsion energy profiles, and electrostatic charge distributions of the ions under study. Validation of the force field consisted of comparison with experimental crystal structure and liquid density data.     

An ab initio force field for the cofactors of bacterial photosynthesis

This article presents a new ab initio force field for the cofactors of bacterial photosynthesis, namely quinones and bacteriochlorophylls. The parameters has been designed to be suitable for molecular dynamics simulations of photosynthetic proteins by being compatible with the AMBER force field. To our knowledge, this is the first force field for photosynthetic cofactors based on a reliable set of ab initio density functional reference data for methyl bacteriochlorophyll a, methyl bacteriopheophytin a, and of a derivative of ubiquinone. Indeed, the new molecular mechanics force field is able to reproduce very well not only the experimental and ab initio structural properties and the vibrational spectra of the molecules, but also the eigenvectors of the molecular normal modes. For this reason it might also be helpful to understand vibrational spectroscopy results obtained on reaction center proteins.     

Using spectroscopic data on imidazolium cation conformations to test a molecular force field for ionic liquids

A molecular force field for the computer simulation of ionic liquids is evaluated a posteriori by confrontation against Raman spectroscopic data, published after the force field had been formulated. Specifically, the terms in the force field describing the conformational aspects of dialkylimidazolium cations, which were specifically developed for these compounds using high level ab initio calculations, are those affecting the distribution of conformers in simulated ionic liquids. Those distributions are compared with analyses of the liquid-phase Raman spectra, and the features of a series of dihedral torsions along the alkyl side chains in 1-alkyl-3-methylimidazolium cations in several ionic liquids are discussed.     

New force field for molecular simulation of guanidinium-based ionic liquids

An all-atom force field was proposed for a new class of room temperature ionic liquids (RTILs), N,N,N',N'-tetramethylguanidinium (TMG) RTILs. The model is based on the AMBER force field with modifications on several parameters. The refinements include (1) fitting the vibration frequencies for obtaining force coefficients of bonds and angles against the data obtained by ab initio calculations and/or by experiments and (2) fitting the torsion energy profiles of dihedral angles for obtaining torsion parameters against the data obtained by ab initio calculations. To validate the force field, molecular dynamics (MD) simulations at different temperatures were performed for five kinds of RTILs, where TMG acts as a cation and formate, lactate, perchlorate, trifluoroacetate, and trifluoromethylsulfonate act as anions. The predicted densities were in good agreement with the experimental data. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) were investigated to depict the microscopic structures of the RTILs.     

Ab initio parametrized polarizable force field for rutile-type SnO<Subscript>2</Subscript>

We report a new, polarizable classical force field for the rutile-type phase of SnO₂, casserite. This force field has been parametrized using results from ab initio (density functional theory) calculations as a basis for fitting. The force field was found to provide structural, dynamical and thermodynamic properties of tin oxide that compare well with both ab initio and experimental results at ambient and high pressures.     

Toward accurate solvation dynamics of lanthanides and actinides in water using polarizable force fields: from gas-phase energetics to hydration free energies

In this contribution, we focused on the use of polarizable force fields to model the structural, energetic, and thermodynamical properties of lanthanides and actinides in water. In a first part, we chose the particular case of the Th(IV) cation to demonstrate the capabilities of the AMOEBA polarizable force field to reproduce both reference ab initio gas-phase energetics and experimental data including coordination numbers and radial distribution functions. Using such model, we predicted the first polarizable force field estimate of Th(IV) solvation free energy, which accounts for −1,638 kcal/mol. In addition, we proposed in a second part of this work a full extension of the SIBFA (Sum of Interaction Between Fragments Ab initio computed) polarizable potential to lanthanides (La(III) and Lu(III)) and to actinides (Th(IV)) in water. We demonstrate its capabilities to reproduce all ab initio contributions as extracted from energy decomposition analysis computations, including many-body charge transfer and discussed its applicability to extended molecular dynamics and its parametrization on high-level post-Hartree–Fock data.     

Nonadiabatic electron wavepacket dynamics of molecules in an intense optical field: an ab initio electronic state study

A theory of quantum electron wavepacket dynamics that nonadiabatically couples with classical nuclear motions in intense optical fields is studied. The formalism is intended to track the laser-driven electron wavepackets in terms of the linear combination of configuration-state functions generated with ab initio molecular orbitals. Beginning with the total quantum Hamiltonian for electrons and nuclei in the vector potential of classical electromagnetic field, we reduce the Hamiltonian into a mixed quantum-classical representation by replacing the quantum nuclear momentum operators with the classical counterparts. This framework gives equations of motion for electron wavepackets in an intense laser field through the time dependent variational principle. On the other hand, a generalization of the Newtonian equations provides a matrix form of forces acting on the nuclei for nonadiabatic dynamics. A mean-field approximation to the force matrix reduces this higher order formalism to the semiclassical Ehrenfest theory in intense optical fields. To bring these theories into a practical quantum chemical package for general molecules, we have implemented the relevant ab initio algorithms in it. Some numerical results in the level of the semiclassical Ehrenfest-type theory with explicit use of the nuclear kinematic (derivative) coupling and the velocity form for the optical interaction are presented.     

Derivation of force field parameters for SnO2-H2O surface systems from plane-wave density functional theory calculations

Plane-wave density functional theory (DFT-PW) calculations were performed on bulk SnO2 (cassiterite) and the (100), (110), (001), and (101) surfaces with and without H2O present. A classical interatomic force field has been developed to describe bulk SnO2 and SnO2-H2O surface interactions. Periodic density functional theory calculations using the program VASP (Kresse et al., 1996) and molecular cluster calculations using Gaussian 03 (Frisch et al., 2003) were used to derive the parametrization of the force field. The program GULP (Gale, 1997) was used to optimize parameters to reproduce experimental and ab initio results. The experimental crystal structure and elastic constants of SnO2 are reproduced reasonably well with the force field. Furthermore, surface atom relaxations and structures of adsorbed H2O molecules agree well between the ab initio and force field predictions. H2O addition above that required to form a monolayer results in consistent structures between the DFT-PW and classical force field results as well.