CC BY
22
УДК 539.19+539.2
DOI: 10.18698/1812-3368-2017-4-75-89
ХАРАКТЕРИСТИКИ ПОВЕРХНОСТИ ЩЕЛОЧНЫХ МЕТАЛЛОВ С УЧЕТОМ ДИСКРЕТНОСТИ КРИСТАЛЛИЧЕСКОЙ РЕШЕТКИ И ФРИДЕЛЕВСКИХ ОСЦИЛЛЯЦИЙ ЭЛЕКТРОННОЙ ПЛОТНОСТИ
В.Л. Глушков О.С. Еркович
glushkov-vl@bmstu.ru erkovic@bmstu.ru
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация
Аннотация
Проведен расчет влияния дискретности кристаллической решетки на характеристики поверхности щелочных металлов в рамках метода функционалов плотности. Электронная плотность в приповерхностной области смоделирована пробными функциями распределения, которые учитывают фриделевские осцилляции. Проведен самосогласованный расчет поверхностной энергии с учетом градиентных поправок на неоднородность электронного газа для кинетической энергии. Определены поверхностные характеристики — работа выхода и высота потенциального барьера. Построены эффективные потенциалы для выбранных металлов. Учет дискретности кристаллической решетки выполнен в рамках теории возмущений. Реальные потенциалы ионов аппроксимированы псевдопотенциалом Ашкрофта. Для анализа учета дискретности все характеристики поверхности металла получены в двух вариантах: с учетом и без учета вклада кристаллической решетки
Ключевые слова
Метод функционалов плотности, электронная плотность, дискретность кристаллической решетки, работа вьхода, фриделевские осцилляции
Поступила в редакцию 26.04.2016 © МГТУ им. Н.Э. Баумана, 2017
Введение. Исследование характеристик поверхности металлов является актуальным научным направлением, которое имеет большое практическое значение при рассмотрении вопросов эмиссионной электроники, вакуумной технологии, оптики и др. Для понимания процессов и явлений, происходящих на поверхности металлов, необходима информация о таких характеристиках, как распределение электронной плотности вблизи поверхности, поверхностная энергия единицы площади поверхности материала, работа выхода электрона и пр. Трудность расчета перечисленных характеристик заключается в том, что электронная плотность вблизи поверхности металла уменьшается от объемного значения до нуля на расстояниях порядка атомных размеров, что создает серьезные трудности при сопоставлении экспериментальных и теоретических результатов. Одним из перспективных направлений в описании сильно неоднородных электронных систем является самосогласованная теория неоднородного электронного газа — метод функционалов плотности (МФП).
Цель настоящей работы — нахождение характеристик поверхности металла (работы выхода, потенциального барьера, поверхностной энергии) с использованием МФП в рамках модели однородного фона для поверхности металла и расчет характеристик поверхности с учетом дискретности кристаллической решетки в рамках теории возмущений, выполненный для разных типов пробных функций распределения электронной плотности.
Использованный в работе подход позволяет представить основные характеристики многочастичных систем в виде однозначных функционалов плотности частиц. Полная энергия основного состояния многочастичной нерелятивистской квантовой системы с фиксированным числом частиц N представляет собой однозначный функционал функции плотности nm (r1, ..., rm), т < N, минимум которого реализуется на функции Пто (ri, ..., r„), соответствующей распределению частиц в основном состоянии системы [1-5]. Использование такого подхода позволяет определить основные характеристики поверхности металлов, не прибегая к поиску ее точной волновой функции, что упрощает применение численных методов. Кроме того, плотность частиц ni(r) можно наблюдать экспериментально [6-8].
В качестве пробных функций распределения были выбраны функции, предложенные в работе [9]. При этом критерием корректности той или иной модели описания поверхности и прогнозируемых ею результатов может служить сопоставление рассчитанных и измеренных значений различных наблюдаемых характеристик металлов.
Методика расчета поверхностных характеристик и электронной плотности у поверхности металла. Для изучения поверхностных характеристик вблизи поверхности металлов часто используют упрощенную модель, в которой решетка заменяется однородным положительным фоном (модель однородного фона или модель «желе» [10, 11]).
В этой модели плотность положительного заряда для кристалла, занимающего полупространство z < 0, имеет вид
Uion (Z) = no0("Z ), (1)
/ ч [1, z > 0;
где ö (z )=|q q — функция Хевисайда.
Ввиду одномерности плотности nion ( z ) электронная плотность также зависит только от координаты z. Кроме того, она должна удовлетворять условиям
lim n(z) = nion (z), (2)
|z
и для электронейтрального кристалла имеет вид
j [n(z)-nion (z)]dz = 0. (3)
—ад
Определение электронного распределения вблизи поверхности металла может быть получено путем решения вариационной задачи о нахождении ми-
нимума функционала E[n] полной энергии неоднородной системы электронов в основном состоянии [11, 12]. При этом в рассмотрение вводятся два семейства пробных функций: n = n(a;, z) — электронная плотность и ф = ф(а;, z) — электростатический потенциал. Эти семейства функций связаны между собой уравнением Пуассона, обеспечивающим самосогласование для кулоновского потенциала. Выбранные пробные функции должны удовлетворять условиям (2) и (1), быть непрерывными на границе металл-среда, приводя к условию непрерывности нормальной составляющей электростатической индукции.
В качестве пробной функции электронной плотности вблизи поверхности раздела металл-вакуум часто используют функцию, представленную без обсуждения и обоснования в работе [11] Смитом
По (г) = По
1 -1 еРг je (-г )+1 е"Рг 0(г )
(4)
где По — объемная электронная плотность; р — вариационный параметр.
Полная энергия электронного газа полуограниченной системы бесконечна, поэтому для поиска вариационного параметра рассматривают эквивалентную задачу о минимизации функционала поверхностной энергии системы. Поверхностная энергия в рамках модели «желе» рассчитывается как разность полной энергии системы, когда электроны распределены в соответствии с функцией (4), и полной энергии системы, когда их распределение задано функцией (1):
ао [п (г )] = <СТ [П (г )] + ахс [П (г )] + аи [П (г)]. (5)
Первое слагаемое в уравнении (5) представляет собой кинетическую часть поверхностной энергии электронного газа:
ак [П (г)] = СТ£ [П(г)] + ^ [П(г)]
ш 3 2 5
где [п (г)]= | —(3л2 )- п (г)з йг — плотность кинетической энергии элек-
—ш 10
тронного газа плотностью п в приближении Томаса — Ферми [11-14] (для упрощения вида формул использована атомная система единиц, е = те = Ь =1); ш УП(г)2
СТ (п )= |-т^г йг — поправка на неоднородность электронного газа для
—ш 72П (г)
плотности кинетической энергии (формула Вейцзеккера — Киржница [15, 16]).
Второе слагаемое в уравнении (5) является обменно-корреляционным вкладом в поверхностную энергию [17-20]:
стхс (п) = — { -Г—3 П(г)з йг — { 0,056п(г)з1 йг. —ш 4 ^ —ш 0,079 + п(г)з
Третье слагаемое в уравнении (5) представляет собой вклад электростатической энергии
ад 1
ае5 (и )= { -ф(г )и (г )ёг,
-ад 2
где ф(г) — потенциал, создаваемый распределением электронов и фоном положительного заряда (возможен учет потенциала внешнего поля).
Проводя численное интегрирование, получаем значение вариационного параметра Р и значения характеристик поверхности металла, которые могут быть сопоставлены с результатами эксперимента [6, 11]. В настоящей работе подобные расчеты были выполнены для следующих пробных функций плотности:
У „ „ Л
« (г ) = по П2 (г) = по
1 -
1 1
2 (1 - г )
е(-г) +
1 1
2 (1 + г)а
-е(г)
1—
a2L
(b2LZ -1)2
е(-г ) + (1 - a2L )e ~a2RZ е(г)
«з (г ) = по
1 -
a3L cos (bзLг - C3L)
(bзLг - 1)2COS (C3L )
е(-г ) + (1 - a3L )e-аз*г е(г)
(6)
(7)
(8)
Здесь й1, в21, Язь, Сзь — вариационные коэффициенты; коэффициенты а2к, азя находятся из условия электронейтральности (3) электронной плотности на границе раздела металл-вакуум; коэффициент Ь2Ь = Ьзь = 2к^ [16].
Как уже было отмечено выше, обоснование подобного выбора пробных функций представлено в работе [9]. Значения поверхностных характеристик, таких как поверхностная энергия, работа выхода и высота потенциального барьера, рассчитанных для пробных функций плотности (6)-(8), будут сопоставляться со значениями тех же характеристик, полученных для функции (4).
Учет дискретности кристаллической решетки проведен в рамках теории возмущений. Реальные потенциалы ионов аппроксимировались псевдопотенциалами [12, 18]. В этом случае поверхностная энергия металла может быть представлена в виде
а[и (г )] = ао [и (г )] + ай [и (г )] + ай [и (г)], (9)
где а [и (г)] — поверхностная энергия (5), в рамках модели «желе»; ае, [и (г)] — вклад от разности электростатических взаимодействий электронов с дискретной решеткой и с однородным фоном (1); а и [и (г)] — вклад в сдвиг поверхностной энергии, обусловленный электростатической энергией конфигурации положительных зарядов [1з].
Для вычисления вклада а и [и (г)] было использовано выражение, полученное в работе [16]:
ац [и(г)] = а2и0,
где а — безразмерная константа, зависящая от структуры кристаллической решетки и ориентации кристаллографической плоскости по отношению к поверхности металла; Z — валентность металла.
При расчете вклада от электрон-ионного взаимодействия ое{ [п (г)] использован подход, представленный в работах [13, 21]:
ш
[п(г)]= | 5у(г)[п(г) — Пш (г)].
—ш
Здесь 5у(г) — усредненная по кристаллографическим плоскостям сумма ионных псевдопотенциалов за вычетом потенциала полубесконечного однородного фона ф+(г),
5v(г) = V (z).
(10)
В качестве потенциала, описывающего электростатическое поле отдельного иона в металле, был выбран псевдопотенциал Ашкрофта
0, \т — ц\<тс;
■(' ) =
-Z
г, r -Щ < тс,
(11)
r - т<\
где тс — параметр обрезания, который рассчитывается путем минимизации объемной энергии электронного газа в металле [12].
Для расчета 8у (г) выражение (10) удобно разбить на два слагаемых:
5v(z) = 8v1 (z) + 5у2 (z) = Z\ Vi (г)-
Z
-Z
Ьф+(г ).
, г — Ш ' \\г—г'\,
Вычисление 5у1 (г) и 5у2 (г) представляет собой несложную электростатическую задачу [16]. Выражение для 8у1 (г) с использованием псевдопотенциала (11) имеет вид
'0, г > 0;
8v1 (z ) = ■
u1 (z + md), -(m + 1)d<z<-md,
где и1 (г) = —2лп0^г + йб|—г — йЦЦ ; т = 0, 1, 2, .
тываемых от поверхности металла.
В случае гс < й/2 вклад 8у2 (г) может быть записан в виде
номера плоскостей, отсчи-
5V2 (z) =
0,
z > 0;
u2 (z + md), -(m + 1)d<z<-md.
Здесь u2 (z) = -2^n0d^rc -
d V Гс - d
х+— х+—
2 J 1 2
; m = 0, 1, 2, . — номера плоско-
стей, отсчитываемых от поверхности металла.
Если d / 2 < Гс < 3 / 2d, то для 5v2 (г) получается выражение
0, г>| Гс--|;
5V2 (г) =
U2 (г) + U2 (г + d), -d< г < ^Гс -;
u2 (г + md) + u2 (г + (m- 1)d) + u2 (г + (m + 1)d), -(m + 1)d < г < md.
Работа выхода электрона с металлической поверхности является одной из основных наблюдаемых характеристик поверхности металла. Большое количество экспериментальных данных по значениям работы выхода обусловливает выбор этого параметра в качестве критерия применимости пробных функций при расчете поверхностных характеристик [22-25]. Работа выхода А0 определяется как минимальная работа, которую необходимо совершить для удаления электрона из металла при Т = 0. Для ее расчета в рамках модели «желе» можно использовать подход [11, 12]
Л = Дф-Цср,
(12)
где Лф — высота электростатического барьера; Лф = ф(да)-ф; Цср — объемный химический потенциал, отсчитываемый от среднего электростатического по. ds( и ) тенциала, усредненный по объему металла, цср = —.
dn
В соответствии с (12) выражение для работы выхода может быть записано
как
Ло = -ф(-оо)-
(3л2 )3 3 0,056п03 + 0,0059пр ( 3 ^3 \
( 1 Y
0,079 + п03
V /
п03.
(13)
Внутри металла каждый электрон понижает свою энергию за счет образования обменно-корреляционной «дырки». Поэтому в значение полного эффективного потенциала внутри металла кроме кулоновского барьера вводят поправки на обменную и корреляционную энергии [10]:
1 2 1 , ч , ч ( 3 Лз , ,1 0,056и (г)з + 0,0059и(г)з
Уе/ (г ) = ф(г )-|-| и (г )3--—-
0,079 + п (г )3
Эффективная высота потенциального барьера определена по выражению
h = Vef (о)- Vef (-о) . (14)
2
Экспериментальные значения работы выхода для реальных металлов зависят от ориентации граней кристаллической решетки [22, 23]. Поэтому при расчете работы выхода (13) и высоты потенциального барьера (14) необходимо ввести поправки на электрон-ионное взаимодействие при учете дискретного распределения заряда ионов. Согласно работам [13, 23], данная поправка может быть представлена в виде
да
| 5у(г)иа(г)dz 5А = —^-,
да
| иа(г)dг
—да
, ч dn(г)
где иа (г) = —.
dг
Учитывая выражение (13), получаем окончательные формулы для работы выхода и высоты потенциального барьера:
А = А0 +5А; Н = Н0 +5А.
Результаты расчетов. Значения вариационных параметров, минимизирующих функционал а [и (г)] и а0 [и (г)], были получены путем численных расчетов. Значения исходных параметров, используемых для расчета поверхностных характеристик, приведены в табл. 1. Численное интегрирование выражений (5) и (9) позволило определить значения вариационных параметров для пробных функций всех трех типов (6)-(8).
Таблица 1
Значения исходных параметров, используемых для расчета поверхностных характеристик
Средняя Межплоскостное Константа а, отвечающая Радиус
Металл Грань плотность расстояние за ориентацию кристал- обрезания
п0, а. е. d, а. е. лической решетки Гс, а. е.
110 7,13 0,00563
K 100 0,00196 5,04 0,03100 2,311
111 2,91 0,06295
110 5,72 0,00563
Na 100 0,00377 4,05 0,03100 1,736
111 2,34 0,06295
110 4,66 0,00563
Li 100 0,00692 3,30 0,03100 1,299
111 1,90 0,06295
Значения полученных вариационных коэффициентов представлены в табл. 2 и 3. Коэффициенты со штрихом рассчитаны с учетом дискретности кристаллической решетки, а коэффициенты без штриха — в рамках модели «желе» и представлены в таблицах для оценки влияния дискретности. Для выполнения условия гладкости функции и3(г), проводилось сглаживание этой функции в области г = 0.
Таблица 2
Значения вариационных коэффициентов функций n1(z) и n2(z)
Элемент Г Коэффициенты rn(z) Коэффициенты m(z)
a' a'iL a2 R
a1 a2L a2 R
110 2,41 0,38 1,37
K 100 2,54 1,98 0,37 0,43 1,43 1,11
111 1,82 0,50 0,84
110 2,48 0,40 1,54
Na 100 2,52 2,10 0,40 0,43 1,54 1,36
111 1,95 0,49 1,07
110 2,56 0,41 1,79
Li 100 2,50 2,23 0,41 0,47 1,79 1,40
111 2,08 0,47 1,40
Таблица 3
Значения вариационных коэффициентов функций n0(z) и n3(z)
Элемент Грань Коэффициенты m(z) Коэффициенты rn(z)
P P' a3L a'L C3L c'l a3R a3R
K 110 100 111 1,32 1,23 0,76 0,67 0,56 0,41 0,36 0,38 —0,30 —1,00 —1,30 —1,30 1,27 1,23 0,86 0,79
Na 110 100 111 1,28 1,24 0,86 0,76 0,57 0,44 0,44 0,47 —0,50 —1,00 —1,20 —1,20 1,27 1,33 0,96 0,85
Li 110 100 111 1,24 1,31 0,97 0,86 0,51 0,54 0,46 0,51 — 1,10 —0,90 —1,20 —1,20 1,08 1,24 1,10 0,90
Результаты расчетов поверхностной энергии в рамках модели однородного фона с учетом дискретности кристаллической решетки, а также сравнение полученных результатов с экспериментальными значениями приведены в табл. 4. Этот расчет проводился только для трех щелочных металлов, так как модель «желе» хорошо зарекомендовала себя для описания поверхностной энергии металлов с достаточно малой концентрацией свободных электронов (п < 0,008) [10].
Таблица 4
Значения поверхностной энергии
Элемент Грань ac-105 aJ-105 0 -105 d -105 02-105 a2-105 О -105 o3-105 dexp • 105
110 5,64 6,67 6,70 6,01
K 100 4,44 7,68 5,70 9,83 6,08 9,27 4,20 7,00 9,66
111 12,49 14,78 14,32 12,39
110 9,74 12,03 13,42 10,28
Na 100 7,11 13,77 9,78 17,85 11,42 18,64 6,77 13,05 15,46
111 21,56 26,52 27,62 21,72
110 14,88 19,70 26,19 15,51
Li 100 8,50 22,10 14,13 29,98 20,65 36,11 8,33 21,77 32,86
111 34,55 44,44 51,52 35,09
Как было отмечено в работе [10], работа выхода и другие наблюдаемые характеристики поверхности металла зависят от вариационной производной 5с[и]/5и функционала поверхностной энергии с[и], а не от значения с[и]. Этим объясняется тот факт, что полученная поверхностная энергия только по порядку величины совпадает с ее экспериментальным значением, в то время как работа выхода и высота потенциального барьера существенно лучше соответствуют результатам эксперимента. Учет дискретности кристаллической решетки значительно улучшает результаты расчета поверхностной энергии. Результаты расчета значений работы выхода и высоты потенциального барьера, а также сравнение их с экспериментальными значениями представлены в табл. 5 и 6.
Таблица 5
Значения работы выхода и высоты потенциального барьера для функций п1(г) и п2(х)
Элемент Грань Функция rn(z) Функция m(z) Эксперименталь-
ные значения
Л AI h1 h1 A2 A2 h2 h2 Aexp hexp
110 3,03 5,07 3,26 5,31
K 100 2,91 3,39 4,95 5,44 3,24 3,18 5,29 5,23 2,22 4,12
111 3,28 5,32 3,03 5,08
110 3,31 6,46 3,48 6,64
Na 100 3,21 3,91 6,37 7,06 3,46 3,37 6,61 6,52 2,35 4,85
111 3,92 7,07 3,18 6,33
110 3,50 8,24 3,45 8,19
Li 100 3,57 4,35 8,31 9,09 3,37 3,36 8,11 8,29 2,38 6,08
111 4,59 9,32 3,05 7,78
Таблица 6
Значения работы выхода и значения высоты потенциального барьера для функций т(г) и п3(г)
Элемент Грань Функция n0(z) Функция n3(z) Эксперименталь-
ные значения
Ac A0 h„ h'n A3 A3 h3 h3 Aexp hexp
110 2,83 4,87 2,64 4,69
K 100 2,76 3,12 4,80 5,16 2,53 2,67 4,58 4,71 2,22 4,12
111 3,10 5,14 2,62 4,66
110 3,03 6,19 2,78 5,94
Na 100 2,93 3,48 6,09 6,64 2,62 2,75 5,78 5,90 2,35 4,85
111 3,56 6,72 2,68 5,82
110 3,07 7,80 2,58 7,31
Li 100 3,10 3,73 7,84 8,46 2,89 2,53 7,63 7,27 2,38 6,08
111 3,96 8,69 2,58 7,31
Значения поверхностных характеристик, рассчитанные для функций По(г) и Пъ(х), наиболее точно соответствуют экспериментальным значениям. Распределения электронной плотности функций По(г) и пз(г) при по = 0,00377 (Ыа) для грани 110 изображены на части а рисунка, эффективный потенциал, полученный с использованием соответствующих функций, — на части б рисунка.
а б
Распределение электронной плотности (а) и эффективного потенциала (б): 1 — функция п0(х); 2 — функция пз(х); 3 — функция п0(г); 4 — функция п3(х)
Заключение. Исследована возможность применения МФП в рамках модели «желе» с учетом дискретности кристаллической решетки с использованием новых пробных функций распределения к описанию характеристик поверхности ряда щелочных металлов. Проведен самосогласованный расчет поверхностной энергии с учетом градиентных поправок на неоднородность электронного газа в приповерхностной области для кинетической энергии. Установлен оптимальный вид этих функций, рассчитаны характеризующие их параметры. Получены значения характеристик поверхности раздела металл-вакуум — работа выхода и высота потенциального барьера для различных граней монокристаллов щелочных металлов. Построены эффективные потенциалы для выбранных металлов. Показано, что в отличие от экспоненциальной пробные функции распределения дают приемлемые результаты, но наиболее хорошо согласие расчетных характеристик с экспериментальными значениями обеспечивает функция, учитывающая фриделевские осцилляции.
Проведенные расчеты поверхностных характеристик показали, что учет дискретности кристаллической решетки позволяет получить хорошее согласие расчетных значений с имеющимися в литературе экспериментальными данными для щелочных металлов. Представленный метод расчета влияния кристаллической решетки дает оптимальные результаты для функции плотности, учитывающей фриделевские осцилляции.
Результаты настоящей работы могут быть использованы при описании как чистой поверхности металлов, так и процессов взаимодействия излучения с ве-
ществом, в частности, при описании элементарных процессов, происходящих на границе раздела металл-плазма, а также при описании воздействия электромагнитного излучения высокой мощности на поверхность металла.
ЛИТЕРАТУРА
1. Kvaal S., Helgaker T. Ground-state densities from the Rayleigh — Ritz variation principle and from density-functional theory // J. Chem. Phys. 2015. Vol. 143. No. 18. DOI: 10.1063/1.4934797 URL: http://aip.scitation.org/doi/10.1063/1.4934797
2. Burke K. Perspective on density functional theory // J. Chem. Phys. 2012. Vol. 136. No. 15. P. 10. DOI: 10.1063/1.4704546 URL: http://aip.scitation.org/doi/10.1063/1.4704546
3. Kryachko E.S., Ludena E.V. Density functional theory: Foundation reviewed // Phys. Rep. — Sect. Phys. Lett. 2014. Vol. 544. Nc. 2. P. 123-239.
4. Глушков В.Л., Еркович О.С. Метод многочастичных функционалов плотности в описании двухкомпонентных систем // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2011. № 1. С. 33-39.
5. Еркович О.С. Описание свойств электронного газа вблизи поверхности металла методом многочастичных функционалов плотности с учетом периодической структуры металла // Инженерный журнал: Наука и инновации. 2013. Вып. 8.
DOI: 10.18698/2308-6033-2013-8-864
URL: http://engjournal.ru/catalog/ fundamentals/physics/864.html
6. Saravanan R., Rani M.P. Metal and alloy bonding: An experimental analysis. Charge density in metals and alloys. London: Springer, 2012. 151 p.
7. Suryanarayana C. Experimental techniques in materials and mechanics. CRC Press, Taylor & Francis Group, 2011. 430 p.
8. Stalke D. Meaningful structural descriptors from charge density // Chem. Eur. J. 2011. Vol. 17. No. 34. P. 9264-9278. DOI: 10.1002/chem.201100615
URL: http://onlinelibrary.wiley.com/doi/10.1002/chem.201100615/full
9. Глушков В.Л., Еркович О.С. Энергетические характеристики поверхностей щелочных металлов с учетом фриделевских осцилляций электронной плотности // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2014. Вып. 4. С. 9-18.
10. Партенский М.Б. Самосогласованная электронная теория металлической поверхности // УФН. 1979. Т. 128. № 1. С. 69-106. DOI: 10.3367/UFNr.0128.197905c.0069
11. Smith J.R. Self-consistent theory of electron work functions and surface potential characteristics for selected metals // Phys. Rev. 1969. Vol. 181. No. 2. P. 522-529.
DOI: 10.1103/PhysRev.181.522
URL: https://journals.aps.org/pr/abstract/10.1103/PhysRev.181.522
12. Вакилов А.Н., Мамонова М.В., Матвеев В.В., Прудников В.В. Теоретические модели и методы в физике поверхности. Омск: Изд-во ОмГУ, 2005. 212 с.
13. Kiejna A., Wojciechowski K.F. Metal surface electron physics. Oxford: Alden Press, 1996. 303 p.
14. Zhao Y., Truhlar D.G. Density functionals with broad applicability in chemistry // Acc. Chem. Res. 2008. Vol. 41. No. 2. P. 157-167. DOI: 10.1021/ar700111a
URL: http://pubs.acs.org/doi/abs/10.1021/ar700111a
15. Jacob C.R., Neugebauer J. Subsystem density-functional theory // WIREs: Computational Molecular Science. 2014. Vol. 4. No. 4. P. 325-326. DOI: 10.1002/wcms.1175
URL: http://onlinelibrary.wiley.com/doi/10.1002/wcms.1175/full
16. Lang N.D., Kohn W. Theory of metal surfaces: Charge density and surface energy // Phys. Rev. B. 1970. Vol. l. No. 12. P. 4555-4568. DOI: 10.1103/PhysRevB.1.4555
URL: https://journals.aps.org/prb/abstract/10.1103/PhysRevB.L4555
17. Giustino F. Materials modelling using density functional theory: Properties and predictions. Oxford University Press, 2014. 303 p.
18. Sholl D., Steckel J.A. Density functional theory: A practical introduction. Wiley, 2009. 252 p.
19. Cohen A.J., Mori-Sanchez P., Yang W.T. Challenges for density functional theory // Chem. Rev. 2012. Vol. 112. No. 1. P. 289-320. DOI: 10.1021/cr200107z
URL: http://pubs.acs.org/doi/abs/10.1021/cr200107z
20. Pevarati R., Truhlar D.G. M11-L: A local density functional that provides improved accuracy for electronic structure calculations in chemistry and physics // J. Phys. Chem. Lett. 2012. Vol. 3. No. 1. P. 117-124. DOI: 10.1021/jz201525m
URL: http://pubs.acs.org/doi/abs/10.1021/jz201525m
21. Lang N.D., Kohn W. Theory of metal surfaces: Work function // Phys. Rev. B. 1971. Vol. 3. No. 4. P. 1215-1223. DOI: 10.1103/PhysRevB.3.1215
URL: https://journals.aps.org/prb/abstract/10.1103/PhysRevB.3.1215
22. Cramer C.J., Truhlar D.G. Density functional theory for transition metals and transition metal chemistry // Phys. Chem. Chem. Phys. 2009. Vol. 11. No. 46. P. 10757-10816.
DOI: 10.1039/b907148b
23. Матвеев А.В. Обобщенная модель поверхностной сегрегации с учетом диэлектрических свойств среды: щелочные металлы и сплавы // Вестник Омского университета. 2013. № 4. С. 91-101.
24. Konig T., Simon G.H., Rust H.P., Heyde M. Work function measurements of thin oxide films on metals-MgO on Ag(001) // J. Phys. Chem. C. 2009. Vol. 113. No. 26. P. 11301-11305.
DOI: 10.1021/jp901226q URL: http://pubs.acs.org/doi/abs/10.1021/jp901226q
25. Schimka L., Gaudoin R., Klimes J., Marsman M., Kresse G. Lattice constants and cohesive energies of alkali, alkaline-earth, and transition metals: Random phase approximation and density functional theory results // Phys. Rev. B. 2013. Vol. 87. No. 21. P. 214102-1-214102-8.
DOI: 10.1103/PhysRevB.87.214102
URL: https://journals.aps.org/prb/abstract/10.1103/PhysRevB.87.214102
Глушков Владимир Леонидович — ассистент кафедры «Физика» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5, стр. 1).
Еркович Ольга Станиславовна — канд. физ.-мат. наук, доцент кафедры «Физика» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5, стр. 1).
Просьба ссылаться на эту статью следующим образом:
Глушков В.Л., Еркович О.С. Характеристики поверхности щелочных металлов с учетом дискретности кристаллической решетки и фриделевских осцилляций электронной плотности // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2017. № 4. C. 75-89. DOI: 10.18698/1812-3368-2017-4-75-89
SURFACE CHARACTERISTICS OF ALKALI METALS WITH THE DISCRETE LATTICE AND FRIEDEL OSCILLATIONS OF THE ELECTRON DENSITY
V.L. Glushkov O.S. Erkovich
glushkov-vl@bmstu.ru erkovic@bmstu.ru
Bauman Moscow State Technical University, Moscow, Russian Federation
Abstract
The study makes an attempt to describe the effect of the crystal lattice on the surface characteristics of alkali metals in the framework of the density functional method. In the first approximation the real potential of the crystal lattice has been replaced by the potential of a homogeneous positive background (the jellium model). We calculated the difference between the discrete-ions potential and the potential of jel-lium in the framework of the perturbation theory. To calculate the impact of electron-ion interaction on the energy characteristics of the surface, the real potentials of ions were approximated by Ashcroft pseudopotential. The electron density in the inhomogeneous electron gas near the surface was approximated by trial distribution functions which take into account Friedel oscillations and lattice effects. Furthermore, we carried out self-consistent calculation of the surface energy taking into account the gradient expansion for the kinetic energy. In this paper, we defined surface characteristics, such as the work function and the potential barrier height and constructed effective potentials for alkaly metals. All the characteristics of the metal surface were calculated within the density functional method
Keywords
Density functional method, electron density, discrete crystal lattice, work function, Friedel oscillations
REFERENCES
[1] Kvaal S., Helgaker T. Ground-state densities from the Rayleigh — Ritz variation principle and from density-functional theory. J. Chem. Phys., 2015, vol. 143, no. 18. DOI: 10.1063/1.4934797 Available at: http://aip.scitation.org/doi/10.1063/L4934797
[2] Burke K. Perspective on density functional theory. J. Chem. Phys., 2012, vol. 136, no. 15. DOI: 10.1063/1.4704546 Available at: http://aip.scitation.org/doi/10.1063/L4704546
[3] Kryachko E.S., Ludena E.V. Density functional theory: Foundation reviewed. Phys. Rep. — Sect. Phys. Lett., 2014, vol. 544, no. 2, pp. 123-239.
[4] Glushkov V.L., Erkovich O.S. Method of many-particle density functionals for description of two-component system. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2011, no. 1, pp. 33-39 (in Russ.).
[5] Erkovich O.S. Application of the multi-particle density functional approach to the properties description of the electron gas of the metal surface with consideration of the periodic structure of metal. Inzhenernyy zhurnal: Nauka i innovatsii [Engineering Journal: Science and Innovation], 2013, no. 8 (in Russ.). DOI: 10.18698/2308-6033-2013-8-864
Available at: http://engjournal.ru/ catalog/fundamentals/physics/864.html
[6] Saravanan R., Rani M.P. Metal and alloy bonding: An experimental analysis. Charge density in metals and alloys. London, Springer, 2012. 151 p.
[7] Suryanarayana C. Experimental techniques in materials and mechanics. CRC Press, Taylor & Francis Group, 2011. 430 p.
[8] Stalke D. Meaningful structural descriptors from charge density. Chem. Eur. J., 2011, vol. 17, no. 34, pp. 9264-9278. DOI: 10.1002/chem.201100615
Available at: http://onlinelibrary.wiley.com/doi/10.1002/chem.201100615/full
[9] Glushkov V.L., Erkovich O.S. Energy response of alcaline metals surfaces taking into account Friedel oscillation of electron density. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki, 2014, no. 4, pp. 9-18 (in Russ.).
[10] Partenskiy M.B. Self-consistent electron theory of a metallic surface. Sov. Phys. Usp., 1979, vol. 22, no. 5, pp. 330-351. DOI: 10.1070/PU1979v022n05ABEH005498
Available at: http://iopscience.iop.org/article/10.1070/PU1979v022n05ABEH005498/meta
[11] Smith J.R. Self-consistent theory of electron work functions and surface potential characteristics for selected metals. Phys. Rev., 1969, vol. 181, no. 2, pp. 522-529.
DOI: 10.1103/PhysRev.181.522
Available at: https://journals.aps.org/pr/abstract/10.1103/PhysRev.181.522
[12] Vakilov A.N., Mamonova M.V., Matveev V.V., Prudnikov V.V. Teoreticheskie modeli i metody v fizike poverkhnosti [Theoretical models and methods in surface physics]. Omsk, OmSU Publ., 2005. 212 p.
[13] Kiejna A., Wojciechowski K.F. Metal surface electron physics. Oxford, Alden Press, 1996. 303 p.
[14] Zhao Y., Truhlar D.G. Density functionals with broad applicability in chemistry. Acc. Chem. Res, 2008, vol. 41, no. 2, pp. 157-167. DOI: 10.1021/ar700111a
Available at: http://pubs.acs.org/doi/abs/10.1021/ar700111a
[15] Jacob C.R., Neugebauer J. Subsystem density-functional theory. WIREs: Computational Molecular Science, 2014, vol. 4, no. 4, pp. 325-326. DOI: 10.1002/wcms.1175
Available at: http://onlinelibrary.wiley.com/doi/10.1002/wcms. 1175/full
[16] Lang N.D., Kohn W. Theory of metal surfaces: Charge density and surface energy. Phys. Rev. B., 1970, vol. l, no. 12, pp. 4555-4568. DOI: 10.1103/PhysRevB.1.4555
Available at: https://journals.aps.org/prb/abstract/10.1103/PhysRevB.1.4555
[17] Giustino F. Materials modelling using density functional theory: Properties and predictions. Oxford University Press, 2014. 303 p.
[18] Sholl D., Steckel J.A. Density functional theory: A practical introduction. Wiley, 2009. 252 p.
[19] Cohen A.J., Mori-Sanchez P., Yang W.T. Challenges for density functional theory. Chem. Rev., 2012, vol. 112, no. 1, pp. 289-320. DOI: 10.1021/cr200107z
Available at: http://pubs.acs.org/doi/abs/10.1021/cr200107z
[20] Pevarati R., Truhlar D.G. M11-L: A local density functional that provides improved accuracy for electronic structure calculations in chemistry and physics. J. Phys. Chem. Lett, 2012, vol. 3, no. 1, pp. 117-124. DOI: 10.1021/jz201525m
Available at: http://pubs.acs.org/doi/abs/10.1021/jz201525m
[21] Lang N.D., Kohn W. Theory of metal surfaces: Work function. Phys. Rev. B., 1971, vol. 3, no. 4, pp. 1215-1223. DOI: 10.1103/PhysRevB.3.1215
Available at: https://journals.aps.org/prb/abstract/10.1103/PhysRevB.3.1215
[22] Cramer C.J., Truhlar D.G. Density functional theory for transition metals and transition metal chemistry. Phys. Chem. Chem. Phys, 2009, vol. 11, no. 46, pp. 10757-10816.
DOI: 10.1039/b907148b
[23] Matveev A.V. Generalized model of the surface segregation with inclusion of dielectric properties of the medium: Alkali metals and alloys. Vestnik Omskogo universiteta [Herald of Omsk University], 2013, no. 4, pp. 91-101 (in Russ.).
[24] Konig T., Simon G.H., Rust H.P., Heyde M. Work function measurements of thin oxide films on metals-MgO on Ag(001). J. Phys. Chem. C., 2009, vol. 113, no. 26, pp. 11301-11305.
DOI: 10.1021/jp901226q Available at: http://pubs.acs.org/doi/abs/10.1021/jp901226q
[25] Schimka L., Gaudoin R., Klimes J., Marsman M., Kresse G. Lattice constants and cohesive energies of alkali, alkaline-earth, and transition metals: Random phase approximation and density functional theory results. Phys. Rev. B., 2013, vol. 87, no. 21, pp. 214102-1-214102-8.
DOI: 10.1103/PhysRevB.87.214102
Available at: https://journals.aps.org/prb/abstract/10.1103/PhysRevB.87.214102
Glushkov V.L. — Assistant of Physics Department, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Erkovich O.S. — Cand. Sc. (Phys.-Math.), Assoc. Professor of Physics Department, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Please cite this article in English as:
Glushkov V.L., Erkovich O.S. Surface Characteristics of Alkali Metals with the Discrete Lattice and Friedel Oscillations of the Electron Density. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2017, no. 4, pp. 75-89. DOI: 10.18698/1812-3368-2017-4-75-89
