In the quest towards coarse-grained potentials and new water models, we present an extension of the force matching technique to parameterize an all-atom force field for rigid water. The methodology presented here allows to improve the matching procedure by first optimizing the weighting exponents present in the objective function. A new gauge for unambiguously evaluating the quality of the fit has been introduced; it is based on the root mean square difference of the distributions of target properties between reference data and fitted potentials.
View Article and Find Full Text PDFWe report non-equilibrium molecular dynamics simulations (NEMD) of water under temperature gradients using a modified version of the central force model (MCFM). This model is very accurate in predicting the equation of state of water for a wide range of pressures and temperatures. We investigate the polarization response of water to thermal gradients, an effect that has been recently predicted using Non-Equilibrium Thermodynamics (NET) theory and computer simulations, as a function of the thermal gradient strength.
View Article and Find Full Text PDFMolecular dynamics simulations have been carried out to investigate structural and dynamical characteristics of NaCl aqueous solutions confined within silica nanopores in contact with a "bulk-like" reservoir. Two types of pores, with diameters intermediate between 20 Å and 37.5 Å, were investigated: The first one corresponded to hydrophobic cavities, in which the prevailing wall-solution interactions were of the Lennard-Jones type.
View Article and Find Full Text PDFRecently, the use of polarizable force fields in Molecular Dynamics simulations has been gaining importance, since they allow a better description of heterogeneous systems compared to simple point charges force fields. Among the various techniques developed in the last years the one based on polarizable point dipoles represents one of the most used. In this paper, we review the basic technical issues of the method, illustrating the way to implement intramolecular and intermolecular damping of the electrostatic interactions, either with and without the Ewald summation method.
View Article and Find Full Text PDF