[gmx-users] question about dielectric constant calculation

Lei Zhou zhoumadison at gmail.com
Thu Mar 30 22:48:09 CEST 2006

Dear GMX-users,

I am trying to calculate the dieletric constant based on MD trajectories.
Two programs are used, g_dipoles and g_dielectric. I have the following

1) In g_dipoles, what number should I use for the flag  -epsilonRF? I used
PME duing MD simulation with the epsilon-r as 1 and epsilon-surface as 0.
2) What other flags are critical for calculation, like -P,  Order of
Legendre polynomial for ACF?
3) Based on the manual, the epsilon calculated by g_dielectric is
frenquency-related. How about the epsilon output from g_dipoles and the
relationship to difference frenquencies?
4) I used the output from g_dipoles, Mtot.xvg, as an input for g_dielectric.
Here is the command used: g_dielectric -f Mtot.xvg, however, I got the
following error messages after frame 49: Fatal error: nparm = 0 in file
expfit.c, line 462.

I am really grateful for any help with those questions. Thank you.

Lei Zhou
Columbia University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20060330/66366906/attachment.html>

More information about the gromacs.org_gmx-users mailing list