к.т.н.
Таджикский технический университет имени академика М.С.Осими
Ассистент
Нематов Дилшод Давлатшоевич,Таджикский технический университет имени академика М.С. Осими; Холмуродов Холмирзо Тагойкулович, Государственный университет «Дубна», г. Дубна, Московская область, Россия
УДК 538.9:538.94
Введение
Вычисление в рамках ТФП, являются одним из самых мощных инструментов исследования, результаты которых совпадают с экспериментом с высокой точностью. В настоящей работе в рамках ТФП [1-2] с использованием пакета WIEN2k [3] нами проведены квантово–механические расчеты с целью исследования электронной структуры наночастиц и тонких пленок ZrO2:Ti, представлены результаты по определению ширины запрещенной зоны и плотность электронного состояния (DOS; densityofstate) - величина, определяющая количество энергетических уровней в интервале энергий на единицу объёма в трёхмерном случае и на единицу площади — в двумерном случае).
Научная новизна
- в рамках ТФП проведены квантово-химические расчёты над ZrO2:Ti, исследованы их энергетические, геометрические и структурные свойства;
- проанализирован вклад энергетических уровней отдельных атомов плотность состояния в элементарных ячейках ZrO2.
- с помощью вычисления в рамках ТФП определена запрещенная зона ZrO2:Ti.
Материалы и методы
Создание оптимального материала требует усиления или подавления каких-либо свойств в существующем материале и невозможно без использования современных экспериментальных и теоретических подходов. Чтобы получить желаемый результат, мы должны построить полезную вычислительную модель, которая тщательно описывает электронные свойства материалов. Однако при этом необходимо изучить атомарную структуру веществ для получении материалов с заданными свойствами. В результате предоставляется возможность контролирования и надстройки свойств таких материалов.
В рамках ТФП можно исследовать структурные, и электронные свойства материалов. Метод ТФП основывается на предположение о том, что важнейшие свойства системы взаимодействующих частиц могут быть выражены с помощью функционала электронной плотности. ТФП-расчеты выполняются на основе данных о кристаллической решетки вещества в виде .сif-файлах. Данные о кристаллической решетке ZrO2 были взяты из базы данных сайта Materials Project [4].
Структурные характеристики ZrO2
На основе заданной геометрии нанокристалла ZrO2, соответственно их вышеприведенных постоянных и пространственной группы, нами выполнены оптимизации геометрии решеток, вычислены координаты всех атомов с учетом их симметрии в элементарной ячейке. Расстояния между атомами в кристаллической решётке представлены на табл. 1.
Связь | Растояние (Å) |
Zr-Zr | 3.56 |
O-O | 2.46 |
Zr-O | 2.23 |
Следуя методикам расчета по WIEN2k и ТФП мы имеем атомы, окруженные сферой Маттина (MuffТin; MT), где общая площадь делится на две зоны. Одна область включает внутреннюю сферу, а вторая часть - пространства имеющая форму промежуточной области. Радиусы сфер МТ равны 2.20 a0 для Zr и 1.99 a0 для O, где a0 - радиус Бора. Рис. 2 демонстрирует разные виды и положений кристаллической решетки ZrO2.
Рис. 2. Трехмерные структуры кристаллической решетки тонких пленок ZrO2
Этих структуры были визуализированы с помощью утилита WESTA с использование данные кристаллическое решетке.
Результаты и обсуждения
Согласно ТФП потенциал обменной корреляции аппроксимируется функционалом электронной плотности. Трёхмерные формы плотности электронов для одного и того же сорта атомов в элементарной ячейке сильно отличаются друг от друга из-за того, что они находятся в разной химической связи и их координатные точки являются разными.
На рис. 3 и 4 соответственно для выбранных плоскостей представлены карты плотности заряда в двумерном виде и электронное распределение атомов верхнего слоя ZrO2.
Рис. 3. Плотность электронов в элементарной ячейке ZrO2
Рис. 4. Распределение электроннов верхнего слоя ZrO2
В элементарной ячейке ZrO2 атомы расположены близко друг к другу, поэтому на энергетические уровни внешних орбитальных электронов также оказывают влияние электроны с внешних энергетических уровней соседнего атома. Электроны из самых внешних орбит атома могут испытывать силу притяжения от ближайшего соседнего атома. В этой связи электронные энергии не будут находиться на одном уровне; энергетические уровни электронов будут меняться до значения, которое может быть выше или ниже исходного энергетического уровня электрона. Для характеристики этих состояний в квантовой физике используется параметр энергетической зоны. Рассчитанная диаграмма энергетических зон ZrO2 в первой зоне Бриллюэна (ЗБ) с указанием точек (W, L, Λ, Г, Δ, X, Z, …, K) и направлений высокой симметрии показана на рис. 5. За счёт симметрии зоны Бриллюэна и возникает запрещенная энергетическая зона. Распределение энергетических зон в системе ZrO2 представлено на рис. 5.
Рис. 5. Диаграмма распределения энергетических зон ZrO2
Плотность состояния в энергетической зоне равна числу разрешённых состояний электронов на каждый интервал энергии. На рисунке 6 показана общая плотность электронных состояний для ZrO2. Пунктирные линии соответствуют уровням Ферми.
Рис. 6. Общая плотность электронных состояний ZrO2
На рис. 7 приведенырезультаты расчетов спин-орбитальных и спин-поляризованных связей, построены графики общей плотности электронных состояний верхней и нижней спинов.
Рисунок 8 показывает образования энергетических состояний за счет отдельных атомов в элементарной ячейке ZrO2.
Рис. 8. Плотность состояний отдельных атомов в элементарной ячейке ZrO2
Как показывает рисунок 7, в образование энергетических зоны имеется вклад широкого набора свободных электронных состояний всех сортов атомов кристалла ZrO2. Ширина запрещенной зоны в рамках расчётов GGA равна 3,3 эВ. Это может служить основанием полагать, что при благоприятных условиях (например, с повышением температуры) может быть осуществлён заброс электронов в эту зону из нижележащих энергетических зон. В этом случае, кристалл может приобрести свойство, характерное для широкозонных полупроводников. В недавней работе [5], результаты исследования электронной структуры и упругих свойств ZrO2 на основе расчетов из первых принципов весьма близки к нашим.
Далее нами было изучено электронных свойств ZrO2:Ti. Приведенные ниже результаты демонстрируют эффект допирования атома Ti (замещение чужих ионов в исходную кристаллическую решетку) для тонкой пленки ZrO2. Рис. 8 и 9 демонстрируют изменения геометрии и энергетических состояний нанокристалла ZrO2 при замещение одного атома Zr на атома Ti.