Dear Gromacs Users,
I’m trying to run a simulation of protein near the graphene sheets in order to study stretching and folding of proteins.
I used these lines in my script to freez the graphene sheets.
; Number of walls, type, atom types, densities and box-z scale factor for Ewald
nwall = 2
wall-type = 12-6
pbc = xy
ewald-geometry = 3dc ; Only Ewald sum compatible with slab geometry,
wall-ewald-zfac = 3 ; Extra empty space for slab geometry calculation (reduc
wall-r-linpot = 0
wall_atomtype = MW MW
freezegrps = CR1
freezedim = Y Y Y
However, when I freeze the graphene sheets in 3 dimensions,I will receive the segmentation error.
I freezed them in two dimensions, so I can run it in NVT and consequently in NPT ensemble without error , but after 2ns of simulation in NPT, I found out some strange dispersion of graphene sheets Carbon atoms (like dash line).
I attached the photo of last snapshot and marked the dash line of atoms.
Do you have any idea what is wrong in this system?
Thank you in advance for your advice,