scholarly journals Prediction of experimental properties of CO2: improving actual force fields

2019 ◽  
Vol 25 (6) ◽  
Author(s):  
Raúl Fuentes-Azcatl ◽  
Hector Domínguez
Keyword(s):  
Author(s):  
John W. Coleman

In the design engineering of high performance electromagnetic lenses, the direct conversion of electron optical design data into drawings for reliable hardware is oftentimes difficult, especially in terms of how to mount parts to each other, how to tolerance dimensions, and how to specify finishes. An answer to this is in the use of magnetostatic analytics, corresponding to boundary conditions for the optical design. With such models, the magnetostatic force on a test pole along the axis may be examined, and in this way one may obtain priority listings for holding dimensions, relieving stresses, etc..The development of magnetostatic models most easily proceeds from the derivation of scalar potentials of separate geometric elements. These potentials can then be conbined at will because of the superposition characteristic of conservative force fields.


2018 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Alejandro Lara ◽  
David L. Mobley ◽  
Toon Vestraelen ◽  
Adelio R Matamala ◽  
...  

<div>Computer simulations of bio-molecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in bio-molecular systems and are therein described by atomic point charges.</div><div>In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute's electron density computed with an implicit solvent model and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the MBIS atomic charge method, including the solvent polarization, with a root mean square error of 2.0 kcal mol<sup>-1</sup> for the 613 organic molecules studied. The largest deviation was observed for phosphor-containing molecules and the molecules with amide, ester and amine functional groups.</div>


2019 ◽  
Author(s):  
Pier Paolo Poier ◽  
Louis Lagardere ◽  
Jean-Philip Piquemal ◽  
Frank Jensen

<div> <div> <div> <p>We extend the framework for polarizable force fields to include the case where the electrostatic multipoles are not determined by a variational minimization of the electrostatic energy. Such models formally require that the polarization response is calculated for all possible geometrical perturbations in order to obtain the energy gradient required for performing molecular dynamics simulations. </p><div> <div> <div> <p>By making use of a Lagrange formalism, however, this computational demanding task can be re- placed by solving a single equation similar to that for determining the electrostatic variables themselves. Using the recently proposed bond capacity model that describes molecular polarization at the charge-only level, we show that the energy gradient for non-variational energy models with periodic boundary conditions can be calculated with a computational effort similar to that for variational polarization models. The possibility of separating the equation for calculating the electrostatic variables from the energy expression depending on these variables without a large computational penalty provides flexibility in the design of new force fields. </p><div><div><div> </div> </div> </div> <p> </p><div> <div> <div> <p>variables themselves. Using the recently proposed bond capacity model that describes molecular polarization at the charge-only level, we show that the energy gradient for non-variational energy models with periodic boundary conditions can be calculated with a computational effort similar to that for variational polarization models. The possibility of separating the equation for calculating the electrostatic variables from the energy expression depending on these variables without a large computational penalty provides flexibility in the design of new force fields. </p> </div> </div> </div> </div> </div> </div> </div> </div> </div>


Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole

<div><div><div><p>Modern molecular mechanics force fields are widely used for modelling the dynamics and interactions of small organic molecules using libraries of transferable force field parameters. For molecules outside the training set, parameters may be missing or inaccurate, and in these cases, it may be preferable to derive molecule-specific parameters. Here we present an intuitive parameter derivation toolkit, QUBEKit (QUantum mechanical BEspoke Kit), which enables the automated generation of system-specific small molecule force field parameters directly from quantum mechanics. QUBEKit is written in python and combines the latest QM parameter derivation methodologies with a novel method for deriving the positions and charges of off-center virtual sites. As a proof of concept, we have re-derived a complete set of parameters for 109 small organic molecules, and assessed the accuracy by comparing computed liquid properties with experiment. QUBEKit gives highly competitive results when compared to standard transferable force fields, with mean unsigned errors of 0.024 g/cm3, 0.79 kcal/mol and 1.17 kcal/mol for the liquid density, heat of vaporization and free energy of hydration respectively. This indicates that the derived parameters are suitable for molecular modelling applications, including computer-aided drug design.</p></div></div></div>


Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole

<div><div><div><p>Modern molecular mechanics force fields are widely used for modelling the dynamics and interactions of small organic molecules using libraries of transferable force field parameters. For molecules outside the training set, parameters may be missing or inaccurate, and in these cases, it may be preferable to derive molecule-specific parameters. Here we present an intuitive parameter derivation toolkit, QUBEKit (QUantum mechanical BEspoke Kit), which enables the automated generation of system-specific small molecule force field parameters directly from quantum mechanics. QUBEKit is written in python and combines the latest QM parameter derivation methodologies with a novel method for deriving the positions and charges of off-center virtual sites. As a proof of concept, we have re-derived a complete set of parameters for 109 small organic molecules, and assessed the accuracy by comparing computed liquid properties with experiment. QUBEKit gives highly competitive results when compared to standard transferable force fields, with mean unsigned errors of 0.024 g/cm3, 0.79 kcal/mol and 1.17 kcal/mol for the liquid density, heat of vaporization and free energy of hydration respectively. This indicates that the derived parameters are suitable for molecular modelling applications, including computer-aided drug design.</p></div></div></div>


2019 ◽  
Author(s):  
Kateryna Goloviznina ◽  
José N. Canongia Lopes ◽  
Margarida Costa Gomes ◽  
Agilio Padua

A general, transferable polarisable force field for molecular simulation of ionic liquids and their mixtures with molecular compounds is developed. This polarisable model is derived from the widely used CL\&P fixed-charge force field that describes most families of ionic liquids, in a form compatible with OPLS-AA, one of the major force fields for organic compounds. Models for ionic liquids with fixed, integer ionic charges lead to pathologically slow dynamics, a problem that is corrected when polarisation effects are included explicitly. In the model proposed here, Drude induced dipoles are used with parameters determined from atomic polarisabilities. The CL\&P force field is modified upon inclusion of the Drude dipoles, to avoid double-counting of polarisation effects. This modification is based on first-principles calculations of the dispersion and induction contributions to the van der Waals interactions, using symmetry-adapted perturbation theory (SAPT) for a set of dimers composed of positive, negative and neutral fragments representative of a wide variety of ionic liquids. The fragment approach provides transferability, allowing the representation of a multitude of cation and anion families, including different functional groups, without need to re-parametrise. Because SAPT calculations are expensive an alternative predictive scheme was devised, requiring only molecular properties with a clear physical meaning, namely dipole moments and atomic polarisabilities. The new polarisable force field, CL\&Pol, describes a broad set set of ionic liquids and their mixtures with molecular compounds, and is validated by comparisons with experimental data on density, ion diffusion coefficients and viscosity. The approaches proposed here can also be applied to the conversion of other fixed-charged force fields into polarisable versions.<br>


Sign in / Sign up

Export Citation Format

Share Document