Dear all,
I’m trying to perform uniaxial test for GO using reaxFF potential in lammps.
I devided my sheet into 3 group and want to perform the force just on top.
After making sure I reached equilibrium I used these lines to perform deformation:
fix mysf2 carbon_top setforce 0 0 0
fix mysf1 carbon_bot setforce 0 0 0
velocity carbon_top set 0 0.1 0
velocity carbon_bot set 0 0 0
but my problem is that I don’t see the atoms moving in the top.
can I ask for help?