Justin,<br><br>I've built my system in accordance to the first way from your Biphgastic system tutorial.<br><br>I've defined system with slight biger ( on 1 nm in each dimension) dimensions that I needed and place maximym CCl4 molecules in that box by<br>
<br><font face="Arial"><font size="3"><pre>genbox -ci ccl4.gro -nmol 900 -box 9.6 7.5 4 -o new_box.gro</pre></font></font>So I've supposed that I've defined maximum density for such box. How else I could specify starting density of my box?<br>
<br>2) I've tried use 1.4 VDV cutoff and crash of my system happened almost alfter first 0.5 ns of simulation. By the way in the KALP tutorial I've found Cutoff = 1.2. Why in that case you've used such uncommon values ?<br>
<br>James <br><br><div class="gmail_quote">2012/2/5 Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><br>
<br>
James Starlight wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Justin,<br>
<br>
<br>
Larger cutoffs would not only make the problem worse (read the error<br>
message carefully and consider the minimum image convention), but it<br>
would also potentially break the validity of the force field model.<br>
The vdW cutoff for Gromos96 should be 1.4 nm, not 0.9 nm, in any case.<br>
<br>
<br>
This sounds controversially alitle :)<br>
<br>
Today I've tried such simulation with the 1.2 cutoffs and obtain the same error on the 3rd ns ( density was 1.35 )<br>
</blockquote>
<br></div>
Still wrong. Set rvdw = 1.4 and rcoulomb = 0.9 for Gromos96. Don't play haphazardly with cutoffs. Use what the force field dictates unless you can demonstrate that what you are doing is superior. The original value of rcoulomb was 0.8 with reaction field, but 0.9 is fine in the case of PME.<div>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
From the error message I've understood that the problem was due to the Z- box vector wich was decreased up to the 3 nm ( this is smaller rhan 2x 1.2 of R_list cutoff ). So in accordance to that message I understood that cutof must be > 1.5 in the case of my system didnt it?<br>
<br>
</blockquote>
<br></div>
No, absolutely not. You don't *want* the cutoffs to be equal to twice the box vector, they must always be *less* than this amount to avoid spurious PBC interactions due to minimum image violations. Are you familiar with this concept and its implications?<br>
<br>
<a href="http://www.gromacs.org/Documentation/Terminology/Minimum_Image_Convention" target="_blank">http://www.gromacs.org/<u></u>Documentation/Terminology/<u></u>Minimum_Image_Convention</a><br>
<br>
You said in your first message:<div><br>
<br>
"2) Initially I've tried to use membrane mimicking dimensions 8.6 6.5 3.0 but during equilibration my box was shrinked to the 8.0 6.0 2.0"<br>
<br></div>
So your problem comes from the box shrinking to 2.0 nm in the z-dimension, which is less than twice your longest cutoff. Your box is compressing in all directions, which means your initial configuration is too diffuse and the application of NPT conditions is squishing it together to arrive at the density the force field model predicts. Build a system with an initial density closer to the target value.<div>
<div><br>
<br>
-Justin<br>
<br>
-- <br>
==============================<u></u>==========<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.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<br>
-- <br>
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/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/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/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>