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