Dear Prof. Spoel, <br>  <br>  Many thanks for your reply.<br>  <br>  Since the fourier_n* values are associated with the grid sizing, will  these values affect my results drastically? do you suggest resubmitting  these runs?<br>  <br>  regards,<br>  Pri...<br>  <br><br><b><i>David van der Spoel &lt;spoel@xray.bmc.uu.se&gt;</i></b> wrote:<blockquote class="replbq" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;">  priyanka srivastava wrote:<br>&gt; Hie,<br>&gt; <br>&gt; Many thanks for your reply.<br>&gt; <br>&gt; Yes I have used semi-isotropic coupling.<br>&gt; <br>&gt; Another doubt is why only a value of 4.5e-5?<br>&gt; <br>&gt; I have used PME and have mentioned rcoulomb in the<br>&gt; .mdp file too, which as is pointed out by you, has no<br>&gt; meaning!! But does this mean that my mdp file is<br>&gt; wrong? Because if PME does not use rcoulomb it shud<br>&gt; simply ignore the value.<br>&gt; Actually I do not want to use twin range
 cutoff. So, I<br>&gt; have given rvdw = rlist. Is this alright?<br>&gt; Portion of my mdp file looks like this:<br>&gt; <br>&gt; ; Electrostatics<br>&gt; coulombtype             =  PME<br>&gt; fourier_nx              =  6.4<br>&gt; fourier_ny              =  5.0<br>&gt; fourier_nz              =  8.1<br>these numbers should be integers<br><br><br>&gt; pme_order               =  4<br>&gt; rcoulomb                =  1.5<br>&gt; rvdw                    =  1.5<br>&gt; <br>&gt; Also since I have not specified "vdwtype" in the mdp<br>&gt; file, so it automatically takes Cut-off.<br>&gt; <br>&gt; regards,<br>&gt; Pri...<br>&gt; <br>&gt; <br>&gt; --- Martin H�fling <martin.hoefling @gmx.de=""> wrote:<br>&gt; <br>&gt;&gt; Am Mittwoch, 20. Juni 2007 schrieb priyanka<br>&gt;&gt; srivastava:<br>&gt;&gt;<br>&gt;&gt;&gt; I am currently involved in doing a lipid-peptide<br>&gt;&gt;&gt; simulation under NPAT conditions. The way I have<br>&gt;&gt;&gt; applied NPAT condition is as
 follows:<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; ref_p                   =  0 1<br>&gt;&gt;&gt; compressibility         =  0.0 4.5e-5<br>&gt;&gt; with semiisotropic?<br>&gt;&gt;<br>&gt;&gt;&gt; Although while performing analysis it shows that<br>&gt;&gt; the x<br>&gt;&gt;&gt; and y dimensions are constant but still I am not<br>&gt;&gt; sure<br>&gt;&gt;&gt; about the way I have applied NPAT. Could somone<br>&gt;&gt; please<br>&gt;&gt;&gt; tell me the basis behind it?<br>&gt;&gt; If i am remembering correctly, this is what above<br>&gt;&gt; values should do:<br>&gt;&gt;<br>&gt;&gt; first dimension is x-y second one is the z<br>&gt;&gt; direction. Reference pressure is 1 <br>&gt;&gt; for both (obviously makes sense) whereas<br>&gt;&gt; compressibility set to 0 for x-y <br>&gt;&gt; plane prevents box changes in x and y direction. So<br>&gt;&gt; only z-direction remains <br>&gt;&gt; for adjusting pressure.<br>&gt;&gt;<br>&gt;&gt;&gt; Another thing is I have used vdwtype:
 cutoff.<br>&gt;&gt;&gt; But my intention is not to use the twin range<br>&gt;&gt; cutoff<br>&gt;&gt;&gt; at all!! In my case rlist = rcoulomb = rvdw, i.e.<br>&gt;&gt; all<br>&gt;&gt;&gt; three are same!! Is this alright? (the coulombtype<br>&gt;&gt;&gt;         =  PME)<br>&gt;&gt; Sorry, didn't get that part of your question. Can<br>&gt;&gt; you post relevant parts of <br>&gt;&gt; your mdp and specify it further? Chapter 7 says that<br>&gt;&gt; parameters for PME are <br>&gt;&gt; fourierspacing and pme_order, rcoulomb should not be<br>&gt;&gt; used then.<br>&gt;&gt;<br>&gt;&gt; Cheers<br>&gt;&gt;  Martin<br>&gt;&gt; _______________________________________________<br>&gt;&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt;&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<br>&gt;&gt; Please search the archive at<br>&gt;&gt; http://www.gromacs.org/search before posting!<br>&gt;&gt; Please don't post (un)subscribe requests to the<br>&gt;&gt; list. Use the
 <br>&gt;&gt; www interface or send it to<br>&gt;&gt; gmx-users-request@gromacs.org.<br>&gt;&gt; Can't post? Read<br>&gt;&gt; http://www.gromacs.org/mailing_lists/users.php<br>&gt;&gt;<br>&gt; <br>&gt; <br>&gt; <br>&gt;        <br>&gt; ____________________________________________________________________________________<br>&gt; Pinpoint customers who are looking for what you sell. <br>&gt; http://searchmarketing.yahoo.com/<br>&gt; _______________________________________________<br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search 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<br><br><br>--
 <br>David.<br>________________________________________________________________________<br>David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,<br>Dept. of Cell and Molecular Biology, Uppsala University.<br>Husargatan 3, Box 596,   75124 Uppsala, Sweden<br>phone: 46 18 471 4205  fax: 46 18 511 755<br>spoel@xray.bmc.uu.se spoel@gromacs.org   http://folding.bmc.uu.se<br>++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br>_______________________________________________<br>gmx-users mailing list    gmx-users@gromacs.org<br>http://www.gromacs.org/mailman/listinfo/gmx-users<br>Please search the archive at http://www.gromacs.org/search before posting!<br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to gmx-users-request@gromacs.org.<br>Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br></martin.hoefling></blockquote><br><p>&#32;
      <hr size=1><a href="http://us.rd.yahoo.com/evt=48226/*http://new.toolbar.yahoo.com/toolbar/features/norton/index.php">Get the free Yahoo! toolbar</a> and rest assured with the added security of spyware protection.