[lammps-users] setting the parameters in pair_style and pair_coeff for ellipse simulation

Dear all,

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.

I suggest you read the LAMMPS doc page, or look in
pair_gayberne.cpp if it's unclear. Mike Brown (wmbrown at sandia.gov)
did the implementation so he's the expert if you're still
confused.

Steve