Расплавы, 2020, № 3, стр. 275-281

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

О. В. Гришенкова a***, А. В. Косов a, Ю. П. Зайков a, В. А. Исаев a

a Институт высокотемпературной электрохимии УрО РАН
Екатеринбург, Россия

* E-mail: olagris@mail.ru
** E-mail: o.grishenkova@ihte.uran.ru

Поступила в редакцию 02.02.2020
После доработки 11.02.2020
Принята к публикации 18.02.2020

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

Аннотация

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

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

ВВЕДЕНИЕ

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

Кинетика начальных стадий образования бинарного сплава на индифферентном электроде при электроосаждении в потенциостатических условиях была проанализирована в рамках различных подходов и приближений. В работе [5] на основе атомистической теории электрохимического сплавообразования [6, 7] был предложен метод определения размера и состава критических зародышей сплава по зависимостям стационарной скорости нуклеации от потенциала и концентрации. Модель, описывающая лимитируемый диффузией и разрядом рост зародыша сплава при высоких перенапряжениях на базе классической теории электрокристаллизации, была разработана в [8]. Аналитические выражения для тока при мгновенной и прогрессирующей нуклеации в условиях кинетического или смешанного контроля роста изолированных зародышей были представлены в статье [9]. Díaz-Morales с соавторами [10] описали кинетику нуклеации и диффузионно-контролируемого роста биметаллической фазы, получили аналитические выражения для определения плотности активных центров на электроде и скорости нуклеации, используя модель Scharifker–Mostany [11].

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

ТЕОРИЯ

Диффузия к единственному полусферическому зародышу при постоянном перенапряжении η может быть описана уравнением второго закона Фика в сферических координатах

