Eric Shamay wrote:
In researching the shell model used to create "polarizable" ions I came across several of the references of the drude oscillator approach and the papers of van der Spoel (for the SW water model) and Roux (for the parameters of the halides and alkalis). However, perhaps due to my own oversight, I did not walk away from those publications having a firm understanding of just how one goes about developing the parameters necessary to implement the shell model, even for simple monovalent ions (i.e. sodium).

I am attempting to simulate an aqueous sodium chloride system using a polarizable (shell) model for the individual sodium and chloride ions, but eventually I will also need data for nitrate and sulfate ions. Is there a more clear reference that illustrates the process (or lists the results) of developing the parameters (force constant [k], charge of the shell particle [q_s], etc.) for using the shell model? I'd like to be able to come up with my own values, but a published set will work as well. As a side note, I'm coming from having used the amber point-polarizable model that assigns an isotropic alpha to each particle. I can calculate the polarizability tensor using DFT methods, but still don't know how to translate that into a force constant or shell-particle charge.

The equations are in the SW paper. However in case of a singly charged atom like Na+ you have one extra degree of freedom. In case of molecular ions you can use e.g. RESP calculations to get the partial charges and use these to determine the spring constant.

Note that you will also need to update Van der Waals parameters.


thank you,
~Eric Shamay


------------------------------------------------------------------------

_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


--
David van der Spoel, Ph.D., Professor of Biology
Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
sp...@xray.bmc.uu.se    sp...@gromacs.org   http://folding.bmc.uu.se
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to