Фазы и амплитуды рассеяния электронов для задач определения локальной атомной структуры методом EXAFS-спектроскопии

Размещено на http://www.allbest.ru/

Фазы и амплитуды рассеяния электронов для задач определения локальной атомной структуры методом EXAFS-спектроскопии

Л.А. Авакян, Н.В. Булат

Аннотация

В работе представлены результаты разработки и тестирования программного комплекса HFampl, предназначенной для расчёта фаз и амплитуд рассеяния из первых принципов с использованием точного обменного потенциала Хартри-Фока (ХФ). Разработка с использованием компиляторов с открытым исходным кодом позволяет запускать программу в любых совместимых операционных системах (Windows, Linux, Mac). Тестирование комплекса на металлах, оксидах и железосодержащих водных сульфатных растворах показало хорошую применимость рассчитываемых фаз и амплитуд рассеяния для проведения высокоточного структурного анализа. Разработанное программное обеспечение предназначено для использования при проведении научных исследований атомного строения неупорядоченных, аморфных структур, нанообъектов, а также в учебном процессе для изучения современных методов структурного анализа сложных, многокомпонентных соединений.

Ключевые слова: рассеяние фотоэлектрона, протяженная тонкая структура рентгеновских спектров поглощения, метод Хартри-Фока, межатомное расстояние, обработка экспериментальных данных

Введение

Метод EXAFS-спектроскопии применяется для изучения структуры на атомном уровне самых разных объектов: активных центров катализаторов [1, 2], фотосинтезирующих систем [3], металлсодержащих растворов [4], аморфных материалов [5] и многих других. Успешное проведение структурного анализа неупорядоченных, аморфных материалов и нанокомпозитов с помощью спектров рентгеновского поглощения требует наличия разработанных алгоритмов и программного обеспечения для расчета надежных значений фаз и амплитуд рассеяния электронов атомами изучаемого вещества [6]. Важным условием для расчета рассеивающих свойств атомов в веществе является использование по возможности более точной модели потенциала, которая, в то же время, должна быть достаточно универсальной для обеспечения трансферабельности фаз и амплитуд рассеяния электронов на атомах в разных материалах.

Существующие подходы анализа EXAFS [7-9] используют фазы и амплитуды рассеяния, полученные в рамках самосогласованного межатомного потенциала с обменной частью, аппроксимируемой потенциалом Хедина-Лунквиста. Использование процедуры самосогласования приводит к проявлению зависимости не только от радиуса координационной сферы, но и от координационных чисел и симметрии расположения атомов, что приводит к необходимости пересчета рассеивающих характеристик в процессе анализа. В некоторых случаях [9] при проведении структурного анализа с фазами и амплитудами рассеяния, полученными с обменным потенциалом Хедина-Лунквиста, значения фактора неупругих потерь фотоионизации S₀² (спектроскоп-фактора) оказываются больше 1, что, вообще говоря, является не физическим результатом.

Альтернативный подход был сформулирован в работе [10], в которой изучалось формирование EXAFS спектров в оксидах титана. В этом подходе не использовалась процедура самосогласования межатомного потенциала, а в качестве обменного потенциала использовался точный Хартри-Фоковский обмен, что позволяло получать более разумные значения S₀². В настоящей работе расширен круг объектов, для которых показана применимость подхода расчета фаз и амплитуд рассеяния. Было проведено усовершенствование программного обеспечения, а также сформированы базисные конфигурации большинства атомов периодической системы. Используемый в работе метод построения потенциала, в котором движется фотоэлектрон, основан на том, что потенциал многоатомной системы можно представить в виде суммы локального кулоновского потенциала и нелокального обменного оператора [11]. Задача построения потенциала для получения трансферабельных фаз рассеяния решается при использовании двух приближений. Первое состоит в переходе к MT-аппроксимации [12] для потенциала многоатомной системы. Второе приближение заключается в использовании дополнительного предположения относительно обменной части ХФ-потенциала, состоящего в том, что его можно рассчитывать независимо в разных атомных центрах. Результаты расчетов как околопороговой (XANES), так и протяженной (EXAFS) областей спектров рентгеновского поглощения атомов в кристаллических и аморфных соединениях, обладающих различными физическими свойствами [13-15], убедительно свидетельствуют об адекватности такого приближения.

