Микроэлектроника, 2020, T. 49, № 1, стр. 46-49

Математическое моделирование электронного строения, синтез и применение наночастиц IrSi

А. Г. Гасанов a, Э. А. Керимов a b*, С. Н. Мусаева c

a Военная академия вооруженных сил
AZ 1065 Баку, ул. Академика Ш. Мехтиева,136, Азербайджан

b Государственная академия нефти и промышленности
AZ 1010 Баку, Пр. Азадлыг, 16/21, Азербайджан

c Азербайджанский технический университет
1073 Баку, пр. Г. Джавида, 25, Азербайджан

* E-mail: E_Kerimov.fizik@mail.ru

Поступила в редакцию 24.04.2019
После доработки 08.07.2019
Принята к публикации 16.08.2019

Полный текст (PDF)

Аннотация

Построены теоретические модели наночастиц IrSi и одним из вариантов метода молекулярных орбиталов исследовано электронное строение наночастиц IrSi полуэмпирическим Расширенным Методом Хюккеля (РМХ). На основе теоретических моделей вычислены значения орбитальных энергий, потенциал ионизации, полной электронной энергии, модули Юнга и твердости наночастиц IrSi. Являясь электрофильными, полупроводниковыми и мягкими материалами они могут применяться в военном деле в будущем как фотоприемные приборы для различных целей.

Ключевые слова: математическое моделирование, квант-механические методы, фотоприемник, орбиталы, полупроводник

ВВЕДЕНИЕ

Диоды Шоттки – перспективный класс полупроводниковых приборов, которые ввиду более высокого быстродействия и малых падений напряжения при прямом смещении обладают преимуществами перед диодами, имеющими p-n переходы, и это делает их привлекательными в качестве элементной базы силовой электроники.

В качестве барьеров Шоттки интерес представляют силициды платиновой группы, имеющие большую высоту потенциального барьера к кремнию, что обуславливает низкие токи утечки и снижает температурную зависимость электрических параметров приборов.

Тонкие пленки силицида иридия получены твердофазной реакцией пленки металлического иридия с кремниевой подложкой. Пленки иридия толщиной 0.1 мкм наносились на подложки монокристаллического кремния (111). Нанесение пленок производилось резистивным испарением высокочистого иридия при давлении 6 × 10–4 Па и температуре 200°С на установке вакуумного напыления УВН РЭ.Э-60. Затем, 20 мин осуществлялся отжиг структур при температуре 450°С.

ТЕОРЕТИЧЕСКИЙ МЕТОД

Для изучения различных свойств наночастиц IrSi, огромную роль имеет применение математического моделирование [13] и исследование их квантмеханическими методами. Для этого сначала рассмотрим основные параметры наночастиц IrSi. Основными параметрами будем считать количество и размеры атомов.

Допустим, что наночастица IrSi имеет форму полого шара (рис. 1). В этом случае количество атомов определяется из формулы:

(1)
$n = \frac{{{{R}^{3}} - {{r}^{3}}}}{{r_{h}^{3}}}.$
Рис. 1.

Теоретические модели молекулы (a) IrSi и наночастицы (б) (IrSi)n.

В этой формуле r = R – 2rh.

Из формулы (1) вычислим количество атомов в наночастице (IrSi)n. Имея формулу шара размер одного соединения IrSi может определятся как rhrIr + rSi (рис. 1). Формуле rIr и rSi ковалентные радиусы соответственно атомов Ir и Si. Для (IrSi)nr = R – 2rh. R – радиус сферы внутри которой сдерживается наночастица. R = D/2 и D – размер наночастицы. А rh это радиус сферы внутри которой сдерживается соединение IrSi.

При rIr = 0.127 нм и rSi = 0.111 нм rh = 0.238 нм и при D = 1.2 нм R = 0.6 нм. В этом случае для наночастицы (IrSi)n получаем n ≈ 16. А общее количество всех атомов Ir и Si будет N = 32. В итоге в зависимости от n можно создать различные модели наночастиц (IrSi)n. Эти модели представлены на рис. 2.

Рис. 2.

Различные модели наночастиц IrSi.

