CHARMM and MMTSB tool set

mmtsb installation:

./install.sh

mener.c:215:14: error: ‘gets’ was not declared in this scope
   gets(datdir);

add #include <iostream> to mener.c

 

Bsic introduction of CHARMM is here:

https://greenlab2.ams.sunysb.edu/index.php/CHARMM#Compiling_CHARMM

MMTSB tool set is here:

http://www.mmtsb.org/

Most probable errors in charmm:http://simulation.haotui.com/viewthread.php?tid=29157

There are many types of warnings that appear in CHARMM. Here are the possible solutions for some of them.
1. "Element negative. No operation" from the SCALAR modulee.g.
CHARMM>    scalar sca4 sum sca5  ! + <y^2>

CHARMM>    scalar sca4 sum sca2  ! - <y>^2

CHARMM>    scalar sca4 sum sca6  ! + <z^2>

CHARMM>    scalar sca4 sum sca3  ! - <z>^2

CHARMM>    scalar sca4 sqrt
***** LEVEL  0 WARNING FROM <SCALAR> *****
***** Element negative. No operation
******************************************
BOMLEV (  0) IS REACHED - TERMINATING. WRNLEV IS  5
[[ Solution ]]:  This is most likely due to the fact that some atoms' coordinates are missing somehow. Try to find out these atoms and fill their internal coordinates.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2. "zero selection specified" when using HBUILDe.g.
CHARMM>    !Rebuild hydrogen coordinates
CHARMM>    HBUIld SELE ALL END
SELRPN>    667 atoms have been selected out of    667

***** LEVEL  4 WARNING FROM <GTHBUI> *****
***** Zero selection specified
******************************************
BOMLEV (  0) IS NOT REACHED. WRNLEV IS  5
[[ Solution ]]:  change the selection to hydrogen atoms only, i.e.
CHARMM>    !Rebuild hydrogen coordinates
CHARMM>    HBUIld SELE hydrogen END
SELRPN>    344 atoms have been selected out of    667
PRNHBD: CUToff Hydrogen Bond  distance =    0.5000   Angle =   90.0000
CuT switching ON HB dist. =     3.5000  OFf HB dist. =    4.0000
CuT switching ON Hb Angle =    50.0000  OFf Hb Angle =   70.0000
ACCEptor antecedents included
All hydrogen bonds for each hydrogen will be found
Hydrogen bonds between excluded atoms will be kept
[[ Reason ]]:  (from CHARMM hbuild.doc)
"<atom-selection> specify the hydrogens to be (re-)constructed."
"By default (if no selection is specified) these are all unknown hydrogens and lone pairs (this is equivalent to a selection "SELEction (LONE .OR. HYDRogen) .AND..NOT INITial")."
[[ Remarks ]]:  What is said in hbuild.doc may be true. However, if one does not make any selections after the HBUILD keyword, the same warning will still appear.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3. "Failed to allocate memory for natomingp, gptoalst" from gbsw.src
e.g.
***** LEVEL -1 WARNING FROM <GB_lookup>gbsw.src *****
***** Failed to allocate memory for natomingp,gptoalst
******************************************
BOMLEV (  0) IS REACHED - TERMINATING. WRNLEV IS  5
[[  Solution ]];   Most likely this is due to the missing coordinates for some  atoms. Try to find  these atoms and  fill their coordinates.

4. "Zero length string being converted to 0" from PHMD
e.g.
CHARMM>    label PERLDONE

CHARMM>    phmd par 23 wri 25 ph 7.0 npri 100 mass 10 barr 1.75 bartau 2 temp 300 lam sele .not. ( resn TIP3 .or. RESN CLA .or. RESN SOD ) end

PHMD> Continous constant pH molecular dynamics

SELRPN>     36 atoms have been selected out of   2163
TITLE> * STATES FILE FOR HYBRID-PHMD
TITLE> * SYNTAX:
TITLE> * RESNAME MODEL_PKA PARA PARB
TITLE> *  ATOM_TYPE CHARGE(1) CHARGE(2) [RAD(1) RAD(2)]
TITLE> *
WARNING from DECODF -- Zero length string being converted to 0.
PHMD> Grp #  1:ARG  (ARG    1) pKa= 12.50 Prot
PHMD> Grp #  1:ARG  (ARG    1) HH12
[[ Solution ]]:   double check the phmd titration parameter file. Most likely its format is not correct, e.g. leaving out the barrier potential.
results:
CHARMM>    phmd par 23 wri 25 ph 7.0 npri 100 mass 10 barr 1.75 bartau 2 temp 300 lam sele .not. ( resn TIP3 .or. RESN CLA .or. RESN SOD ) end

PHMD> Continous constant pH molecular dynamics

