ФИЗИКА МЕТАЛЛОВ И МЕТАЛЛОВЕДЕНИЕ, 2015, том 116, № 7, с. 709-713
СТРУКТУРА, ФАЗОВЫЕ ПРЕВРАЩЕНИЯ И ДИФФУЗИЯ
УДК 669.2122:536.77
МЕЖАТОМНЫЕ ВЗАИМОДЕЙСТВИЯ И ТЕРМОДИНАМИЧЕСКИЕ ПАРАМЕТРЫ В МАЛОКОНЦЕНТРИРОВАННЫХ ТВЕРДЫХ РАСТВОРАХ
СИСТЕМЫ Ag-Au
© 2015 г. В. И. Богданов, Л. А. Большое, Е. А. Корнейчук, В. А. Попов, С. К. Корнейчук, Д. А. Баданин
Вологодский государственный технический университет 160000 Вологда, ул. Ленина, 15 e-mail: korn62@mail.ru Поступила в редакцию 22.07.2014 г.; в окончательном варианте — 24.12.2014 г.
Важное значение в теории растворов имеют термодинамические параметры взаимодействия и эн-тальпийные параметры — коэффициенты разложения парциальных избыточных термодинамических функций в ряд по концентрациям растворенных компонентов. В приближении парного взаимодействия между атомами примеси в растворе рассматриваемые параметры могут быть вычислены теоретически с использованием методов теории функционала плотности в электронной теории сплавов. В качестве примера были выбраны твердые растворы замещения Au в Ag, образованные атомами из близких по химическим свойствам компонентов, деформационные взаимодействия в которых должны быть малы, и не нужно учитывать сложные магнитные вклады в парные потенциалы. Полная энергия малоконцентрированного раствора Au в Ag и вклады от химических и деформационных взаимодействий в потенциалы парного взаимодействия рассчитаны до седьмой координационной сферы включительно. Получено вполне удовлетворительное согласие с термодинамическими параметрами, полученными из экспериментальных данных.
Ключевые слова: химические и деформационные взаимодействия, твердые растворы замещения Аи в Ag, термодинамические параметры.
DOI: 10.7868/S0015323015070025
ВВЕДЕНИЕ
В последние два десятилетия термодинамические свойства сплавов успешно исследуются, исходя из первых принципов ("ab-initio"), на основе теории функционала плотности [1], которая позволяет найти полную энергию электронной системы в основном состоянии для фиксированных позиций атомов. Ряд методов решения уравнений Кона—Шэма [2] основан на решении уравнения Шредингера, использующего в качестве базисных функций плоские волны, по аналогии с методом псевдопотенциала и с методом PAW (projector augmented waves) [3]. Именно последние наиболее часто применяются в исследовании фазовых равновесий в сплавах из-за эффективного, точного и удобного осуществления в пакете программ VASP (Vienna Ab-initio Simulation Package) [4].
С помощью этого пакета программ, применяя метод суперячеек, появилась возможность учитывать деформационное взаимодействие в твердых растворах. Начальная концепция деформа-
ционного взаимодействия основана на модели Кривоглаза—Канзаки [5], развитой в работах Хачатуряна [6], в которой каждый примесный атом вызывает статические смещения окружающих атомов растворителя. Локальные релаксации этих атомов минимизируют возникшую упругую энергию. Так в [7] указанные первопринципные методы применены для расчетов статических смещений атомов и вкладов от химического и деформационного взаимодействия в эффективные кластерные взаимодействия в сплаве Си083 Мп0.17. Показано, что рассчитанные параметры ближнего порядка в сплаве хорошо воспроизводят результаты экспериментов по рассеянию нейтронов. Подобные же расчеты выполнены для анализа упорядоченных структур в сплавах систем №—Р1 [8] и Аи—Рё [9], опытные данные о которых получены методом диффузного рассеяния рентгеновских лучей.
В [10, 11] нами оценены термодинамические параметры взаимодействия и энтальпийные параметры для ряда твердых растворов замещения.
710
БОГДАНОВ и др.
При этом использованы экспериментальные [12] значения избыточных термодинамических функций в бинарных растворах. В данной работе на основе результатов расчетов методами теории функционала плотности, определены термодинамический параметр взаимодействия и энтальпийный параметр первого порядка в твердых растворах замещения системы Л§—Ли. Выбраны твердые растворы на основе благородных металлов, поскольку в них не нужно учитывать довольно сложные магнитные вклады в эффективные кластерные взаимодействия. Кроме того, из феноменологической теории [5, 6] следует, что в этих твердых растворах, образованных атомами из близких по свойствам компонентов, деформационные взаимодействия малы и быстро затухают. Это существенно при расчетах термодинамических параметров — необходимо, чтобы соответствующие ряды, учитывающие взаимодействия атомов примеси в твердых растворах для различных координационных сфер (КС), сходились.
I. ДЕТАЛИ ТЕОРЕТИЧЕСКИХ РАСЧЕТОВ
1.1. Полные энергии. Полные энергии суперя-чеек из 108 атомов 3 х 3 х 3(х 4) чистых металлов Ag, Au и твердых растворов замещения были рассчитаны методом PAW, реализованном в пакете программ VASP. Для учета эффектов обмена и корреляции в электронном газе обычно используют приближения локальной плотности (Local Density Approximation) LDA [13] и градиентное приближение (Generalized gradient Approximation) GGA [14]. Тщательный анализ роли этих двух приближений в первопринципных исследованиях чистых ГЦК благородных металлов показал [15], что LDA следует отдать предпочтение при расчете свойств основного состояния в Ag и, особенно, в Au. Поэтому расчеты полных энергий чистых Ag и Au и их твердых растворов были выполнены, учитывая эффекты обмена и корреляции в электронном газе в приближении локальной плотности. Расчет полной энергии в релаксированном состоянии производился с помощью итерационной процедуры до тех пор, пока силы Геллмана-Фейн-мана не становились меньше 0.01 эВ/А. Объем и геометрия суперячейки не менялись. Метод Монк-хорста—Пака [16] использован для выбора специальных k-точек в зоне Бриллюэна. Электронные заселенности определялись в соответствии со схемой [17]. Метод PAW был также применен и для расчета энергии локальных релаксаций и энергий деформационного взаимодействия.
1.2. Атомный конфигурационный гамильтониан и эффективные кластерные взаимодействия. Будем исходить, как и в [18] (см. также [19]) из того, что конфигурационная энергетика твердых растворов
замещения может быть изучена в модели Изинга. Конфигурационный гамильтониан имеет вид:
н=1 xxvp^axj
p tj^p
1ZZ V^bfixj8xk +...,
(1)
t
где I, ], к — индексы, определяющие номер узла,
УР2) и V((3) — эффективные двух- и трехчастичные взаимодействия атомов примеси замещения в координационных ячейках р и дx¡ = XI - (х^ — отклонение чисел заполнения х1 от их среднего значения (х^ = х; XI принимает значение 1, если узел I занят атомом примеси и 0, если он занят атомом растворителя. Как в [18], будем предполагать, что
химические V{2)~x и деформационные Vp2)~8 взаимодействия доминируют, и что результирующие эффективные взаимодействия определяются их
суммой Vр = Vр + Vр .
Химический вклад в
потенциалы парного взаимодействия связан с изменением конфигурационной энергии, обусловленным перераспределением атомов в идеальной (неискаженной) решетке. Статические смещения атомов растворителя, понижающие полную энергию твердого раствора, дают дополнительный вклад в потенциалы межатомного взаимодействия, рассматриваемый как деформационное взаимодействие. В [20] обсуждались два метода,
которыми могут быть рассчитаны Vp2hx и Vp2)~д. Если парные взаимодействия преобладают, то наиболее прямой путь для вычисления химического вклада в эти взаимодействия при малых концентрациях — это расчет разности полных энергий двух больших суперячеек: в одной атомы примеси (в нашем случае атомы Ли) находятся в позициях, отстоящих друг от друга на расстоянии, равном радиусу рассматриваемой координационной сферы, а в другой — эти атомы достаточно далеко расположены друг от друга, так что их взаимодействием можно пренебречь (метод I). Решеточные позиции должны быть фиксированы в обоих расчетах. Если же размеры суперячейки недостаточны, чтобы обеспечить изоляцию атомов примеси, тогда эффективные кластерные взаимодействия могут быть получены в эквивалентных расчетах трех суперячеек:
VP2)-x = E%p + Eо - 2Ei,
(2)
Здесь Е0, Ех и Е2,р — полные энергии суперячеек (одинакового размера) без атомов примеси, с одним атомом примеси и двумя атомами примеси, один из которых находится в р-той КС другого (метод II).
Следует отметить два недостатка подхода, основанного на применении формулы (2) для расчета эффективных кластерных взаимодействий. Первый связан с тем, что данный подход определяет взаимодействия для растворов в пределе малых концентраций, в то время как взаимодействия, определяемые из (2), основаны на расчете энергий систем, имеющих различные составы. Однако в пределе малых концентраций, когда размер суперячейки достаточно велик, можно пренебречь небольшой разностью концентраций. Второй очевидный недостаток связан с тем, что рассчитываемые взаимодействия атомов примеси, расположенных в различных КС, ограничены естественно размерами суперячейки.
Деформационный вклад в эффективные взаимодействия был рассчитан как разность энергий релаксации:
V?-8 = (< - ЕХр) - 2 (ЕГ1 - Е), (3)
где Е2,р и Е — энергия суперячейки с двумя атомами примеси, расположенными на расстоянии, соответствующем ^-той КС, и с одним атомом примеси после релаксации кристаллической решетки.
Поскольку на 5 и 7 КС для выбранной нами су-перячейки взаимодействия учитываются дважды, то в формулах (2) и (3) для этих КС результат необходимо брать в два раза меньшим.
2. РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
В качестве теста применительно к чистым металлам (Лё и Аи) были рассчитаны параметры решетки а и объемные модули В в основном состоянии в приближении LDA. Для Лё мы нашли а = 4.02 А, В = 1.35 х 1011 Н/м2; для Ли: а = 4.08 А, В = 2.56 х 1011 Н/м2. Параметры решетки в основном состоянии, рассчитанные в этих металлах подобным же методом [21], соответствуют нашим результатам: а = 4.02 А (Лё) и 4.08 А (Ли). Экспериментальные данные при комнатной температуре [22] для Лё: а = 4.09 А, В = 1.09 х 1011 Н/м2; для Ли: а = 4.08 А, В = 1.80 х 1011 Н/м2. Сравнение теории с экспериментом показывает, что приближение LDA несколько недооценивает экспериментальные параметры решетки и переоцен
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.