На основе созданных моделей можно произвести теоретические вычисления и изучит свойства наночастиц. В частности такие вычисления производятся с помощью метода молекулярных орбиталей. Также известно, что Расширенный Метод Хюккеля (РМХ) простой полуэмпрический вариант метода молекулярных орбиталей (МО). В методе МО считается, что электроны в молекуле в эффективных полях созданных ядрами и электронами двигаются независимо от других электронов. Состояние электрона в молекуле описывается одноэлектронной волновой функцией называемой молекулярный орбитал [5]. Эти функции являются многоцентравими функциями. Так как в формулах таких функций присутствуют расстояние электрона от различных ядер атомов. Существуют различные варианты нахождения молекулярных орбиталей. Одним из таких вариантов является метод нахождения Ui молекулярных орбиталей в качестве линейных комбинаций атомных орбиталей атомов входящих в состав молекул (метод MOLCAO):

(2)
${{U}_{i}} = \sum\limits_{q = 1}^m {{{c}_{{qi}}}{{\chi }_{q}}} ,$

в формуле сqi – неизвестные коэффициенты, χq – атомные орбиталы, выбранные как базисные функции. В статье использованы функции Гаусса в качестве атомных орбиталей [5]. Неизвестные коэффициенты сqi находим из следующих систем уравнений:

(3)
$\sum\limits_q {\left( {{{H}_{{pq}}} - {{\varepsilon }_{i}}{{S}_{{pq}}}} \right)} {{c}_{{qi}}} = 0.$

В формуле (2) обозначены:

(4)
${{H}_{{pq}}} = \int {\chi _{p}^{*}} {{\hat {H}}_{{{\text{ef}}}}}{{\chi }_{q}}dV,$
(5)
${{S}_{{pq}}} = \int {\chi _{p}^{*}{{\chi }_{q}}dV} .$

Spq – интеграл перекрытие между атомными орбиталами χq и χp. ${{\hat {H}}_{{{\text{ef}}}}}$ – оператор Гамильтона для одного электрона движущимся на некотором эффективном поле независимо от других электронов:

(6)
${{\hat {H}}_{{{\text{ef}}}}} = - \frac{1}{2}{{\nabla }^{2}} + U(r).$

Так как U(r) не имеет четкого выражения, величины Hpq не могут вычисляться точно. Они оцениваются с помощью некоторых экспериментальных параметров. В Расширенном Методе Хюккеля диагональные элементы матрицы Hpq считается равными ионизационного потенциала соответствующих валентных состояний атомов. А не диагональные элементы матрицы Hpq вычисляются приближением [5]:

(7)
${{H}_{{pq}}} = 0.5K{{S}_{{pq}}}({{H}_{{pp}}} + {{H}_{{qq}}}).$

K – находим или из сравнения экспериментальных фактов или же из условия минимума энергии. На работе K = 1.75. Уравнения (2) система линейных однородных уравнений. Решая эти уравнения с общепризнанными правилами, находим численное значения для параметров εi и cqi. Используя численные значения εi, определяются электронная энергия и потенциал ионизации наночастиц IrSi и исследуются механические, электрические, магнитные и др. свойства этих наночастиц.

РАСЧЕТ И АНАЛИЗ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ

В предлагаемой работе значения орбитальных энергий, потенциала ионизации, полной электронной энергии наночастиц IrSi рассчитаны с помощью полуэмпирического метода РМХ. РМХ сохраняя простого варианта молекулярных орбиталей, метода Хюккеля, предполагает вычисления интегралов перекрытия и одновременно рассмотрит π и σ – электронов.

Для каждого объекта начиная с самого низшего уровня энергии, электроны находятся на уровнях которые расположены по две. Определена энергия наивысшего молекулярного орбитала в котором расположены электроны как – εHOMO и энергия низшего пустого молекулярного орбитала как – εLUMO. Потенциал ионизации определен из формулы: Ip = – εHOMO, значение запрещенной зоны Eg = εLUMO – εHOMO, а твердость с помощью формулы η = ½Eg [14]. Длина волны фотона, которого излучает материал может вычислена из формулы

$\lambda = \frac{{ch}}{{1.6{{E}_{g}}}} \times {{10}^{{28}}}\,\,{\text{нм}}{\text{.}}$
Здесь h – постоянная Планка, с – скорость света в вакууме. При вычислении λ пользовались значениями эВ для Eg. При η < 1 эВ материал считается мягким, а при η > 1 эВ он твердый.

Стабильность материала получаем из формулы

$\Delta E = E - \sum\limits_A {{{E}_{A}}} .$
Здесь Е – полная энергия системы, ЕА – энергия А – числа атомов входящих в систему, ΔЕ – параметр характеризирующий стабильность системы. При ΔЕ > 0 система считается нестабильным, а при ΔЕ < 0 она стабильна. Полученные результаты даны в табл. 1.

Таблица 1.  

