Расплавы, 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
Аннотация
Предложен подход к математическому описанию электрохимического образования и роста зародышей бинарного сплава на поверхности индифферентного электрода в условиях постоянного пересыщения (перенапряжения), рассматривающий общий поток двух сортов осаждаемых ионов в объеме электролита к поверхности полусферических зародышей новой фазы. Усреднение профиля концентраций проводили с учетом мольно-долевого соотношения компонентов. Исследованы случаи мгновенной и прогрессирующей нуклеации и диффузионно-контролируемого роста. Перекрытие диффузионных зон зародышей проанализировано в рамках Scharifker–Hills (SH) модели. Получены уравнения для расчета кажущегося коэффициента диффузии осаждаемых ионов, количества зародышей на единице поверхности электрода, стационарной скорости нуклеации.
ВВЕДЕНИЕ
Сплавы обладают более широким спектром свойств, чем индивидуальные металлы. Постоянно растущая потребность в новых конструкционных и функциональных материалах с повышенными эксплуатационными характеристиками для нужд современной техники, электроники, энергетики требует лучшего понимания основных закономерностей электрохимического образования сплавов для эффективного управления процессом на наноуровне. Кроме того, проблема совместного выделения металлических примесей при электроосаждении индивидуального металла также отчасти связана со сплавообразованием. Это обуславливает особый интерес к теории электрокристаллизации двух или нескольких металлов [1–4].
Кинетика начальных стадий образования бинарного сплава на индифферентном электроде при электроосаждении в потенциостатических условиях была проанализирована в рамках различных подходов и приближений. В работе [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]$(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].$(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],$Теперь допустим, что в электролите присутствуют не один, а два сорта осаждаемых ионов ${{{\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],$(7)
${{r}_{0}} = {{(2\tilde {D}{{c}_{0}}\tilde {\upsilon }\tau )}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}},$При множественной нуклеации вокруг каждого зародыша формируется полусферическая диффузионная зона, размер которой увеличивается по мере роста зародыша. Это приводит к постепенному перекрытию диффузионных полей и изменению условий массопереноса. Задача о перекрытии 3D зон может быть решена только в рамках некоторых приближений [11, 16–24]. Наиболее простой вариант решения был предложен в модели 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}}}}}),$(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),$(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),$(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).$Транзиенты плотности тока 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}}}.$(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],$(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}}.$(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],$(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}}.$(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}}}}}},$ВЫВОДЫ
Исследованы начальные стадии электрохимического образования бинарного сплава на индифферентном электроде в потенциостатических условиях. В рамках подхода, описывающего общий диффузионный поток двух сортов осаждаемых ионов в объеме электролита к поверхности полусферических зародышей сплава, найдено выражение для плотности тока роста. Рассмотрены случаи мгновенной и прогрессирующей нуклеации. При учете перекрытия диффузионных зон зародышей использовано приближение модели Scharifker–Hills. Получены уравнения для расчета кажущегося коэффициента диффузии осаждаемых ионов $\tilde {D},$ плотности и насыщенной плотности зародышей на поверхности электрода N0 и Ns, стационарной скорости нуклеации J0.
Исследование выполнено при поддержке гранта Российского научного фонда (проект № 16-13-00061).
Список литературы
Гамбург Ю.Д. Электрохимическая кристаллизации металлов и сплавов. М.: Янус-К, 1997.
Gamburg Yu.D., Zangari G. Theory and Practice of Metal Electrodeposition. N.Y.: Springer, 2011.
Zangari G. Electrodeposition of Alloys and Compounds in the Era of Microelectronics and Energy Conversion Technology // Coatings. 2015. 5. P. 195–218.
Landolt D. Fundamental aspects of alloy plating // Plat. Surf. Finish. 2001. 88. P. 70–79.
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.
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.
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.
Milchev A., Lacmann R. Electrochemical growth of single metal and alloy clusters – Part II. Potentiostatic conditions. Theory // Electrochim. Acta. 1995. 40. P. 1475–1478.
Milchev A. Electrochemical alloy formation – theory of progressive and instantaneous nucleation without overlap // Electrochim. Acta. 1997. 42. P. 1533–1536.
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.
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.
Bard A., Faulkner L. Electrochemical methods: fundamentals and applications. New York: Wiley, 2001.
Галюс З. Теоретические основы электрохимического анализа. М.: Мир, 1974.
Исаев В.А. Электрохимическое фазообразование. Екатеринбург: УрО РАН, 2007.
Fletcher S. Electrochemical deposition of hemispherical nuclei under diffusion control. Some theoretical considerations // J. Chem. Soc. Faraday Trans. I. 1983. 79. P. 467–479.
Gunawardena G., Hills G.J., Montenegro I., Scharifker B. Electrochemical nucleation // J. Electroanal. Chem. 1982. 138. P. 225–239.
Scharifker B.R., Hills G.J. Theoretical and experimental studies of multiple nucleation // Electrochim. Acta. 1983. 28. P. 879–889.
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.
Mirkin M.V., Nilov A.P. Three-dimensional nucleation and growth under controlled potential // J. Electroanal. Chem. 1990. 283. P. 35–51.
Heerman L., Tarallo A. Theory of the chronoamperometric transient for electrochemical nucleation with diffusion-controlled growth // J. Electroanal. Chem. 1999. 470. P. 70–76.
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.
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.
Исаев В.А., Гришенкова О.В., Зайков Ю.П. Анализ геометрико-вероятностных моделей электрокристаллизации // Расплавы. 2016. № 5. С. 355–370.
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.
Колмогоров А.Н. К статистической теории кристаллизации металлов // Изв. АН СССР. Сер. мат. 1937. № 3. С. 355–359.
Беленький В.З. Геометрико-вероятностные модели кристаллизации. Феноменологический подход. М: Наука, 1980.
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.
Дополнительные материалы отсутствуют.