Dear users,<br><br><b>Berk Hess says:<br>[gmx-users] a) rlist vs rvdw/rcoulomb size confusion, and b) reduced units<br>Mon Jul 16 14:17:12 CEST 2007</b><br><pre>There are three options in Gromacs.
The option you want is rcoulomb < rlist and rvdw < rlist.
This works and gives the most accurate and also the most costly integration.
A faster and very commonly used setting is: rlist=rcoulomb=rvdw.
With PME the Coulomb interactions are very small at the cut-off,
as are the LJ interactions.
So with a small sacrifice in integration accuracy one can gain a lot
of performance, also because analytical LJ is cheaper than tabulated.
The last option is rcoulomb > rlist and/or rvdw > rlist.
Then the energies and forces beyond rlist are only updated every nslist
steps.
This gives less integration accuracy but can give a lot of "interaction"
accuracy
at a small computational cost.<br><br></pre><b>Justin A. Lemkul says:<br>[gmx-users] Twin-range cut-off<br>Tue Sep 13 01:25:28 CEST 2011</b><br><pre><br>A twin-range cutoff just means that your short-range cutoffs aren't all the same
value, such that they form two interaction zones. Within the shortest, the
neighbor list is updated every step. Between the shortest and longest cutoffs,
the neighbor list is updated every nstlist steps. For instance:
rlist = 0.9
rcoulomb = 0.9
rvdw = 1.4
are common settings for Gromos96 force fields (in conjunction with PME). Thus
there are two interaction zones - the first is if two atoms (or charge groups,
depending on the algorithm) are within 0.9 nm, and the second is if the two
interacting species are beyond 0.9 nm but within 1.4 nm of each other.<br></pre>rcoulomb: distance for Coulomb cut-off (nm)<br>rvdw:distance for LJ or Buckingham cut-off (nm)<br>nstlist: neighbor list update frequency<br>
rlist: cut-off distance of the short-range neighbor<br>Twin range cutoff consists of rcoulomb and rvdw, isnt it?<br><br><span id="result_box" class="short_text" lang="en"><span class="hps">Let's imagine</span> <span class="hps">a sphere</span></span>(two concentric spheres)<span id="result_box" class="short_text" lang="en"><span class="hps">:<br>
</span></span><span id="result_box" class="short_text" lang="en"><span class="hps">radius of the </span></span><span id="result_box" class="short_text" lang="en"><span class="hps">inside</span> small </span><span id="result_box" class="short_text" lang="en"><span class="hps">sphere</span></span>=rvdw<br>
<span id="result_box" class="short_text" lang="en"><span class="hps">radius of the</span></span><span id="result_box" class="short_text" lang="en"> big </span><span id="result_box" class="short_text" lang="en"><span class="hps">sphere</span></span>=rcoulomb<br>
distance between two of our nested spheres:rlist<br>is this approach correct?<br> <div id="gt-res-content" class="almost_half_cell"><div dir="ltr" style="zoom:1"><span id="result_box" class="short_text" lang="en"><span class="hps">I could not understand</span> <span class="hps">the fourierspacing and</span> <span class="hps">rlist.</span></span></div>
</div><br><br>Thanks in advance<br><br><br><br><div class="gmail_quote">21 Mart 2012 20:53 tarihinde ahmet yıldırım <span dir="ltr"><<a href="mailto:ahmedo047@gmail.com">ahmedo047@gmail.com</a>></span> yazdı:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear users,<br><br>I have two configuration as the following related to Neigborsearching, Electrostatics and vdw options. I checked the literature<span lang="en"><span>:<br>Generally the rlist, rcoulomb and rvdw have used as the following.<br>
rlist=1<br>rcoulomb=0.8<br>rvdw=1.4<br></span></span><br><b>Is there much difference between the following two options in the calculation/the results? Is there one significant difference between the two options. If yes, then what is it? What is relationship between rlist, nstlist and rvdw/rcoulomb?</b><br>
<br>Furthermore,<br>fourierspacing = 0.16 <br>or<br>fourierspacing = 0.12<br><b>difference between these two options?</b><br><br><u><b>1.choice</b></u><br>.....<br>; NEIGHBORSEARCHING PARAMETERS<br>nstlist = 5 <br>
ns-type = Grid <br>pbc = xyz <br><b>rlist = 1.0 </b> <br><br>; OPTIONS FOR ELECTROSTATICS AND VDW<br>coulombtype = PME <br>
pme_order = 4 <br>fourierspacing = 0.16 <br><b>rcoulomb = 1.0</b><br>vdw-type = Cut-off<br><b>rvdw = 1.0</b><br>...<br><br><u><b>2.choice</b></u><br>
..<br>; NEIGHBORSEARCHING PARAMETERS<br>nstlist = 5 <br>ns-type = Grid <br>pbc = xyz <br><b>rlist = 0.9 </b> <br>
<br>; OPTIONS FOR ELECTROSTATICS AND VDW<br>coulombtype = PME <br>pme_order = 4 <br>fourierspacing = 0.16 <br><b>rcoulomb = 0.9 </b> <br>vdw-type = Cut-off <br>
<b>rvdw = 1.4 </b><br clear="all">...<br><br>Thanks in advance<span class="HOEnZb"><font color="#888888"><br>-- <br>Ahmet Yıldırım<br><br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br>Ahmet Yıldırım<br><br>