<div dir="ltr">Hi,<div><br></div><div>Thanks for the report. A lot of this stuff was thrown in the air during some recent bug fixing before 4.6.6, so possibly something got broken. It would be interesting to know what you observe with 4.6.6, and with earlier 4.6.x also.</div>
<div><br></div><div>Mark</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jul 23, 2014 at 5:36 AM, Igor Leontyev <span dir="ltr"><<a href="mailto:ileontyev@ucdavis.edu" target="_blank">ileontyev@ucdavis.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear gmx-developers,<br>
Testing with gmx-5.0 different schemes of Long-Range electrostatics for NMA liquid I noticed unreasonable behavior with "coulombtype=Switch". Electrostatic energies calculated for the same configuration (i.e. at 0 MD step) with PME and Cut-off coulombtype reasonably agree between each other. But once I use coulombtype=Switch, the electrostatic energy jumps by the order of magnitude which is suspicious. See the results bellow. Any idea what can be wrong ?<br>
<br>
Thank you,<br>
Igor<br>
<br>
<br>
<br>
<coulombtype=PME, rcoulomb=1.3><br>
Coulomb-14 LJ (SR) Disper. corr. Coulomb(SR)+Recipr. Potential<br>
-1.29147e+04 -3.28852e+03 -1.28735e+02 -3.54258e+03 -7.86555e+03<br>
<br>
<coulombtype=Cut-off, rcoulomb=1.3><br>
Coulomb-14 LJ (SR) Disper. corr. Coulomb (SR) Potential<br>
-1.29147e+04 -3.28852e+03 -1.28735e+02 -3.53416e+03 -7.85714e+03<br>
<br>
<coulombtype=Switch, rcoulomb=1.3, rcoulomb-switch=1.2><br>
Coulomb-14 LJ (SR) Disper. corr. Coulomb (SR) Potential<br>
-1.29147e+04 -3.28852e+03 -1.28735e+02 -5.25858e+04 -5.69088e+04<br>
<br>
==============================<u></u>==============================<u></u>=<br>
Long-Range block in mdp-file ==============================<u></u>==============================<u></u>= <br>
nstlist = 5 ; update pairlist<br>
rlist = 1.3 ; cut-off for ns<br>
cutoff-scheme = group;<br>
coulombtype = Switch/PME/Cut-off ;<br>
rcoulomb = 1.3 ; cut-off for coulomb<br>
rcoulomb_switch = 1.2 ; switch for coulomb<br>
vdwtype = Cut-off<br>
rvdw = 1.3 ; cut-off for vdw<br>
rvdw_switch = 1.0 ; where to start switching vdw<br>
dispcorr = EnerPres ; LR dispersion corr for En and Pres<br>
;=============================<u></u>==============================<u></u>============<span class="HOEnZb"><font color="#888888"><br>
-- <br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/GMX-<u></u>developers_List</a> before posting!<br>
<br>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" target="_blank">https://maillist.sys.kth.se/<u></u>mailman/listinfo/gromacs.org_<u></u>gmx-developers</a> or send a mail to <a href="mailto:gmx-developers-request@gromacs.org" target="_blank">gmx-developers-request@<u></u>gromacs.org</a>.<br>
</font></span></blockquote></div><br></div>