Loading…

Force Field Parametrization from the Hirshfeld Molecular Electronic Density

The Hirshfeld charges are linearly increased to reproduce the experimental dielectric constant of 10 polar solvents having values between 13 (pyridine) and 182 (N-methylformamide). The OPLS/AA force field is used to obtain the new parameters. The surface tension and liquid density are also target pr...

Full description

Saved in:
Bibliographic Details
Published in:Journal of chemical theory and computation 2018-11, Vol.14 (11), p.5949-5958
Main Authors: Pérez de la Luz, Alexander, Aguilar-Pineda, Jorge Alberto, Méndez-Bermúdez, José Guillermo, Alejandre, José
Format: Article
Language:English
Subjects:
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The Hirshfeld charges are linearly increased to reproduce the experimental dielectric constant of 10 polar solvents having values between 13 (pyridine) and 182 (N-methylformamide). The OPLS/AA force field is used to obtain the new parameters. The surface tension and liquid density are also target properties to determine the new nonbonding parameters. The charge scaling factor is between 1.2 and 1.3. In addition, properties that were not used in the parametrization procedure, such as the heat of vaporization, self-diffusion coefficient, shear viscosity, isothermal compressibility, and volumetric expansion coefficient are obtained. Binary mixtures of amide/water and amide/amide are also studied. The original parameters of OPLS/AA, CGenFF, and GAFF force fields are evaluated. The TIP4P/ε force field is used to simulate water. The results from this work with the new parameters, for both pure components and binary mixtures, are in better agreement with experimental data than those obtained with the original values for most of the calculated properties. The maximum density of N-methylformamide in aqueous solutions is correctly predicted only with the new parameters. The high value of the dielectric constant of acetamide, formamide, and N-methylformamide is discussed in terms of the chain formation from the hydrogen bond interactions.
ISSN:1549-9618
1549-9626
DOI:10.1021/acs.jctc.8b00554