<br><br>----- Original Message -----<br>From: Gaurav Goel &lt;gauravgoeluta@gmail.com&gt;<br>Date: Thursday, June 17, 2010 5:21<br>Subject: [gmx-users] how small rcoulomb can be<br>To: Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<br><br>&gt; For my NPT/NVT simulations I'm using a BKS potential for dispersion<br>&gt; interactions (VdWtype=USER)&nbsp; which has a cut-off of 0.55 <br>&gt; nm. I also<br>&gt; set rcoulomb=rlist=0.55 nm.<br>&gt; <br>&gt; My first question is regarding the choice of rcoulomb=0.55nm. Any<br>&gt; thoughts on whether this rcoulomb is large enough to ascertain that<br>&gt; real space electrostatics interactions are fairly small (in comparison<br>&gt; to total coulombic interactions)?<br><br>That's around half the size of typical rcoulomb values used for PME, so you'll need to adjust ewald_rtol and fourier_spacing to a considerable degree in order that the electrostatics your PME model is comparably accurate. If this model has been used with PME in the past, then you will have some guidance in its literature (though GROMACS does not make it easy to choose beta directly). This will mean the reciprocal-space processors are doing lots more work than normal, and you'll want to set -npme to be half of -np, or maybe zero. Even so, your parallel scaling will be noticeably worse than normal PME.<br>&nbsp;<br>Mark<br><br>&gt; My second question is in regards to doing long (~15-100 ns) NVE<br>&gt; simulations. It appears using coulombtype=PME-Switch will be the best<br>&gt; option to accurately conserve energy. Any thoughts on how to choose<br>&gt; rcoulomb_switch and rcoulomb?<br>&gt; <br>&gt; Thanks,<br>&gt; Gaurav<br>&gt; -- <br>&gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search <br>&gt; before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php