My problem is setting the parameters in pair_style and pair_coeff for ellipse simulation.
Usual format GB( see/sss, eee/ess, n,m) is needed to be changed to the format in Lammps.
If I want to build an ellipse bulks with three radiuses a=b<c, then I set the following command:
shape 1 1 1 2 #name a b c
pair_style gayberne 1.0 3.0 1.0 4.0 # type g n m cutoff
pair_coeff 1 1 1.0 1.0 1 1 0.2 0 0 0 2.5 # name_i name_j e s sia sib sic sja sjb sjc cutoff
I know the setting is very important. However I am not sure about it. Could you just check it in general idea? whether the n,m inGB( see/sss, eee/ess, n,m) is equal to that used in lammps?
Thanks a million.