Communication: Modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions
- Department of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115 (United States)
- Computation Institute, The University of Chicago, Chicago, Illinois 60637 (United States)
We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley “bracelet” and “rod” test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, “Charge asymmetries in hydration of polar solutes,” J. Phys. Chem. B 112, 2405–2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.
- OSTI ID:
- 22436530
- Journal Information:
- Journal of Chemical Physics, Vol. 141, Issue 13; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Instantaneous, parameter-free methods to define a solute’s hydration shell
First Principles Simulations fo the Supercritical Behavior of Ore Forming Fluids