Thanks a lot, Berk. This is exactly the answer I was hoping for, 
meaning it's clear for me now.

Jeroen


> Date: Thu, 12 Apr 2007 13:42:42 +0200
> From: "Berk Hess" <[EMAIL PROTECTED]>
> Subject: RE: [gmx-users] soft-core potential in combination with PME
>       (sorry, again)
> To: [EMAIL PROTECTED]
> Message-ID: <[EMAIL PROTECTED]>
> Content-Type: text/plain; format=flowed
> 
> I don't know about Anwars paper.
> But I can say that their is no fundamental difference between PME
> and RF with soft-core. With the GROMACS setup they converge
> to exactly the same energies at any lambda value when you increase the RF 
> cut-off.
> 
> I don't understand what David Mobley meant exactly.
> There is no Coulomb singularity with soft-core.
> Maybe one could have an unfortunate situation where the LJ
> is already very soft, but the Coulomb not very soft,
> which could lead to instabilities.
> But I have never encountered this.
> 
> Berk.
> 
> 
> >From: "Jeroen van Bemmelen" <[EMAIL PROTECTED]>
> >Reply-To: Discussion list for GROMACS users <[EMAIL PROTECTED]>
> >To: [EMAIL PROTECTED]
> >Subject: RE: [gmx-users] soft-core potential in combination with PME 
> >(sorry,again)
> >Date: Thu, 12 Apr 2007 12:38:56 +0200
> >
> >Hi Berk (and others),
> >
> >Thanks a lot for your quick answer. And sorry to bother you again
> >with this long email, but I still feel I'm missing something really
> >essential here.
> >
> >I agree that in essence the convergence and the singularity are two
> >different problems. And if I understand you correctly, my first
> >conclusion (that it's completely fine to use soft-core together with
> >PME, that it takes care of the singularity problem, both for the VDW
> >and the electrostatic part at the same time, and that it does not
> >lead to additional problems) was right. So numerical instabilities
> >(simulation crash) during free energy calculations resulting from
> >singularities (overlapping atoms) should not be an issue with ANY
> >type of electrostatics method in GROMACS (even PME), just as long as
> >you use the soft-core potential.
> >
> >But in that case I really don't understand why Anwar claims in his
> >aforementioned paper that the electrostatics singularity problem is
> >only solved by the standard soft-core potential for the simple
> >truncation and reaction field schemes, and not for PME. And why David
> >Mobley mentions numerical instabilities caused by "charge-charge
> >overlap" (e.g. in http://www.gromacs.org/pipermail/gmx-users/2006-
> >March/020785.html and other mails) while using soft-core. Are the
> >mentioned numerical instabilities (simulation blowing up) resulting
> >from charge overlap not a direct consequence of the singularity at
> >r=0? And is the soft-core potential not supposed to take care of
> >exactly THIS problem, by effectively limiting the forces as r-->0? At
> >least that's how I understand it from the Beutler paper and from the
> >manual.
> >
> >So I guess somebody is wrong here, most likely I myself, and I hope
> >you or someone else can clarify this a bit further.
> >
> >Thanks,
> >Jeroen
_______________________________________________
gmx-users mailing list    [EMAIL PROTECTED]
http://www.gromacs.org/mailman/listinfo/gmx-users
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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to