<br>Hello,<br><br>I am sorry, I did mean 5 nm box.<br><br>I did do the energy minimization but somehow overlooked the error message at that stage.<br><br>Things seems to be working now.<br><br>Thanks for the time.<br><br>Shivangi<br>
<br><div class="gmail_quote">On Tue, Mar 29, 2011 at 2:06 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div class="im"><br>
<br>
shivangi nangia wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hello gmx-users,<br>
<br>
I have a 5 ang box with equal number of water and methanol molecules (1800), a neutral peptide, 8 hydroxide ions and 10 Li+ ions. The system is overall +2 charged.<br>
<br>
</blockquote>
<br></div>
If your box truly is 5 A (0.5 nm), grompp should have failed based on your cutoffs. I assume you have a 5-nm box?<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
I am trying to to do a NVT equilibration which runs into the following error:<br>
<br>
<br>
t = 0.000 ps: Water molecule starting at atom 11816 can not be settled.<br>
Check for bad contacts and/or reduce the timestep.<br>
<br>
</blockquote>
<br></div>
Did you do energy minimization first? If so, what was the outcome? Were the energies and forces acceptable? What force field are you using? The generic information about your problem can be found at:<br>
<br>
<a href="http://www.gromacs.org/Documentation/Terminology/Blowing_Up" target="_blank">http://www.gromacs.org/Documentation/Terminology/Blowing_Up</a><br>
<br>
-Justin<div><div></div><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
the nvt.mdp is:<br>
<br>
title = hist NVT equilibration<br>
define = -DPOSRES ; position restrain the protein<br>
; Run parameters<br>
integrator = md ; leap-frog integrator<br>
nsteps = 50000 ; 2 * 50000 = 100 ps<br>
dt = 0.002 ; 2 fs<br>
; Output control<br>
nstxout = 100 ; save coordinates every 0.2 ps<br>
nstvout = 100 ; save velocities every 0.2 ps<br>
nstenergy = 100 ; save energies every 0.2 ps<br>
nstlog = 100 ; update log file every 0.2 ps<br>
; Bond parameters<br>
continuation = no ; first dynamics run<br>
constraint_algorithm = lincs ; holonomic constraints<br>
constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained<br>
lincs_iter = 1 ; accuracy of LINCS<br>
lincs_order = 4 ; also related to accuracy<br>
; Neighborsearching<br>
ns_type = grid ; search neighboring grid cells<br>
nstlist = 5 ; 10 fs<br>
rlist = 1.0 ; short-range neighborlist cutoff (in nm)<br>
rcoulomb = 1.0 ; short-range electrostatic cutoff (in nm)<br>
rvdw = 1.0 ; short-range van der Waals cutoff (in nm)<br>
; Electrostatics<br>
coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics<br>
pme_order = 4 ; cubic interpolation<br>
fourierspacing = 0.16 ; grid spacing for FFT<br>
; Temperature coupling is on<br>
tcoupl = V-rescale ; modified Berendsen thermostat<br>
tc-grps = Protein Non-Protein ; two coupling groups - more accurate<br>
tau_t = 0.1 0.1 ; time constant, in ps<br>
ref_t = 250 250 ; reference temperature, one for each group, in K<br>
; Pressure coupling is off<br>
pcoupl = no ; no pressure coupling in NVT<br>
; Periodic boundary conditions<br>
pbc = xyz ; 3-D PBC<br>
; Dispersion correction<br>
DispCorr = EnerPres ; account for cut-off vdW scheme<br>
; Velocity generation<br>
gen_vel = yes ; assign velocities from Maxwell distribution<br>
gen_temp = 250 ; temperature for Maxwell distribution<br>
gen_seed = -1 ; generate a random seed<br>
<br>
I did try to reduce the timestep from 2 fs to 1 fs but ran into the same error.<br>
<br>
<br>
<br>
Searching through a few archives I got a hint it is something do to with water's density per cubic nm.<br>
Is that whats going on here?<br>
If yes, then how do I randomly delete some water molecules from my system?<br>
<br>
Please help.<br>
<br>
Thanks,<br>
Shivangi<br>
<br>
<br>
<br>
<br>
<br>
</blockquote>
<br></div></div>
-- <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> | <a href="tel:%28540%29%20231-9080" target="_blank">(540) 231-9080</a><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><font color="#888888">
-- <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/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/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/Support/Mailing_Lists</a><br>
</font></blockquote></div><br>