Hello Justin,<br><br>Thanks for your comments and also I really liked your reasoning for getting no answer from people that are not awake :)))) You are right ;)<br><br>Actually, I have no charged particles in my system. All I have is carbon and hydrogen. I read through the manual meticulously about nonbonded interaction terms. I would like you to correct me if I am wrong. See the summary below (my understanding) please:<br>
<br>Treating nonbonded interactions for a hydrocarbon system:<br><br>1-Short range: <br><br>vdw: shift<br>electrostatics: shift<br><br>so the settings would be :<br>rvdw=1<br>rcoulomb=1<br>rlist=1.1<br>rcoulomb-switch=0 <br>
(when this is zero means shift function is added to actual function for all r?)<br>
<br>I have a question here: in manual I read setting this last parameter nonzero means it acts as a switch function. I thought if it is nonzero shift acts just at a nonzero-r between zero and r_c (rcoulomb or rvdw)<br><br>
2-Long range <br><br>vdw: Dispcorr:Enerpres<br>electrostatics: PME <br><br>There are 2 issues from last post (quotes) I am unsure about. <br><br>"in md.mdp file I have only coulombtype=PME" <br>
<br>
"This contradicts your statement above. You said you had only uncharged particles, but clearly this is not the case."<br><br>What you mean is Basically that I dont need to consider PME for hydrocarbon system? <br>
<br>and if you mean I am using wrong setting:<br>
<br>I do not know how else I can treat LR electrostatics in my system. I read in the manual that PME almost for all cases performs better than Ewald and also PPPM cant do P coupling so I picked PME :) <br><br>Also you said:<br>
"A method like PME is significantly more accurate, but you can also
achieve reasonable results using switch, shift, etc. Read about the
relevant mechanics of the force field you're using for what might be
appropriate."<br><br>From what your statement above I realize that PME and shift or switch func's are different alternatives. To be honest I am a little confused about the methods for treating nonbonded interactions. Thats why I wrote that little summary above so that you can correct me. This could be useful for other beginners as well. My understanding from reading the manual is that PME is for LR and switch and shift functions for SR. How could they be alternatives?<br>
<br>Thank you,<br><br><div class="gmail_quote">On 9 July 2010 20:15, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div><br>
<br>
jojo J wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear experts,<br>
<br>
I posted this message before but I think you did not see it since I noticed messages after this have been replied. I have a hydrocarbon <br>
</blockquote>
<br></div>
Messages do not receive replies in any sort of order. Probably no one has had any relevant comments yet, or they live in a time zone where they are not awake :)<div><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
system with no charged particles. Today In one of the GROMACS tutorials I saw the following setting in em.mdp file:<br>
fourierspacing = 0.12<br>
;fourier_nx = 0<br>
;fourier_ny = 0<br>
;fourier_nz = 0<br>
;pme_order = 4<br>
;ewald_rtol = 1e-5<br>
;optimize_fft = yes<br>
<br>
has been used for PME . Sofar I did not even include coulombtype=PME in my em.mdp file. What I had was:<br>
constraints = all-bonds<br>
integrator = steep<br>
dt = 0.002 ; ps !<br>
nsteps = 200<br>
nstlist = 10<br>
ns_type = grid<br>
rlist = 1.0<br>
rcoulomb = 1.0<br>
rvdw = 1.0<br>
;<br>
; Energy minimizing stuff<br>
;<br>
emtol = 1000.0<br>
<br>
Can you please guide me whether it is necessary to include these lines in em file for system having only partial charges like hydrocarbon (alkane). Also, in md.mdp file I have only coulombtype=PME and not the <br>
</blockquote>
<br></div>
This contradicts your statement above. You said you had only uncharged particles, but clearly this is not the case.<div><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
above Ewald settings. Do you recommend to include those settings? <br>
</blockquote>
<br></div>
For energy minimization, the effects of plain cutoffs may or may not be significant. Using plain cutoffs for actual simulation is a bad idea. The artifacts are well-documented and for modern simulations any good reviewer should raise an immediate concern if plain cutoffs were used and not thoroughly justified. A method like PME is significantly more accurate, but you can also achieve reasonable results using switch, shift, etc. Read about the relevant mechanics of the force field you're using for what might be appropriate.<div>
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Without the mentioned lines I am getting reasonabel results (potential energy, kinetic, ,, also T coupling works perfectly). If I had better to unclude them what is the best value for fourierspacing?<br>
<br>
</blockquote>
<br></div>
I don't think there's ever been any sort of systematic study on the effects of Fourier grid spacing. The trade-off is between speed and accuracy. Larger values speed up calculations, but make PME less accurate. Smaller values have the opposite effect. You could do a small test for yourself (0.12, 0.16, 0.18...) to see if there are any adverse effects.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Thank you,<br>
<br>
</blockquote><div>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
-- <br></div><div><div></div><div>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br><input type="hidden"><input type="hidden"><div>
</div>