(1)
$\frac{{\partial c(r,t)}}{{\partial t}} = D\left[ {\frac{{{{\partial }^{2}}c(r,t)}}{{\partial {{r}^{2}}}} + \frac{2}{r}\frac{{\partial c(r,t)}}{{\partial r}}} \right]$
с начальным и граничными условиями $c(r,0) = {{c}_{0}},$ $\mathop {\lim }\limits_{r \to \infty } c(r,t) = {{c}_{0}},$ $c({{r}_{0}},t) = {{c}_{{{\text{sr}}}}},$ где c(r, t) – концентрация осаждаемых ионов; r – расстояние от центра зародыша; t – время; D – коэффициент диффузии; c0 и csr – концентрация осаждаемых ионов в объеме электролита и у поверхности зародыша, соответственно; r0 – радиус зародыша. Решение этого дифференциального уравнения с использованием преобразования Лапласа дает профиль концентрации [12, 13]:
(2)
$с(r,t) = ({{c}_{0}} - {{c}_{{{\text{sr}}}}})\left[ {1 - \frac{{{{r}_{0}}}}{r}{\text{erfc}}\frac{{r - {{r}_{0}}}}{{2{{{(Dt)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}}} \right].$
Отсюда для плотности тока роста зародыша, $i = zeD{{({{\partial c} \mathord{\left/ {\vphantom {{\partial c} {\partial r}}} \right. \kern-0em} {\partial r}})}_{{r = {{r}_{0}}}}},$ имеем
(3)
$i = zeD({{c}_{0}} - {{c}_{{{\text{sr}}}}})\left[ {\frac{1}{{{{r}_{0}}}} + \frac{1}{{{{{(\pi Dt)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}}} \right],$
где z – валентность осаждаемых ионов; e – элементарный электрический заряд. Учитывая зависимость csr от перенапряжения, а также тот факт, что диффузию к зародышу малого размера, характерного для начальных стадий электрокристаллизации, с хорошим приближением можно считать стационарной, запишем уравнение (3) в виде [14]
(4)
$i = \frac{{zeD{{c}_{0}}}}{{{{r}_{0}}}}\left[ {1 - \exp ( - f{\eta })} \right],$
где $f = {{ze} \mathord{\left/ {\vphantom {{ze} {kT}}} \right. \kern-0em} {kT}};$ k – постоянная Больцмана; T – температура. Поскольку обычно процессы нуклеации на чужеродной подложке происходят при достаточно высоких перенапряжениях, формулы (3) и (4) можно упростить, полагая ${{c}_{0}} - {{c}_{{{\text{sr}}}}} \approx {{c}_{0}}$ и $1 - \exp ( - f{\eta }) \approx 1.$

Теперь допустим, что в электролите присутствуют не один, а два сорта осаждаемых ионов ${{{\text{A}}}^{{{{z}_{{\text{A}}}} + }}}$ и ${{{\text{B}}}^{{{{z}_{{\text{B}}}} + }}}$ с коэффициентами диффузии DA и DB и концентрациями cA и cB, соответственно, в результате восстановления которых на поверхности индифферентного катода образуются зародыши бинарного сплава AB. Если рост зародыша сплава контролируется диффузией, то для каждого из компонентов сплава можно записать выражение, подобное уравнению (1), т.е. $\frac{{\partial {{c}_{{\text{A}}}}}}{{\partial t}}$ = ${{D}_{{\text{A}}}}\left[ {\frac{{{{\partial }^{2}}{{c}_{{\text{A}}}}}}{{\partial {{r}^{2}}}} + \frac{2}{r}\frac{{\partial {{c}_{{\text{A}}}}}}{{\partial r}}} \right]$ и $\frac{{\partial {{c}_{{\text{B}}}}}}{{\partial t}}$ = = ${{D}_{{\text{B}}}}\left[ {\frac{{{{\partial }^{2}}{{c}_{{\text{B}}}}}}{{\partial {{r}^{2}}}} + \frac{2}{r}\frac{{\partial {{c}_{{\text{B}}}}}}{{\partial r}}} \right].$ Предположим, аналогично [10], что ${{c}_{{\text{A}}}} = \gamma {{c}_{{\text{B}}}} = c;$ тогда при суммировании этих уравнений получим

(5)
$\frac{{\partial c}}{{\partial t}} = \tilde {D}\left[ {\frac{{{{\partial }^{2}}c}}{{\partial {{r}^{2}}}} + \frac{2}{r}\frac{{\partial c}}{{\partial r}}} \right],$
где $\tilde {D}$ = ${{(\gamma {{D}_{{\text{A}}}} + {{D}_{{\text{B}}}})} \mathord{\left/ {\vphantom {{(\gamma {{D}_{{\text{A}}}} + {{D}_{{\text{B}}}})} {(\gamma + 1)}}} \right. \kern-0em} {(\gamma + 1)}}$ – коэффициент диффузии “псевдоионов” ${\text{A}}{{{\text{B}}}^{{\tilde {z} + }}}$, зависящий от соотношения концентраций осаждаемых ионов. Заряд этих частиц также будет зависеть от доли ионов ${{{\text{A}}}^{{{{z}_{{\text{A}}}} + }}}$ и ${{{\text{B}}}^{{{{z}_{{\text{B}}}} + }}}$ в объеме электролита, т.е. ${{x}_{{\text{A}}}} = {{{{c}_{0}}_{{\text{A}}}} \mathord{\left/ {\vphantom {{{{c}_{0}}_{{\text{A}}}} {({{c}_{0}}_{{\text{A}}} + {{c}_{0}}_{{\text{B}}})}}} \right. \kern-0em} {({{c}_{0}}_{{\text{A}}} + {{c}_{0}}_{{\text{B}}})}}$ и ${{x}_{{\text{B}}}}$ = = ${{{{c}_{0}}_{{\text{B}}}} \mathord{\left/ {\vphantom {{{{c}_{0}}_{{\text{B}}}} {({{c}_{0}}_{{\text{A}}} + {{c}_{0}}_{{\text{B}}})}}} \right. \kern-0em} {({{c}_{0}}_{{\text{A}}} + {{c}_{0}}_{{\text{B}}})}}.$ Тогда $\tilde {z} = {{z}_{{\text{A}}}}{{x}_{{\text{A}}}} + {{z}_{{\text{B}}}}{{x}_{{\text{B}}}}$ или $\tilde {z}$ = ${{(\gamma {{z}_{{\text{A}}}} + {{z}_{{\text{B}}}})} \mathord{\left/ {\vphantom {{(\gamma {{z}_{{\text{A}}}} + {{z}_{{\text{B}}}})} {(\gamma + 1)}}} \right. \kern-0em} {(\gamma + 1)}}$ при условии, что ${{c}_{0}}_{{\text{A}}} = \gamma {{c}_{0}}_{{\text{B}}} = {{c}_{0}}.$ При таком подходе вместо уравнения (4) получим
(6)
$i = \frac{{\tilde {z}e\tilde {D}{{c}_{0}}}}{{{{r}_{0}}}}.$
Радиус r0 полусферического кластера, растущего в условиях диффузионного контроля при достаточно высоких перенапряжениях, когда ${{c}_{{{\text{sr}}}}} \to 0,$ равен ${{r}_{0}} = {{(2D{{c}_{0}}\upsilon \tau )}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$ [15], а в нашем случае
(7)
${{r}_{0}} = {{(2\tilde {D}{{c}_{0}}\tilde {\upsilon }\tau )}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}},$
где $\tilde {\upsilon } = {{\upsilon }_{{\text{A}}}}{{x}_{{\text{A}}}} + {{\upsilon }_{{\text{B}}}}{{x}_{{\text{B}}}}$ = ${{(\gamma {{\upsilon }_{{\text{A}}}} + {{\upsilon }_{{\text{B}}}})} \mathord{\left/ {\vphantom {{(\gamma {{\upsilon }_{{\text{A}}}} + {{\upsilon }_{{\text{B}}}})} {(\gamma + 1)}}} \right. \kern-0em} {(\gamma + 1)}};$ υA и υB – объемы атомов А и В, соответственно (считаем, что объем новой фазы аддитивно складывается из объемов отдельных компонентов с учетом их мольных долей); τ – время жизни кластера.

При множественной нуклеации вокруг каждого зародыша формируется полусферическая диффузионная зона, размер которой увеличивается по мере роста зародыша. Это приводит к постепенному перекрытию диффузионных полей и изменению условий массопереноса. Задача о перекрытии 3D зон может быть решена только в рамках некоторых приближений [11, 1624]. Наиболее простой вариант решения был предложен в модели Scharifker–Hills (SH) [16, 17]: 3D поля были заменены эквивалентными плоскими диффузионными зонами, степень перекрытия (θ) которых обуславливала переход от сферической диффузии к линейной. Для определения θ была использована теория кристаллизации Kolmogorov–Johnson–Mehl–Avrami [25, 26]. Применив аналогичный подход к формированию зародышей сплава AB, получим

(8)
$i = \tilde {z}e{{c}_{0}}{{({{\tilde {D}} \mathord{\left/ {\vphantom {{\tilde {D}} {\pi t}}} \right. \kern-0em} {\pi t}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}\theta ,\,\,\,\,\theta = 1 - \exp ( - {{\theta }_{{{\text{ex}}}}}),$
где θex – доля поверхности электрода, покрытая плоскими диффузионными зонами радиусом rd${{(\tilde {D}\tau )}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$ без учета их перекрытия. Тогда
(9)
${{\theta }_{{{\text{ex}}}}} = \pi r_{{\text{d}}}^{2}{{N}_{0}}\,\,\,\,\left( {{\text{мгновенная}}\,\,{\text{нуклеация}}} \right),$
(10)
${{\theta }_{{{\text{ex}}}}} = \pi {{J}_{0}}\int\limits_0^t {r_{{{\text{d }}}}^{2}{\text{d}}u} \,\,\,\,\left( {{\text{прогрессирующая}}\,\,{\text{нуклеация}}} \right),$
где N0 – максимальная плотность зародышей на электроде (равна количеству активных центров на единице поверхности электрода); J0 – стационарная скорость нуклеации. Из уравнений (8)–(10) имеем
(11)
$i = \tilde {z}e{{c}_{0}}{{({{\tilde {D}} \mathord{\left/ {\vphantom {{\tilde {D}} {\pi t}}} \right. \kern-0em} {\pi t}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}[1 - \exp ( - \pi {\kern 1pt} \tilde {k}\tilde {D}{{N}_{0}}t)]\,\,\,\,\left( {{\text{мгновенная}}\,\,{\text{нуклеация}}} \right),$
(12)
$i = \tilde {z}e{{c}_{0}}{{({{\tilde {D}} \mathord{\left/ {\vphantom {{\tilde {D}} {\pi t}}} \right. \kern-0em} {\pi t}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}[1 - \exp ({{ - 2\pi {\kern 1pt} \tilde {k}\tilde {D}{{J}_{0}}{{t}^{2}}} \mathord{\left/ {\vphantom {{ - 2\pi {\kern 1pt} \tilde {k}\tilde {D}{{J}_{0}}{{t}^{2}}} 3}} \right. \kern-0em} 3})]\,\,\,\,\left( {{\text{прогрессирующая}}\,\,{\text{нуклеация}}} \right),$
где $\tilde {k} = {{(8\pi {{c}_{0}}\tilde {\upsilon })}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$ – коэффициент пропорциональности, найденный из условия, что плотности тока, описываемые уравнениями (11) и (12) при малых t, когда θ ≈ θex, должны быть равны плотностям тока роста изолированных кластеров при соответствующих режимах нуклеации:
(13)
$i = \pi {\kern 1pt} \tilde {z}e{{\tilde {\upsilon }}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}{{(2\tilde {D}{{c}_{0}})}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}}{{N}_{0}}{{t}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}\,\,\,\,\left( {{\text{мгновенная}}\,\,{\text{нуклеация}}} \right),$
(14)
$i = 2\pi \tilde {z}e{\kern 1pt} {{\tilde {\upsilon }}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}{{(2\tilde {D}{{c}_{0}})}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}}{{{{J}_{0}}{{t}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}}} \mathord{\left/ {\vphantom {{{{J}_{0}}{{t}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}}} 3}} \right. \kern-0em} 3}\,\,\,\,\left( {{\text{прогрессирующая}}\,\,{\text{нуклеация}}} \right).$
При больших t (θ ≈ 1), выражения (11) и (12) сводятся к $i = \tilde {z}e{{c}_{0}}{{({{\tilde {D}} \mathord{\left/ {\vphantom {{\tilde {D}} {\pi t}}} \right. \kern-0em} {\pi t}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}},$ т.е. к виду, подобному уравнению Котрелла, справедливому для линейной диффузии к плоскому электроду. Таким образом, этот подход позволяет правильно представить временны́ е зависимости тока (потенциостатические транзиенты тока) в предельных случаях.

Транзиенты плотности тока i(t), описываемые уравнениями (11) и (12), имеют максимум (im, tm), координаты которого могут быть найдены из условия di/dt = 0. Для мгновенной нуклеации имеем:

(15)
${{t}_{{\text{m}}}} = {{1.2564} \mathord{\left/ {\vphantom {{1.2564} {\pi \tilde {k}{{N}_{0}}\tilde {D}}}} \right. \kern-0em} {\pi \tilde {k}{{N}_{0}}\tilde {D}}}.$
Подставляя (15) в уравнение (11), получим
(16)
${{i}_{{\text{m}}}} = \tilde {z}e{{c}_{0}}{{({{\tilde {D}} \mathord{\left/ {\vphantom {{\tilde {D}} {\pi {{t}_{{\text{m}}}}}}} \right. \kern-0em} {\pi {{t}_{{\text{m}}}}}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}\left[ {1 - \exp ( - 1.2564)} \right],$
(17)
$i_{{\text{m}}}^{2}{{t}_{{\text{m}}}} = 0.1629{{(\tilde {z}e{{c}_{0}})}^{2}}\tilde {D},$
(18)
$\frac{{{{i}^{2}}}}{{i_{{\text{m}}}^{2}}} = 1.9543\frac{{{{t}_{{\text{m}}}}}}{t}{{\left[ {1 - \exp \left( { - 1.2564\frac{t}{{{{t}_{{\text{m}}}}}}} \right)} \right]}^{2}}.$
Подобным образом, с учетом уравнения (12), для прогрессирующей нуклеации имеем:
(19)
${{t}_{{\text{m}}}} = {{({{3.5050} \mathord{\left/ {\vphantom {{3.5050} {\pi \tilde {k}{{J}_{0}}\tilde {D}}}} \right. \kern-0em} {\pi \tilde {k}{{J}_{0}}\tilde {D}}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}},$
(20)
${{i}_{{\text{m}}}} = \tilde {z}e{{c}_{0}}{{({{\tilde {D}} \mathord{\left/ {\vphantom {{\tilde {D}} {\pi {{t}_{{\text{m}}}}}}} \right. \kern-0em} {\pi {{t}_{{\text{m}}}}}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}\left[ {1 - \exp ( - 2.3367)} \right],$
(21)
$i_{{\text{m}}}^{2}{{t}_{{\text{m}}}} = 0.2598{{(\tilde {z}e{{c}_{0}})}^{2}}\tilde {D},$
(22)
$\frac{{{{i}^{2}}}}{{i_{{\text{m}}}^{2}}} = 1.2254\frac{{{{t}_{{\text{m}}}}}}{t}{{\left[ {1 - \exp \left( { - 2.3367\frac{{{{t}^{2}}}}{{t_{{\text{m}}}^{2}}}} \right)} \right]}^{2}}.$
Формулы (18) и (22) совпадают с классическими зависимостями (i/im)2 от (t/tm) SH модели [16, 17], полученными для электроосаждения индивидуального металла; они позволяют определить механизм нуклеации в случае диффузионного контроля роста, сравнивая соответствующие экспериментальные и теоретические безразмерные зависимости. Выражения (15), (17) и (21) подобны уравнениям SH модели, но содержат комбинированные значения $\tilde {z},$ $\tilde {D}$ и $\tilde {k}.$ Кажущийся коэффициент диффузии $\tilde {D}$ = = ${{(\gamma {{D}_{{\text{A}}}} + {{D}_{{\text{B}}}})} \mathord{\left/ {\vphantom {{(\gamma {{D}_{{\text{A}}}} + {{D}_{{\text{B}}}})} {(\gamma + 1)}}} \right. \kern-0em} {(\gamma + 1)}}$ можно оценить по уравнениям (17) или (21) и затем рассчитать DA при известном значении DB или наоборот. Выражения для определения параметров нуклеации в случае образования зародышей сплава найдем из уравнений (15) и (17), (19) и (21) при мгновенной и прогрессирующей нуклеации соответственно:
(23)
${{N}_{{\text{0}}}} = 0.0130{{({{\tilde {z}e} \mathord{\left/ {\vphantom {{\tilde {z}e} {{{i}_{{\text{m}}}}{{t}_{{\text{m}}}}}}} \right. \kern-0em} {{{i}_{{\text{m}}}}{{t}_{{\text{m}}}}}})}^{2}}{{c_{0}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}} \mathord{\left/ {\vphantom {{c_{0}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}} {{{{\tilde {\upsilon }}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}}} \right. \kern-0em} {{{{\tilde {\upsilon }}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}},$
(24)
${{J}_{{\text{0}}}} = 0.0578{{({{\tilde {z}e} \mathord{\left/ {\vphantom {{\tilde {z}e} {{{i}_{{\text{m}}}}{{t}_{{\text{m}}}}}}} \right. \kern-0em} {{{i}_{{\text{m}}}}{{t}_{{\text{m}}}}}})}^{2}}{{c_{0}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}} \mathord{\left/ {\vphantom {{c_{0}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}} {({{{\tilde {\upsilon }}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}{{t}_{{\text{m}}}})}}} \right. \kern-0em} {({{{\tilde {\upsilon }}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}{{t}_{{\text{m}}}})}},$
(25)
${{N}_{{\text{s}}}} = 0.0335{{({{\tilde {z}e} \mathord{\left/ {\vphantom {{\tilde {z}e} {{{i}_{{\text{m}}}}{{t}_{{\text{m}}}}}}} \right. \kern-0em} {{{i}_{{\text{m}}}}{{t}_{{\text{m}}}}}})}^{2}}{{c_{0}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}} \mathord{\left/ {\vphantom {{c_{0}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}} {{{{\tilde {\upsilon }}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}}} \right. \kern-0em} {{{{\tilde {\upsilon }}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}},$
где Ns – насыщенная плотность зародышей на электроде при прогрессирующей нуклеации, ${{N}_{{\text{s}}}} = {{({{{{J}_{0}}} \mathord{\left/ {\vphantom {{{{J}_{0}}} {2k{\kern 1pt} 'D}}} \right. \kern-0em} {2k{\kern 1pt} 'D}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$ [17, 27]. Значения $\tilde {D},$ N0, J0 могут быть также оценены из анализа начальных и конечных участков транзиентов по формулам (13), (14) и уравнению $i = \tilde {z}e{{c}_{0}}{{({{\tilde {D}} \mathord{\left/ {\vphantom {{\tilde {D}} {\pi t}}} \right. \kern-0em} {\pi t}})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}},$ соответственно.

ВЫВОДЫ

Исследованы начальные стадии электрохимического образования бинарного сплава на индифферентном электроде в потенциостатических условиях. В рамках подхода, описывающего общий диффузионный поток двух сортов осаждаемых ионов в объеме электролита к поверхности полусферических зародышей сплава, найдено выражение для плотности тока роста. Рассмотрены случаи мгновенной и прогрессирующей нуклеации. При учете перекрытия диффузионных зон зародышей использовано приближение модели Scharifker–Hills. Получены уравнения для расчета кажущегося коэффициента диффузии осаждаемых ионов $\tilde {D},$ плотности и насыщенной плотности зародышей на поверхности электрода N0 и Ns, стационарной скорости нуклеации J0.

Исследование выполнено при поддержке гранта Российского научного фонда (проект № 16-13-00061).

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

  1. Гамбург Ю.Д. Электрохимическая кристаллизации металлов и сплавов. М.: Янус-К, 1997.

  2. Gamburg Yu.D., Zangari G. Theory and Practice of Metal Electrodeposition. N.Y.: Springer, 2011.

  3. Zangari G. Electrodeposition of Alloys and Compounds in the Era of Microelectronics and Energy Conversion Technology // Coatings. 2015. 5. P. 195–218.

  4. Landolt D. Fundamental aspects of alloy plating // Plat. Surf. Finish. 2001. 88. P. 70–79.

  5. Milchev A., Michailova E., Zapryanova T. Initial stages of electrochemical alloy formation: size and composition of critical nuclei // Electrochem. Comm. 2004. 6. P. 713–718.

  6. Milchev A., Lacmann R. On the nucleation theory of electrochemical alloy formation I. Overvoltage dependence of the stationary nucleation rate // J. Cryst. Growth. 1991. 110. P. 919–924.

  7. Milchev A., Lacmann R. On the nucleation theory of electrochemical alloy formation II. Concentration dependence of the stationary nucleation rate // J. Cryst. Growth. 1991. 110. P. 925–929.

  8. Milchev A., Lacmann R. Electrochemical growth of single metal and alloy clusters – Part II. Potentiostatic conditions. Theory // Electrochim. Acta. 1995. 40. P. 1475–1478.

  9. Milchev A. Electrochemical alloy formation – theory of progressive and instantaneous nucleation without overlap // Electrochim. Acta. 1997. 42. P. 1533–1536.

  10. Díaz-Morales O., Mostany J., Borrás C., Scharifker B.R. Current transient study of the kinetics of nucleation and diffusion-controlled growth of bimetallic phases // J. Solid State Electrochem. 2013. 17. P. 345–351.

  11. Scharifker B.R., Mostany J. Three-dimensional nucleation with diffusion controlled growth: Part I. Number density of active sites and nucleation rates per site // J. Electroanal. Chem. 1984. 177. P. 13–23.

  12. Bard A., Faulkner L. Electrochemical methods: fundamentals and applications. New York: Wiley, 2001.

  13. Галюс З. Теоретические основы электрохимического анализа. М.: Мир, 1974.

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

  15. Fletcher S. Electrochemical deposition of hemispherical nuclei under diffusion control. Some theoretical considerations // J. Chem. Soc. Faraday Trans. I. 1983. 79. P. 467–479.

  16. Gunawardena G., Hills G.J., Montenegro I., Scharifker B. Electrochemical nucleation // J. Electroanal. Chem. 1982. 138. P. 225–239.

  17. Scharifker B.R., Hills G.J. Theoretical and experimental studies of multiple nucleation // Electrochim. Acta. 1983. 28. P. 879–889.

  18. Sluyters-Rehbach M., Wijenberg J.H.O.J., Bosco E., Sluyters J.H. The theory of chronoamperometry for the investigation of electrocrystallization. Mathematical description and analysis in the case of diffusion-controlled growth // J. Electroanal. Chem. 1987. 236. P. l–20.

  19. Mirkin M.V., Nilov A.P. Three-dimensional nucleation and growth under controlled potential // J. Electroanal. Chem. 1990. 283. P. 35–51.

  20. Heerman L., Tarallo A. Theory of the chronoamperometric transient for electrochemical nucleation with diffusion-controlled growth // J. Electroanal. Chem. 1999. 470. P. 70–76.

  21. D’Ajello P.C.T., Munford M.L., Pasa A.A. Transient equations for multiple nucleation on solid electrodes: A stochastic description // J. Chem. Phys. 1999. 111. P. 4267–4273.

  22. Matthijs E., Langerock S., Michailova E., Heerman L. The potentiostatic transient for 3D nucleation with diffusion-controlled growth: theory and experiment for progressive nucleation // J. Electroanal. Chem. 2004. 570. P. 123–133.

  23. Исаев В.А., Гришенкова О.В., Зайков Ю.П. Анализ геометрико-вероятностных моделей электрокристаллизации // Расплавы. 2016. № 5. С. 355–370.

  24. Isaev V.A., Grishenkova O.V., Zaykov Yu.P. On the theory of 3D multiple nucleation with kinetic controlled growth // J. Electroanal. Chem. 2018. 818. P. 265–269.

  25. Колмогоров А.Н. К статистической теории кристаллизации металлов // Изв. АН СССР. Сер. мат. 1937. № 3. С. 355–359.

  26. Беленький В.З. Геометрико-вероятностные модели кристаллизации. Феноменологический подход. М: Наука, 1980.

  27. Isaev V.A., Zaykov Yu.P., Grishenkova O.V., Kosov A.V., Semerikova O.L. Analysis of Potentiostatic Current Transients for Multiple Nucleation with Diffusion and Kinetic Controlled Growth // J. Electrochem. Soc. 2019. 166. P. D851–D856.

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