<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Hi, I'm new to gromacs and I can't figure out what I'm doing wrong with parallel runs. I'm running gromacs 3.3.2 on a 215 residue protein: 1rz4.pdb. The 2 missing residues were added with modloop followed by energy minimization. The protein was put in a dodecahedron box and solvated with a layer of water at least 1.2 nm from the protein and 14 Na+ ions were added to bring the total charge of the system to zero. Energy minimization was done to remove close contacts.<BR>
<BR>
When I run on a single processor my MD runs seem fine (heatup position restrained MD, constant temperature position restrained MD, and now the production MD run is almost 86% complete). But when I run on a cluster some I get strange fatal errors, maybe tripped by boolean checks during writing of energy info to file. During position restrained MD I got:<BR>fatal error in function do_enx (line 343 of enxio.c) <BR>during the call:<BR>gmx_fatal(FARGS,"could not write energies")<BR>The fatal error happened after printing Coulomb1-4 energy in the log file, the next line of headings was NOT printed.<BR>
<BR>
During production run of MD I got:<BR>Source code file: stat.c, line: 283<BR>Fatal error: XTC error<BR>
happening after writing these lines in the log:<BR>
Step Time Lambda<BR> 257600 515.20001 0.00000<BR>
<BR>
These errors do not appear when I run equilibriation and production MD on a single processor on my own machine, but I’d like to take advantage of the cluster to do runs faster and free up my own machine.<BR>
<BR>
Here’s my .mdp file for the production run:<BR>integrator = md<BR>nsteps = 500000<BR>dt = 0.002<BR>nstlist = 10<BR>rlist = 1.2<BR>coulombtype = pme<BR>rcoulomb = 1.2<BR>vdw-type = cut-off<BR>rvdw = 1.2<BR>tcoupl = Berendsen<BR>tc-grps = protein non-protein<BR>tau-t = 0.1 0.1<BR>ref-t = 300 300<BR>pcoupl = Berendsen<BR>tau-p = 1.0<BR>compressibility = 5e-5<BR>ref-p = 1.0<BR>nstxout = 1000<BR>nstvout = 1000<BR>nstxtcout = 100<BR>nstenergy = 100<BR>constraints = allbonds<BR>; Generate velocites is on at 300 K.<BR>gen_vel = yes<BR>
<BR>
The same .top, .gro, .mdp and .tpr files are used on both single processor and parallel runs, with the GROMOS96 53a6 force field. Any help you can offer is much appreciated!<BR>
<BR>
Thanks,<BR>Patricia Francis-Lyon<BR><BR><br /><hr />Get the power of Windows + Web with the new Windows Live. <a href='http://www.windowslive.com?ocid=TXT_TAGHM_Wave2_powerofwindows_012008' target='_new'>Get it now!</a></body>
</html>