SELRPN>     36 atoms have been selected out of   2163
TITLE> * STATES FILE FOR HYBRID-PHMD
TITLE> * SYNTAX:
TITLE> * RESNAME MODEL_PKA PARA PARB BARR
TITLE> *  ATOM_TYPE CHARGE(1) CHARGE(2) [RAD(1) RAD(2)]
TITLE> *
PHMD> Grp #  1:ARG  (ARG    1) pKa= 12.50 Prot
PHMD> Grp #  1:ARG  (ARG    1) HH12
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5. "extraneous characters" from NBOND settings
e.g.
CHARMM>    label PERLDONE

CHARMM>    nbond elec atom cdie shift vdw vswitch ctonnb 10.0 ctofnb 12.0 cutnb 14.0 cutim 14.0 inbfrq -1 imgfrq -
1 ewald pmew fftx 64 ffty 64 fftz 64  kappa .34 spline order 6
<PME> Total heap storage needed =     1034887
Fill Ewald table: Number of points=     10000 EWXmax=    4.250000
fill erfc table: linear inter has rms error = 0.979220D-08 maximum error = 0.218740D-07
fill erfc table: cubic spline has rms error = 0.360331D-11 maximum error = 0.108438D-10

NONBOND OPTION FLAGS:
ELEC     VDW      ATOMs    CDIElec  SHIFt    VATOm    VSWItch
BYGRoup  NOEXtnd  EWALd
CUTNB  = 14.000 CTEXNB =999.000 CTONNB = 10.000 CTOFNB = 12.000
WMIN   =  1.500 WRNMXD =  0.500 E14FAC =  1.000 EPS    =  1.000
NBXMOD =      5
PME EWALD OPTIONS: KAPPA  =  0.340  QCOR =  0.000  Bspline order = 6
FFTX=  64  FFTY=  64  FFTZ=  64
Using Pub FFT
Real-to-Complex FFT
There are        0 atom  pairs and        0 atom  exclusions.
There are        0 group pairs and        0 group exclusions.
<MAKINB> with mode   5 found  27483 exclusions and   2472 interactions(1-4)
<MAKGRP> found    733 group exclusions.
Generating nonbond list with Exclusion mode = 5
== PRIMARY == SPACE FOR  7400454 ATOM PAIRS AND        0 GROUP PAIRS
VEHEAP> Expanding heap size by        11436032 words.
== PRIMARY == SPACE FOR 11100701 ATOM PAIRS AND        0 GROUP PAIRS
VEHEAP> Expanding heap size by        16990208 words.
== PRIMARY == SPACE FOR 16651072 ATOM PAIRS AND        0 GROUP PAIRS

General atom nonbond list generation found:
11314422 ATOM PAIRS WERE FOUND FOR ATOM LIST
641446 GROUP PAIRS REQUIRED ATOM SEARCHES

**** Warning ****  The following extraneous characters
were found while command processing in CHARMM
CUTIM 14.0  IMGFRQ -1
[[ Comments ]]: CUTIM should not appear in NBOND settings . The actual value of CUTIm should be calculated by "crystal build" based on the information on the water box dimensions (see below).
[[ Solution ]]: Thus CUTIM should be removed from the NBOND settings and IMGFRQ can be kept. After this modification, no more warning is issued.
CHARMM>    label PERLDONE

CHARMM>    crystal defined octa 70.488 70.488 70.488 109.4712206344907 109.4712206344907 109.4712206344907
Crystal Parameters : Crystal Type = OCTA
A     =   70.48800 B    =   70.48800 C     =   70.48800
Alpha =  109.47122 Beta =  109.47122 Gamma =  109.47122

CHARMM>    label PERLDONE

CHARMM>    crystal build cutoff 35.244 noper 0
XBUILD> Building all transformations with a minimum atom-atom
contact distance of less than   35.24 Angstroms.

Range of Grid Search for Transformation     1 :
Lattice Vector A    -3 TO     3
Lattice Vector B    -3 TO     3
Lattice Vector C    -2 TO     2

Note: the actual value of CUTIm is 35.244 in this case.
NSTEP =      500    NSAVC =      500    NSAVV =        0
ISCALE =        0   ISCVEL =        0   IASORS =        1
IASVEL =        1   ICHECW =        1   NTRFRQ =      500
IHTFRQ =        0   IEQFRQ =       50   NPRINT =      500
INBFRQ =       -1   IHBFRQ =        0   IPRFRQ =      500
ILBFRQ =       50   IMGFRQ =       -1    ISEED =          2023018546
Note: the value for IMGFRQ is "-1".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%

1. the default setting in aarex.pl is gb=1, so I edit the aarex.pl and change this value to 0.
2. if for some reason the simulation stoped, and you just type "aarex.pl" would not work on condition replica exchange. But if you rerun the simulation by typing the whole command, the simulation will start from beginning.
3. if you rerun the whole command after the previous one completed, the simulation will continue, but lose the previous energy files.