Hi, Lammps users
I have some question about the material’s geometry after thermal relaxation with NPT.
I’m working on single layer graphene on Cu film.
I have used AIREBO potential and EAM potential for graphene and Cu respectively.
LJ potential was used for C-Cu interaction.
I have currently checked my script with different distance between single graphene layer and Cu top surface.
When I have 0.8nm for the distance and ran the relaxation, its geometry was fine.
However, its geometry was strange when I have 0.5nm for the distance.
It means that some Cu atoms were on top of graphene and graphene is in between two Cu layer as attached picture.
Also, graphene and Cu structure were bent at 0.5nm distance but there was not any bending at 0.8nm distance.
Here is my lammps scripts.
I’d appreciate it if you give me some comments.
This is the control script for LAMMPS