Hi everybody
I am wondering how to calculate torque when rotating a group,and I got a few questions.
1)What’s the difference between displace_atoms and fix move.
2) When I use displace_atoms + setforce 0.0 0.0 0.0 + velocity set 0.0 0.0 0.0, the graphene do rotate but the torque output by torque() is 0. I kwon the reason is that “setforce 0.0 0.0 0.0” would zero the torque and force, and some one said I should use fix rigid and turn off the torque and force. Then I used it but the graphene didn’t rotate, so I hope you may point out what’s wrong with my code.
Below is my code.
thanks a lot
Sydney
#LAMMPS input file for graphene
units metal
boundary p p p
atom_style atomic
read_data data.Graphene
pair_style airebo 0.0 0 0
pair_coeff * * CH.airebo C H
#pair_style tersoff
#pair_coeff * * BN_albe B N
#neighbor 0.2 bin
#neigh_modify every 20 delay 100 check yes
#mass 1 10.81100000
#mass 2 14.00670000
mass 1 12.0107
mass 2 1.00
group g_C type 1
group g_H type 2