Search

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. Four rigid water models have been parameterized; the matching procedure has been used to assess the role of the ghost atom in TIP4P-like models and of electrostatic damping. In the former case, burying the negative charge inside the molecule allows to fit better the torques. In the latter, since short-range interactions are damped, a better fit of the forces is obtained. Overall, the best performing model is the one with a ghost atom and with electrostatic damping. The approach shown in this paper is of general validity and could be applied to any matching algorithm and to any level of coarse graining, also for non-rigid molecules. [ABSTRACT FROM AUTHOR]

Molecular 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. In addition, we also examined the behavior of solutions trapped within hydrophilic cavities, in which a set of unsaturated O-sites at the wall were transformed in polar silanol Si-OH groups. In all cases, the overall concentrations of the trapped electrolytes exhibited important reductions that, in the case of the narrowest pores, attained 50% of the bulk value. Local concentrations within the pores also showed important fluctuations. In hydrophobic cavities, the close vicinity of the pore wall was coated exclusively by the solvent, whereas in hydrophilic pores, selective adsorption of Na+ ions was also observed. Mass and charge transport were also investigated. Individual diffusion coefficients did not present large modifications from what is perceived in the bulk; contrasting, the electrical conductivity exhibited important reductions. The qualitative differences are rationalized in terms of simple geometrical considerations. [ABSTRACT FROM AUTHOR]

Recently, 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. We also show how to reduce the computational overhead for evaluating the dipoles, introducing to the state-of-the-art methods: the extended Lagrangian method and the always stable predictor corrector method. Finally we discuss the importance of screening the electrostatic interactions at short range, defending this technique against simpler approximations usually made. We compare results of density functional theory and classical force field-based Molecular Dynamics simulations of chloride in water. [ABSTRACT FROM AUTHOR]