[gmx-users] atom types
Gavin Melaugh
gmelaugh01 at qub.ac.uk
Mon Aug 1 16:25:54 CEST 2011
Hi all
A very quick question. I have an atom-type labelled CH in the atom-types
with a particular charge, and in the atom list I assign some of these
specific atoms with zero charge as below. When I generate 1,4
interactions using gen_pairs =yes, what charge for the CH type does it
use? Does gromacs assign the CH with the different charge as a new atom
type.
;Parameter level
[defaults]
; nbfunc comb-rule gen-pairs fudgeLJ fudgeQQ
1 3 yes 0.5 0.5
[atomtypes]
;type mass charge ptype sigma(nm)
epsilon(kjmol-1)
CB 12.011000 0.000000 A 0.355000 0.292880
CA 12.011000 -0.115000 A 0.355000 0.292880
HC 1.008000 0.115000 A 0.242000 0.125520
CU 13.019000 0.265000 A 0.350000 0.334720
NU 14.007000 -0.597000 A 0.325000 0.711280
CH 13.019000 0.332000 A 0.385000 0.334720
C3 15.035000 0.000000 A 0.391000 0.669440
C2 14.027000 0.000000 A 0.390500 0.493712
;Molecular level
[moleculetype]
; name nrexcl
isotridecylcage 3
[atoms]
.................
72 CH 1 CGE CH 24 0.3320 13.0190
73 C2 1 CGE C2 25 0.0000 14.0270
74 C2 1 CGE C2 25 0.0000 14.0270
75 C2 1 CGE C2 25 0.0000 14.0270
76 C2 1 CGE C2 26 0.0000 14.0270
77 C2 1 CGE C2 26 0.0000 14.0270
78 C2 1 CGE C2 26 0.0000 14.0270
79 C2 1 CGE C2 27 0.0000 14.0270
80 C2 1 CGE C2 27 0.0000 14.0270
81 C2 1 CGE C2 27 0.0000 14.0270
82 C2 1 CGE C2 28 0.0000 14.0270
83 CH 1 CGE CH 28 0.0000 13.0190
84 C3 1 CGE C3 29 0.0000 15.0350
Many Thanks
Gavin
More information about the gromacs.org_gmx-users
mailing list