Электрохимия, 2019, T. 55, № 2, стр. 250-252

Образование адатомов и кристаллических зародышей при электрокристаллизации металлов в гальваностатических условиях

Ю. Д. Гамбург *

Институт физической химии и электрохимии им. А.Н. Фрумкина РАН
119071 Москва, Ленинский просп., 31, Россия

* E-mail: gamb@list.ru

Поступила в редакцию 14.12.2017
После доработки 06.06.2018
Принята к публикации 25.07.2018

Полный текст (PDF)

Аннотация

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

Ключевые слова: металлы, электроосаждение, нуклеация, гальваностатический режим

1. ВВЕДЕНИЕ

Образование кристаллических зародышей при электроосаждении металлов в гальваностатических условиях i = const было предметом ряда публикаций как теоретических, так и экспериментальных [13]. Сначала было установлено, что при включении относительно небольшого катодного тока происходит смещение потенциала в отрицательном направлении. Смещение сначала возрастает, затем проходит через экстремум и далее снижается. Это, как полагали, связано с тем, что первоначально происходит накопление адатомов выделяемого металла на чужеродной поверхности, после чего при достаточно высоком достигнутом перенапряжении начинается образование кристаллических зародышей, которые поглощают часть адатомов. Момент максимума перенапряжения считали соответствующим началу нуклеации. В дальнейшем было достигнуто понимание того, что нуклеация начинается несколько раньше достижения максимума.

Большая часть более поздних работ в данной области выполнялась уже в потенциостатических условиях, существенно упрощающих анализ происходящих процессов. Тем не менее, интерес к процессам при i = const сохранился, так как основная часть процессов в технике проводится именно гальваностатически. Так, в [3] было установлено, что при невысокой плотности тока (сравнимой с плотностью тока обмена) вначале наблюдается почти линейное возрастание перенапряжения, связанное с линейным характером зависимости поверхностной концентрации адатомов Cs от времени. При этом предполагалось, что текущая величина электродного потенциала целиком определяется концентрацией адатомов на поверхности и связана с ним уравнением Нернста, которое в данном случае для одноэлектронного процесса принимает вид

(1)
$\eta = {{\left( {{{RT} \mathord{\left/ {\vphantom {{RT} F}} \right. \kern-0em} F}} \right){\text{ln}}{{C}_{{\text{s}}}}} \mathord{\left/ {\vphantom {{\left( {{{RT} \mathord{\left/ {\vphantom {{RT} F}} \right. \kern-0em} F}} \right){\text{ln}}{{C}_{{\text{s}}}}} {C_{{\text{s}}}^{0}}}} \right. \kern-0em} {C_{{\text{s}}}^{0}}},$
что при невысоком отношении ${{{{C}_{{\text{s}}}}} \mathord{\left/ {\vphantom {{{{C}_{{\text{s}}}}} {C_{{\text{s}}}^{0}}}} \right. \kern-0em} {C_{{\text{s}}}^{0}}}$ переходит в
(2)
$\eta = {{iRTt} \mathord{\left/ {\vphantom {{iRTt} {\left( {{{F}^{{\text{2}}}}C_{{\text{s}}}^{0}} \right)}}} \right. \kern-0em} {\left( {{{F}^{{\text{2}}}}C_{{\text{s}}}^{0}} \right)}},$
поскольку текущая концентрация адатомов Cs составляет $C_{{\text{s}}}^{0} + {{it} \mathord{\left/ {\vphantom {{it} F}} \right. \kern-0em} F}$ [моль/см2]. Далее подобный вид зависимости перенапряжения от времени нарушается в момент начала образования зародышей. Необходимо отметить еще несколько сравнительно недавних публикаций по гальваностатической нуклеации [46].

Ниже мы покажем, что допущение о линейной зависимости η(t), как и допущение о справедливости уравнения Нернста, в данном случае не соответствует действительности, противоречит эксперименту и приводит к существенным ошибкам при определении параметров процесса по экспериментальным данным.

2. КИНЕТИКА НАКОПЛЕНИЯ АДАТОМОВ НА ПОВЕРХНОСТИ И СООТВЕТСТВУЮЩИЕ ИЗМЕНЕНИЯ ПЕРЕНАПРЯЖЕНИЯ

