From: JC Gumbart (gumbart_at_ks.uiuc.edu)
Date: Sun May 13 2012 - 09:27:26 CDT
Probably the atom type needs to start with H to get flagged properly. I
don’t know if the check for rigid bonds also relies on typing, but better
safe than sorry.
From: owner-namd-l_at_ks.uiuc.edu [mailto:owner-namd-l_at_ks.uiuc.edu] On Behalf
Of ypg
Sent: Sunday, May 13, 2012 9:14 AM
To: JC Gumbart
Cc: namd-l
Subject: RE: namd-l: five 'ERROR' in the log file of NAMD2.9
Oh, I think that is the point, because I have defined five new type 'H'
atoms in my parameter file for convenience.
Thank you so much for your help!!
BTW, would such atoms cause any problems in simulation? Since my timestep is
2fs, I think they should be constrained like other hydrogen atoms in
principle.
2012-05-13
_____
ypg
_____
发件人: JC Gumbart
发送时间: 2012-05-13 21:23:48
收件人: 'ypg'; 'namd-l'
抄送:
主题: RE: namd-l: five 'ERROR' in the log file of NAMD2.9
I grepped the source code and found it has to do with GBIS (specifically in
Molecule.C)
if ( (atoms[i].mass<1.5) != (lcpoParamType[i]==0) ) {
CkPrintf("ERROR\n");
So apparently you have some atoms that aren’t hydrogen (lcpo type of 0) but
have a mass of less than 1.5?
From: owner-namd-l_at_ks.uiuc.edu [mailto:owner-namd-l_at_ks.uiuc.edu] On Behalf
Of ypg
Sent: Sunday, May 13, 2012 2:15 AM
To: namd-l
Subject: namd-l: five 'ERROR' in the log file of NAMD2.9
Dear NAMD experts,
There are five strange ERROR symbol appear in my log file produced by NAMD2.
9 just like this:
Info: PARAMETER file: CHARMM format!
Info: PARAMETERS
/gpfsTMP/wangwei/ypg/NAMD/ff/par_all27_prot_lipid_CRO2.txt
Info: USING ARITHMETIC MEAN TO COMBINE L-J SIGMA PARAMETERS
Info: BINARY COORDINATES YFP+_eq3.restart.coor
Info: SUMMARY OF PARAMETERS:
Info: 194 BONDS
Info: 473 ANGLES
Info: 601 DIHEDRAL
Info: 51 IMPROPER
Info: 6 CROSSTERM
Info: 123 VDW
Info: 0 VDW_PAIRS
Info: 0 NBTHOLE_PAIRS
ERROR
ERROR
ERROR
ERROR
ERROR
Info: TIME FOR READING PSF FILE: 0.243933
Info: TIME FOR READING PDB FILE: 0.0932858
These ERRORs would not appear if I use NAMD2.7.
It seems that there is something wrong in my force field file, and I indeed
have modified my parameter file to define a new type residue. However, these
ERRORs would not terminate my simulations and I can still obtain plausible
MD results. Also, ERRORs only appear when my system involve the new type
residue.
Now I'm really worrying about the validity of my simulated results. I will
be highly appreciate for any help!
ypg
_____
2012-05-13
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:32 CST