Dear all,
There are several parameters need to be set for ellipse as follows:
GB(sigma_ee/sigma_ss,epsilon_ss/epsilon_ee,upsilon,mu) (3,5,3,1)1
shape 1 3 1 1
pair_style gayberne 1.0 3.0 1.0 4
pair_coeff 1 1 1 1 1 1 0.2 1 1 0.2 4
This time the parameter for GB(sigma_ee/sigma_ss,epsilon_ss/epsilon_ee,upsilon,mu) (3,5,3,1) will be needed. But I am totally confused with the introduction from Lammps manual. It is the relationship between “shape” and “pair_coeff”.
When a longest axis is set to x-axis as the “shape” command set, it means the original ellipse orients to x-axis. Then the “end to end “ parameter would be set in the x-axis. But whether “pair_coeff” know which direction is the orientation? Different potential curves were gotten with “shape 1 3 1 and “shape 1 even if the “pair_coeff 1 1 1 1 1 1 0.2 1 1 0.2 were same. So it looks like that “pair_coeff” does not know which direction is the orientation. The curve from“shape 1 would be more similar to article report however the value was just half to that in the article.
Then how should I set these parameters? Please help me out! Thanks a million.