Метою даної роботи було підвищення стійкості симуляції на основі методу Нумерова при розрахунку хвильових функцій електронів в металевому кластері та забезпечення стабільності самоузгоджених розрахунків енергетичних характеристик металевих кластерів, шляхом обмеження зміни електростатичного потенціалу. Об'єктом дослідження є методика розрахунку хвильових функцій електронів методом Нумерова та знаходження власних значень енергетичних рівнів методом пристрілки. На етапі моделювання, з метою підвищення ефективності моделі було використано метод функціоналу
густини (DFT) у поєднанні з формалізмом Кона-Шема, що дозволило спростити складну задачу про
взаємодію електронів в полі заряджених іонів та отримати модель, в якій незалежні електрони рухаються в деякому ефективному потенціалі. Для розрахунку обмінних і кореляційних енергій ми використовували модель стабільного желе (SJM) та наближення локальної густини (LDA). Зміна обмеження профілю електростатичного потенціалу здійснювалася шляхом введення коефіцієнтів, що визначають внески попереднього та поточного профілів електростатичного потенціалу до результуючого. Розроблено модель металевого кластера із центрованою моновакансією та методику для розрахунку його параметрів. Для забезпечення збіжності та стійкості методу Нумерова запропоновано методику двостороннього обчислення з "зшиванням" хвильової функції в емпірично обраній точці. Розроблено методику симуляції з оптимальним кроком, яку реалізовано в програмному коді для розрахунку енергетичних характеристик металевих нанокластерів, що містять моновакансії. Моделювання дозволило отримати профілі електронної густини та ефективного потенціалу для заряджених та нейтральних кластерів. Результати розрахунків порівнювалися з даними експерименту, а також із ab initio обчисленнями. Розроблені підходи та методики моделювання можуть бути рекомендовані для аналізу низькорозмірних металевих систем, систем із шаруватою структурою в тому числі.The goal of this work was to increase the stability of the simulation based on the Numerov method during the calculation of the electron wave functions in a metal cluster and to ensure the stability of selfconsistent calculations of the energy characteristics of metal clusters by limitation of the electrostatic potential changes. The object of the study is the method for calculation of the electron wave functions and energy eigenvalues by the Numerov and Shooting methods. At the modeling stage, in order to increase the efficiency of the model, the density functional theory (DFT) was used in conjunction with the Kohn-Shem formalism, which allowed us to simplify the complex problem of the electron interaction in the field of charged ions and to obtain the model of independent electrons moving in some effective potential. We used the stabilized jellium model (SJM) and local density approximation (LDA) for the exchange and correlation energies. The change of the electrostatic potential profile limitation is used by adding coefficients that determine the contributions of the previous and current electrostatic potential profile to the resultant one. Models of a metal cluster with a centered monovacancy and approaches for calculating its parameters were developed. To ensure the convergence and stability of the Numerov method, methods of two-side calculation with "cross-linking" of the wave function at the empirically selected point were proposed. The method of simulation with the optimal step is developed and implemented in the program code for calculating the energy characteristics of metal nanoclusters containing the monovacancy. Simulation allowed to obtain the electron density and effective potential profiles for charged and neutral clusters. The results of the calculations were compared with experimental data, as well as with ab initio computations. Developed approaches and simulation techniques can be recommended for the analysis of low-dimensional metal systems, including systems with a layered structure
