It could be an error in protein structure ( because 1 - 9 interaction ), and we usually overcome similar issues by a pre -energy minimization before submitting to gromacs.<br><br>i m carryin a protein explicit solvent simulation &nbsp;using gromacs &nbsp;
3.3<br>for intially minimizing the hydrogen i did some small mini and equi in the<br>vaccum by maintaing the restraints<br>then i added ions and water but in the very first minimization i m getting<br>the following error<br>
<br><br>1-4 interaction between 1 and 9 at distance 8.454 which is larger than the<br>1-4 table size 1.000 nm<br>These are ignored for the rest of the simulation<br>This usually means your system is exploding,<br><br><br>
in the log file im geting :<br><br>4 inconsistant shifts check ur topology<br><br>following is my input file :<br><br>title &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1min in water<br>define &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= -DPOSRES<br>integrator &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= steep<br>
dt &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 0.001<br>constraints &nbsp; &nbsp; &nbsp; &nbsp; = none<br>nsteps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 500<br>nstlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10<br>rcoulomb &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1<br>rlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1<br>rvdw &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1.4<br>ns_type &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = grid
<br>coulombtype &nbsp; &nbsp; &nbsp; &nbsp; = PME<br>pbc &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = xyz<br>gen_vel &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = no<br>fourierspacing &nbsp; &nbsp; &nbsp;= 0.12<br>pme_order &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 6<br>ewald_rtol &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1e-5<br>optimize_fft &nbsp; &nbsp; &nbsp; &nbsp;= yes<br>comm_mode &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = angular
<br>emtol &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 100.0<br>emstep &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 0.01<br><br><br>after adding the ions i have modified the topology file,<br>plz tel me what is the exact reason for this error or what are the &nbsp;other<br>things i should check in relation to the &nbsp;error
<br><br><br>thanks<br><br><br>Gurpreet<br>IIT MADRAS<br clear="all"><br>-- <br>M.N. MANOJ<br>PROJECT LEADER - BIOINFORMATICS<br>VITTAL MALLYA SCIENTIFIC RESEARCH FOUNDATION<br>BANGALORE<br>INDIA