Una,<br>
<br>
Sorry. I was in a hurry and made a typo. With PME, you don't have to
have rlist>rcoulomb. PME, of course, is only for coulomb
interactions, so it has no bearing on how big rvdw has to be. <br>
<br>
It turns out that you have to have rvdw>rlist with vdwtype=cutoff,
as GROMACS is telling you. But I'm not sure that cutoff vdw
interactions is really the best bet, as that means the potential
changes discontinuously. We tend to prefer switched vdw, which has the
additional side-effect that you don't have to have rvdw>rlist.<br>
<br>
David<br><br><div><span class="gmail_quote">On 1/27/06, <b class="gmail_sendername">Una Bjarnadottir</b> <<a href="mailto:una.bjarnadottir@ucd.ie">una.bjarnadottir@ucd.ie</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear all,<br><br>Thank you David for your reply,<br><br>Now when I try to use rvdw = 1.0 and rlist 1.1 and rcoulomb = 1.0 I get<br>the error; With vdwtype = Cut-off,rvdw must be >= rlist.<br>I don't have any rvdw-switch or rcoulomb-switch
<br>I'm using version 3.2.1 and I thought it had been fixed according to the<br>discussion on the list.<br><br>Is there some trich that one has to do to have rlist>rvdw?<br><br>Best wishes to everyone, Una<br><br>> With PME, you don't have to have rlist<rvdw. You might, however, want to
<br>> have rlist slightly larger than rvdw. You don't have to, but since the<br>> neighbor list is only updated every nstlist steps, if rlist is the same as<br>> rvdw, you can run into the situation where an atom has moved inside the vdw
<br>> cutoff in a couple of steps, but your calculation doesn't realize it has<br>> moved inside the cutoff until after nstlist steps. You can solve the<br>> problem, of course, by using nstlist=1, but that could slow things down a
<br>> lot, so it may be better to simply use, say, rlist=1.1 and rvdw=1.0 so atoms<br>> appear in the neighbor list before getting counted for vdw interactions.<br>><br>> David<br>><br>><br>> On 1/24/06, Una Bjarnadottir <
<a href="mailto:una.bjarnadottir@ucd.ie">una.bjarnadottir@ucd.ie</a>> wrote:<br>><br>>>Dear all,<br>>><br>>>Just to have the content of this discussion clear<br>>><a href="http://www.gromacs.org/pipermail/gmx-users/2004-August/011683.html">
http://www.gromacs.org/pipermail/gmx-users/2004-August/011683.html</a><br>>>and all other discussion on this matter.<br>>><br>>>When using coulombtype = PME and vdw-type = Cut-off is it alright to<br>>>have the;
<br>>><br>>>rlist=rvdw=rcoulomb = 1<br>>>rvdw-switch=rcoulomb-swithc = 0<br>>><br>>>Or does it have to be;<br>>><br>>>rlist=rcoulmb = 0.9 and rvdw = 1.4<br>>><br>>>Is the 2nd thing only for cut-off but not for PME?
<br>>><br>>>Thanks in advance for any replies, this is confusing me a lot.<br>>><br>>>Una<br><br><br>_______________________________________________<br>gmx-users mailing list<br><a href="mailto:gmx-users@gromacs.org">
gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to
<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br><br><br></blockquote></div><br>