Hi,<br><br><i><span class="w">quote from a paper:</span><span class="st"> </span>A twin range cut-off for van der Waals (0.9/1.4 nm) and a smooth particle mesh Ewald algorithm for Coulomb interactions (switching distance of 0.9 nm) were used.</i><br>
<br>i.e.<span id="result_box" class="short_text" lang="en"><span class="hps">this means</span></span>:<br><br>; NEIGHBORSEARCHING PARAMETERS<br>.....         <br>rlist                     = 0.9<br><br>; OPTIONS FOR ELECTROSTATICS AND VDW<br>
coulombtype           = PME          <br>rcoulomb                 = 0.9    <br>vdw-type                 = Cut-off<br>rvdw                       = 1.4   <br><br>Am I wrong? <br>A twin range cut-off for van der Waals (0.9/1.4 nm) means rlist/rvdw. isnt it?<i><br>
</i><br>Thanks in advance<i><br><br></i><div class="gmail_quote">22 Mart 2012 01:32 tarihinde Oliver Stueker <span dir="ltr">&lt;<a href="mailto:ostueker@gmail.com">ostueker@gmail.com</a>&gt;</span> yazdı:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Ahmet,<br>
<div class="im"><br>
&gt; Let&#39;s imagine a sphere(two concentric spheres):<br>
&gt; radius of the inside small sphere=rvdw<br>
&gt; radius of the big sphere=rcoulomb<br>
&gt; distance between two of our nested spheres:rlist<br>
&gt; is this approach correct?<br>
<br>
</div>No. I suggest you read section 4.6.3 (and  probably also 4.6.2) in the<br>
Gromacs Manual.<br>
The Difference between rcoulomb (or rvdw) and rlist is a buffer-zone<br>
for the fact that the<br>
neighbor-lists are only updated every nstlist steps (often nstlist = 5).<br>
<br>
Oliver<br>
<br>
2012/3/21 ahmet yıldırım &lt;<a href="mailto:ahmedo047@gmail.com">ahmedo047@gmail.com</a>&gt;:<br>
<div class="HOEnZb"><div class="h5">&gt; Dear users,<br>
&gt;<br>
&gt; Berk Hess says:<br>
&gt; [gmx-users] a) rlist vs rvdw/rcoulomb size confusion, and b) reduced units<br>
&gt; Mon Jul 16 14:17:12 CEST 2007<br>
&gt;<br>
&gt; There are three options in Gromacs.<br>
&gt;<br>
&gt; The option you want is rcoulomb &lt; rlist and rvdw &lt; rlist.<br>
&gt; This works and gives the most accurate and also the most costly integration.<br>
&gt;<br>
&gt; A faster and very commonly used setting is: rlist=rcoulomb=rvdw.<br>
&gt; With PME the Coulomb interactions are very small at the cut-off,<br>
&gt; as are the LJ interactions.<br>
&gt; So with a small sacrifice in integration accuracy one can gain a lot<br>
&gt; of performance, also because analytical LJ is cheaper than tabulated.<br>
&gt;<br>
&gt; The last option is rcoulomb &gt; rlist and/or rvdw &gt; rlist.<br>
&gt; Then the energies and forces beyond rlist are only updated every nslist<br>
&gt; steps.<br>
&gt; This gives less integration accuracy but can give a lot of &quot;interaction&quot;<br>
&gt; accuracy<br>
&gt; at a small computational cost.<br>
&gt;<br>
&gt; Justin A. Lemkul says:<br>
&gt; [gmx-users] Twin-range cut-off<br>
&gt; Tue Sep 13 01:25:28 CEST 2011<br>
&gt;<br>
&gt;<br>
&gt; A twin-range cutoff just means that your short-range cutoffs aren&#39;t all the<br>
&gt; same<br>
&gt; value, such that they form two interaction zones.  Within the shortest, the<br>
&gt; neighbor list is updated every step.  Between the shortest and longest<br>
&gt; cutoffs,<br>
&gt; the neighbor list is updated every nstlist steps.  For instance:<br>
&gt;<br>
&gt; rlist = 0.9<br>
&gt; rcoulomb = 0.9<br>
&gt; rvdw = 1.4<br>
&gt;<br>
&gt; are common settings for Gromos96 force fields (in conjunction with PME).<br>
&gt; Thus<br>
&gt; there are two interaction zones - the first is if two atoms (or charge<br>
&gt; groups,<br>
&gt; depending on the algorithm) are within 0.9 nm, and the second is if the two<br>
&gt; interacting species are beyond 0.9 nm but within 1.4 nm of each other.<br>
&gt;<br>
&gt; rcoulomb: distance for Coulomb cut-off (nm)<br>
&gt; rvdw:distance for LJ or Buckingham cut-off (nm)<br>
&gt; nstlist: neighbor list update frequency<br>
&gt; rlist: cut-off distance of the short-range neighbor<br>
&gt; Twin range cutoff consists of rcoulomb and rvdw, isnt it?<br>
&gt;<br>
&gt; Let&#39;s imagine a sphere(two concentric spheres):<br>
&gt; radius of the inside small sphere=rvdw<br>
&gt; radius of the big sphere=rcoulomb<br>
&gt; distance between two of our nested spheres:rlist<br>
&gt; is this approach correct?<br>
&gt; I could not understand the fourierspacing and rlist.<br>
&gt;<br>
&gt;<br>
&gt; Thanks in advance<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; 21 Mart 2012 20:53 tarihinde ahmet yıldırım &lt;<a href="mailto:ahmedo047@gmail.com">ahmedo047@gmail.com</a>&gt; yazdı:<br>
&gt;<br>
&gt;&gt; Dear users,<br>
&gt;&gt;<br>
&gt;&gt; I have two configuration as the following related to Neigborsearching,<br>
&gt;&gt; Electrostatics and vdw options. I checked the literature:<br>
&gt;&gt; Generally the rlist, rcoulomb and rvdw have used as the following.<br>
&gt;&gt; rlist=1<br>
&gt;&gt; rcoulomb=0.8<br>
&gt;&gt; rvdw=1.4<br>
&gt;&gt;<br>
&gt;&gt; Is there much difference between the following two options in the<br>
&gt;&gt; calculation/the results? Is there one significant difference between the two<br>
&gt;&gt; options. If yes, then what is it? What is relationship between rlist,<br>
&gt;&gt; nstlist and rvdw/rcoulomb?<br>
&gt;&gt;<br>
&gt;&gt; Furthermore,<br>
&gt;&gt; fourierspacing         = 0.16<br>
&gt;&gt; or<br>
&gt;&gt; fourierspacing         = 0.12<br>
&gt;&gt; difference between these two options?<br>
&gt;&gt;<br>
&gt;&gt; 1.choice<br>
&gt;&gt; .....<br>
&gt;&gt; ; NEIGHBORSEARCHING PARAMETERS<br>
&gt;&gt; nstlist                 = 5<br>
&gt;&gt; ns-type               = Grid<br>
&gt;&gt; pbc                    = xyz<br>
&gt;&gt; rlist                    = 1.0<br>
&gt;&gt;<br>
&gt;&gt; ; OPTIONS FOR ELECTROSTATICS AND VDW<br>
&gt;&gt; coulombtype              = PME<br>
&gt;&gt; pme_order                = 4<br>
&gt;&gt; fourierspacing          = 0.16<br>
&gt;&gt; rcoulomb                 = 1.0<br>
&gt;&gt; vdw-type                 = Cut-off<br>
&gt;&gt; rvdw                       = 1.0<br>
&gt;&gt; ...<br>
&gt;&gt;<br>
&gt;&gt; 2.choice<br>
&gt;&gt; ..<br>
&gt;&gt; ; NEIGHBORSEARCHING PARAMETERS<br>
&gt;&gt; nstlist                  = 5<br>
&gt;&gt; ns-type                = Grid<br>
&gt;&gt; pbc                      = xyz<br>
&gt;&gt; rlist                     = 0.9<br>
&gt;&gt;<br>
&gt;&gt; ; OPTIONS FOR ELECTROSTATICS AND VDW<br>
&gt;&gt; coulombtype           = PME<br>
&gt;&gt; pme_order              = 4<br>
&gt;&gt; fourierspacing         = 0.16<br>
&gt;&gt; rcoulomb                 = 0.9<br>
&gt;&gt; vdw-type                 = Cut-off<br>
&gt;&gt; rvdw                       = 1.4<br>
&gt;&gt; ...<br>
&gt;&gt;<br>
&gt;&gt; Thanks in advance<br>
&gt;&gt; --<br>
&gt;&gt; Ahmet Yıldırım<br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; Ahmet Yıldırım<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don&#39;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>
Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br>Ahmet Yıldırım<br><br>