Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.
Citation: P. Vashishta, R.K. Kalia, A. Nakano, and J.P. Rino (2007), "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, 101(10), 103515. DOI: 10.1063/1.2724570.
Abstract: An effective interatomic interaction potential for SiC is proposed. The potential consists of two-body and three-body covalent interactions. The two-body potential includes steric repulsions due to atomic sizes, Coulomb interactions resulting from charge transfer between atoms, charge-induced dipole-interactions due to the electronic polarizability of ions, and induced dipole-dipole (van der Waals) interactions. The covalent characters of the Si–C–Si and C–Si–C bonds are described by the three-body potential. The proposed three-body interaction potential is a modification of the Stillinger-Weber form proposed to describe Si. Using the molecular dynamics method, the interaction potential is used to study structural, elastic, and dynamical properties of crystalline (3C), amorphous, and liquid states of SiC for several densities and temperatures. The structural energy for cubic (3C) structure has the lowest energy, followed by the wurtzite (2H) and rock-salt (RS) structures. The pressure for the structural transformation from 3C-to-RS from the common tangent is found to be 90 GPa. For 3C-SiC, our computed elastic constants (C11, C12, and C44), melting temperature, vibrational density-of-states, and specific heat agree well with the experiments. Predictions are made for the elastic constant as a function of density for the crystalline and amorphous phase. Structural correlations, such as pair distribution function and neutron and x-ray static structure factors are calculated for the amorphous and liquid state.