[gmx-users] NPT equilibration in KALP15-DPPC

Shima Arasteh shima_arasteh2001 at yahoo.com
Sun Aug 26 14:54:30 CEST 2012



 Hi,

For NPT equilibration of KALP15 in DPPC following the Justin's tutorial, when I run the grompp command, I get 4 notes. Is it necessary to pay attention to the notes ? Are they notable?


NOTE 1 [file npt.mdp]:
  You are using a cut-off for VdW interactions with NVE, for good energy
  conservation use vdwtype = Shift (possibly with DispCorr)


NOTE 2 [file npt.mdp]:
  You are using a cut-off for electrostatics with NVE, for good energy
  conservation use coulombtype = PME-Switch or Reaction-Field-zero

Generated 837 of the 2346 non-bonded parameter combinations
Excluding 3 bonded neighbours molecule type 'Protein'
Excluding 3 bonded neighbours molecule type 'DPPC'
Excluding 2 bonded neighbours molecule type 'SOL'
Excluding 1 bonded neighbours molecule type 'CL'

NOTE 3 [file topol.top, line 929]:
  The bond in molecule-type Protein between atoms 4 N and 5 H has an
  estimated oscillational period of 1.0e-02 ps, which is less than 10 times
  the time step of 1.0e-03 ps.
  Maybe you forgot to change the constraints mdp option.

Number of degrees of freedom in T-Coupling group rest is 33897.00

NOTE 4 [file npt.mdp]:
  You are using a plain Coulomb cut-off, which might produce artifacts.
  You might want to consider using PME electrostatics.


Reading Coordinates, Velocities and Box size from old trajectory
Will read whole trajectory
Last frame         -1 time  100.000   
Using frame at t = 100 ps
Starting time for run is 0 ps
Largest charge group radii for Van der Waals: 0.239, 0.234 nm
Largest charge group radii for Coulomb:       0.239, 0.234 nm
This run will generate roughly 1 Mb of data

There were 4 notes

Thanks for your suggestions in advance.

Sincerely,
Shima



More information about the gromacs.org_gmx-users mailing list