На основе первопринципного метода функционала плотности выполнена оптимизация геометрии кристаллов BeXAs₂ (X = Si, Ge, Sn), определены геометрические параметры примитивных ячеек и вычислены фононные спектры в центре зоны Бриллюэна. Интерполирование фононных спектров по всей зоне Бриллюэна приведено с использованием феноменологической модели Китинга. Частоты, вычисленные из первых принципов, использовались при нахождении параметров Китинга путем программной оптимизации. Численные значения феноменологических параметров находятся в соответствии с вычисленными из первых принципов характеристиками химической связи. Исследован парциальный состав фононных спектров. Качественно структура акустических и оптических ветвей в основном определяется соотношением масс ионов исследуемых соединений.