Пусть концентрация ионов-реагентов в растворе составляет С0, плотность тока обмена i0, а равновесная концентрация атомов $C_{{\text{s}}}^{0},$ причем к моменту включения тока эта равновесная концентрация уже установилась. После включения небольшого катодного тока (например, i = i0/2), поверхностная концентрация адатомов, как уже сказано выше, увеличивается согласно ${{C}_{{\text{s}}}} = {{C_{{\text{s}}}^{0} + it} \mathord{\left/ {\vphantom {{C_{{\text{s}}}^{0} + it} F}} \right. \kern-0em} F}$ (процесс опять-таки для простоты считаем одноэлектронным).

Величина плотности тока i связана с величиной текущего перенапряжения η; поскольку ток мы приняли небольшим, то допустима линеаризация компонентов суммарной зависимости

(3)
$\begin{gathered} i\left( t \right) = {{i}_{{\text{c}}}}\left( t \right) + {{i}_{{\text{a}}}}\left( t \right) = {{i}_{0}}\left\{ {{\text{exp}}[\alpha f\eta \left( t \right)]--} \right. \\ \left. { - \,\,\left( {{{{{C}_{{\text{s}}}}} \mathord{\left/ {\vphantom {{{{C}_{{\text{s}}}}} {C_{{\text{s}}}^{0}}}} \right. \kern-0em} {C_{{\text{s}}}^{0}}}} \right){\text{exp}}[--({\text{1}}--\alpha )f\eta \left( t \right)]} \right\} \\ \end{gathered} $
(катодный ток мы приняли положительным, как и перенапряжение).

В таком случае

(4)
${{i}_{{\text{c}}}}\left( t \right) = {{i}_{0}}\left[ {1 + \alpha f\eta \left( t \right)} \right],$
(5)
${{i}_{{\text{a}}}}\left( t \right) = --{{i}_{0}}\left( {{{{\text{1}} + it} \mathord{\left/ {\vphantom {{{\text{1}} + it} {C_{{\text{s}}}^{0}F}}} \right. \kern-0em} {C_{{\text{s}}}^{0}F}}} \right)\left[ {{\text{1}}--\left( {{\text{1}}--\alpha } \right)f\eta \left( t \right)} \right],$
где учтено накопление адатомов на поверхности, увеличивающее анодную составляющую плотности тока.

Отсюда, опуская обозначение (t), получаем

(6)
$\begin{gathered} i = {{i}_{{\text{c}}}} + {{i}_{{\text{a}}}}\quad = {{i}_{0}}\left[ {{{f\eta --it} \mathord{\left/ {\vphantom {{f\eta --it} {C_{{\text{s}}}^{0}F}}} \right. \kern-0em} {C_{{\text{s}}}^{0}F}} + } \right. \\ \left. { + \,\,{{itf\eta \left( {{\text{1}}--\alpha } \right)} \mathord{\left/ {\vphantom {{itf\eta \left( {{\text{1}}--\alpha } \right)} {C_{{\text{s}}}^{0}F}}} \right. \kern-0em} {C_{{\text{s}}}^{0}F}}} \right] = {{i}_{0}}\left[ {{{f\eta --it} \mathord{\left/ {\vphantom {{f\eta --it} {C_{{\text{s}}}^{0}F{\text{ }}}}} \right. \kern-0em} {C_{{\text{s}}}^{0}F{\text{ }}}} + } \right. \\ \left. { + \,\,{{it\eta \left( {{\text{1}}--\alpha } \right)} \mathord{\left/ {\vphantom {{it\eta \left( {{\text{1}}--\alpha } \right)} {C_{{\text{s}}}^{0}RT}}} \right. \kern-0em} {C_{{\text{s}}}^{0}RT}}} \right], \\ \end{gathered} $
где f = F/RT, α – коэффициент переноса для катодного процесса.

Выражая отсюда η через остальные параметры, получаем

(7)
$\eta = {{\left[ {iRT\left( {{{C_{{\text{s}}}^{0}F} \mathord{\left/ {\vphantom {{C_{{\text{s}}}^{0}F} {{{i}_{0}} + t}}} \right. \kern-0em} {{{i}_{0}} + t}}} \right)} \right]} \mathord{\left/ {\vphantom {{\left[ {iRT\left( {{{C_{{\text{s}}}^{0}F} \mathord{\left/ {\vphantom {{C_{{\text{s}}}^{0}F} {{{i}_{0}} + t}}} \right. \kern-0em} {{{i}_{0}} + t}}} \right)} \right]} {\left[ {F\left( {C_{{\text{s}}}^{0}F + i\left( {{\text{1}}--\alpha } \right)t} \right)} \right]}}} \right. \kern-0em} {\left[ {F\left( {C_{{\text{s}}}^{0}F + i\left( {{\text{1}}--\alpha } \right)t} \right)} \right]}}.$

