влияние адсорбции атомов щелочных металлов на величину

advertisement
ФИЗИКА
Вестник Омского университета, 2004. №1. С. 31–36.
c Омский государственный университет
°
УДК 539.612.001
ВЛИЯНИЕ АДСОРБЦИИ АТОМОВ ЩЕЛОЧНЫХ
МЕТАЛЛОВ НА ВЕЛИЧИНУ РАБОТЫ ВЫХОДА
ЭЛЕКТРОНА С ПОВЕРХНОСТИ МЕТАЛЛОВ
А.В. Матвеев, М.В. Мамонова, В.В. Прудников
Омский государственный университет, кафедра теоретической физики
644077, Омск, пр. Мира, 55а 1
Получена 5 ноября 2003 г.
The results of recent density functional theory calculations of work function changes induced by
the activated adsorption of alkali metal atoms on metallic surfaces are discussed.
Адсорбция щелочных металлов уже давно занимает центральное место среди объектов теоретического исследования свойств адсорбции на
поверхности кристаллов. Интерес к этим системам обусловлен, главным образом, технологическим значением адсорбатов из щелочных металлов для получения эффективных электродов с
низкой работой выхода и осуществления гетерогенного катализа с улучшенной активностью
и селективностью катализаторов. Современные
экспериментальные исследования распределения
атомов щелочных металлов на поверхностном
слое подложки выявили новые адсорбционные
фазы и, в частности, структуры с заместительной (или активированной) адсорбцией, при которой поверхность подложки испытывает значительную реконструкцию [1]. В связи с этим,
возникает необходимость теоретического описания таких структур и расчета работы выхода
электрона с поверхности подложки. Настоящая
работа посвящена объяснению этих последних
данных экспериментальных исследований на основе предложенной нами модели активированной
адсорбции [2]. Опуская детальное описание методики расчета работы выхода электрона с поверхности металлов в рамках теории функционала
плотности, которое можно найти в [3], здесь мы
приведем лишь основные расчетные формулы.
Рассмотрим полубесконечный металл, ограниченный плоской поверхностью и занимающий область z < −D . Пленка адсорбата занимает область D < z < D + h . Между пленкой и поверхностью вводится вакуумный зазор шириной
2D . В результате процессов взаимного перераспределения атомов адсорбата и атомов подлож1
e-mail: prudnikv@univer.omsk.su
ки электронная плотность поверхностной области подложки будет отличаться от ее объемной
электронной плотности n1 , а сама поверхность
может подвергаться различным реконструкциям.
При этом значительную реконструкцию испытывает обычно первый приповерхностный слой подложки [1]. В связи с этим, первый приповерхностный слой подложки выделен в отдельную область
со средней плотностью заряда n2 и толщиной l .
При заместительной адсорбции в первом приповерхностном слое и пленке адсорбата образуется твердый раствор из атомов подложки и
адсорбата. При рассмотрении бинарных растворов мы использовали следующую терминологию.
Обозначим адсорбат символом A, субстрат – S.
Тогда для состава бинарного раствора пленки
справедлива формула: Ap S1−p , где индекс p обозначает относительную долю атомов адсорбата.
Электронную плотность пленки можно представить в виде
n3 (Θ, p) = Z3 ns3 /h,
(1)
где заряд ионов пленки Z3 определяется выражением
Z3 = pZ a + (1 − p)Z s ,
(2)
ns3 – поверхностная концентрация атомов в пленке, h – толщина пленки,
h = pda + (1 − p)ds ,
(3)
da , ds – расстояния между наиболее плотноупакованными плоскостями в кристаллах щелочного
металла и подложки соответственно. Аналогичные формулы использовались для бинарного раствора в первом приповерхностном слое (относительная доля атомов адсорбата в этом случае обозначалась индексом p/ : Ap/ S1−p/ ) с электронной
32
А.В. Матвеев, М.В. Мамонова, В.В. Прудников
плотностью приповерхностного слоя и его толщиной, задаваемых соотношениями:
n2 (Θ, p/ ) = Z2 ns2 /l,
(4)
l = p/ da + (1 − p/ )d1 ,
(5)
где ns2 – поверхностная концентрация атомов в
слое, Z2 – заряд ионов слоя определяется аналогично выражению (2) при замене p на p/ .
Предполагая равномерное распределение атомов по поверхности с повторением ее симметрии,
параметр заполнения Θ можно выразить через
параметры бинарных растворов p и p/ ,
Θ = pq + p/ q / ,
(6)
где параметры q и q / определяются выражениями:
q = ns3 /ns1 , q / = ns2 /ns1
(7)
и характеризуют степень заполнения атомами
пленки и приповерхностного слоя соответственно. Параметры Θ, q и q / связаны друг с другом
соотношением
Θ = q + q / − 1.
(8)
Степень заполнения Θ адатомами поверхности
задается таким образом, чтобы для Θ = 1 число
адатомов равнялось числу атомов в поверхностной грани нереконструированной подложки.
Мы получили следующее выражение для величины работы выхода электронов с поверхности
подложки с адсорбированной моноатомной пленкой щелочного металла [3]:
a
s
Φ(Θ) = Φ0 − µ + Φm
ei + Φei + Φei .
(9)
Здесь Φ0 – дипольный потенциальный барьер на
поверхности подложки:
Φ0 =
4πn2
,
β2
(10)
β – вариационный параметр модели, значения
которого определялись из условия минимальности полной межфазной энергии системы (см. [2]).
Выражение для химического потенциала µ c учетом обменно-корреляционных и псевдопотенциальных поправок имеет вид
µ = 0, 5(3π 2 n1 )2/3 −
2/3
−
0,056n1
2/3
−
1/3
1/3 2
)
¡ 4πn ¢1/3
1
3
1
π
+0,0059n1
(0,079+n1
−0, 4Z1
¡ 3n ¢1/3
−
(11)
2
+ 4πn1 rc1
,
где rc1 – параметр псевдопотенциала Ашкрофта
для ионов подложки [3].
s
a
В формуле (9) Φm
ei , Φei , Φei представляют поправки на электрон-ионное взаимодействие при
учете дискретного распределения заряда ионов в
узлах кристаллической решетки подложки, приповерхностного слоя и адсорбированной пленки
соответственно. Для них мы получили следующие выражения [3]:
4πn1
Φm
ei = − β 2 (1 + β(D + h))×
−β(2D+l+h) −βd1 /2
×e
e
{βd1 ch(βrc1 )−2 sh(βd1 /2)}
2−(1+β(D+h))e−β(2D+l+h+d1 )
(12)
– для металлической подложки,
2
Φsei = − 4πn
β 2 (1 + β(D + h))×
×e
e
{βl ch(βrc2 )−2 sh(βl/2)}
2−(1+β(D+h))e−β(2D+l+h+d1 )
−β(2D+h) −βl/2
(13)
– для приповерхностного слоя,
3
Φaei = − 4πn
β 2 (1 + β(D + h))×
{βh ch(βrc3 )−2 sh(βh/2)}
× e2−(1+β(D+h))e
−β(2D+l+h+d1 )
−βh/2
(14)
– для пленки адсорбата.
На рис. 1 приведены графики зависимости
изменения работы выхода электрона ∆Φ(Θ) =
Φ(Θ) − Φm , обусловленного адсорбцией атомов
натрия, от состава p раствора в пленке (A p S 1−p )
при параметре покрытия Θ = 0, 3 для плотноупакованной грани (111) подложки из алюминия. На
рисунке видно, что для всех значений параметров
p и q , характеризующих состав бинарного раствора и распределение атомов в адсорбированной
пленке, происходит понижение величины работы
выхода электрона по сравнению с работой выхода чистой поверхности Al(111). Минимальное
уменьшение работы выхода на 0,34 эВ соответствует состоянию системы, при котором бинарные растворы в пленке и приповерхностной области характеризуются следующими значениями
параметров: q = 0, 6, p = 0, 5 (q / = 0, 7, p/ = 0)
– все атомы натрия расположены в пленке. Максимальное уменьшение работы выхода на 2,39 эВ
соответствует состоянию системы, при котором
бинарные растворы в пленке и приповерхностной области характеризуются значениями параметров: q = 1, p = 0 (q / = 0, 3, p/ = 1) – все
атомы натрия расположены в приповерхностном
слое, а атомы алюминия из этого слоя смещены в пленку. Состоянию системы с монослойным
покрытием атомами натрия (неактивированная
адсорбция, q = 0, 3, p = 1) соответствует промежуточное значение работы выхода: ∆Φ = −1, 44
эВ.
Согласно экспериментальным данным [1,4],
при параметре покрытия Θ ≈ 0, 3 в системе
Na/Al(111) наблюдается
образование структур
√
√
с симметрией ( 3 × 3)R300 , где на каждые два
Влияние адсорбции атомов щелочных металлов...
Рис. 1. Рассчитанная зависимость изменения работы
выхода электрона ∆Φ от состава p раствора в пленке
при параметре покрытия Θ = 0, 3 для системы
Na p Al 1−p /Al(111). Приведены графики для следующих
значений параметра q : 1 – q = 0, 3 ; 2 – q = 0, 4 ;
3 – q = 0, 5 ; 4 – q = 0, 6 ; 5 – q = 0, 7 ; 6 – q = 0, 8 ;
7 – q = 0, 9 ; 8 - q = 1
атома алюминия приходится один атом натрия. В
нашей модели этой структуре можно сопоставить
значения параметров q = 1, p ≈ 0, 3 (q / = 0, 3,
p/ = 0). На рис. 1 этим значениям параметров соответствует ∆Φ = −1, 46 эВ, которое находится в
достаточно хорошем соответствии с экспериментальной величиной ∆Φ = −1, 7 эВ [4].
На рис. 2 приведены графики зависимости изменения работы выхода электрона ∆Φ(Θ), обусловленного адсорбцией атомов натрия, от состава p раствора в пленке (A p S 1−p ) при параметре
покрытия Θ = 0, 5 для плотноупакованной грани (111) подложки из алюминия. Минимальное
уменьшение работы выхода на 0,07 эВ соответствует состоянию системы со значениями параметров: q = 0, 8, p = 0, 6 (q / = 0, 7, p/ ≈ 0) –
все атомы натрия расположены в пленке. Максимальное уменьшение работы выхода на 1,84 эВ
соответствует состоянию системы с q = 1, p = 0
(q / = 0, 5, p/ = 1) – все атомы натрия расположены в приповерхностном слое, а все атомы
алюминия из этого слоя смещены в пленку. Для
монослойного покрытия ( q = 0, 5, p = 1) расчет
дает значение ∆Φ = −1, 06 эВ.
Согласно экспериментальным данным [4], при
параметре покрытия Θ ≈ 0, 5 в системе Na/Al(111)
наблюдается образование структур с симметрией
2 × 1, где так же, как и при Θ ≈ 0, 3, на каждые
два атома алюминия приходится один атом натрия. Но в этом случае данной структуре соответствуют значения параметров q = 1, p ≈ 0, 3 при
q / = 0, 5, p/ ≈ 0, 4. Согласно рис. 2, этим значениям параметров соответствует ∆Φ = −1, 15 эВ,
находящееся в хорошем соответствии с экспериментальной величиной ∆Φ = −1, 25 эВ [4].
33
Рис. 2. Рассчитанная зависимость изменения работы
выхода электрона ∆Φ от состава p раствора в пленке
при параметре покрытия Θ = 0, 5 для системы
Na p Al 1−p /Al(111). Приведены графики для следующих
значений параметра q : 1 – q = 0, 5 ; 2 – q = 0, 6 ; 3 –
q = 0, 7 ; 4 – q = 0, 8 ; 5 – q = 0, 9 ; 6 – q = 1
[1] Renee D. Diehl, Ronan McGrath. Structural studies
of alkali metal adsorption and coadsorption on metal
surfaces // Surf. Sci. Rep. 1996. № 23. P. 43–171.
[2] Матвеев А.В., Мамонова М.В., Прудников В.В.
Модель активированной адсорбции атомов щелочных металлов на металлических поверхностях
// Вестн. Ом. ун-та. 2003. № 4. С. 27–29.
[3] Мамонова М.В., Матвеев А.В., Прудников В.В.
Многопараметрическая модель адсорбции атомов
щелочных металлов на металлических поверхностях // ФММ. 2002. 94. № 5. С. 16–25.
[4] Nagao T., Iizuka Y., Umeuchi M., Shimazaki T.,
Oshima C. Vibrations of alkali-metal atoms
chemisorbed on the Al(111) surface // Surf. Sci.
1995. 329. P. 269–275.
Download