David:
Thank you. This really resolves a lot of our "head-scratching" on TIP4P model effort.
-keith
-----Original Message-----
From: owner-amber.scripps.edu [mailto:owner-amber.scripps.edu] On Behalf Of David A. Case
Sent: Saturday, November 29, 2008 11:21 AM
To: amber.scripps.edu
Subject: Re: AMBER: Questions on TIP4P model?
On Fri, Nov 28, 2008, Chin, Keith B wrote:
>
> 1. In amber amber9/dat/solvents/tip4p dir, there is a parm.top file of
> tip4p water. The atomic mass of the extra point, EPW, is 3.0 instead of
> 0.0. For what purpose was the EPW point assigned a value of 3.0?
As the 00_README files indicate, these are pretty old files. The mass of
extra points is ignored; I suspect the 3.0 value comes from some earlier times
(these files were made by Jim Caldwell in 2002) when extra points were handled
in a different fashion.
> 3. Most importantly. In xleap, I created sand parmtop and parmcrd file
> using frcmod.tip4pew and frcmod.tip model files.
> For tip4pew model the partial charges in the parmtop file are:
> O H1 H2 EPW
> 0.00000000E+00 9.55249411E+00 9.55249411E+00 -1.91049882E+01
>
> For tip4p model the partial charges in the parmtop file are:
> O H1 H2 EPW
> 0.00000000E+00 9.47559600E+00 9.47559600E+00 -1.89511920E+01
>
As Carlos explained, charges in a prmtop file are in internal units. If you
divide the EPW charge value of 18.951192 by 18.2223, you get -1.04, in units
of the electron charge.
TIP4P/EW is recommended over TIP4P itself for most Amber simulations. Be sure
you really want TIP4P itself if you do simulations with it.
....dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Fri Dec 05 2008 - 18:12:45 PST
