I want to come back to the question of the NMA in the Gromacs :)<br><br>I&#39;ve found in manual possible algorithms of this analysis-  I must calculate Hessian matrix via Md-run and then calculate modes with g_nmens program<br>
<br><br>1- I&#39;ve performed CG minimization of my initianl structure<br><br>2- I&#39;ve loadet to the gromp .tpr file of my minimized structure<br><br>3-  I&#39;ve used mdrun and obtain error <br>Fatal error:<br>Constraints present with Normal Mode Analysis, this combination is not supported<br>
<br>So I suppose that the problem with my .mdp file wich consist of the parametres for the NMA of Lyzosyme in water :)<br><br>integrator    = nm        ; Normal Mode Analysis<br><br>; Parameters describing how to find the neighbors of each atom and how to calculate the interactions<br>
nstlist        = 1        ; Frequency to update the neighbor list and long range forces<br>ns_type        = grid        ; Method to determine neighbor list (simple, grid)<br>rlist        = 1.0        ; Cut-off for making neighbor list (short range forces)<br>
coulombtype    = PME        ; Treatment of long range electrostatic interactions<br>rcoulomb    = 1.0        ; Short-range electrostatic cut-off<br>rvdw        = 1.0        ; Short-range Van der Waals cut-off<br>pbc        = xyz         ; Periodic Boundary Conditions (yes/no)<br>
<br><br>How I should edit this configure file for my NMA?  <br>If I edit the above parametries for rvdw etc in accardance to the <a href="http://www.gromacs.org/Documentation/How-tos/Normal_Mode_Analysis?highlight=nm">http://www.gromacs.org/Documentation/How-tos/Normal_Mode_Analysis?highlight=nm</a> I&#39;ve obtained error that<br>
<br>ERROR 1 [file nma.mdp]:<br>  With coulombtype = PME, rcoulomb must be equal to rlist<br>  If you want optimal energy conservation or exact integration use PME-Switch<br><br><br>ERROR 2 [file nma.mdp]:<br>  With vdwtype = Cut-off, rvdw must be &gt;= rlist<br>
<br>What are the most suitable coulombtype for the NMA and combination of the other options?<br><br><br>James<br><br> <br>