Dear gromacs users,<br><br>In section 4.6.3 of gromacs 4 manual it says ``... the interactions between pairs that do not fall within rlist but do fall within max(rcoulomb,rvdw) are computed during neighbor search (NS), and the forces and energy are stored separately, and added to short-range forces at every time step between successive NS.&#39;&#39;<br>
<br>Now I find that if I use rlist (= rvdw) = rcoulomb then I only have Coulomb (SR) values in the md.log file but if I use rlist (= rvdw) &lt; rcoulomb then I get both Coulomb (SR) and Coulomb (LR) values in the md.log. In both cases I have assigned vdwtype=user and also rcoulomb=user for a tabulated potential that I use.  <br>
<br>Now, my question is how in the second case when rlist &lt; rcoulomb I am getting the Coulomb (LR) values in md.log? Particularly, I want to know how the Coulomb (LR) is calculated? Is it calculated in the same way as Coulomb (SR) and simply denoted as  Coulomb (LR), since rlist &lt; rcoulomb? Or, in this case it is suggested that I should always use rlist=rcoulomb.<br>
<br>Any suggestion is appreciated and thanks in advance.<br><br>-Swarnendu<br>