Численные методы моделирования свойств нанокристаллов. Нагорнов Ю.С. - 43 стр.

UptoLike

Составители: 

43
Здесь использовались параметры оболочечной модели и
потенциалы (2), приведенные в работе [13]. Константы A
ij
и B
ij
этих
потенциалов восстанавливались совместно с параметрами оболочечной
модели по диэлектрическим и упругим свойствам кристалла UO
2
, тогда как
константы дисперсионного притяжения рассчитывались по формуле
Слетера-Кирквуда (3).
Ядро (
C
), заряд
-
Q
C
S
), заряд
-
Q
S
Жесткость связи
-
K
Заряд иона:
Q = Q
C
+ Q
S
.
Поляризуемость
иона:
a
= (Q
S
)
2
/K.
В поле единичной
напряженности:
X = Q
S
/K.
поле
поле
X
Энергия поляризации иона:
2
K
×
X
2
U
Pol.
=
Рис. 2.2 Оболочечная модель иона
Таблица 2.1
Параметры модели кристалла UO
2
[13]
Ион Заряд
иона Q, e
Заряд оболочки
Q
S
, e
Жесткость иона
k, эВ×Å
-2
U
4+
+4 6.54 98.23
O
2
-
-2 -4.40 296.80
Потенциал A, эВ В, Å
-
1
С
6
, эВ× Å
6
U
4+
-U
4+
18600 3.6417 32.64
U
4+
-O
2
-
2494.2 2.9306 40.16
O
2
-
-O
2
-
108.0 2.6316 56.06
1.4. Граничные условия
1.4.1. Метод Мотта-Литтлетона
Размер модельного кристалла ограничен эффективностью
используемого алгоритма вычисления сил и быстродействием ЭВМ. На
границах моделируемого объема должны выполняться условия,
связывающие его с бесконечным кристаллом.
Многими авторами (см., например, [14]) для постановки граничных
условий при моделировании ионных кристаллов в описанных выше
приближениях используется метод Мотта-Литтлетона (Mott-Littleton