[gmx-users] calculating dielectric constant

[David van der Spoel]

Qiao Baofu wrote:
> Hi all,
> 
> I have some question about calculating dielectric constant:
> 1. From the manual, it seems that there are two methods to calculate the 
> dielectric constant: a. use g_dipoles -temp 400 -f -s directly, b) use 
> g_dipoles -nonormalize -corr total  -c dipcorr.xvg, then g_dielectric -f 
> dipcorr.xvg -o epsw.xvg. (Please correct me if the commandlines are 
> wrong). Thus I want to know the difference between them.

g_dielectric does another analyis. You want to run g_dipoles either with 
the -enx flag or with a trajectory.

> 2. If I want to calculate the dielectric constant, must I use the 
> Reaction Field in the couloumbtype? Now I have some results from PME. 
> Can I use these results directly to calculate the dielectric constant?

PME is fine. You have to supply the right dielectric constant to the 
program (though the default of 0 == infinity) is what most people use.

> 
> Thanks very much in advances!
> 
> 
> 
> Sincerely yours,
> Baofu Qiao
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> 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.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
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
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++