Dear Users,

In the documentation page of GayBerne potential there is a link to a pdf file containing further information about the formulation. I have found some discrepancies between the file and the citation it included (Evaraers), in the definition of parameters \chi, **B**_{12} and **E**_1 (**E**_2).

In the paper,

(1) **E**_1=diag(\epsilon_a1^(-1/\mu),\epsilon_b1^(-1/\mu),\epsilon_c1^(-1/\mu)),

(2) **B**_12=**A**_1 **E**_1 **A**_1^T+**A**_2 **E**_2 **A**_2^T, and

(3) \chi=[**r**_12 **B**_12 **r**_12^T]^\mu.

In lammps explanation,

(1) **E**_1=diag(\epsilon_a1,\epsilon_b1,\epsilon_c1),

(2) **B**_12=**A**_1 **E**_1^2 **A**_1^T+**A**_2 **E**_2^2 **A**_2^T, and

(3) \chi=[**r**_12 **B**_12 **r**_12^T]^\mu

In the former, the well depth is proportional to the parameters \epsilon_i, however in the latter, it is inversely proportional to some exponential of them. The latter formulation cannot be found anywhere else and is incompatible to what is being done in the source code according to my understanding. Could anyone please check it up?

Thanks a lot!

Regards,

Patrick Hejian Zhu