Based on the Keating interatomic potential an analytical derivation of force constant matrices of two- and three body interactions in crystals is performed. Using the derive d force constant matrices the in-plane phonon ...