Кинетика образования и начальных стадий роста новой фазы при электроосаждении из расплавленных солей тема автореферата и диссертации по химии, 02.00.05 ВАК РФ
Гришенкова, Ольга Владимировна
АВТОР
|
||||
кандидата химических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Екатеринбург
МЕСТО ЗАЩИТЫ
|
||||
2012
ГОД ЗАЩИТЫ
|
|
02.00.05
КОД ВАК РФ
|
||
|
На правах рукописи
Гришенкова Ольга Владимировна
КИНЕТИКА ОБРАЗОВАНИЯ И НАЧАЛЬНЫХ СТАДИЙ РОСТА НОВОЙ ФАЗЫ ПРИ ЭЛЕКТРООСАЖДЕНИИ ИЗ РАСПЛАВЛЕННЫХ СОЛЕЙ
Специальность 02.00.05 - Электрохимия
АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата химических наук
1 5 мдр 20 ¡2
Екатеринбург - 2012
005013865
Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте высокотемпературной электрохимии Уральского отделения Российской академии наук
Научный руководитель: доктор химических наук, старший научный сотрудник Исаев Владимир Александрович
Официальные оппоненты: Останина Татьяна Николаевна,
Ведущая организация: ФГБОУ ВПО «Кабардино-Балкарский
государственный университет им. Х.М. Бербекова», г. Нальчик
Защита состоится 28 марта 2012 г. в 14 часов на заседавши диссертационного совета Д 004.002.01 в Институте высокотемпературной электрохимии УрО РАН по адресу: г. Екатеринбург, ул. Академическая, 20, конференц-зал.
Ваши отзывы в двух экземплярах, подписанные и заверенные гербовой печатью с датой подписания, просим высылать по адресу: 620990, г. Екатеринбург, ул. Академическая, 20. Институт высокотемпературной электрохимии УрО РАН. Ученому секретарю диссертационного совета Кулик Нине Павловне, E-mail: N.P.Kulik@ihte.uran.ru. Факс: (343)3745992
С диссертацией можно ознакомиться в Центральной библиотеке УрО РАН
Автореферат разослан « февраля 2012 г.
Ученый секретарь диссертационного совета,
доктор химических наук, профессор, ФГАОУ ВПО «Уральский федеральный университет имени первого Президента России Б.Н. Ельцина», профессор
Елшина Людмила Августовна, кандидат химических наук, старший научный сотрудник, Институт высокотемпературной электрохимии УрО РАН, старший научный сотрудник
кандидат химических наук
Кулик Н.П.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы
Изучение начальных стадий фазообразования в электрохимических системах является актуальной задачей фундаментальной науки, так как позволяет определить механизмы и лимитирующие стадии зарождения и роста новой фазы, получить прямую информацию о кинетике перехода частиц через межфазную границу электролит / кластер.
Критические кластеры новой фазы при электрокристаллизации имеют наноразмеры, поэтому исследование электрохимического зародышеобразования дает возможность установить закономерности образования и роста нанокластеров, оценить размерные эффекты. В настоящее время в связи с широким внедрением нанотехнологий актуальность данной тематики особенно возросла.
Результаты исследования начальных стадий электрокристаллизации в расплавленных солях могут быть использованы и при решении прикладных задач, связанных с определением условий формирования электродных осадков и покрытий с заданной структурой.
Работа выполнена в соответствии с планами НИР Института высокотемпературной электрохимии УрО РАН "Разработка научных основ безотходных экологически чистых, энерго- и ресурсосберегающих высокотемпературных электрохимических технологий получения металлов и неорганических соединений в солевых расплавах" (№ гос. регистрации 01.2.007 01883), "Физико-химические основы электрохимических технологий рафинирования и получения металлов и их соединений в расплавах солей" (№ гос. регистрации 01.2.010 00812), а также в рамках проектов: РФФИ-"УРАЛ" № 04-0396124 р2004урал_а "Электрохимическое образование волокнистых осадков металлов с заданной каталитической активностью из расплавов и водных растворов"; совместного проекта УрО РАН-СО РАН № 81 "Исследование процессов перехода веществ и систем между различными агрегатными состояниями" (2006-2008 г.г.); программ Президиума РАН "Разработка методов
получения химических веществ и создание новых материалов: создание новых конструкционных и функциональных материалов на базе нанотехнологий" по теме "Исследование условий получения нанопорошков и нановолокон из серебра (модельный металл) электролизом растворов и низкотемпературных серебросодержащих солевых расплавов" (2006-2008 г.г.) и "Основы фундаментальных исследований нанотехнологий и наноматериалов" по проекту 09-П-3-1026 "Получение нанокристаллических осадков электролизом солевых расплавов" (2009-2011 г.г.).
Цель работы
Определение основных закономерностей образования и роста кластеров новой фазы в области их нано- и микроразмеров при электроосаждении из расплавленных солей в потенцио- и гальваностатических условиях.
Задачи
1. Теоретически изучить закономерности зарождения кластеров новой фазы при постоянном пересыщении (в потенциостатических условиях).
2. Теоретически исследовать кинетику роста кластеров на начальных этапах гальваностатического электроосаждения из расплавленных солей в модельной системе (осаждаемый металл-А§; расплав- А£1МОз, А^КОз-КЫОз-ЫаШ3; катод - 1г, Р1).
3. Разработать компьютерную модель этого процесса, сравнить результаты расчета с экспериментальными данными, оценить возможность применения предложенной модели для анализа начальных стадий фазообразования в расплавленных солях.
4. Использовать компьютерный эксперимент для определения механизма, лимитирующих стадий и основных параметров электрокристаллизации в исследуемых системах.
5. Разработать математическую модель образования и роста нитевидных кристаллов при электроосаждении из расплавленных солей, провести компьютерное моделирование процесса. !
Научная новизна
1. Получено новое решение основного кинетического уравнения, описывающего процесс зародышеобразования (уравнения Фоккера-Планка), для случая электрохимического потенциостатического фазообразования в расплавленных солях. Найдены выражения для скорости зародышеобразования, учитывающие лимитирующие стадии образования критических кластеров новой фазы при электроосаждении из расплава индивидуальной соли и расплава, содержащего избыток фонового электролита.
2. Разработана компьютерная модель единичного и множественного зародышеобразования в гальваностатических условиях электроосаждения из расплавленных солей.
3. Впервые получено количественное подтверждение адекватности предложенной модели гальваностатического фазообразования в расплавленных солях. Оценены основные параметры электрокристаллизации в исследуемых системах: удельная емкость двойного электрического слоя, концентрация адатомов (мономеров) на электроде, плотность тока обмена на границе кристалл / расплав, коэффициент диффузии осаждаемых ионов в объеме расплава.
4. Разработана модель роста нитевидных кристаллов новой фазы, формирующихся на катоде при гальваностатическом электроосаждении из расплавленных солей.
Научная ценность и практическая значимость
Результаты работы и сделанные в ней обобщения представляют теоретический интерес для развития представлений о закономерностях начальных стадий электрохимического фазообразования. Аналитические методы решения и предложенный метод компьютерного моделирования могут быть использованы для определения механизма и лимитирующих стадий процесса, а также основных параметров электрокристаллизации металлов в расплавленных солях и, следовательно, будут полезны при выборе оптимальных условий получения осадков с заданными свойствами.
На защиту выкосятся:
1. Решение основного кинетического уравнения (уравнения Фоккера-Планка), описывающего процесс зародышеобразования в потенциостатических условиях электроосаждения из расплавленных солей. Уравнения для расчета скорости зародышеобразования, учитывающие лимитирующие стадии образования критических кластеров новой фазы при электрокристаллизации.
2. Результаты компьютерного моделирования образования и начальных стадий роста кристаллов серебра при гальваностатическом электролизе солевых расплавов.
3. Оценка основных параметров электрокристаллизации в исследуемых системах: удельной емкости двойного электрического слоя, концентрации адатомов на электроде, плотности тока обмена на границе кристалл / расплав, коэффициента диффузии осаждаемых ионов в объеме расплава.
4. Математическая модель образования и роста на электроде нитевидных кристаллов, результаты компьютерной реализации данной модели.
Личный вклад автора
Теоретическое исследование фазообразования при постоянном и переменном пересыщении в системе, вывод уравнений, компьютерное моделирование, в том числе разработка оригинальных программ расчета, обработка экспериментальных данных, анализ, выводы и оформление результатов выполнены лично автором. Постановка задач осуществлялась научным руководителем д.х.н. Исаевым В.А. Экспериментальное исследование электрохимического зарождения и роста кристаллов серебра было выполнено к.х.н. Семериковой О.Л.
Апробация работы
Основные результаты и положения исследований были представлены на российских и международных конференциях и семинарах: "Проблемы электрокристаллизации металлов" и "Современные аспекты электрокристаллизации металлов" (Екатеринбург, 2000 и 2005); X Кольском семинаре по электрохимии редких металлов (Апатиты, 2000); XII, XIII, XIV, XV
Российской конференциях по физической химии расплавленных и твердых электролитов (Нальчик-Екатеринбург, 2001, 2004, 2007, 2010); EUCIIEM Molten Salts Conferences (Oxford, UK, 2002 и Piechowice, Poland, 2004); втором семинаре CO РАН-УрО РЛН "Новые неорганические материалы и химическая термодинамика" (Екатеринбург, 2002); конференции, посвященной 80-летию кафедры ТЗХП УГТУ-УПИ (Екатеринбург, 2003); Всероссийской конференции "Химия твердого тела и функциональные материалы-2004" (Екатеринбург, 2004); 55th Annual Meeting of ISE (Thessaloniki, Greece, 2004); VI Международной научной конференции "Кинетика и механизм кристаллизации. Самоорганизация при фазообразовании" (Иваново, 2010).
Публикации
По результатам исследований опубликовано 13 статей в российских и международных журналах из списка ВАК и 18 тезисов докладов в материалах российских и международных конференций.
Структура и объем диссертации
Диссертационная работа состоит из введения, трех глав, обобщающих выводов и списка использованных источников. Диссертация изложена на 138 страницах, содержит 34 рисунка, 2 таблицы; в списке использованных источников 174 работы отечественных и зарубежных авторов.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность работы, сформулированы цели и задачи исследования, научная новизна и практическая значимость полученных результатов, перечислены положения, выносимые на защиту.
В первой главе рассмотрены основные аспекты энергетики и кинетики электрохимического зародышеобразования. Сделан обзор литературы по вопросу определения работы образования кластеров (зародышей) в рамках классического и атомистического подходов, обсуждены границы их корректного использования. Показано, что закономерности электрохимического фазообразования в расплавленных солях обычно адекватно описываются классической теорией
зарождения. Обсуждены механизмы нуклеации, реализующиеся в различных условиях электрокристаллизации и факторы, снижающие работу образования кластеров новой фазы на электроде.
Проведен кинетический анализ процесса зародышеобразования. Показано, что с позиций классической теории дискретный подход Беккера-Дёринга и континуальная модель нуклеации (уравнение Фоккера-Планка) для стационарной скорости зарождения Зст дают одинаковый результат. Сравнены зависимости от пересыщения (перенапряжения) в классической и атомистической теориях. Описаны основные экспериментальные методики исследования стационарного электрохимического зарождения и способы их интерпретации с привлечением различных модельных представлений для определения работы и констант зародышеобразования, размера критического кластера.
Теоретически исследована проблема нестационарного электрохимического зарождения. Решено уравнение Фоккера-Планка для случая замедленного перераспределения кластеров по размерам на начальных стадиях потенциостатического электроосаждения из расплавленных солей. Исходное уравнение путем замены переменных было сведено к уравнению Фика с начальным, и граничными условиями, соответствующими случаю диффузии в конечном слое. Получено выражение для нестационарной скорости зародышеобразования ./(/).
Показано, что лимитирующие стадии образования критических кластеров новой фазы при электрокристаллизации могут быть определены по вкладам стадий перехода и доставки в кинетический коэффициент уравнения для Д?) -коэффициент диффузии в пространстве размеров зародышей Вкр. Установлено, что при электроосаждении из индивидуального расплава нитрата серебра лимитирующей стадией зародышеобразования является разряд ионов на поверхности кластера. В случае электроосаждения серебра из нитратного расплава, содержащего избыток фонового электролита, процесс зарождения контролируется диффузией осаждаемых ионоег в объеме расплава и их разрядом на поверхности кластера.
Полученное решение уравнения Фоккера-Планка и анализ лимитирующих стадий позволили найти выражения для расчета стационарной скорости зарождения Jcm и индукционного времени г при электроосаждении из индивидуальных расплавов (лимитирующая стадия - разряд):
где и - объем, приходящийся на одну частицу в новой фазе; а - межфазное натяжение границы кластер/расплав; я - валентность иона; е - заряд электрона; к - постоянная Больцмана; Т - температура; ¡0 - плотность тока обмена на межфазной границе кластер/расплав; Га - концентрация одиночных адатомов, находящихся на активных центрах электродной поверхности при нулевом перенапряжении; / = ге/кТ\ г] - перенапряжение; К2 - константа, определяемая работой зародышеобразования; с0 - концентрация осаждаемых ионов в глубине электролита; П - коэффициент диффузии ионов в расплаве.
Показано, что продолжительность индукционного периода, определяемая в эксперименте по потенциостатическому электроосаждению из расплавленных солей, вероятнее всего, связана с активированием поверхности электрода.
В главе 2 представлены результаты теоретического исследования и компьютерного моделирования роста единичных кристаллов серебра при электроосаждении из расплавленных солей в гальваностатических условиях и их сравнение с экспериментальными данными.
Характер процесса образования и роста новой фазы в гальваностатических условиях электроосаждения значительно усложняется по сравнению с фазообразованием в потенциостатическом режиме. Изменение пересыщения в
(1)
и расплавов с избытком фонового электролита (разряд + диффузия)
(2)
системе приводит к нестационарности процесса. Заряжение и разряжение двойного электрического слоя, изменение концентрации адатомов на электроде, изменение условий доставки вещества к поверхности катода и растущих кластеров, перераспределение тока на электроде при появлении сверхкритических кластеров - факторы, существенно затрудняющие теоретическое описание гальваностатического фазообразования. Кроме того нужно учитывать, что кинетика зарождения и кинетика роста кластеров взаимосвязаны: потребление тока образующимися кластерами ведет к изменению перенапряжения, а значит и скорости появления новых зародышей.
Представлен литературный обзор основных результатов теоретического и экспериментатьного изучения электрохимического фазообразования в гальваностатических условиях. Обсуждены возможности использования различных вариантов компьютерных моделей для получения информации о зародышеобразовании при изменяющемся со временем пересыщении. Показано, что метод компьютерного моделирования, основанный на численном решении системы дифференциальных уравнений, описывающих изменение параметров электрохимической системы в ходе образования и роста новой фазы, позволяет установить основные закономерности процесса.
Проанализированы особенности кинетики преодоления ионами межфазной границы электролит/полусферический кластер новой фазы. Показано, что влияние кривизны поверхности кластера на кинетику разряда наиболее ощутимо при размерах, близких к критическому гкр = 2аи/гет], то есть в области наноразмеров:
где г'*, г'г - плотности катодного и анодного токов для кластера, состоящего из g атомов; а, /3 - коэффициенты переноса заряда, а + /? = 1; г - радиус кластера. Кривизна поверхности кластера сказывается .и на стадии доставки ионов к межфазнсй границе, что проявляется в изменении закономерностей диффузии или
£=г0ехраг/?7
(3)
и
миграции, а также в существовании зависимости концентрации осаждаемых ионов вблизи поверхности кластера и величины омического падения потенциала в слое электролита, прилегающем к нему, от размера кластера.
Для малых кластеров (10"7ч-10'4 см) плотности токов роста
сопоставимы по величине с плотностью тока обмена на границе кластер/расплав, что позволяет выделить сопротивление разряда на фоне сопротивления, связанного с доставкой ионов к зародышу, а значит оценить ц и определить лимитирующие стадии роста.
Представлены результаты исследования роста единственного кластера новой фазы (Лё) на индифферентной подложке при электроосаждении из нитратных расплавов с различной концентрацией осаждаемых ионов (расплав индивидуальной соли и расплав, содержащий избыток фонового электролита) в гальваностатических условиях.
Моделирование электрохимического роста кластера, аппроксимированного полусферой, проводили в рамках классической теории нуклеации. Считали, что за время электроосаждения не происходит встраивания одиночных адатомов в электрод; образование и рост кластера осуществляются путем прямого присоединения осаждаемых ионов из объема электролита к поверхности кластера. Математическая модель процесса учитывала баланс тока на электроде, изменение концентрации адатомов при изменении перенапряжения, особенности кинетики разряда на границах расплав/кластер, расплав/электрод и закономерности стадии массопереноса в электролите.
В расплаве индивидуальной соли доставка вещества к электроду и растущему кластеру реализуется за счет миграции ионов в электрическом поле. Процесс роста единичного кластера в этом случае описывается системой уравнений
¿■д _ ¿-2я-г2;г/5 ■ • 71~71Р
Сд+ге/Г0ехр/г}' гег' ** ~ (\/Ла) + (г/к)'
где Г - время осаждения; г, гг - плотности внешнего тока и тока, идущего на рост кластера; 5 - площадь электрода; Сд - удельная емкость двойного
электрического слоя; Г0 - равновесное (начальное) значение концентрации адатомов; rjp - фазовое перенапряжение, при котором зародыш радиуса г находится в неустойчивом равновесии с расплавом, rjp = 2avlzer; к- удельная
электропроводность расплава.
Система уравнений (4) решалась численно методом конечных разностей в MathCad Professional с постоянным, малым шагом по времени А/, обеспечивавшим необходимую точность расчета. Разработанная программа расчета позволяла получить массив значений искомых величин в каждый момент времени и построить точные временные зависимости перенапряжения, тока роста и радиуса кластера. Появление критического зародыша в системе осуществлялось при достижении определенного перенапряжения.
Результаты расчетов, проведенных с использованием типичных параметров электроосаждения серебра из нитратного расплава, показали качественное соответствие расчетных зависимостей ф) (рис. 1) экспериментальным данным, позволили описать ход процесса электроосаждения и выявить ряд особенностей роста кластера.
Установлено, что наиболее существенный рост кластера происходит в начальный период - от момента введения критического кластера (стрелка на рис. 1) до максимального перенапряжения, причем увеличение размеров кластера примерно на 2 порядка происходит за время ~104с. Причина этого - очень высокая плотность тока, идущего на рост наноразмерного зародыша.
т}, мВ 40
Рис. 1. Зависимость перенапряжения от времени (/ = 5-10"4 А-см"2, z'o = 103 А-см"2, Т= 543 К, и= 2-Ю"23 см3,5 = 7.85-10"3 см2,
30
20
10
0
2 3 НО2, с
3
Сд = 8-Ю"5 Ф-см"2, Г0 = 1.2-1013 см"2, сг= 10"5 Дж-см'2, к= 0.88 Ом"'-см'1). Стрелкой указан момент введения критического кластера. R0JRa при ц,
шах
равно 0.9, в точках 1 - 3.1; 2 - 4; 3 - 4.6.
С ростом кластера происходит изменение соотношения сопротивлений стадий разряда Ла =1/(2яг2/'о) и доставки Лом =1/2жгк (рис. 1), что приводит к смене лимитирующей стадии роста. Скорость роста крупных зародышей определяется миграцией ионов в электролите (см. точки на рис. 1).
Установлено, что введение критического кластера при различных перенапряжениях (20, 30, 40 мВ) не меняет общий вид зависимостей ф), г{{) и характер процессов, протекающих на электроде. Размер растущего кластера после спада перенапряжения не зависит от момента появления критического зародыша, так как все осаждаемое вещество после пика перенапряжения, связанного с фазообразованием, сосредоточено в новой фазе.
Количественное подтверждение адекватности модели было получено при сравнении расчетов с результатами эксперимента по росту единичного кристалла серебра (гальваностатический метод; расплав - AgNOз; Т= 500 -г- 525 К; / = 6.38-10'5 -г 2.55-10"4 А-см"2; атмосфера - воздух; контейнер -кварцевая пробирка; рабочий электрод - отполированный алмазной пастой торцевой микрокатод из иридиевой проволоки диаметром 0.5 мм, заплавленной в пирексовое стекло; анод и электрод сравнения - серебряные проволоки диаметром 1 мм; предварительная анодная обработка при потенциале от +1 до +20 мВ в течение 1-г10 мин; кривые ф) регистрировали с помощью запоминающего осциллографа).
По начальным участкам полученных экспериментальных зависимостей, используя уравнение баланса количества электричества на электроде до появления кластера, определяли равновесную концентрацию адатомов и удельную емкость двойного электрического слоя. Затем численно решали систему (4) с найденными значениями Го ((5 + 33)-1012 см"2) и Сд (15 -г 54 мкФ-см"2) при различной плотности тока обмена. Сравнивая расчетные и опытные кривые ф), находили значение ц, при котором совпадение было максимальным. Учитывая случайный характер процесса зародышеобразования, расчет проводился для каждой снятой кривой отдельно, значения Сд и Го не усреднялись.
Типичные теоретические и экспериментальные кривые i](t) и /-(/) приведены на рис. 2.
г), мВ 30
б H О2 с
6 /Юте
Рис. 2. Зависимости перенапряжения (а, 6) и радиуса (в) кластера от времени. Начальная точка расчета: а) ПтспУ, б) (? = 0, Г\ = 0). Кластер вводили при Г = 1 - экспериментальная кривая п(()\ 2-5 - расчет 77(0 при г0, А'См'2: 2 - со, 3 - 1000,4 - 500, 5 - 100; 6-7 - расчет г(і), соответствующий: 6 - кривой 3(а), 7 - кривой 3(6); точки рассчитаны по кривой 1. Г= 510 К, Сд = 5-Ю"5 Ф-см-2, Г0 = 2.38-1013 см-2, і = 2.55-10"4 А-см-2,1тах = 2.66-10"2 с, цтса = 30.8 мВ, 1 Ом"'-см"', = 1.96-10'3 см2, и = 2-Ю"23 см3, <х= 10"5 Дж-см"2.
Отклонение расчетной зависимости ф) от опытной при возрастающем перенапряжении (рис. 26), вероятно связано с осаждением примесных ионов. Во всех случаях максимальное совпадение кривых наблюдалось при ц > 1000 А'см'2, для найденной низшей границы плотности тока обмена частота присоединения ионов к критическому кластеру В = Іжг^/ге сос тавляет ~ 108 с"1.
Аналогичное исследование было проведено в расплаве, содержащем избыток фонового электролита, где доставка вещества к поверхности электрода и кластерам, осуществляется диффузией осаждаемых ионов в объеме расплава. Так как в данной системе влияние /0 проявляется лишь в самые начальные моменты роста (при г ~ 10'7 см), плотность тока роста определяется выражением і = гес0Р[ 1 - ехр/(?7 -т])]/г. (5)
Число образовавшихся кристаллов N в данном случае всегда было больше одного, поэтому в компьютерном эксперименте одновременно (при t=t:/iax) вводили соответствующее число сверхкритических кластеров, суммируя их токи роста при расчете drj/dt. При вычислениях использовали найденные из экспериментальных кривых значения Г0 ((1 * 4)-1013 см"2) и Сд (24 85 мкФ-см"2); коэффициент диффузии ионов Ag+ варьировали в пределах (0.7-^2.2)-10"5 см2-с"'. Типичные зависимости /?(0 и КО приведены на рис. 3.
Рис. 3. Зависимости перенапряжения (а) и размера кластера (б) от времени. 1-3 -расчет 7?(0 и КО с точки (и*> Птах) при Д см2-с"': 1 - 1-Ю'5, 2 - 1.5-10"5,3 - 2-Ю"5;
4 - экспериментальная кривая ф); точки - расчет по кривой 4. Т= 520 К, Сд = 3.43-10"5 Ф-см"2, Г0 = 3.35-1013 см"2, г = 6.38-10"5 А-см"2, /0 = °° (диффузионный режим), 1тах = 3.92-10"2 с, цпах = 12.3 мВ, с0 = 6.68-Ю20 см"3, ^ = 1.96-10° см2, и= 2-Ю"23 см3, о-= 10'5 Дж-см"2, 2.
Более резкий спад теоретических кривых г}{{) по сравнению с экспериментальными (рис. За) говорит о замедленности роста реальных кристалликов вследствие, по-видимому, пассивации поверхности растущих зародышей продуктами термического разложения расплава и примесями.
В главе 3 рассмотрена кинетика множественного зародышеобразования в расплавах. Исследованы зарождение и начальные стадии роста кластеров, имеющих форму полусферы, и нитевидных кристаллов.
Моделирование процесса образования кластеров новой фазы проводили на основе уравнения Фольмера, то есть считали, что число кластеров Аг, существующих на электроде к данному моменту времени /, равно
мо^М' 1 к
о
еХР|- ,Л2
м, (6)
К] и К2 - константы зародышеобразования, найденные из эксперимента по изучению статистики образования первого зародыша при потенциостатическом злектроосаждении. Для моделирования процесса роста группы невзаимодействующих полусферических кластеров в индивидуальном расплаве использовали систему (4), учитывая, что при расчете перенапряжения необходимо суммировать токи роста всех сверхкритических кластеров. В расплаве с избытком фонового электролита плотность тока роста рассчитывали по уравнению (5).
Процедура вычислений была аналогична описанной в главе 2, за исключением того, что критерием образования кластеров было достижение N целочисленного значения. При этом в систему вводилось соответствующее число кластеров, размер которых несколько превышал критический для текущего значения з], начинался расчет их радиуса, тока роста, фазового перенапряжения. Начальные условия расчета ДО) = Г0, 77(0) = О, Лг(0) = 0, г(0) = 0; вычисления проводили с использованием типичных параметров элекгрокристаллизации серебра из нитратных расплавов.
Расчеты показали, что при плотностях постоянного внешнего тока /> 5-10"4 А-см"2 на катоде происходит процесс множественного зарождения и роста, что соответствует экспериментальным результатам и подтверждает возможность использования уравнения (6) для определения числа образующихся зародышей в гальваностатических условиях.
Установлено, что количество кристаллов, образовавшихся в расплаве индивидуальной соли значительно меньше, чем в присутствии избытка фонового электролита (рис. 4). Это связано со снижением токов роста при уменьшении концентрации осаждаемых ионов в электролите, более медленным падением перенапряжения, удлинением периода зарождения. Повышение плотности внешнего тока ускоряет процесс фазообразования, приводит к росту Лтя и N как в индивидуальном расплаве, так и в электролите с избытком фона.
При медленном снижении г/ в компьютерном эксперименте наблюдалось растворение небольшого числа (единиц) кластеров, образовавшихся в конце процесса зарождения, при убывающем перенапряжении, за счет этого ускорялся рост раньше появившихся крупных зародышей.
Рис. 4. Зависимости перенапряжения от времени при /0 = 103 Л-см"2: а) индивидуальный расплав; б) расплав с избытком фона. А-см"2:1- Ю"3; Н-5-10"4; III- 10"4.
Числа у кривых - количество образовавшихся кластеров, в скобках - N до г]тах\ стрелка - момент возникновения первого кластера.
В ряде случаев электроосаждения на катоде образуются и растут кристаллы в виде тонких длинных нитей. В работе приведен литературный обзор, касающийся механизмов образования нитевидных кристаллов при электрокристаллизации и вариантов моделирования их роста.
Предложена следующая модель нитевидного кристалла: боковая поверхность имеет форму цилиндра, а вершина - полусферы. Удлинение кристалла (рост в направлении, перпендикулярном поверхности электрода) происходит в результате прямого встраивания разряжающихся ионов в места роста на вершине кристалла. В ходе электроосаждения боковая поверхность кристалла, а также радиус г полусферической и цилиндрической его частей не изменяются. На межфазной границе расплав/электрод реализуются условия незаторможенного перехода, происходит накопление адатомов. Баланс тока на электроде при изменении перенапряжения з системе достигается за счет протекания процессов заряжения-разряжения двойного электрического слоя,
накопления одиночных адатомов на электроде или их реионизации в объем расплава, и расходования тока на рост всех существующих к данному моменту времени нитевидных кристаллов.
В индивидуальном расплаве подвод вещества к электроду и растущему кристаллу осуществляется миграцией осаждаемых ионов в объеме электролита к поверхности электрода и вершине нити. Скорости образования и роста определяются величиной тока обмена на границе расплав/вершина кристалла и омическим сопротивлением слоя расплава вблизи нее, процесс зарождения и роста нитевидных кристаллов описывается системой уравнений
Ж Сд+ге/Г0ещ>/т) { К V) ^
. _т]-{2аи/гег) ¿1 _ 2 о .
1н'ШМ7/к)' 7ГТе1"
где /„ - плотность тока, идущего на рост нити; / - длина цилиндрической части растущего кристалла. При расчетах методом конечных разностей использовали следующие' начальные условия ДО) = Д>, .'/(0) = 0, МО) = 0, /(0) = 0 и параметры электрокристаллизации серебра из нитратного расплава: Т= 523 К, С# = 8-10'5 Ф-см"2, Го= 1.2-1013 см"2, и= 2-10'23 см3, к= 1 Ом"1-см"1, а~-1.5-10"4 Дж-см"2, 107см"2-с"', К2 = 10"2 В2, в = 7.85-10'3 см2. Введение соответствующего числа зародышей (полусфер с радиусом г=10"3см) осуществлялось при достижении N целочисленного значения, с этого же момента начинался расчет их длины и г'„. Расчетные временные зависимости перенапряжения и длины первого кристалла приведены на рис. 5.
Компьютерный эксперимент показал, что после включения постоянного внешнего тока происходит заряжение двойного электрического слоя, увеличение концентрации мономеров, рост перенапряжения. Образование кристаллов начинается при концентрации адатомов, составляющей 1.8-4.7 % от концентрации атомов в монослое гоютноупахованной грани серебра. Если ток роста не превышает внешний, то наблюдается дальнейшее увеличение г], Г, N. При т]„ш суммарный ток роста кристаллов равен или превышает внешний ток за счет
разряжения двойного электрического слоя и реионизации адатомов в объем расплава. На конечных стадиях процесса весь ток идет на рост кристаллов, концентрация адатомов близка к стационарному значению.
б ------------
.11 \ \4 / \ V 5-- 187(83)
1 \6---
Щ22\ 6(3)
5 10 15 !■ 102, с
- 4
5
5 10 15 МО2,с
, 5. Временные зависимости перенапряжения и длины первого кристалла: а) ¡о = °о; г, А-см"2: 1 -10'1; 2 - 10"2; 3 - 10"3; б) г = 10"3 А-см"2, /о, А-см"2:4 - 100; 5 - 500; 6 - оо.
Увеличение плотности тока приводит к росту г}тах и N (рис. 5а); при г = 10"4 А-см"2 образуется только один нитевидный кристалл. С ростом плотности тока обмена на границе расплав/вершина нити происходит увеличение тока, расходуемого на ее рост, что вызывает быстрое падение потенциала и уменьшение скорости появления кристаллов (рис. 56).
При гальваностатическом электроосаждении из расплава, содержащего избыток фонового электролита, доставка вещества к электроду и нитевидным кристаллам происходит за счет диффузии в объеме расплава. Скорость зарождения и роста нитей зависит от величины тока обмена на границе расплав/кластер новой фазы, концентрации и коэффициента диффузии
осаждаемых ионов в расплаве. При используемом в модели значении радиуса вершины 10"5 см, рост кристалла лимитирует диффузия, уравнение для плотности тока роста нитей в системе (7) заменяется выражением (5). Расчетные зависимости и /(?) при £) = 2-Ю'5 см2-с"' приведены на рис. 6.
0 10 20 30 МО2,с 0 1 0 2 0 30 /.102,с
Рис. 6. Временные зависимости перенапряжения и длины первого кристалла: а) с0 = 6-Ю19 см"3, А-см'2: 1 - Ю-3; 2 - 5-Ю"4; 3 - 10"4; б) I = 10'4 А-см"2, с0, см'3: 4 - 6-Ю20; 5 - 6-Ю15; 6 - 6-Ю18.
Количество нитей, образовавшихся при электроосаждении из расплава индивидуальной соли намного меньше, а их длина больше, чем в случае расплава с избытком фона, что объясняется большими токами роста в индивидуальном расплаве. Повышение плотности внешнего тока вызывает увеличение т]тах и N (рис. 6а). Рост концентрации осаждаемых ионов в объеме электролита приводит к увеличению тока роста нити, более резкому падению перенапряжения и уменьшению скорости зарождения (рис. 66). Полученные результаты качественно согласуются с экспериментом по росту волокнистых осадков серебра в нитратных расплавах.
Выводы
1. Найдено новое решение основного кинетического уравнения, описывающего процесс зародышеобразования (уравнения Фоккера-Планка)_для электрохимического потенциостатического зарождения новой фазы в расплавленных солях. Получены выражения для скорости зарождения, индукционного периода зарождения, коэффициента диффузии в пространстве размеров зародышей.
2. Установлено, что образование нанокластеров серебра при электроосаждении из расплава AgNOз лимитирует разряд серебра на поверхности кластера. В нитратном расплаве AgNOз-KNOз-NaNOз, содержащем избыток фонового электролита, процесс зародышеобразования контролируется разрядом и диффузией осаждаемых ионов в объеме расплава к поверхности нанокластера.
3. Впервые получено количественное соответствие результатов компьютерного моделирования и эксперимента для начальных стадий роста (от нано- до микроразмеров) кластеров серебра при электроосаждении на индифферентном иридиевом электроде из расплавов нитратов серебра.
4. Показано, что наиболее быстрое увеличение размера кластера происходит в начальные моменты его роста вследствие высокой плотности тока роста нанокластера, влияние кривизны поверхности кластера на кинетику разряда сказывается только в области размеров, близких к критическому (в области наноразмеров). На рост нанокластеров серебра в исследованных расплавах оказывает влияние разряд ионов серебра на поверхности кластера, в то время как скорость роста микрокристаллов определяется стадией доставки осаждаемых ионов в электролите.
5. Оценены основные параметры электрокристаллизации серебра в нитратных расплавах: удельная емкость двойного электрического слоя, концентрация адагомов на электроде, плотность тока обмена на границе расплав/кристалл серебра, частота присоединения частиц к критическому кластеру новой фазы, коэффициент диффузии ионов серебра в объеме расплава.
Найдено, что начальная (при пулевом перенапряжении) концентрация одиночных адатомов (мономеров) на электроде составляет 0.3 - 2.3 %, образование первого сверхкритического кластера происходит при концентрации адатомов 0.4 - 4 % относительно концентрации атомов в монослое плотноупакованной грани серебра.
6. Разработана математическая модель образования и роста нитевидных кристаллов новой фазы при электроосаждении из расплавленных солей, учитывающая баланс тока на электроде и изменение пересыщения в системе.
Получены временные зависимости перенапряжения, числа и длины нитевидных кристаллов. Показано, что повышение плотности постоянного внешнего тока удлиняет период зародышеобразования. Рост концентрации осаждаемых ионов в объеме электролита и увеличение плотности тока обмена на границе кристалл/расплав приводят к уменьшению числа образующихся кристаллов.
Основное содержание диссертации опубликовано в работах:
1. Исаев В.А., Гришенкова О.В. Токи обмена в процессах электрохимического фазообразования // Расплавы. 1999. № 5. С. 84-89.
2. Isaev V.A., Grishenkova O.V. Kinetics of electrochemical nucleation and growth // Electrochemistry Communication. 2001. V. 3. № 9. P. 500-504.
3. Исаев B.A., Гришенкова О.В. Кинетика разряда ионов в электрохимическом фазообразовании//Расплавы. 2002. №3. С. 58-63.
4. Isaev V.A., Grishenkova O.V. Exchange current densities at electrochemical phase formation from molten salts // J. Mining and Metallurgy. 2003.39 (1-2) B. P. 167-175.
5. Исаев B.A., Гришенкова О.В. Кинетика электроосаждения металлов // Вестник УГТУ-УПИ. Серия химическая. Теория и практика электрохимических процессов. Екатеринбург: ГОУ ВПО УГТУ-УПИ, 2004. № 14(44). С. 93-97.
6. Исаев В.А., Гришенкова О.В. Образование и рост нитевидных кристаллов в расплавах индивидуальных солей // Расплавы. 2007. № 1. С. 40-45.
7. Исаев В.А., Гришенкова О.В. Кинетика формирования нитевидных кристаллов // Электрохимия. 2008. т. 44. № 6. С. 712-715.
8. Исаев В.А., Гришенкова О.В. Скорость электрохимического зародышеобразования в расплавах // Расплавы. 2008. № 1. С. 80-84.
9. Исаев В.А., Гришенкова О.В. Нестационарное электрохимическое зародышеобразование в расплавах // Расплавы. 2009. N° 3. С. 45-49.
10. Гришенкова О.В., Исаев В.А. Образование нано- и микрокристаллов при электроосаждении. I. Индивидуальный расплав // Расплавы. 2009. № 4. С. 79-84.
11. Гришенкова О.В., Исаев В.А. Образование нано- и микрокристаллов при электроосаждении. II. Фоновый расплав // Расплавы. 2009. № 5. С. 65-70.
12. Исаев В.А., Гришенкова О.В. Скорость электрохимического образования нанокластеров // Расплавы. 2010. № 2. С. 77-82.
13. Гришенкова О.В., Семерикова О.Л., Исаев В.А. Рост единичного кристалла серебра при электроосаждении из расплава. Нано- и микроразмеры // Расплавы. 2010. №5. С. 56-61.
14. Гришенкова О.В. Электрохимический рост зародышей новой фазы // Проблемы электрокристаллизации металлов: Тез. докл. конф., посвящ. 80-летию со дня рождения акад. А.Н. Барабошкина. Екатеринбург. 14-15 нояб. 2000. С. 16.
15. Исаев В.А., Гришенкова О.В., Вакарин C.B. Кинетика разряда ионов в процессах электроосаждения металлов // Тез. докл. X Кольского семинара по электрохимии редких металлов. Апатиты. 4-7 дек. 2000. С. 33.
16. Гришенкова О.В., Исаев В.А. Моделирование роста новой фазы при электролизе расплавов // Тез. докл. ХП Российск. конф. по физ. химии расплавленных и твердых электролитов. Нальчик. 17-23 сент. 2001. Т. 1. С. 236-237.
17. Исаев В.А., Гришенкова О.В. Кинетика начальных стадий электроосаждения металлов // Там же. С. 239-241.
18. Isaev V.A., Grishenkova O.V. Galvanostatic phase formation at the electrode-molten salt interface // Abstr. EUCHEM Molten Salts Conferences. Oxford. 1-6 Sept. 2002. P33.
19. Гришенкова O.B., Исаев В.А. Кинетика разряда ионов при электрохимическом образовании новой фазы // Новые неорганические материалы и химическая термодинамика: Тез. докл. второго семинара СО РАН-УрО РАН. Екатеринбург. 24-26 сент. 2002. С. 56.
20. Исаев В.А., Гришенкова О.В. Кинетика электроосаждения металлов // Тез. докл. научно-практ. конф., посвящ. 80-летию кафедры ТЭХП УГТУ-УПИ. Екатеринбург. 22-25 окт. 2003. С. 34-36.
21. Isaev V.A., Grishenkova O.V. Electrodeposition from molten salts // Abstr. EUCHEM Molten Salts Conferences. Piechowice. Poland. 20-25 June 2004. P. 138.
22. Исаев В.А., Гришенкова О.В. Кинетика разряда ионов в расплавах // Тез. докл. XIII Российск. конф. по физ. химии и электрохимии расплавленных и твердых электролитов. Екатеринбург. 27 сент.-1 окт. 2004. Т. 1. С. 126-127.
23. Исаев В.А., Мурашова И.Б., Зайков Ю.П., Гришенкова О.В., Вакарин C.B. Рост волокнистых осадков металлов // Там же. С. 160.
24. Исаев В.А., Гришенкова О.В. Моделирование роста наночастиц // Химия твердого тела и функциональные материалы-2004: Тез. докл. Всероссийск. конф. Екатеринбург. 25-28 окт. 2004. С. 167.
25. Isaev V.A., Grishenkova О.У. Computer modeling of galvanostatic nucleation // Electrochemistry: From Nanostructures to Power Plants. I. Abstr. 55th Annual Meeting of ISE. 19-24 Sept. 2004. Thessaloniki. Greece. P. 566.
26. Исаев B.A., Гришенкова O.B., Вакарин C.B. Моделирование роста волокнистых осадков серебра // Современные аспекты электрокристаллизации металлов: Тез. докл. конф., посвящ. 80-летию со дня рождения акад. А.Н. Барабошкина. Екатеринбург. 16-17 нояб. 2005. С. 29.
27. Исаев В .А., Поляков Е.Г., Гришенкова О.В., Вакарин C.B. Подходы к моделированию процесса кристаллизации при металлотермическом получении танталовых порошков // Там же. С. 30-31.
28. Мурашова И.Б., Зайков Ю.П., Петров Л.А., Хохлов С.Л., Малков В.Б., Зотин И.В., Исаев В.А., Гришенкова О.В. Динамика роста и структура дендритного осадка серебра в условиях конвективной подачи вещества // Там же. С. 116-117.
29. Исаев В.А., Гришенкова О.В. Скорость электрохимического зародышеобразования в расплавах // Физическая химия и электрохимия расплавленных и твердых электролитов: Тез. докл. XIV Российск. конф. 2007. Т. 1.С. 78-79.
30. Гришенкова О.В., Семерикова О.Л., Исаев В.А. Образование нанокристаллов серебра при электроосаждении из расплавов // Физическая химия и электрохимия твердых электролитов. Прикладные аспекты высокотемпературной электрохимии: Тез. докл. XV Российск. конф. по физической химии и электрохимии расплавленных и твердых электролитов (с междунар. участием). Нальчик: КБГУ. 13-19 септ. 2010. С. 213-215.
31. Гришенкова О.В., Семерикова О.Л., Исаев В.А. Кинетика начальных стадий роста кластеров серебра при электролизе расплавов // Кинетика и механизм кристаллизации. Самоорганизация при фазообразовании: Тез. докл. VI Междунар. науч. конф. Иваново. Россия. 21-24 септ. 2010. С. 69-70.
Подписано в печать 17.02.2012. Формат 60x84 1/16. Усл. печ. л. 1,5. Тираж 110 экз. Заказ ¿$0?
Отпечатано с готового оригинал-макета
ООО Издательство УМЦ УПИ 620002, Екатеринбург, ул. Гагарина, 35-а
61 12-2/469
Федеральное государственное бюджетное учреждение науки Институт высокотемпературной электрохимии Уральского отделения Российской академии наук
На правах рукописи
Гришенкова Ольга Владимировна
КИНЕТИКА ОБРАЗОВАНИЯ И НАЧАЛЬНЫХ СТАДИЙ РОСТА НОВОЙ ФАЗЫ ПРИ ЭЛЕКТРООСАЖДЕНИИ ИЗ РАСПЛАВЛЕННЫХ СОЛЕЙ
Специальность 02.00.05 - Электрохимия
диссертация на соискание ученой степени кандидата химических наук
Научный руководитель доктор химических наук Исаев Владимир Александрович
Екатеринбург - 2012
Оглавление
Список условных обозначений 4
Введение 6 Глава 1. Энергетика и кинетика электрохимического
зарод ышеобразования......................................................................................................11
1.1 Работа электрохимического зародышеобразования..................................................11
1.2 Кинетический анализ процесса зарождения новой фазы........................................20
1.2.1 Дискретный анализ (подход Беккера-Деринга)................................................20
1.2.2 Континуальная модель. Основное кинетическое уравнение..................24
1.3 Стационарная скорость электрохимического зарождения..............................................27
1.4 Нестационарное электрохимическое зародышеобразование в расплавах (потенциостатические условия)........................................................................35
1.5 Выводы по главе..........................................................................................................................................47
Глава 2. Рост единичных кристаллов серебра при электроосаждении
из расплавов в гальваностатических условиях......................................49
2.1 Основные закономерности гальваностатического фазообразования. Литературный обзор результатов эксперимента и моделирования.... .. 49
2.2 Кинетика разряда ионов на зародыше новой фазы......................................................55
2.3 Компьютерное моделирование роста единичных кристаллов. Сравнительный анализ расчетных и экспериментальных данных..... ... 63
2.3.1 Индивидуальный расплав......................................................................................................63
2.3.2 Расплав, содержащий избыток фонового электролита................................80
2.4 Выводы по главе............................. ................................................................................91
Глава 3. Кинетика множественного зародышеобразования в
расплавах..........................................................................................................................................93
3.1 Моделирование образования и роста нано- и микрокристаллов....................93
3.2. Кинетика формирования нитевидных кристаллов при
электроосаждении......................................................................................................................................105
3.2.1 Индивидуальный расплав......................................................................................................105
3.2.2 Образование нитевидных кристаллов при электроосаждении в присутствии избытка фонового расплава..................................................................113
3.3 Выводы по главе..........................................................................................................................................118
Основные результаты и выводы......................................................................................................120
Список использованных источников..........................................................................................122
Основные обозначения
А, В - кинетический коэффициент и коэффициент диффузии в пространстве размеров в уравнении Фоккера-Планка, с"1
а, Ъ - частоты присоединения частиц к кластеру и отрыва от кластера, с"1
л
Со ~ концентрация осаждаемых ионов в глубине электролита, см"
о
cs - концентрация осаждаемых ионов вблизи поверхности электрода, см"
csr - концентрация осаждаемых ионов вблизи растущего зародыша, см"3
Сд - удельная емкость двойного электрического слоя, Ф-см"2
D - коэффициент диффузии, см2-с-1
Е - потенциал электрода, В
е ~ заряд электрона, е = 1,602-10"19 Кл
/ - коэффициент равный zefkT
G - свободная энергия Гиббса, Дж
g - число атомов в зародыше
/ - ток, А
i ~ плотность тока, А-см"2
¿о - плотность тока обмена, А-см'
if ~~ плотность фарадеевского тока, А-см"
J - скорость (частота) зародышеобразования, см'^с"1
Jem ~~ стационарная скорость зародышеобразования, см'2-с"1
-2 -1
J - плотность потока, см -с
KhK2 - константы зародышеобразования
7 лт 1
к - постоянная Больцмана, к = 1,3 81 • 10" Дж-К"
к+, к~ - константы скорости катодного и анодного процессов
N - концентрация образовавшихся зародышей, см"
q - количество электричества, Кл
г - радиус зародыша, расстояние, см
s - площадь электрода, см"2
sg - площадь поверхности раздела ^-атомный кластер/электролит, см"2
Т ~ температура, К
г - время, с
и - внутренняя энергия, энергия активации, Дж
Ж - работа образования кластера, Дж
2 - функция распределения кластеров по размерам, см"2
2е - равновесная функция распределения кластеров по размерам, см"
2п - фактор неравновесности (Зельдовича)
г - валентность иона
а, ¡5 ~ коэффициенты переноса заряда (а + /3= 1)
Г - концентрация адатомов (мономеров), см"
Га - концентрация адатомов, находящихся на активных центрах электродной поверхности при нулевом перенапряжении, см"2
т-1 -2
Го - равновесная концентрация адатомов, см
Ун ~ краевое натяжение, Дж-см"1
д ~ толщина диффузионного слоя электрода, см
г( ~ перенапряжение, В
к ~ удельная электропроводность электролита, Ом"1-см"1
¡л ~ химический потенциал, Дж
/л - электрохимический потенциал, Дж
<7 ~ межфазное натяжение, Дж-см"
г ~ индукционное время зародышеобразования, с
и ~ объем новой фазы, приходящийся на одну частицу, см"3
(р - электрический потенциал, В
Введение
Актуальность работы
Изучение начальных стадий фазообразования в электрохимических системах является актуальной задачей фундаментальной науки, так как позволяет определить механизмы и лимитирующие стадии зарождения и роста новой фазы, получить прямую информацию о кинетике перехода частиц через межфазную границу электролит / кластер.
Критические кластеры новой фазы при электрокристаллизации имеют наноразмеры, поэтому исследование электрохимического
зародышеобразования дает возможность установить закономерности образования и роста нанокластеров, оценить размерные эффекты. В настоящее время в связи с широким внедрением нанотехнологий актуальность данной тематики особенно возросла.
Результаты исследования начальных стадий электрокристаллизации в расплавленных солях могут быть использованы и при решении прикладных задач, связанных с определением условий формирования электродных осадков и покрытий с заданной структурой.
Работа выполнена в соответствии с планами НИР Института высокотемпературной электрохимии УрО РАН "Разработка научных основ безотходных экологически чистых, энерго- и ресурсосберегающих высокотемпературных электрохимических технологий получения металлов и неорганических соединений в солевых расплавах" (№ гос. регистрации 01.2.007 01883), "Физико-химические основы электрохимических технологий рафинирования и получения металлов и их соединений в расплавах солей" (№ гос. регистрации 01.2.010 00812), а также в рамках проектов: РФФИ-"УРАЛ" № 04-03-96124 р2004урал_а "Электрохимическое образование волокнистых осадков металлов с заданной каталитической активностью из расплавов и водных растворов"; совместного проекта УрО РАН-СО РАН № 81 "Исследование процессов перехода веществ и систем между различными
агрегатными состояниями" (2006-2008 г.г.); программ Президиума РАН "Разработка методов получения химических веществ и создание новых материалов: создание новых конструкционных и функциональных материалов на базе нанотехнологий" по теме "Исследование условий получения нанопорошков и нановолокон из серебра (модельный металл) электролизом растворов и низкотемпературных серебросодержащих солевых расплавов" (2006-2008 г.г.) и "Основы фундаментальных исследований нанотехнологий и наноматериалов" по проекту 09-П-3-1026 "Получение нанокристаллических осадков электролизом солевых расплавов" (20092011 г.г.).
Цель работы
Определение основных закономерностей образования и роста кластеров новой фазы в области их нано- и микроразмеров при электроосаждении из расплавленных солей в потенцио- и гальваностатических условиях.
Задачи
1. Теоретически изучить закономерности зарождения кластеров новой фазы при постоянном пересыщении (в потенциостатических условиях).
2. Теоретически исследовать кинетику роста кластеров на начальных этапах гальваностатического электроосаждения из расплавленной соли в модельной системе (осаждаемый металл - расплав - А§Ж)з, А§Ж)з-КШз-ЫаЖЬ; катод - 1г, К).
3. Разработать компьютерную модель этого процесса, сравнить результаты расчета с экспериментальными данными, оценить возможность применения предложенной модели для анализа начальных стадий фазообразования в расплавленных солях.
4. Использовать компьютерный эксперимент для определения механизма, лимитирующих стадий и основных параметров электрокристаллизации в исследуемых системах.
5. Разработать математическую модель образования и роста нитевидных кристаллов при электроосаждении из расплавленных солей, провести компьютерное моделирование процесса.
Научная новизна
1. Получено новое решение основного кинетического уравнения, описывающего процесс зародышеобразования (уравнения Фоккера-Планка), для случая электрохимического потенциостатического фазообразования в расплавленных солях. Найдены выражения для скорости зародышеобразования, учитывающие лимитирующие стадии образования критических кластеров новой фазы при электроосаждении из расплава индивидуальной соли и расплава, содержащего избыток фонового электролита.
2. Разработана компьютерная модель единичного и множественного зародышеобразования в гальваностатических условиях электроосаждения из расплавленных солей.
3. Впервые получено количественное подтверждение адекватности предложенной модели гальваностатического фазообразования в расплавленных солях. Оценены основные параметры электрокристаллизации в исследуемых системах: удельная емкость двойного электрического слоя, концентрация адатомов (мономеров) на электроде, плотность тока обмена на границе кристалл / расплав, коэффициент диффузии осаждаемых ионов в объеме расплава.
4. Разработана модель роста нитевидных кристаллов новой фазы, формирующихся на катоде при гальваностатическом электроосаждении из расплавленных солей.
Научная ценность и практическая значимость
Результаты работы и сделанные в ней обобщения представляют теоретический интерес для развития представлений о закономерностях начальных стадий электрохимического фазообразования. Аналитические
методы решения и предложенный метод компьютерного моделирования могут быть использованы для определения механизма и лимитирующих стадий процесса, а также основных параметров электрокристаллизации металлов в расплавленных солях и, следовательно, будут полезны при выборе оптимальных условий получения осадков с заданными свойствами.
На защиту выносятся:
1. Решение основного кинетического уравнения (уравнения Фоккера-Планка), описывающего процесс зародышеобразования в потенциостатических условиях электроосаждения из расплавленных солей. Уравнения для расчета скорости зародышеобразования, учитывающие лимитирующие стадии образования критических кластеров новой фазы при электрокристаллизации.
2. Результаты компьютерного моделирования образования и начальных стадий роста кристаллов серебра при гальваностатическом электролизе солевых расплавов.
3. Оценка основных параметров электрокристаллизации в исследуемых системах: удельной емкости двойного электрического слоя, концентрации адатомов на электроде, плотности тока обмена на границе кристалл / расплав, коэффициента диффузии осаждаемых ионов в объеме расплава.
4. Математическая модель образования и роста на электроде нитевидных кристаллов, результаты компьютерной реализации данной модели.
Личный вклад автора
Теоретическое исследование фазообразования при постоянном и переменном пересыщении в системе, вывод уравнений, компьютерное моделирование, в том числе разработка оригинальных программ расчета, обработка экспериментальных данных, анализ, выводы и оформление результатов выполнены лично автором. Постановка задач осуществлялась
научным руководителем д.х.н. Исаевым В. А. Экспериментальное исследование электрохимического зарождения и роста кристаллов серебра было выполнено к.х.н. Семериковой O.JI.
Апробация работы
Основные результаты и положения исследований были представлены на российских и международных конференциях и семинарах: "Проблемы электрокристаллизации металлов" и "Современные аспекты электрокристаллизации металлов" (Екатеринбург, 2000 и 2005); X Кольском семинаре по электрохимии редких металлов (Апатиты, 2000); XII, XIII, XIV, XV Российской конференциях по физической химии расплавленных и твердых электролитов (Нальчик-Екатеринбург, 2001, 2004, 2007, 2010); EUCHEM Molten Salts Conferences (Oxford, UK, 2002 и Piechowice, Poland, 2004); втором семинаре CO РАН-УрО РАН "Новые неорганические материалы и химическая термодинамика" (Екатеринбург, 2002); конференции, посвященной 80-летию кафедры ТЭХП УГТУ-УПИ (Екатеринбург, 2003); Всероссийской конференции "Химия твердого тела и функциональные материалы-2004" (Екатеринбург, 2004); 55th Annual Meeting of ISE (Thessaloniki, Greece, 2004); VI Международной научной конференции "Кинетика и механизм кристаллизации. Самоорганизация при фазообразовании" (Иваново, 2010).
Публикации
По результатам исследований опубликовано 13 статей в российских и международных журналах из списка ВАК и 18 тезисов докладов в материалах российских и международных конференций.
Структура и объем диссертации
Диссертационная работа состоит из введения, трех глав, обобщающих выводов и списка использованных источников. Диссертация изложена на 138 страницах, содержит 34 рисунка, 2 таблицы; в списке использованных источников 174 работы отечественных и зарубежных авторов.
Глава 1. Энергетика и кинетика электрохимического зародышеобразования
1.1 Работа электрохимического зародышеобразования
Первоначально теория фазообразования была разработана для процесса появления капель жидкости в пересыщенном паре [1-3]. Дж. В. Гиббсом было показано, что работа обратимого изотермического образования зародыша, имеющего форму шара, складывается из работы, выигрываемой при возникновении объема новой фазы, и работы, затрачиваемой на формирование ее поверхности, и равна
4 ж г3
Ж = AG = -gAJU + scr =--Д// + 4жг о, (1.1)
3 и
где АС - изменение свободной энергии Гиббса; g - число частиц в капле; А ¡и - изменение химического потенциала при изменении агрегатного состояния при данных давлении и температуре, то есть пересыщение системы; я - поверхность образовавшегося кластера; сг - удельная поверхностная свободная энергия на границе пар/жидкость; г - радиус капли; и - объем, приходящийся на одну частицу в новой фазе.
При заданном пересыщении зависимость работы образования капли от ее радиуса имеет вид кривой с максимумом (рис. 1), соответствующим критическому значению
2 сто А/л
Тогда энергетический барьер, который необходимо преодолеть для образования капли критического размера с площадью поверхности равен
= 16жст3и2
кр ЗА//2 ' 1 ;
КР=^крст. (1.4)
(1-2)
Рисунок 1 - Зависимости работы образования зародыша новой фазы и ее объемной и поверхностной составляющих от размера согласно классической теории зарождения
Критический кластер новой фазы находится в состоянии
W W О 1 U
неустойчивого равновесия с пересыщенной материнской фазой; он с равной вероятностью может терять и приобретать атомы; работа его флуктуационного образования максимальна. В рамках классической теории нуклеации, оперирующей с непрерывными величинами размера кластера и работы его образования, эти три определения эквивалентны [4]. Кластеры с радиусом г < гкр имеют тенденцию к растворению, а для зародышей с г> гкр энергетически выгодно увеличивать свой размер.
С ростом пересыщения в системе кривая W(r) смещается к началу координат, критический размер кластеров и работа их образования уменьшаются.
Основные закономерности гетерогенного зарождения были установлены Дж. В. Гиббсом и М. Фольмером [1, 2]. Работа образования зародыша жидкой фазы, имеющего форму шарового сегмента (рис. 2), на твердой гладкой изотропной подложке определяется уравнением
W = -gAjU + ser + s0(a2(1-5)
где cr, <X|, (7*2 ~ межфазные натяжения границ раздела кластер/среда, подложка/среда и кластер/подложка; s = 2ят2(1 - cosé?), s0 = яг2 sin2 в -
площади поверхностей раздела кластер/среда и кластер/подложка.
<Т
среда
кластер
'фтшщшФшшшйшшшшт^
подяожш
Рисунок 2 - Поперечное сечение зародыша на подложке [5]
Используя соотношение Юнга crl - <г2 =<J cos ^> и учитывая, что объем шарового сегмента V = - 3cos в + cos3 £?), получим [6-8]
3
W = — (2-3 cos 9 + cos3 в)\-^-Аи + 4тгг2а 4 7 3 и
(1.6)
и для кластера критического размера
ЗА//
где Ф = ^-(2-Зсо8# + со83#) - фактор формы, зависящий от угла
смачивания и равный отношению объемов шарового сегмента жидкой фазы на поверхности подложки и шара того же радиуса. Сравнение уравнений (1.1) и (1.6) или (1.3) и (1.7) показывает, что при одинаковом пересыщении работа гетерогенного зарождения всегда меньше работы гомогенного образования зародышей (в пре