Электрохимия, 2019, T. 55, № 5, стр. 637-640
Размер критического зародыша при электрокристаллизации
Ю. Д. Гамбург *
Институт физической химии и электрохимии им. А.Н. Фрумкина РАН
119071 Москва, Ленинский просп., 31, Россия
* E-mail: gamb@list.ru
Поступила в редакцию 14.09.2018
После доработки 26.11.2018
Принята к публикации 12.12.2018
Аннотация
Предложена новая трактовка критического зародыша, основанная на кинетическом анализе нуклеации. Под критическим размером подразумевается такой, при котором поток перехода к следующему классу размеров является наименьшим. Установлена связь критического размера с перенапряжением и с величиной равновесного заполнения поверхности адатомами, либо с величинами удельной поверхностной энергии и энергии адгезии.
ВВЕДЕНИЕ
Размер кристаллических зародышей, образуемых при электрохимической кристаллизации, зависит, как известно [1, 2], прежде всего от перенапряжения η и удельной поверхностной энергии кристалла σ. Объем сферического критического жидкого зародыша, лежащего на абсолютно несмачиваемой поверхности электрода, то есть в отсутствие взаимодействия материалов субстрата и капли, равен
(1)
$V_{{\text{c}}}^{{\hom }} = ({{32\pi } \mathord{\left/ {\vphantom {{32\pi } 3}} \right. \kern-0em} 3}){{\sigma }^{3}}{{V_{{\text{m}}}^{{\text{3}}}} \mathord{\left/ {\vphantom {{V_{{\text{m}}}^{{\text{3}}}} {{{{(zF\eta )}}^{3}}}}} \right. \kern-0em} {{{{(zF\eta )}}^{3}}}},$Физическое взаимодействие материалов зародыша и субстрата друг с другом и с электролитом учитывается путем введения понятия краевого угла (угла смачивания) α, причем объем гетерогенно образованного кластера в виде сферического сегмента равен
(2)
$\begin{gathered} V_{{\text{c}}}^{{{\text{het}}}} = V_{{\text{c}}}^{{\hom }}\left( {0.5--0.75\,{\text{cos}}\alpha + 0.25{\text{co}}{{{\text{s}}}^{3}}\alpha } \right) = \\ = V_{{\text{c}}}^{{\hom }}{\text{Ф }}, \\ \end{gathered} $Угол смачивания непосредственно связан с величиной σ, а также с энергией адгезии металла зародыша к основе σs, а именно σs = σ (1 + cos α). Комбинируя это выражение с (3), получаем
(4)
$\Phi \approx 1 - {{{{\sigma }_{{\text{s}}}}} \mathord{\left/ {\vphantom {{{{\sigma }_{{\text{s}}}}} {2\sigma }}} \right. \kern-0em} {2\sigma }}.$Формула (4) для кластера в виде сферического сегмента является, как мы видим, приближенной. Однако, как строго показано в [2], для твердых цилиндрических зародышей, находящихся на инородном субстрате, она является точной; при этом правда, возникает дополнительный численный коэффициент, в результате чего вместо множителя (32π/3) появляется число 2К того же порядка, зависящее от типа кристаллографической грани, которая граничит с субстратом.
Итак, объем гетерогенно образованного трехмерного критического зародыша равен, как это следует из (1), (2) и (4),
(5)
$V_{{\text{c}}}^{{{\text{het}}}} = {\text{K}}{{\sigma }^{2}}(2\sigma --{{\sigma }_{{\text{s}}}}){{{{v}^{3}}} \mathord{\left/ {\vphantom {{{{v}^{3}}} {{{{(ze\eta )}}^{3}}}}} \right. \kern-0em} {{{{(ze\eta )}}^{3}}}},$(6)
$g = {\text{K}}{{v}^{2}}{{\sigma }^{2}}{{(2\sigma --{{\sigma }_{{\text{s}}}})} \mathord{\left/ {\vphantom {{(2\sigma --{{\sigma }_{{\text{s}}}})} {{{{(ze\eta )}}^{3}}}}} \right. \kern-0em} {{{{(ze\eta )}}^{3}}}},$(7)
$g = 20{{s}^{3}}{{\sigma }^{3}}{{\text{Ф }} \mathord{\left/ {\vphantom {{\text{Ф }} {{{{(ze\eta )}}^{3}}}}} \right. \kern-0em} {{{{(ze\eta )}}^{3}}}}\quad = 10{{s}^{3}}{{\sigma }^{2}}{{(2\sigma --{{\sigma }_{{\text{s}}}})} \mathord{\left/ {\vphantom {{(2\sigma --{{\sigma }_{{\text{s}}}})} {{{{(ze\eta )}}^{3}}}}} \right. \kern-0em} {{{{(ze\eta )}}^{3}}}}.$Из (7) следует, что при возрастании перенапряжения величина g монотонно и быстро уменьшается, и при некотором η = η* становится равной единице. Этот момент имеет важное значение для электроосаждения металлов, так как при более высоких перенапряжениях кристаллизация идет фактически без образования зародышей: всякий вновь образованный атом металла становится зародышем. При более высоких величинах η уравнение (7), как и (1), (2), теряет физический смысл. Приравнивая в (7) величину g единице, получаем в случае Ф = 1 соответствующее перенапряжение η*
(8)
$\eta * \approx 2.7s{\sigma \mathord{\left/ {\vphantom {\sigma {\left( {ze} \right)}}} \right. \kern-0em} {\left( {ze} \right)}}$При двумерной нуклеации, которая имеет место в случае достаточно высокой энергии адгезии зародыша к субстрату, вместо (8) для η* получается сходное выражение:
(9)
$\eta * \approx s{{\left( {3\sigma - {{\sigma }_{{\text{s}}}}} \right)} \mathord{\left/ {\vphantom {{\left( {3\sigma - {{\sigma }_{{\text{s}}}}} \right)} {\left( {ze} \right)}}} \right. \kern-0em} {\left( {ze} \right)}},$Это выражение следует признать лучше соответствующим физическому смыслу, так как кластеры уже из 2–3 частиц всегда являются двумерными, что делает естественным переход к одной частице. Уравнение (9) обладает и тем преимуществом, что в нем в явном виде участвует величина энергии адгезии частицы к поверхности σs.
Тем не менее, приведенные уравнения плохо согласуются с экспериментальными данными. Например, для зародыша серебра на хорошо смачиваемой основе (9) дает сильно завышенную величину η* ≈ 1 В. Поэтому следует искать другие методы определения критического размера зародышевого кластера.
КИНЕТИЧЕСКОЕ ОПРЕДЕЛЕНИЕ КРИТИЧЕСКОГО КЛАСТЕРА
В наших работах [4, 5] было установлено, что при стохастическом процессе разряда ионов на поверхности электрода с образованием адатомов, стационарную скорость J образования кластеров из g адатомов (g-кластеров) можно выразить как
(10)
$\ln J = g\ln {{\theta }_{0}} + \ln \left( {g!} \right) + \left[ {\left( {g + \alpha } \right){{zF\eta } \mathord{\left/ {\vphantom {{zF\eta } {RT}}} \right. \kern-0em} {RT}}} \right] + C,$(11)
${{{\text{d}}\ln J} \mathord{\left/ {\vphantom {{{\text{d}}\ln J} {{\text{d}}\eta }}} \right. \kern-0em} {{\text{d}}\eta }} = \left( {g + \alpha } \right){{zF\eta } \mathord{\left/ {\vphantom {{zF\eta } {RT}}} \right. \kern-0em} {RT}},$Формула (11) полностью совпадает с выражением, которое из совершенно других соображений выводится в атомистической теории нуклеации [1].
Но в [4, 5] не было приведено никаких объяснений, почему именно g*, а не какая-либо иная величина g соответствует критическому зародышу. Было отмечено только, что именно такие кластеры по какой-то причине оказываются способными к преимущественному дальнейшему росту. Следует подчеркнуть, что при выводе приведенных формул не использовались никакие понятия и соображения, относящиеся к энергетике нуклеации, удельной поверхностной энергии, равновесной и неравновесной функциям распределения кластеров по размерам, и т.д.
Здесь необходимо обратить внимание на то, что существование критической величины g* в данной концепции вызвано очень простой причиной. Для того, чтобы это стало очевидным, продифференцируем выражение (11) по g при постоянном перенапряжении. Такое дифференцирование, если воспользоваться тем, что dln(g!)/dg ≈ ln g, дает
(12)
${{{\text{d}}\ln J} \mathord{\left/ {\vphantom {{{\text{d}}\ln J} {{\text{d}}g}}} \right. \kern-0em} {{\text{d}}g}} = \ln {{\theta }_{0}} + \ln g + {{zF\eta } \mathord{\left/ {\vphantom {{zF\eta } {RT}}} \right. \kern-0em} {RT}},$(13)
$g* = \left( {{1 \mathord{\left/ {\vphantom {1 {{{\theta }_{0}}}}} \right. \kern-0em} {{{\theta }_{0}}}}} \right)\exp \left( { - {{zF\eta } \mathord{\left/ {\vphantom {{zF\eta } {RT}}} \right. \kern-0em} {RT}}} \right).$При данной величине g величина потока J проходит через минимум. Следовательно, именно переход от g*-клаcтера к (g + 1)-кластеру является лимитирующим, т.е. определяет скорость всего дальнейшего процесса роста (которая в стационарном состоянии равна скорости этого перехода). Это фактически совпадает с определением критического зародыша: именно начиная с этого размера начинается устойчивый рост, скорость которого может быть (в принципе) выше скорости указанного перехода.
Сопоставление формулы (13) с имеющимися немногочисленными экспериментальными данными о размере критических зародышей показывает следующее. В [6], например, показано (независимым методом), что для критического зародыша серебра на подложке из стеклоуглерода при перенапряжении 25 мВ g = 4. Поскольку в этом случае (–zFη/RT) ≈ 1, то θ0 ≈ 0.09, что является вполне правдоподобной величиной. К сожалению, приводимые в [6] данные не позволяют определить истинную величину θ0, поэтому такое сопоставление не вполне корректно.
Итак, формула (16) определяет размер критического зародыша, соответствующего заданному перенапряжению. Она позволяет также найти величину того достаточно большого перенапряжения η*, при котором критический зародыш будет состоять из одного атома: (1/θ0)exp(–zFη*/RT) = 1, откуда
(14)
$\eta * = - \left( {{{RT} \mathord{\left/ {\vphantom {{RT} {zF}}} \right. \kern-0em} {zF}}} \right)\ln {{\theta }_{0}} = - \left( {{{{{k}_{{\text{B}}}}T} \mathord{\left/ {\vphantom {{{{k}_{{\text{B}}}}T} {ze}}} \right. \kern-0em} {ze}}} \right)\ln {{\theta }_{0}}.$Для условий [6] эта формула дает перенапряжение 0.1 В, что соответствует данным указанной работы.
При развитом здесь подходе искомое перенапряжение зависит только от величины равновесного заполнения θ0, а не от удельных энергий σ и σs, как в уравнении (9). Указанные величины, однако, связаны между собой уравнением изотермы адсорбции адатомов выделяемого металла на поверхности: θ0 = f(σ, σs). При этом качественно очевидно, что чем больше энергия адгезии σs, тем выше θ0. Поэтому, согласно уравнению (14), при увеличении энергии адгезии уменьшается η*. То же самое следует и из “классического” уравнения (9).
Форма уравнения θ0 = f(σ, σs) неизвестна. Если постулировать, что
(15)
$\ln {{\theta }_{0}} = s{{\left( {3\sigma - {{\sigma }_{{\text{s}}}}} \right)} \mathord{\left/ {\vphantom {{\left( {3\sigma - {{\sigma }_{{\text{s}}}}} \right)} {{{k}_{{\text{B}}}}T}}} \right. \kern-0em} {{{k}_{{\text{B}}}}T}},$НАХОЖДЕНИЕ РАВНОВЕСНОГО ЗАПОЛНЕНИЯ
Величину θ0 для подстановки в (13) можно определить экспериментально, например, с помощью контроля потенциала электрода после включения постоянного тока. Действительно, после включения катодного тока поверхностная концентрация адатомов (вплоть до момента начала нуклеации) увеличивается согласно
(18)
${{C}_{{\text{s}}}} = C_{{\text{s}}}^{0} + {{it} \mathord{\left/ {\vphantom {{it} {zF}}} \right. \kern-0em} {zF}},$(19)
$\begin{gathered} i\left( t \right) = {{i}_{{\text{c}}}}\left( t \right) + {{i}_{{\text{a}}}}\left( t \right) = {{i}_{0}}\left\{ {\exp \left( {{{\alpha zF\eta \left( t \right)} \mathord{\left/ {\vphantom {{\alpha zF\eta \left( t \right)} {RT}}} \right. \kern-0em} {RT}}} \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)\exp \left[ { - \left( {1 - \alpha } \right)} \right]{{zF\eta \left( t \right)} \mathord{\left/ {\vphantom {{zF\eta \left( t \right)} {RT}}} \right. \kern-0em} {RT}}} \right\} \\ \end{gathered} $(20)
$\begin{gathered} \left( {{{it} \mathord{\left/ {\vphantom {{it} {C_{{\text{s}}}^{0}F}}} \right. \kern-0em} {C_{{\text{s}}}^{0}F}}} \right)\exp \left[ { - \left( {1 - \alpha } \right){{zF\eta } \mathord{\left/ {\vphantom {{zF\eta } {RT}}} \right. \kern-0em} {RT}}} \right] = \\ = - {i \mathord{\left/ {\vphantom {i {{{i}_{0}} + }}} \right. \kern-0em} {{{i}_{0}} + }}\exp \left( {{{\alpha zF\eta } \mathord{\left/ {\vphantom {{\alpha zF\eta } {RT}}} \right. \kern-0em} {RT}}} \right) - \\ - \exp \left[ { - \left( {1 - \alpha } \right){{zF\eta } \mathord{\left/ {\vphantom {{zF\eta } {RT}}} \right. \kern-0em} {RT}}} \right], \\ \end{gathered} $Если применяемая при таких измерениях плотность тока существенно ниже плотности тока обмена, уравнение (20) можно линеаризовать, что сильно упрощает анализ получаемых данных.
Необходимо еще раз подчеркнуть, что уравнение (20) соблюдается только при накоплении адатомов на поверхности вплоть до начала роста кластеров. В этот момент данная зависимость начинает искажаться, что дает дополнительную информацию о концентрации адатомов на поверхности в данный момент, а также о длительности процесса накопления адатомов на поверхности.
Величина $C_{{\text{s}}}^{0}$ непосредственно связана с равновесным заполнением θ0. Для вычисления последнего нужно только знать полное количество мест для посадки адатомов на единице площади поверхности N: тогда равновесное заполнение определится как
ЗАКЛЮЧЕНИЕ
Развита концепция, связывающая число атомов в критическом зародыше g с минимальным значением потока g → g + 1, т.е. с лимитирующей стадией процесса роста кластера. При этом число g зависит как от перенапряжения, так и от равновесной степени заполнения поверхности электрода адатомами осаждаемого металла.
Список литературы
Milchev, A., Electrocrystallization: fundamentals of nucleation and growth. Kluwer Academic Publishing, 2002.
Gamburg, Yu.D. and Zangari, G., Theory and Practice of Metal Electrodeposition. Springer, 2011 [Гамбург, Ю.Д., Зангари, Дж. Теория и практика электроосаждения металлов. М., Бином, 2015.]
Дамаскин, Б.Б., Петрий, О.А., Цирлина, Г.А. Электрохимия., М., Химия, Колос С, 2008.
Гамбург, Ю.Д. Кинетическая модель электрохимической нуклеации. Электрохимия. 2009. Т. 45. С. [Gamburg, Yu.D., Kinetic Model of Electrochemical Nucleation, Russ. J. of Electrochemistry, 2009, vol. 45, p. 1397.]
Gamburg, Yu.D., Some novel efforts to describe the nucleation and growth at electrocrystallization, J. of Solid State Electrochemistry, 2012, vol. 17, pp. 353–359.
Milchev, A., Vassileva, E., and Kertov, V., Electrolytic nucleation of silver on a glassy carbon electrode, J. of Electroanal. Chem., 1980, vol. 107, pp. 323–345.
Дополнительные материалы отсутствуют.
Инструменты
Электрохимия