Схема работы программы HFampl

Программный комплекс состоит из нескольких компонент, взаимосвязь между которыми проиллюстрирована на Рисунке 1. В него входят: программа расчета волновых функций и логарифмических производных в произвольном энергетическом интервале ELWF; программа расчета характеристик рассеивания для пар атомов XAFSS в заданном интервале энергий; программа формирования выходных файлов в различных форматах DataMaker. Сложность взаимодействия между отдельными частями компенсируется наличием графической оболочкой HFampl, настраивающей и запускающей отдельные модули в нужном порядке. Такая архитектура позволяет проводить расчет фаз и амплитуд рассеяния в полуавтоматическом режиме без изменения стандартных значений параметров.

Рис. 1. Схема взаимодействий между отдельными модулями программного комплекса HFampl.

В первом модуле, а именно, в программе ELWF, расчет волновых функций электронов и соответствующих логарифмических производных осуществляется для изолированного атома, что соответствует ситуации внутри МТ-сферы, после перехода к системе атомов.

Входными параметрами, передаваемые через текстовые a-файлы для каждого энергетического интервала, являются: 1) тип обменно-корреляционного взаимодействия; 2) Радиус MT-сферы (R_MT), на границе которой рассчитывается логарифмическая производная; 3) выбор величины среднего межатомного потенциала - т.н. МТ-нуля (E_MT); 4) электронные конфигурации поглощающего и рассеивающих атомов в веществе.

В случае произвольно выбранных значений радиусы атомных сфер (R_MT) и EMT возможно возникновение скачка межатомного потенциала (Рис. 2) на поверхности МТ-сфер, что приводит к появлению не физического надбарьерного рассеяния, сильно искажающего структуру спектра в области энергий в несколько десятков электрон-вольт за краем поглощения. Чтобы избежать этого эффекта следует выбирать R_MT так, чтобы, с одной стороны, обеспечить достаточное перекрытие потенциалов соседних атомов, а с другой - сделать одинаковыми величины скачков на поверхностях всех МТ-сфер. В таком случае, подобрав единым образом энергию МТ-нуля (E_мт) относительно вакуумного уровня (V.L.) атомов (Рис. 2), можно устранить скачки МТ-потенциала на поверхностях всех атомных сфер.

Рис. 2. Схематическое изображение энергетического скачка в потенциале атома, устраняемого в разработанном программном обеспечении выбором величин E_MT = -V(R_MT) при расчете фаз рассеяния фотоэлектронных волн.

При выборе электронных конфигураций поглощающего и окружающих его атомов, для расчета волновых функций и фаз рассеяния электрона учтен перенос экранирующего электронного заряда на поглощающий атом, содержащий остовную дырку, обусловленный процессами сверхатомной релаксации. В металлах и ковалентных соединениях такой процесс ведет к увеличению на единицу числа валентных электронов на ионизируемом атоме с остовной дыркой. Для атомов окружения используются конфигурации нейтральных атомов.

Программный пакет HFampl содержит предварительно сформированные входные файлы ELWF для большинства элементов периодической системы, так что пользователю достаточно выбрать нужный элемент из списка (Рис. 3а).

Рис. 3. Экранные снимки программы-оболочки HFamp: а) окно атомного расчёта; б) главное окно - сборка файла фаз и амплитуд рассеяния.

В основе используемой модели MT-потенциала с ХФ-обменом между фотоэлектроном и остовными электронами рассеивающих атомов лежит идея о том, что вычисленные на радиусах MT-сфер логарифмические производные радиальной волновой функции, которые определяют энергетическое поведение фаз рассеяния электрона, сдвигаются как целое по энергии, практически не меняя своей формы при переходе от изолированного атома к атому в молекуле или твёрдом теле. Такой подход открывает возможность определения функций из ХФ-расчёта для изолированного атома, а полученные функции используются затем при построении кристаллического MT-потенциала для расчета фаз рассеяния, которые определяются из условия “сшивки” на радиусе MT-сферы атомных волновых функций с волновыми функциями в постоянном потенциале области между сфер. На основе выбранных приближений, параметров изучаемой системы и электронных конфигураций атомов рассчитываются фазы и амплитуды рассеяния электронов на атомах изучаемого соединения (Рис. 3б).