Рассчитанные значения энергетических параметров для наночастиц (IrSi)n

N Объект εHOMO εLUMO Полная энергия
E (a.е.)
Параметр стабильности
ΔE (a.е.)
Потенциал ионизации
Ip (эВ)
Запрещенная зона Eg (эВ) Параметр твердости
η (эВ)
Длина волны, λ (нм)
1 IrSi –11.600
068
–9.0912 –6.28111
9798
–0.30825
4595
11.600 2.508868 1.254 495.493
2 (IrSi)2 –11.842
59
–10.367
675
–12.7475
3132
–0.80180
0912
11.843 1.474915 0.737 842.845
3 (IrSi)3 –10.979
948
–10.458
483
–19.0221
7159
–1.10357
5985
10.980 0.521465 0.261 2383.909
4 (IrSi)4 –11.492
561
–10.769
836
–25.4492
9043
–1.55782
9613
11.493 0.722725 0.361 1720.053
5 (IrSi)5 –10.746
609
–10.601
305
–31.7904
8026
–1.92615
4248
10.747 0.145304 0.073 8555.339

В табл. 1 по расчетным значениям стабилизированный энергии можно вычислит модуль Юнга и другие механические параметры: G – модуль скольжения, K – объемно эластический модуль и HN – модуль твердости для наночастиц (IrSi)5 [6, 7]. Для этого можно использовать следующую формулу:

(8)
$Y = \frac{F}{S},\,\,\,\,F \approx ~\,\,\frac{{\Delta E}}{r},\,\,\,\,G = \frac{Y}{{2 \times \left( {1 + \nu } \right)}},\,\,\,\,~K = \frac{Y}{{3 \times \left( {1 - 2\nu } \right)}}~,\,\,\,~HN = GA{{e}^{{ - ВT}}}.$

В формуле (8):

F = 2.866 × 10–8 Н сила распада наночастицы;

S = 2.176195715539529 × 10–9 м2 поперечная площадь для (IrSi)5;

r = 2.926 × 10–10 м длина связи между атомами входящими в наночастицу;

– ν = 0.225 коэффициент Пуассона, A = 0.807, B = 2.204 × 10–3 и T = 300 ${\text{K}}$.

В табл. 2 показаны результаты расчетов.

Таблица 2.  

Расчетные значения механических параметров для (IrSi)5

Модуль скольжения (ГПа) Объемно эластический модуль (ГПа) Модуль Юнга (ГПа) Модуль твердости (ГПа)
1 53.76 79.83 131.7 22.4

ЗАКЛЮЧЕНИЕ

Одним из вариантов метода молекулярных орбиталов исследовано электронное строение наночастиц IrSi полуэмпирическим Расширенным Методом Хюккеля . На основе компьютерных расчетов вычислены орбитальные энергии, потенциал ионизации, полной электронной энергии, модули Юнга и твердости. Выводы показывают, что эти материалы могут быть использованы в будущем в военном промышленности в качестве фотоприемных устройств.

Список литературы

  1. Гасанов А.Г., Байрамов А.А. Моделирование электронной структуры композитного материала графен-поливинилиденфторид // Физика твердого тела. 2019. Т. 61. Вып. 1s. 194–200.

  2. Gasanov A.G., Bayramov A.A.,Hashimov E.G. Mathematical modeling of the electron structure of SiO2 nanoparticle. FİZİKA Volume XXIII, Number 1, Section: En, April. 2017, s. 34–39.

  3. Pashaev F.G., Gasanov A.G., Mahmood A.T. The Study of Gold Nanoparticles in basis of Slater Functions // J. Nano. Adv. Mat. 2014. V. 2. № 1. P. 35–41.

  4. Gasanov A.G., Pashaev F.G. The Computer Program for the Study of Nanoparticles in Basis of Slater Atomic Orbitals. Romanian // Journal of Information Science And Technology. 2016. V. 19. № 4. P. 331–337.

  5. Minkin V.I., Simkin B.Y., Minyaev R.M. Theory ofstructure of molecule. Rostov at Don: Phoenix. 1997, 558 p.

  6. Черняева Т.П., Грицина В.М., Михайлов Е.А., Остапов А.В. Корреляция между упругостью и другими свойствами циркония. Научно-тех. комп. “Ядерный топливный цикл” Нац. науч. Центра “Харьковский физ.-тех. институт”, Харьков, Украина, 2008.

  7. Опенов Л.А., Подливаев А.И. // ФТТ. 2017. Т. 59. Вып. 6.

Дополнительные материалы отсутствуют.