Sorry for sending out the incomplete one by mis-pressed the send button.

Dear Lammps users,

I am writing to make sure that my understanding in pair_coeff is correct.

My confusion lies in here:

For better understanding,I listed these values in a table form. (The information below is taken from Lammps default input file of ellipse)

type i type j epsilon sigma epsilon_i_a epsilon_i_b epsilon__c epsilon_J_a epsilon_J_b epsilon_J_c

pair_coeff 1 1 3.0 1.0 1 1 1 1 1 1

pair_coeff 1 2 3.0 1.0 1 1 1 0 0 0

pair_coeff 2 2 1.0 1.0 1 1 0. 2 0 0 0

According to the Lammps description of values of these epsilons, where it mentioned ’ **the typical way to define the epsilon_i and epsilon_j coefficients is to list their values in “pair_coeff I J” commands when I = J, but set them to 0.0 when I != J**.', and then here comes my question.

If say pair_coeff for type 1 to 1 is supposed to be listed (b/c list values when I=J),especially epsilon_J_(a,b,c), then why for type 2 to 2 (I=J) the pair_coeffs are not listed for epsilon_J(a,b,c)? According to my understanding after reading the manual, the epsilon J_a, epsilon J_b, epsilon J_c should be assigned values as well.

Can someone kindly explain to me how are these values supposed to be assigned?

Thanks in advance for any ideas!

Sincerely,

LQ