Электрохимия, 2019, T. 55, № 5, стр. 637-640

Размер критического зародыша при электрокристаллизации

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

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

* E-mail: gamb@list.ru

Поступила в редакцию 14.09.2018
После доработки 26.11.2018
Принята к публикации 12.12.2018

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

Аннотация

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

Ключевые слова: электрокристаллизация, кристаллические зародыши, кластеры, перенапряжение

ВВЕДЕНИЕ

Размер кристаллических зародышей, образуемых при электрохимической кристаллизации, зависит, как известно [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}}}},$
где Vm – мольный объем, zF – количество электричества, переносимое при образовании моля металла. Эта формула непосредственно следует из формулы (10.6.7), приведенной в учебнике [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} $
где всегда Ф ≤ 1. В линейном приближении этому соответствует зависимость
(3)
$\Phi \approx 0.5\left( {1 - \cos \alpha } \right),$
которая получается из (2) заменой 0.25cos3α на 0.25cos α.

Угол смачивания непосредственно связан с величиной σ, а также с энергией адгезии металла зародыша к основе σ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}}}},$
(в последнем выражении мы перешли от мольных единиц к атомным: $\text{v}$ это объем, приходящийся на один атом, е – заряд электрона). Теперь, поделив обе части (5) на $\text{v},$ получаем выражение для количества атомов g в критическом зародыше:
(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}}}},$
причем, например, для зародыша, контактирующего с субстратом гранью (111), K = $8\sqrt 3 $. Если пренебречь кристаллографическими особенностями, то для цилиндрического зародыша К = 4π, что мало отличается от 8√3. Далее мы можем для удобства дальнейших расчетов перейти от объема $v$ к площади поверхности соприкосновения одного атома зародыша с основой s (для плотноупакованной структуры, лежащей на поверхности плотноупакованной гранью, s ≈ 1.09$v$2/3, т.е. ${{v}^{2}}$ ≈ ≈ 0.77s3), и тогда с точностью до единицы в коэффициенте

(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)}}$
(при Ф = 0.5 коэффициент снижается c 2.7 до 2.15, что указывает на роль σs).

При двумерной нуклеации, которая имеет место в случае достаточно высокой энергии адгезии зародыша к субстрату, вместо (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)}},$
где коэффициенты при σ и σS (как и сами эти величины) несколько зависят от кристаллографических факторов [1].

Это выражение следует признать лучше соответствующим физическому смыслу, так как кластеры уже из 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,$
где θ0 – равновесная степень заполнения поверхности адатомами, а остальные обозначения стандартны. При этом J определяется как общее количество частиц, присоединяемых в единицу времени ко всем n-кластерам, находящимся на единице поверхности. В [4] отмечалось, что с помощью соотношения, которое следует из (10), а именно
(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}},$
можно по экспериментальным измерениям скорости нуклеации (при нескольких близких потенциалах) определить величину g для критического зародыша, т.е. g* (в данном интервале потенциалов), а если эти измерения выполнить при нескольких более существенно различающихся потенциалах, то получить и другую информацию.

Формула (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}},$
а поскольку θ0$ \ll $ 1, т.е. ln θ0$ \ll $ 0, то при некотором g = g* данная производная проходит через нуль: ln g* + zFη/RT + ln θ0 = 0, или
(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}},$
то (14) принимает вид, совпадающий с (9); при этом формула (15) позволяет записать уравнение для числа атомов g* в критическом зародыше, находящемся на подложке, в очень наглядном виде
(16)
${{k}_{{\text{B}}}}T\ln g* \approx - ze\eta + s\left( {3\sigma - {{\sigma }_{{\text{s}}}}} \right),$
который, однако, представляется неполным и дает сильно завышенную величину g*. По нашему мнению, здесь требуют уточнения коэффициенты. Перенапряжение также следует исправить на некоторую величину, т.е. полное уравнение должно иметь вид
(17)
${{k}_{{\text{B}}}}T\ln g* \approx - ze\left( {\eta + \Delta \eta } \right) + s\left( {A\sigma - B{{\sigma }_{{\text{s}}}}} \right),$
где параметры Δη, A и B неизвестны. Поэтому форма (13) является предпочтительной.

НАХОЖДЕНИЕ РАВНОВЕСНОГО ЗАПОЛНЕНИЯ

Величину θ0 для подстановки в (13) можно определить экспериментально, например, с помощью контроля потенциала электрода после включения постоянного тока. Действительно, после включения катодного тока поверхностная концентрация адатомов (вплоть до момента начала нуклеации) увеличивается согласно

(18)
${{C}_{{\text{s}}}} = C_{{\text{s}}}^{0} + {{it} \mathord{\left/ {\vphantom {{it} {zF}}} \right. \kern-0em} {zF}},$
причем величина плотности тока i связана с величиной текущего перенапряжения η(t) по уравнению Батлера–Фольмера как
(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} $
(катодный ток мы приняли положительным, как и перенапряжение). С учетом (17) это уравнение принимает вид
(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} $
что позволяет, используя достаточный массив экспериментально найденных величин η при соответствующих t и нескольких различных плотностях тока численно определить параметры электрохимического процесса величин i0, $C_{{\text{s}}}^{0}$ и α.

Если применяемая при таких измерениях плотность тока существенно ниже плотности тока обмена, уравнение (20) можно линеаризовать, что сильно упрощает анализ получаемых данных.

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

Величина $C_{{\text{s}}}^{0}$ непосредственно связана с равновесным заполнением θ0. Для вычисления последнего нужно только знать полное количество мест для посадки адатомов на единице площади поверхности N: тогда равновесное заполнение определится как

(21)
${{\theta }_{0}} = {{C_{{\text{s}}}^{0}{{N}_{{\text{A}}}}} \mathord{\left/ {\vphantom {{C_{{\text{s}}}^{0}{{N}_{{\text{A}}}}} N}} \right. \kern-0em} N}\,\,({{N}_{{\text{A}}}}--{\text{ч и с л о А в о г а д р о }}).$

ЗАКЛЮЧЕНИЕ

Развита концепция, связывающая число атомов в критическом зародыше g с минимальным значением потока gg + 1, т.е. с лимитирующей стадией процесса роста кластера. При этом число g зависит как от перенапряжения, так и от равновесной степени заполнения поверхности электрода адатомами осаждаемого металла.

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

  1. Milchev, A., Electrocrystallization: fundamentals of nucleation and growth. Kluwer Academic Publishing, 2002.

  2. Gamburg, Yu.D. and Zangari, G., Theory and Practice of Metal Electrodeposition. Springer, 2011 [Гамбург, Ю.Д., Зангари, Дж. Теория и практика электроосаждения металлов. М., Бином, 2015.]

  3. Дамаскин, Б.Б., Петрий, О.А., Цирлина, Г.А. Электрохимия., М., Химия, Колос С, 2008.

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

  5. 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.

  6. 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.

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