Тестирование программного обеспечения

Тестирование программного комплекса осуществлялось путем выполнения Фурье-анализа экспериментальных EXAFS спектров соединений с известной структурой, такие как металлические медь, серебро и никель; объемный оксид цинка (ZnO); а также водные сульфатные растворы железа с различной степенью окисления (Fe²⁺ и Fe³⁺). Особый интерес представляют именно водные растворы, поскольку, с одной стороны, они наиболее важны с точки зрения практического применения в медицине и биологии, а с другой стороны их экспериментальное изучение затруднительно и ограничивается рентгеновской спектроскопией [16, 17] и ЯМР методами. Наноразмерный оксид цинка является перспективным материалом для детекторов газов [18]. Экспериментальные спектры были получены из открытых источников (Cu, ZnO), а также измерены на российских и европейских синхротронных центрах (Ni, Ag, Fe-aqua). Предварительная обработка большинства спектров и выделение осциллирующей EXAFS функции ч(k) проводилась стандартным образом, с помощью программы Athena [19]. В случае Fe K-XAFS спектров энергетический интервал ограничен 120 эВ над краем поглощения, что делает стандартный подход неприменимым. В этом случае использовался оригинальный подход, описанный в работе [20]. Фурье-преобразование и последующий фитинг теоретических и экспериментальных кривых проводился с помощью программы Larch [21], с использованием фаз и амплитуд рассеяния рассчитанных в HFampl для пар атомов, находящихся на расстоянии, отличающимся от дифракционного не более чем на 0.1 Е.

Рис. 4. Сопоставление модулей экспериментальных (сплошные кривые) и теоретических кривых (пунктир), полученных в результате Фурье-анализа Zn K-EXAFS оксида цинка: а) Фурье-образы |F(R)| функции ч(k); б) функции ч(k) с k-весом 1.

Таблица 1. Значения структурных параметров первой координационной сферы атомов, полученных в результате Фурье-анализа EXAFS спектров модельных соединений.

*Данные соответствуют кристаллической фазе

Согласие теоретических и экспериментальных кривых, полученных в результате фитинга, иллюстрируется на Рисунке 4 для случая серебряной фольги. Результаты Фурье-анализа, представленные в Таблице 1, демонстрируют хорошее согласие полученных данных с дифракционными, что позволяет сделать вывод об адекватности сделанных приближений и применимости получаемых фаз амплитуд рассеяния для анализа структуры материалов.

Преимущество выбранного потенциала ХФ иллюстрируется на примере анализа K-спектра поглощения меди в металлической меди. В Таблице 2 структурные параметры первой координационной сферы поглощающего атома меди в кристалле, полученные с помощью фаз и амплитуд рассеяния, рассчитанных на основе волновых функций электронов для разных форм обменно-корреляционного потенциала, сравниваются с известными дифракционными данными. Из представленных данных следует, что использование ХФ амплитуд позволяет получить более правильное значение амплитудных параметров (S₀², N) по сравнению с традиционно используемым потенциалом Хедина-Лундквиста.

рассеяние электрон аморфный железосодержащий

Таблица 2. Значения структурных параметров первой координационной сферы атома меди в металлической меди.

Заключение

Разработан программный комплекс, позволяющий рассчитывать рассеивающие характеристики фотоэлектрона на атомах в соединении. Тестирование на металлах, оксидах и железосодержащих водных растворах показало применимость рассчитываемых фаз и амплитуд рассеяния для проведения высокоточного структурного анализа. Ошибка в определении межатомных расстояний составляет порядка 1 % при разумных физических значениях спектроскоп фактора (0 < S₀² < 1) и параметра Дебая-Валлера (у² ~ 0.005 Е).