Необходимо повторить, что это справедливо только для небольших перенапряжений (в пределах линеаризации экспонент). Ясно, что (7) не совпадает с (2) даже при бесконечно большом токе обмена; в последнем случае

(8)
$\eta = {{\left( {iRTt} \right)} \mathord{\left/ {\vphantom {{\left( {iRTt} \right)} {\left[ {F\left( {C_{{\text{s}}}^{0}F + i\left( {{\text{1}}--\alpha } \right)t} \right)} \right]}}} \right. \kern-0em} {\left[ {F\left( {C_{{\text{s}}}^{0}F + i\left( {{\text{1}}--\alpha } \right)t} \right)} \right]}},$
то есть перенапряжение оказывается в любой момент времени меньшим и, кроме того, зависит от времени нелинейно. Это связано с тем, что в реальности ситуация на поверхности не является равновесной. Имеющиеся экспериментальные данные действительно показывают, что линейность зависимости η(t) не соблюдается даже в самый начальный период измерений после включения тока.

Продифференцировав (8) по времени и найдя производную при t → 0, получаем следующее выражение для величины $C_{{\text{s}}}^{0}{\text{:}}$

(9)
$C_{{\text{s}}}^{0} = {{iRT\left[ {{{{\text{1}}--i\left( {{\text{1}}--\alpha } \right)} \mathord{\left/ {\vphantom {{{\text{1}}--i\left( {{\text{1}}--\alpha } \right)} {{{i}_{0}}}}} \right. \kern-0em} {{{i}_{0}}}}} \right]} \mathord{\left/ {\vphantom {{iRT\left[ {{{{\text{1}}--i\left( {{\text{1}}--\alpha } \right)} \mathord{\left/ {\vphantom {{{\text{1}}--i\left( {{\text{1}}--\alpha } \right)} {{{i}_{0}}}}} \right. \kern-0em} {{{i}_{0}}}}} \right]} {\left[ {{{F}^{{\text{2}}}}{{{\left( {{{{\text{d}}\eta } \mathord{\left/ {\vphantom {{{\text{d}}\eta } {{\text{d}}t}}} \right. \kern-0em} {{\text{d}}t}}} \right)}}_{0}}} \right]}}} \right. \kern-0em} {\left[ {{{F}^{{\text{2}}}}{{{\left( {{{{\text{d}}\eta } \mathord{\left/ {\vphantom {{{\text{d}}\eta } {{\text{d}}t}}} \right. \kern-0em} {{\text{d}}t}}} \right)}}_{0}}} \right]}},$
что также отличается от известной формулы $C_{{\text{s}}}^{0} = {{iRT} \mathord{\left/ {\vphantom {{iRT} {[{{F}^{{\text{2}}}}{{{({{{\text{d}}\eta } \mathord{\left/ {\vphantom {{{\text{d}}\eta } {{\text{d}}t}}} \right. \kern-0em} {{\text{d}}t}})}}_{0}}]}}} \right. \kern-0em} {[{{F}^{{\text{2}}}}{{{({{{\text{d}}\eta } \mathord{\left/ {\vphantom {{{\text{d}}\eta } {{\text{d}}t}}} \right. \kern-0em} {{\text{d}}t}})}}_{0}}]}},$ хотя и совпадает с ней при i0 → ∞. Это дает возможность определить равновесную величину заполнения поверхности адатомами.

Формула (8) показывает, что, откладывая η–1 как функцию t–1 можно определить коэффициент переноса, а также – не прибегая к (9), величину $C_{{\text{s}}}^{0},$ так как угол наклона этой зависимости составляет ${{{{F}^{{\text{2}}}}C_{{\text{s}}}^{0}} \mathord{\left/ {\vphantom {{{{F}^{{\text{2}}}}C_{{\text{s}}}^{0}} {iR}}} \right. \kern-0em} {iR}},$ а начальная ордината (1 – α)/RT. Кроме того, такое построение позволяет определить момент начала нуклеации – это происходит, когда функция начинает отклоняться от линейности.

Если отказаться от линеаризации, которая была нами предпринята только с целью более наглядного представления результатов, то можно распространить данное рассуждение на более широкий диапазон перенапряжений. В этом случае обработка экспериментальных данных потребует несколько более сложного численного анализа. Уравнение для определения величин i0, $C_{{\text{s}}}^{0}$ и α при этом принимает вид

