Dear LAMMPS users,
Just out of curiosity has anyone implemented or made use of the SiC potential reported in:
P. Vashishta, R. K. Kalia, A. Nakano, and J. P. Rino, “Interaction potential for silicon carbide: A molecular dynamics study of elastic constants and vibrational density of states for crystalline and amorphous silicon carbide,” Journal of Applied Physics, vol. 101, no. 10, pp. 103515–103515–12, May 2007.
I wanted to check before I attempt to do so. Any help or suggestions would be great!
Not that I am aware of.
[[email protected] potentials] grep \-i vash \*
\[[email protected]\.\.\.2929\.\.\. potentials\]
From what I recall, that SiC potential is similar to Stillinger-Weber,
but with a different two-body term. So, you may be able to create it
by combining Stillinger-Weber with one or more two-body potentials
using pair_style hybrid/overlay. More elegantly, you could create a
new pair_style that is a C++ child class of PairSW.
Thank you for the quick response. Yes , the potential is essentially a two-body pairwise term with a modified Stillinger-Weber term.
Both the two-body and SW potentials are not conventional so hybrid/overlay won’t work. It looks like (as you stated) the best way to do this is by creating a new pair_style.