Литература

1. Deng D., Chen X., Yu L., Wu X., Liu Q., Liu Y., Yang H., Tian H., Hu Y., Du P., Si R., Wang J., Cui X., Li H., Xiao J., Xu T., Deng J., Yang F., Duchesne P.N., Zhang P., Zhou J., Sun L., Li J., Pan X., Bao X. // Science Advances. 2015. V.1. P.e1500462

2. Ding K., Gulec A., Johnson A.M., Schweitzer N.M., Stucky G.D., Marks L.D., Stair P.C. // Science American. 2015. V. 350. pp.189-192.

3. Sun L. // Science. 2015.V.348. pp.635-636

4. Sadeghi O., Zakharov L.N., Nyman M. // Science. 2015. V.347 P.1359-1362

5. Gordon L.M., Cohen M.J., MacRenaris K.W., Pasteris J.D., Seda T., Joester D. // Science. 2015. V.347. pp.746-750.

6. Koningsberger D., Prins R. X-ray absorption: principles, applications, techniques of EXAFS, SEXAFS, and XANES. New York:Wiley, 1988. 688 p.

7. Ankudinov A.L., Ravel B., Rehr J.J., Conradson S.D. // Physical Review B. 1998. V.58. pp. 7565-7576.

8. Di Cicco A., Filipponi A. GNXAS Extended Suite of Programs for Advanced X-ray Absorption Data-analysis: Methodology and Practice. Task, 2009. 124 p.

9. Rehr J.J., Kas J.J., Vila F.D., Prange M.P., Jorissen K. // Physical Chemistry Chemical Physics 2010. V.12. pp. 5503-5513.

10. Bugaev L., Vedrinskii R., Levin I. // Physica B: Condensed Matter . 1989. V.158. pp. 378-382.

11. Марч Н., Янг У., Сампантхар С. Проблема многих тел в квантовой механике. M.: Мир, 1969. 496 c.

12. Займан Д., Бонч-Брюевич В. Принципы теории твердого тела: Пер. с англ. M.: Мир, 1974. 472 c.

13. Bugaev L.A., Gegusin I.I., Datsyuk V.N., Novakovich A.A., Vedrinskii R.V. // Physica Status Solidi (B) 1986. V.133. pp. 195-202.

14. Bugaev L.A., Shuvaeva V.A., Zhuchkov K.N., Rusakova E.B., Vedrinskii R.V. // Journal of Synchrotron Radiation 2001. V.8. pp. 308-310.

15. Bugaev L., Farges F., Rusakova E., Sokolenko A., Latokha Y., Avakyan L. // Physica Scripta 2005. V.115. pp. 215-217.

16. Положенцев О.Е., Гуда А.А., Сафонова О.В., Ван Бокховен Д.А. и Солдатов А.В. Методика in-situ исследования катализаторов с помощью спектроскопии рентгеновского поглощения // Инженерный вестник Дона, 2013, №2 URL: ivdon.ru/en/magazine/archive/n2y2013/1649.

17. Кочкина В.К., Положенцев О.Е., Мазалова В.Л. Синтез и спектроскопическое исследование структуры и окислительной способности наночастиц Fe₃O₄ магнитной жидкости // Инженерный вестник Дона, 2015, №1 URL: ivdon.ru/en/magazine/archive/n1y2015/2576.

18. Пушкарев В.И., Лянгузов Н.В., Кайдашев Е.М. Низкотемпературный термический синтез наностержней оксида цинка из паров цинка без катализатора // Инженерный вестник Дона, 2013, №4 URL: ivdon.ru/en/magazine/archive/n4y2013/1884.

19. Ravel B., Newville M. // Journal of Synchrotron Radiation 2005. V.12. pp. 537-541.

20. Bugaev L.A., Bokhoven J.A., Khrapko V.V., Avakyan L.A., Latokha V.J. // Journal of Physical Chemistry B 2009. V.113. pp. 4614-4618.

21. Newville M. // Journal of Physics: Conference Series 2013. V.430. P. 012007.

Размещено на Allbest.ru