<br><br><div class="gmail_quote">On Thu, Mar 1, 2012 at 5:43 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
<br>
Steven Neumann wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
<br>
<br>
On Thu, Mar 1, 2012 at 3:58 PM, Steven Neumann <<a href="mailto:s.neumann08@gmail.com" target="_blank">s.neumann08@gmail.com</a> <mailto:<a href="mailto:s.neumann08@gmail.com" target="_blank">s.neumann08@gmail.com</a>><u></u>> wrote:<br>
<br>
<br>
<br>
On Thu, Mar 1, 2012 at 3:32 PM, Steven Neumann<br></div><div class="im">
<<a href="mailto:s.neumann08@gmail.com" target="_blank">s.neumann08@gmail.com</a> <mailto:<a href="mailto:s.neumann08@gmail.com" target="_blank">s.neumann08@gmail.com</a>><u></u>> wrote:<br>
<br>
<br>
<br>
On Thu, Mar 1, 2012 at 2:55 PM, Justin A. Lemkul<br></div><div><div class="h5">
<<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
Steven Neumann wrote:<br>
<br>
Dear Gmx Users,<br>
I am trying to run nvt simulation (equilibration) of<br>
the protein in implicit solvent. My mdp:<br>
integrator = md ; leap-frog integrator<br>
<br>
nsteps = 1000000 ; 0.0005 * 1000000 = 0.5 ns<br>
<br>
dt = 0.0005 ; 0.5 fs<br>
<br>
; Output control<br>
<br>
nstxout = 10000<br>
<br>
nstxtcout = 10000 ; xtc compressed trajectory output<br>
every 2 ps<br>
<br>
nstenergy = 10000 ; save energies every 2 ps<br>
<br>
nstlog = 10000 ; update log file every 2 ps<br>
<br>
; Bond parameters<br>
<br>
continuation = no<br>
<br>
constraints = none<br>
<br>
; Neighborsearching<br>
<br>
ns_type = simple ; search neighboring grid cells<br>
<br>
nstlist = 0 ; 10 fs<br>
<br>
rlist = 0 ; short-range neighborlist cutoff (in nm)<br>
<br>
; Infinite box with no cutoffs<br>
<br>
pbc = no<br>
<br>
rcoulomb = 0 ; short-range electrostatic cutoff (in nm)<br>
<br>
coulombtype = cut-off<br>
<br>
vdwtype = cut-off<br>
<br>
rvdw = 0 ; short-range van der Waals cutoff (in nm)<br>
<br>
epsilon_rf = 0<br>
<br>
comm_mode = angular ; remove angular and com motions<br>
<br>
; implicit solvent<br>
<br>
implicit_solvent = GBSA<br>
<br>
gb_algorithm = OBC<br>
<br>
gb_epsilon_solvent = 80.0<br>
<br>
sa_surface_tension = 2.25936<br>
<br>
rgbradii = 0<br>
<br>
sa_algorithm = Ace-approximation<br>
<br>
nstgbradii = 1<br>
<br>
; Temperature coupling is on<br>
<br>
Tcoupl = v-rescale<br>
<br>
tau_t = 0.1<br>
<br>
tc_grps = system<br>
<br>
ref_t = 298<br>
<br>
; Velocity generation<br>
<br>
gen_vel = yes ; Velocity generation is on<br>
<br>
gen_temp = 298.0<br>
<br>
gen_seed = -1<br>
<br>
Then after grompp I am trying to run the simulation on<br>
the cluster:<br>
<br>
mdrun -pd -deffnm nvt500ps<br>
<br>
My log file:<br>
<br>
Back Off! I just backed up nvt500ps.log to<br>
./#nvt500ps.log.1#<br>
<br>
Reading file nvt500ps.tpr, VERSION 4.5.4 (single precision)<br>
<br>
Starting 8 threads<br>
<br>
Back Off! I just backed up nvt500ps.trr to<br>
./#nvt500ps.trr.1#<br>
<br>
Back Off! I just backed up nvt500ps.xtc to<br>
./#nvt500ps.xtc.1#<br>
<br>
Back Off! I just backed up nvt500ps.edr to<br>
./#nvt500ps.edr.1#<br>
<br>
starting mdrun 'Protein'<br>
<br>
1000000 steps, 500.0 ps.<br>
<br>
Segmentation fault<br>
<br>
Do you have any clue what is happening?<br>
<br>
<br>
Try running in serial or with a maximum of 2 threads. Your<br>
problem could be related to<br></div></div>
<a href="http://redmine.gromacs.org/__issues/777" target="_blank">http://redmine.gromacs.org/__<u></u>issues/777</a><br>
<<a href="http://redmine.gromacs.org/issues/777" target="_blank">http://redmine.gromacs.org/<u></u>issues/777</a>>. You will need to<div class="im"><br>
upgrade to 4.5.5 (serial should work on 4.5.4).<br>
<br>
-Justin<br>
<br>
Thank you. What do you mean by running in serial? Well... with 2<br>
threads it does not make sense to use implicit solvent. Will<br>
4.5.5 resolve this problem?<br>
<br>
Indeed, it works in serial. Will version 4.5.5 resolve it run it on<br>
e.g. 12 nodes?<br>
<br>
Sorry, the same problem :<br>
<br>
</div></blockquote>
<br>
I think, as stated in the redmine issue cited before, the limit is 2 threads/processors.<br>
<br>
-Justin<br></blockquote><div><br><br>I tried on one thread and on two threads and the error is still the same. Without particle decomposition the problem remains. Any suggestions? <br></div><div><br>Steven <br></div><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
<br>
Initial temperature: 298.398 K<br>
<br>
Started mdrun on node 0 Thu Mar 1 15:56:04 2012<br>
<br>
Step Time Lambda<br>
<br>
0 0.00000 0.00000<br>
<br>
Energies (kJ/mol)<br>
<br>
Bond U-B Proper Dih. Improper Dih. CMAP Dih.<br>
<br>
2.27391e+02 4.27461e+02 1.27733e+03 2.73432e+01 -6.95645e+02<br>
<br>
GB Polarization Nonpolar Sol. LJ-14 Coulomb-14 LJ (SR)<br>
<br>
inf 3.59121e+02 6.85627e+02 3.68572e+04 -7.57140e+02<br>
<br>
Coulomb (SR) Potential Kinetic En. Total Energy Conserved En.<br>
<br>
-3.30643e+04 inf nan nan nan<br>
<br>
Temperature Pressure (bar)<br>
<br>
nan 0.00000e+00<br>
<br>
<br>
-<br></div>
==============================<u></u>__==========<div class="im"><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></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu/" target="_blank">http://vt.edu/</a>> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><br>
<tel:%28540%29%20231-9080><br>
<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank"><u></u>vt.edu/Pages/Personal/justin</a><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>__==========<div class="im"><br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<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><div class="im"><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<br></div>
<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><div class="im"><br>
<<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<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>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>
<<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>><br>
<br>
<br><span class="HOEnZb"><font color="#888888">
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
<br>
-- <br></font></span><div class="HOEnZb"><div class="h5">
==============================<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> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><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>