CC BY
23
УДК: 539.91/37
В. В. Погосов, Е. В. Васютин, А. В. Коротун, А. В. Бабич
ДЕФОРМАЦИОННАЯ ЗАВИСИМОСТЬ РАБОТЫ ВЫХОДА ЭЛЕКТРОНОВ И КОНТАКТНОЙ РАЗНОСТИ ПОТЕНЦИАЛОВ В МЕТАЛЛАХ
Поверхностный стресс и контактная разность потенциалов упруго деформированных граней А1, Си, Аи, N1 и Л кристаллов вычислены с помощью самосогласованного метода Хартри-Фока-Кона-Шема. Расчеты демонстрируют уменьшение/увеличение работы выхода при растяжении/сжатии металлического кристалла. Результаты вычислений указывают, что измерения методом Кельвина контактной разности потенциалов деформированной поверхности соответствуют не изменению работы выхода электронов, а изменению поверхностного значения эффективного потенциала. Полученные величины стресса, работы выхода и контактной разности потенциалов находятся в хорошем согласии с результатами вычислений из первых принципов.
1 Введение
К настоящему времени накопился определенный объем экспериментальных исследований за. „, „ м о,, „ р. 6 о, ы .„ход, М от деформации. По этой зависимости можно судить о параметрах напряженного состояния металла: величине остаточных механических напряжений, дислокационной структуре и пр. Химическая активность, определяемая величиной поверхностной энергии или стресса, также чувствительна к деформации поверхности металла [1].
Прямые измерения, использующие метод Кельвина (метод динамического конденсатора) [2-4], указывают на уменьшение/увеличение контактной разности потенциалов (КРП) упруго растянутого/ сжатого плоского металлического образца. Эти, на первый взгляд, неожиданные результаты означают, что работа выхода увеличивается/уменьшается при одноосном растяжении/сжатии металлического образца. Этот факт противоречит другому факту: работа выхода простых металлов уменьшается с уменьшением концентрации электронов металла (т.е. при переходе Al ^ № ^ Cs в таблице Менделеева).
Общепризнанный метод измерения работы выхода в зависимости от деформации вдоль х-оси (см. рис.1) основан на выражении:
КРП, (1)
т.е. работа выхода как бы увеличивается для растянутого образца. Изменение работы выхода
АШ = Ш (ихх) - Ш (0) было измерено в [2-4] на грани металлического образца,перпендикулярной у -или г -направлениям, ихх - относительная деформация образца вдоль х -оси.
Рис. 1. Качественная схема растянутого образца вдоль х-направления, ихх > 0 - относительная деформация, 5 - площадь грани, >> 8у,8х.
В данной работе вычислены деформационные зависимости поверхностной энергии, работы выхода и КРП для различных плоскостей таких металлов, как А1, И, N1, Си и Аи. Обсуждается проблема корректного определения работы выхода. Точность соотношения (1) тестирована полностью самосогласованными вычислениями. Показано, что использование выражения (1) может привести к неверным результатам в диагностике упруго-напряженной металлической поверхности.
2 Методика вычисления
В рамках метода функционала плотности полная энергия металла является функционалом неоднородной электронной концентрации п(г), которая стремится п(г) ^ п0 к своему объемному (постоянному) значению в объеме металла и быстро убывает за поверхностью в вакуум. Задавая деформацию ихх по оси х, деформация по другим направлениям определяется коэффициентом Пуассона. Равновесный профиль электронов рассчитывается из условия минимума полной энергии кристалла. Для этого используется метод Хартри-Фока-Кона-Шема, учитывающий обменно-корреляцион-ные эффекты в неоднородном электронном газе на
© В. В. Погосов, Е. В. Васютин, А. В. Коротун, А. В. Бабич 2006 г.
- 0219яшВестникяИвигателестроенияя1 4/т006
- 115 -
z
У
х
£
Б
У
и
хх
фоне ионного желе [1,5]. Когда равновесный профиль известен, рассчитывается поверхностная энергия и работа выхода. Диагональная хх - компонента поверхностного стресса для верхней грани образца, на которой в [2-4] измерялась КРП, равна
(2)
где у - удельная поверхностная энергия этой грани. Работа выхода определяется как
(3)
где Vf < 0 - положение дна зоны проводимости в металле (объемное значение эффективного потенциала для электронов проводимости), sF > 0 -энергия Ферми (кинетическая энергия фермиевс-ких электронов).
3 Результаты и их обсуждение
Вначале вычисления значений работы выхода
Wface (uxx) и поверхностной энергии Y face (uxx)
выполнены для ненапряженной металлической поверхности, а затем для напряженной в области упругой деформации: -0,01 < uxx <+0,01 для Ni и -0,03 < uxx < +0,03 для Al, Au, Cu, и Ti, соответственно. Положительная/отрицательная деформация uxx эквивалентна растяжению/сжатию образца в x - направлении. Верхняя грань образца предполагается упакованной как (100) или (110), (111), (0001).
Расчетные значения работы выхода и поверхностной энергии ненапряженных поверхностей Al, Au, Cu, Ni и Ti находятся в согласии с хорошо известными экспериментальными данными и расчетами других авторов (см. [1]). Деформационные зависимости указанных величин линейны по отношению к отрицательным и положительным деформациям, т.е. деформационные градиенты положи-
тельны. Величина компоненты стресса тхх меняется в интервале (1,15 -1,75 ч) у/асе (рис. 2).
Соответствующее изменение в работе выхода равно примерно 1% при максимальных растяжениях. С ростом сжатия (ихх < 0) хвосты электронного профиля и, соответственно, эффективного потенциала становятся более крутыми при убывании в вакуум, зануляясь вдали за поверхностью. При растяжении наблюдается противоположная тенденция этих величин. Полное уменьшение/увеличение работы выхода ж определяется положительным/ отрицательным сдвигом значения эффективного потенциала в объеме металла относительно деформации (пренебрегая деформационной зависимостью
еF(ихх) можно считать АЖ « -АУе^). Наши вычисления имитируют глобальную зависимость работы выхода от электронной концентрации в металлах, т.е. "переход" А1 ^ № ^ Сэ. С другой стороны, выражение (1) дает неверную зависимость | в упругой области, что, на первый взгляд, противоречит экспериментам [2-4].
Экспериментальные наблюдения могут быть объяснены не деформационным изменением дна зоны проводимости, а изменением эффективного потенциала за поверхностью металла на мнимой поверхности изображения, отстоящей на расстоянии г 0 (примерно одного радиуса Бора), т.е.
КРП = АУе^(ихх,го) В работе проведены независимые вычисления АЖ и КРП без использования выражения (1). Рис 3 демонстрирует, с одной стороны, хорошее качественное согласие вычисленных величин КПР (ихх) с экспериментальными данными [2-4], а с другой стороны - зависимость Ж(ихх), обратную той, что следует из (1). Экспериментальные значения КРП взяты из [2,3] для сжатых (I I) поликристаллов А1, Си и Аи
образцов, и из [1, 3] - для растянутых А1
Рис. 2. Расчетные значения производной ¿у/йихх для оценки поверхностного стресса (2). Левая и правая части рисунка соответствуют сжатию (ихх < 0) и растяжению (ихх > 0 ) образца, соответственно
□ -гцк(100) Д- (110) (111) О -гбц(0001) Х- експ
Рис. 3. Расчетные значения работы выхода и контактной разности потенциалов
(ихх = +0,03) и № (ихх = +0,01) образцов.
Выводы
1. По оригинальной вычислительной схеме выполнен расчет деформационных зависимостей поверхностной энергии, стресса и работы выхода А1, N Си, Аи и И. Результаты вычислений показывают, что при одноосной деформации (независимо от ее знака и от индексов кристаллографического направления) наблюдается линейный рост поверхностной энергии. Величина работы выхода убывает линейно с ростом деформации в упругой области.
2. Решена важная, с практической точки зрения, задача адекватной интерпретации результатов измерений деформационного изменения КРП по методу Кельвина. Доказано, что результаты таких измерений демонстрируют изменение потенциала поверхности, а не работы выхода.
Список литературы
1. Шпак А.П., Погосов В.В., Куницкий Ю.А. Введение в физику ультрадисперсных сред - К.: Ака-
демпериодика, 2006. - 424 с.
2. Craig P.P. Direct observation of stress-induced shifts in contact potentials // Physical Review Letters. - 1969. - V. 22, №14. - P. 700-703.
3. Левитин В.В., Лоскутов С.В., Погосов В.В. Влияние деформации и механического напряжения в металлах на работу выхода электронов // Физика металлов и металловедение. - 1990. - №9. - C. 73-79.
4. Li W., Li D.Y. Effects of elastic and plastic deformations on the electron work function of metals during bending tests // Philosophical Magazine. - 2004. - V.84, №35. - P. 3717-3727.
5. Pogosov V.V., Kurbatsky V.P. Density-functional theory of elastically deformed finite metallic sample: work function and surface stress // Журнал экспериментальной и теоретической физики. - 2001. - Т.119, №2. - C. 350-358.
Поступила в редакцию 18.05.2006 г.
Поверхневий стрес i контактна р1зниця потенц1ал1в пружно деформованих граней Al, Cu, Au, Ni i Ti кристалiв обчисленi за допомогою самоузгодженого методу Хартрi-Фока-Кона-Шема. Розрахунки демонструють зменшення^льшення роботи виходу при розтяганнi/стисканнi металевого кристала. Результати обчислень указують, що вим-iрювання методом Кельвiна контактно!'рiзницi потенцiалiв деформовано!' поверхнi в'дпо-вiдають не зм'1н'1 роботи виходу електронiв, а зм'1н'1 поверхневого значення ефективного потенцалу. Отриман i величини стресу, роботи виходу i контактно!'рiзницi потенцiалiв добре узгоджуються з результатами обчислень з перших принципiв.
The surface stress and the contact potential differences of elastically deformed faces of Al, Cu, Au, Ni, and Ti crystals are calculated using the self-consistent Hartry-Fox-Kohn-Sham method. The obtained values of the surface stress are in agreement with the results of the available first-principal calculations. We find that the work function decreases/increases linearly with elongation/ compression of crystals. Our results confirm that the available experimental data for the contact potential difference obtained for the deformed surface by the Kelvin method correspond not to the change of the work function but to the change of the surface potential.
—Отй 9яшВестникяИвигателестроенияя1 4/mD06
- 117 -