(10)
$\begin{gathered} \left( {{{it} \mathord{\left/ {\vphantom {{it} {C_{{\text{s}}}^{0}F}}} \right. \kern-0em} {C_{{\text{s}}}^{0}F}}} \right){\text{exp}}[--({\text{1}}--\alpha )f\eta ] = \\ = {{--i} \mathord{\left/ {\vphantom {{--i} {{{i}_{0}}}}} \right. \kern-0em} {{{i}_{0}}}} + {\text{exp}}(\alpha f\eta )--{\text{exp}}[--({\text{1}}--\alpha )f\eta ]. \\ \end{gathered} $

Это позволяет, используя достаточный массив экспериментально найденных величин η при соответствующих t и заданной плотности тока i (в особенности при нескольких различных плотностях тока), численно определить указанные параметры электрохимического процесса.

3. АНАЛИЗ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ

С помощью приведенных формул мы проанализировали представленные в [2] экспериментальные данные, относящиеся к гальваностатическому осаждению кадмия (с учетом того, что это двухэлектронный процесс). Оказалось, что нуклеация начиналась в этих исследованиях существенно раньше (и, соответственно, при меньших заполнениях), чем предполагали авторы: например, при плотности тока 11 мкА/см2 зародыши возникали не через 100 с после начала накопления адатомов, а уже через 27 с. Найденный по формуле (8) коэффициент переноса составил 0.6, а величина равновесной концентрации адатомов $C_{{\text{s}}}^{0}$оказалась равной 2.4 × 10–10 моль/см2, т.е. около 1.4 × × 1014 атомов на 1 см2. Это примерно в 1.6 раза меньше, чем нашли авторы [2]. Плотность тока обмена по этим данным вычислить не удалось, так как она была как минимум на порядок выше применявшихся плотностей тока (от 11 до 80 мкА/см2). Для ее вычисления необходимо выполнить аналогичные эксперименты при более высоких токах и с миллисекундным разрешением.

Началу процессов нуклеации соответствовали накопленные концентрации $C_{{\text{s}}}^{ * }$ адатомов Cd от 1.5 × 10–9 до 2 × 10–9, что отвечает интервалу перенапряжений от 22 до 26 мВ, в полном согласии с экспериментальными данными [2].

При этом мы полагаем, что непосредственной причиной начала образования зародышей является не достигнутое перенапряжение, а именно достаточно высокая поверхностная концентрация адатомов $C_{{\text{s}}}^{ * }.$ И притом не пересыщение по адатомам по сравнению с начальной величиной, т.е. ${{C_{{\text{s}}}^{ * }} \mathord{\left/ {\vphantom {{C_{{\text{s}}}^{ * }} {C_{{\text{s}}}^{0}}}} \right. \kern-0em} {C_{{\text{s}}}^{0}}},$ а именно сама величина $C_{{\text{s}}}^{ * },$ при которой вероятность нуклеации резко возрастает, в соответствии с нашей моделью, развитой в [7]. Такой подход, в частности, снимает все вопросы, относящиеся к выделению истинной величины кристаллизационного перенапряжения.

Список литературы

  1. Самарцев, А.Г., Евстропьев, К.С. Явления перенапряжения при электролитическом выделении металла. Журн. физической химии. 1934. Т. 5. № 7. С. 854.

  2. Гуцов, И. Кинетика электролитического фазообразования в гальваностатических условиях. Известия Института физикохимии Болгарской АН. 1964. Т. 4. С. 69.

  3. Трофименко, В.В., Коваленко, В.С., Житник, В.П., Лошкарев, Ю.М. Кинетика нестационарного зародышеобразования в гальваностатическом режиме электролиза. Электрохимия. 1983. Т. 19. С. 887

  4. Milchev, A. and Montenegro, M.I., A galvanostatic study of electrochemical Nucleation, J. Electroanal. Chem. 1992, vol. 333, p. 93–102.

  5. Исаев, В.А. Электрохимическое фазообразование. Екатеринбург: УрО РАН, 2007. 123 с.

  6. Isaev, V.A. and Grishenkova, O.V., Galvanostatic phase formation, J. Solid State Electrochemistry, 2014, vol. 18, p. 2383–2386.

  7. Гамбург, Ю.Д. Кинетическая модель электрохимической нуклеации. Электрохимия. 2009. Т. 45. С. 1510. [Gamburg, Yu.D., Kinetic model of electrochemical nucleation, Russ. J. Electrochem., 2009, vol. 45, p. 1397.]

Дополнительные материалы отсутствуют.