[gmx-users] CG Graphene - dihedral problem/parameter

Sajjad Kavyani kavyani.sajjad at gmail.com
Tue Apr 25 10:27:26 CEST 2017

Dear experts,
I want to simulate a graphene sheet with MARTINI force field. (In MARTINI
FF, the banzene is represented with 3 beads)
The problem I encountered with is the dihedral parameter.
I searched papers and found that the dihedral parameters for the sheet are
K=3.1, n=2 and phi =180. just that

The parameter (itp) and the gro I used is in the zip-file below:


It is noteworthy that I tried both dihedral type
3   4.000     0.000     -4.000    0.000     0.000     0.000
1     180.00   2.0      2

but in every EM run, the forces are high.

Could you please help me how to place the parameter in the topology? and
what should order of indices be? Is they correct?


More information about the gromacs.org_gmx-users mailing list