[gmx-users] Polymer Dihedral Angle Definition

[decaiyu]

Dear All,

When I look at the dihedral angle definition (see below) for gromos43a1
force field, most equilibrium angles are either 180 or 0 degrees.
I assume that the second number represent the strength of dihedral angle
term and the third one is the power.
How do I define the gd term involving a tetrahedral structure, such as
C-C-N-(C3). The equilibrium angle here is not 180 or 0 degree. How do I
decide the value for the second and third term?

; ICP(H)[N]  CP[N] PD[N] NP[N]
;
#define gd_1    180.000       5.86          2
; -C-C- 1.4
;
#define gd_2    180.000       7.11          2
; -C-OA- (at ring)      1.7
;

Best Regards,
decaiyu

--
View this message in context: http://gromacs.5086.x6.nabble.com/Polymer-Dihedral-Angle-Definition-tp5013768.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.