Re: [AMBER] Analysis of minimization stage

From: Valentina Romano <valentina.romano.unibas.ch>
Date: Thu, 15 May 2014 13:07:43 +0000

Hi

>Look at the components of the energy at this step. Compare them to the
>individual components of step "0" in your MD run. Report any differences.

I did it. Few components are different:

Min.

   NSTEP ENERGY RMS GMAX NAME NUMBER
   2500 -1.0766E+05 7.5114E-01 7.8713E+01 C6 3838

 BOND = 7906.5852 ANGLE = 529.6831 DIHED = 2310.5740
 VDWAALS = 19126.7051 EEL = -147423.9068 HBOND = 0.0000
 1-4 VDW = 822.2452 1-4 EEL = 8936.0878 RESTRAINT = 134.0700

 EAMBER = -107792.0264

MD

 NSTEP = 0 TIME(PS) = 0.000 TEMP(K) = 438.35 PRESS = 0.0
 Etot = -89978.9428 EKtot = 25598.9152 EPtot = -115577.8580
 BOND = 120.7537 ANGLE = 529.6831 DIHED = 2310.5740
 1-4 NB = 822.2452 1-4 EEL = 8936.0877 VDWAALS = 19126.7051
 EELEC = -147423.9068 EHBOND = 0.0000 RESTRAINT = 0.0000
 Ewald error estimate: 0.1085E-03

I think that the RESTRAINT component should not be 0.0 at the step 0 of MD. Could this cause a problem during the MD?

>What are your starting velocities? You have a temperature of 508K at step
>2 of the MD run. And you already have potential energy that is 10,000
>kcal/mol higher than in the minimization. You need to try to track down
>where such difference are coming from.

The input file for the MD was:

 &cntrl
  imin=0,
  irest=0,
  ntx=1,
  ig=-1,
  ntb=1,
  ntr=1,
  cut=10,
  ntc=2,
  ntf=2,
  tempi=300.0,
  temp0=300.0,
  ntt=3,
  gamma_ln=5.0,
  nstlim=15000, dt=0.001,
  ntpr=1, ntwx=1, ntwr=1000
  restraint_wt = 10.0,
  restraintmask=':1-246.CA,C,O,N | .3834-3843'

Since I want to relax in a NVT ensemble solvent molecules of the structure previously minimized, might be more correct set irest=1 and ntx=7?


Vale





_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu May 15 2014 - 06:30:02 PDT
Custom Search