Ранее была предложена двумерная статистическая модель деформированных тонких пленок с дефектами типа вакансий. Для однородного плосконапряженного состояния решетки составлена система двух интегральных уравнений для потенциалов средних сил, которая решена аналитически в приближении Гаусса. Потенциалы средних сил определяют младшие коррелятивные функции распределения частиц в окрестности узлов плоской деформированной решетки. Знание функций распределения необходимо для статистического расчета термодинамических потенциалов молекулярной системы. Благодаря этому появилась возможность получить аналитическое выражение для свободной энергии деформированных кристаллических молекулярных пленок с дефектами.