[gmx-users] What do parameters in topology file mean ?
大木啓輔
keisuke.ohki at keio.jp
Mon Nov 28 06:22:16 CET 2016
Dear gmx users,
I have questions about topology file parameters.
I tried get topology files with topolbuild. But I was required to fix the *.top file
because H2O parameters were not given.
I could give some parameters in atoms and bonds, but I don't know some meanings of parameters.
topology file has some sections, for example atoms, bonds and angles etc..
I show the topology file below {
[ bonds ]
; ai aj funct b0 kb
1 2 1 0.10900 284512. ; C3- H
1 3 1 0.10900 284512. ; C3- H
: : : : :
6 7 ; HW- OW
: : : : :
[ angles ]
; ai aj ak funct th0 cth
3 1 2 1 107.800 276.1440 ; H- C3- H
4 1 2 1 107.800 276.1440 ; H- C3- H
: : : : : :
43 42 44 ; HW- OW- HW-
: : :
}
In bonds section, what do b0 and kb parameters mean?
I think respective they show distance of bond and energy of bond.
In angle section, what do th0 and cth mean?
Probably th0 means angle degree. But I cannot expect cth meaning.
Then, At the last of topology file {
; Include Position restraint file
; WARNING: Position restraints and distance restraints ought not be done together
#ifdef POSRES
#include "posreHYDRATE_4PT.itp"
#endif
; Include water topology
#include "spce.itp"
#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
; i funct fcx fcy fcz
1 1 1000 1000 1000
#endif
}
What do these mean?
I think that ifdef POSRES is about restrain field.
And #include "spce.itp" means that SPC/E water model is exported to this topology file
that has no parameters of water? So I have no need to give parameters to water?
Could you tell me solution to this problem, please?
Keisuke Ohki
More information about the gromacs.org_gmx-users
mailing list