Разработаны основы статистической теории упругости, которая развивается в рамках двухуровневого молекулярно-статистического описания неоднородных систем. Одновременно используются метод корреляционных функций Боголюбова-Борна-Грина-Кирквуда-Ивона (ББГКИ), метод условных распределений Л. А. Ротта и метод термодинамических функционалов. В результате получено статистическое выражение для функционала свободной энергии, зависящее от поля микроскопического тензора деформации кристаллической решетки с вакансиями и поля чисел заполнения ячеек решетки (узлов). Это выражение учитывает, помимо потенциальной энергии деформации (силового фактора), структурные особенности кристалла с вакансиями (энтропийный фактор).