[gmx-users] g_dipoles gives wrong epsilon for SPC water
yesint4 at yahoo.com
Wed Dec 24 21:08:26 CET 2008
I've found that the dielectric constant of pure SPC water according to g_dipoles is only 3.94, which is obviously wrong!
I simulated the box of ~1500 SPC waters at 325 K and computed epsilon with
g_dipoles -temp 325
That's what I get:
Finite system Kirkwood g factor G_k = 0.190603
Infinite system Kirkwood g factor g_k = 0.143163
Epsilon = 3.94756
After some googling I've found that eps=65 for SPC model. So, what is wrong? Is there any other way of computing epsilon?
----- Original Message ----
From: David van der Spoel <spoel at xray.bmc.uu.se>
To: Discussion list for GROMACS users <gmx-users at gromacs.org>
Sent: Wednesday, December 24, 2008 6:01:38 PM
Subject: Re: [gmx-users] g_dipoles questions
Semen Esilevsky wrote:
> Dear David,
> Thank you for reply! However there are still some question:
>>> Dear All,
>>> I have few questions about g_dipoles, which are not explained in the manual:
>>> 1) What method is used to calculate the epsilon? Is there a reference which can be cited? Is it suitable for water solutions?
>> In principle it is suitable for solutions of neutral molecules only.
> So, it should work for SPC water for example?
> I'm confused by the phrase in the man page "you can compute e.g. the dielectric constant for low dielectric media". What is "low"? Does it mean that the method works for, say, benzene but not for water?
anything of neutral molecules.
>> The references are to Neuman papers from the 1980s.
> Could you please give a little bit more detailed citation? Neuman is a widespread name and 1980s is 10 years... At least a title or some keywords are needed to find these papers.
check my 1998 JCP paper, the ref is given in g_dipoles.
> gmx-users mailing list gmx-users at gromacs.org
> 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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
David van der Spoel, PhD, Professor of Biology
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596, 75124 Uppsala, Sweden
phone: 46 18 471 4205 fax: 46 18 511 755
spoel at xray.bmc.uu.se spoel at gromacs.org http://folding.bmc.uu.se
gmx-users mailing list gmx-users at gromacs.org
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-request at gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
More information about the gromacs.org_gmx-users