English
 
Help Privacy Policy Disclaimer
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

Is the solution activity derivative sufficient to parametrize ion–ion interactions? Ions for tip5p water

MPS-Authors
/persons/resource/persons204558

Satarifard,  Vahid
Andrea Grafmüller, Theorie & Bio-Systeme, Max Planck Institute of Colloids and Interfaces, Max Planck Society;

/persons/resource/persons186287

Kashefolgheta,  Sadra
Ana Vila Verde, Theorie & Bio-Systeme, Max Planck Institute of Colloids and Interfaces, Max Planck Society;

/persons/resource/persons125918

Verde,  Ana Vila
Ana Vila Verde, Theorie & Bio-Systeme, Max Planck Institute of Colloids and Interfaces, Max Planck Society;

/persons/resource/persons121351

Grafmüller,  Andrea
Andrea Grafmüller, Theorie & Bio-Systeme, Max Planck Institute of Colloids and Interfaces, Max Planck Society;

External Resource
No external resources are shared
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available
Citation

Satarifard, V., Kashefolgheta, S., Verde, A. V., & Grafmüller, A. (2017). Is the solution activity derivative sufficient to parametrize ion–ion interactions? Ions for tip5p water. Journal of Chemical Theory and Computation, 13(5), 2112-2122. doi:10.1021/acs.jctc.6b01229.


Cite as: http://hdl.handle.net/11858/00-001M-0000-002D-3D37-2
Abstract
Biomolecular processes involve hydrated ions, and thus molecular simulations of such processes require accurate force-field parameters for these ions. In the best force-fields, both ion-water and anion-cation interactions are explicitly parameterized. First, the ion Lennard-Jones parameters are optimized to reproduce, e.g., single ion solvation free energies; then ion-pair interactions are often adjusted to match experimental activity or activity derivatives. Here we apply this approach to derive optimized parameters for concentrated NaCl, KCl, MgCl2 and CaCl2 salt solutions, to be used with the TIP5P water model. These parameters are of interest because of a number of desirable properties of the TIP5P water model, especially for the simulation of carbohydrates. The results show, that this approach is insufficient, because the activity derivative often reaches a plateau near the target experimental value, for a wide range of parameter values. The plateau emerges from the interconversion between different types of ion pairs, so parameters leading to equally good agreement with the target solution activity or activity derivative yield very different solution structures. To resolve this indetermination, a second target property, such as the experimentally determined ion-ion coordination number, is required to uniquely determine anion-cation interactions. Simulations show that combining activity derivatives and coordination number as experimental target properties to parameterize ion-ion interactions, is a powerful method for reliable ion-water force field parameterization, and gives insight into the concentration of contact or solvent shared ion pairs in a wide range of salt concentrations. For the alkali and halide ions Li+, Rb+, Cs+, F-, Br-, I-, we present ion-water parameters appropriate at infinite dilution only.