В рамках модели стабильного желе выполнены самосогласованные вычисления пространственных распределений электронов, потенциалов, энергий диссоциации, когезии, образования вакансии, прилипания электронов и потенциала ионизации сплошных кластеров AlN, NaN (N ≤ 254) и кластеров, содержащих вакансию (N ≥ 12). Обсужден вклад моновакансии в энергетику кластера, размерные зависимости характе-ристик и их асимптотик. Расчет выполнeн на кластере СКИТ-3 Института Кибернетики им. В.М. Глушкова НАН Украины (Rpeak = 7.4 Tflops). DOI: 10.21883/JTF.2018.02.45405.2258 Введение Одна из точек зрения на плавление твердых тел [1] предполагает скачкообразное увеличение концентрации дырок-вакансий в тройной точке, а также уменьшение энергии образования вакансий с ростом их концентра-ции. В точке плавления концентрация вакансий в метал-лах составляет доли процента. Несмотря на столь малые концентрации, вакансии оказывают большое влияние на свойства твердых тел.Равновесная концентрация вакансий вычисляется при наличии энергии образования вакансии, которая может быть извлечена из спектра аннигиляции, локализован-ных в вакансии позитронов. При такой процедуре кон-центрация вакансий остается свободным параметром.Первоначально в экспериментах установлено, что температура плавления кластеров на подложке и сво-бодных кластеров уменьшается с уменьшением их раз-меров. Интерпретации и моделированию этого мезо-скопического явления посвящено значительное число работ [2][3][4][5][6][7].Современные масс-спектрометрические и калоримет-рические методы, позволяющие детально исследовать процесс предплавления и плавления металлических кла-стеров, состоящих из счетного числа атомов [8][9][10][11], показали, что температура плавления характеризуется осцилляционной размерной зависимостью, а также име-ет размерные аномалии (например, для Al), плохо опи-сываемые в простых моделях [8,10]. Также в процессе плавления диффузия поверхностных вакансий в объем более выгодна для кластеров с незаполненными элек-тронными оболочками, чем для кластеров с магическим числом атомов [9]. Эти факты стимулируют повышенный интерес к рассмотрению фазового перехода из твердого в жидкое состояние, как к конфигурационному возбуж-дению дырок-вакансий в кластерах.Популярна точка зрения [12-14] с позиции термоди-намики: вблизи температуры плавления энергия обра-зования вакансий тем меньше, чем меньше кластер, а концентрация вакансий не зависит от его размера. Од-нако термодинамика не дает указаний на механизм об-разования вакансий, и вопрос о размерной зависимости энергии образования вакансий, их концентрации и связи вакансий с процессом плавления остается открытым.Несмотря на то что энергетические характеристи-ки сплошных металлических кластеров рассчитывались неоднократно в различных моделях ([15-17] и ссылки в них), самосогласованные вычисления энергии обра-зования вакансии в кластерах и влияние на нее кван-тования электронного спектра до сих пор не выполня-лись. Помимо этой задачи, одной из актуальных задач, которые можно сформулировать в связи с плавлением малоразмерных агрегатов,...