Doimiy zaryad (constant-charge, CC) va doimiy potensial (constant-potential, CP) hisoblari elektrokimyoviy elektrod–elektrolit interfeysi energiyasini ikki xil statistik ansambl yondashuvi bilan modellashtiradigan hisoblash usullaridir. CC yondashuvida elektronlar soni doimiy saqlanadi va qoʻllangan elektrod potensialining taʼsiri keyin termodinamik tuzatishlar orqali qoʻshiladi; CP yondashuvida esa elektronlarning elektrokimyoviy potensiali, demak elektrod potensiali doimiy saqlanadi va tizimdagi elektronlar sonining oʻzgarishiga ruxsat beriladi. Ushbu tadqiqot Ni, Pd, Pt, Cu, Ag va Au metallarining (111) yuzalarida *H, *O va *OH adsorbsiyalarini taqqoslab, arzonroq CC yondashuvi qimmatroq CP hisobining energetik natijalarini qaysi sharoitlarda ishonchli qayta hosil qila olishini oʻrganadi.

Tadqiqotning markaziy natijasi shuki, CC va CP oʻrtasidagi farqni adsorbsiya nol zaryad potensialini (potential of zero charge, PZC) qanchalik oʻzgartirishi bilan katta darajada izohlash mumkin. Adsorbsiya tufayli yuzaga kelgan \(|\Delta U_{\mathrm{PZC}}|\) kichik boʻlsa, yashirin erituvchi tutgan doimiy zaryad yondashuvi (CCsolv) CP natijalariga yaqin qolishi mumkin. Manba oʻrganilgan tizimlarda \(|\Delta U_{\mathrm{PZC}}|\leq0{,}40\ \mathrm{V_{SHE}}\) yoki vakuumda hisoblangan sirtga-normal dipol zichligi boʻyicha \(|\Delta P_{\mathrm{norm,vac}}|\leq10{,}09\ \mathrm{me\,Å^{-1}}\) qiymatlarini CCsolvdan foydalanish uchun amaliy mezon sifatida taklif qiladi.

*H adsorbsiyasi 1,00 monolayergacha (ML) barcha oltita metall yuzasida CCvac, CCsolv va CP oʻrtasida juda oʻxshash energiyalarni bergan. *O va ayniqsa *OH qoplanishlarida esa PZC siljishlari kattalashishi, elektrod potensiali adsorbsiya joyining barqarorligini oʻzgartirishi va CC usullari CP natijalaridan farqlanishi mumkin. *OH uchun turli adsorbsiya joylari orasidagi PZC oʻzgarishi \(\Delta\Delta U_{\mathrm{PZC}}\geq0{,}50\ \mathrm{V_{SHE}}\) boʻlganda potensialga bogʻliq adsorbsiya joyi oʻzgarishi muhim boʻlib qoladi.

Tadqiqot shuningdek \(\Delta U_{\mathrm{PZC}}\) bilan adsorbsiya hosil qilgan sirtga-normal dipol zichligi oʻrtasida kuchli chiziqli bogʻlanish topgan. Yashirin erituvchi hisoblaridagi global regressiyaning qiyaligi taxminan 7,91 V Å e⁻¹, R² qiymati taxminan 0,97 va bu bogʻlanishdan olingan samarali dielektrik doimiy 20,69. Vakuumda hisoblangan dipol momenti ham PZC oʻzgarishini yetarlicha yaxshi bashorat qilgani sababli, qimmat CP hisobi haqiqatan zarurmi yoki yoʻqmi, arzonroq hisobdan oldindan baholash mumkin boʻladi.

Bu natijalar barcha elektrokimyoviy tizimlar uchun isbotlangan universal chegaralar sifatida talqin qilinmasligi kerak. Tadqiqot oltita oʻtish metallining (111) yuzalari, *H/*O/*OH adsorbatlari, ishlatilgan DFT funksionali, yashirin erituvchi modeli va koʻrsatilgan hisoblash sharoitlari bilan cheklangan. Mualliflar mezonlarning murakkabroq interfeyslar va reaksiya tarmoqlari uchun umumlashtirilishi kelajakda alohida sinovdan oʻtkazilishi kerakligini taʼkidlaydilar.

Elektrokimyoviy sirt qoplanishi nega muhim?

Elektrokatalizator ish vaqtida “yalangʻoch metall yuza” boʻlib qoladi deb faraz qilish koʻp hollarda yetarli emas. Elektrolitdan keladigan vodorod, gidroksil, kislorod yoki boshqa turlar yuzadagi faol joylarni qoplashi mumkin. Ushbu adsorbsiyalangan turlar reaksiya oraliq mahsulotlari boʻlishi, faol joylarni toʻsishi yoki yangi faol sirt fazalarini hosil qilishi mumkin.

Shu sababli sirt qoplanishi muayyan adsorbat yuzadagi mos adsorbsiya joylarining qanchasini egallashini ifodalaydi. Tadqiqotda qoplanish miqdori monolayer (ML) birliklarida koʻrib chiqiladi. Masalan, 0,25 ML belgilangan sirt joylarining taxminan toʻrtdan biriga mos keladigan model qoplanishni bildiradi.

Elektrodning eng barqaror sirt qoplanishi potensial, pH va elektrolit muhiti bilan oʻzgarishi mumkinligi sababli, katalitik faollikni faqat yalangʻoch yuzaning energiyasiga qarab baholash xato boʻlishi mumkin. Tadqiqotning asosiy maqsadi elektrokatalizatorning haqiqiy ish sharoitlaridagi “dam olish holati” yoki termodinamik jihatdan afzal qoplanish fazasini ishonchliroq hisoblashdir.

Elektrokimyoviy ikki qavat nima?

Elektrod elektrolit bilan kontaktda boʻlganda, yuzadagi elektr zaryadi bilan eritmadagi ionlar bir-birini elektrostatik tarzda tartibga soladi va elektr ikki qavat (electric double layer, EDL) hosil boʻladi. Manba konseptual ravishda Gouy–Chapman–Stern modelidan foydalanadi.

Elektrod yuzasiga yaqin zich ion sohasi Stern qatlamini, undan uzoqdagi tarqoqroq ion taqsimoti esa diffuz qatlamni hosil qiladi. Tashqi Helmholtz tekisligi (outer Helmholtz plane, OHP) ushbu ikki sohani tavsiflashda ishlatiladigan referenslardan biridir. Potensialning eritmada qanchalik tez soʻnishini Debye ekranlanish uzunligi \(L_D\) tavsiflaydi. Kichikroq \(L_D\) elektr maydoni qisqaroq masofada kuchliroq ekranlanishini anglatadi.

*H, *OH va *O belgilari nimani anglatadi?

Sirt ilmida yulduzcha belgisi \(*\) elektrod yuzasidagi boʻsh yoki adsorbsiya qilinishi mumkin boʻlgan joyni ifodalaydi. Shunga koʻra:

\[ \mathrm{H^+_{(aq)} + e^- + * \rightleftharpoons H*} \]

ushbu reaksiya proton va elektron birikib, yuzada adsorbsiyalangan vodorod hosil boʻlishini ifodalaydi.

\[ \mathrm{H_2O_{(l)} + * \rightleftharpoons *OH + H^+_{(aq)} + e^-} \]

suv molekulasidan yuzaga bogʻlangan gidroksil hosil boʻlishini koʻrsatadi.

\[ \mathrm{H_2O_{(l)} + * \rightleftharpoons *O + 2H^+_{(aq)} + 2e^-} \]

esa yuzaga bogʻlangan kislorod turining hosil boʻlishini ifodalaydi.

Bu tenglamalar faqat kimyoviy stexiometriya emas; CHE yondashuvida adsorbsiya erkin energiyalarini elektrod potensialiga bogʻlashning asosini tashkil qiladi.

Doimiy Zaryad va Doimiy Potensial Hisoblari Oʻrtasidagi Asosiy Farq Nima?

Doimiy zaryad yondashuvida hisoblash yacheykasidagi elektronlar soni reaksiya koordinatasi boʻylab doimiy saqlanadi va elektrod potensialining taʼsiri keyin termodinamik tuzatish orqali qoʻshiladi; doimiy potensial yondashuvida esa elektrodning elektron elektrokimyoviy potensiali doimiy saqlanadi, tizim zarurat boʻlsa elektron olishi yoki berishi mumkin va interfeys energiyasi qoʻllangan potensial bilan oʻz-oʻzidan mos ravishda oʻzgaradi.

Kanonik ansambl va Computational Hydrogen Electrode

Tadqiqot doimiy zaryad hisoblarini kanonik ansambl yondashuvi doirasida koʻrib chiqadi. Bu doirada umumiy elektronlar soni doimiy. Elektrokimyoviy potensial taʼsirini keyin qoʻshishda eng keng tarqalgan usullardan biri Computational Hydrogen Electrode (CHE) yondashuvidir.

CHEning asosiy referensi:

\[ \mathrm{H^+_{(aq)}+e^- \rightleftharpoons \frac{1}{2}H_2(g)} \]

tenglamasidir.

0 V RHE va birlik proton faolligida proton–elektron juftining kimyoviy potensiali:

\[ \mu_{\mathrm{H^++e^-}}=\frac{1}{2}\mu_{\mathrm{H_2}} \]

deb olinadi. RHE shkalasida istalgan \(U_{\mathrm{RHE}}\) potensialida:

\[ \mu_{\mathrm{H^++e^-}}= \frac{1}{2}\mu_{\mathrm{H_2}}-eU_{\mathrm{RHE}} \]

boʻladi. Bu yerdagi \(e\) elementar elektr zaryadidir. Potensial oshganda proton–elektron juftlarining kimyoviy potensiali ushbu chiziqli had orqali oʻzgaradi.

RHE, SHE va pH munosabati

Manba Nernst munosabatini quyidagi shaklda beradi:

\[ E_{\mathrm{RHE}} = E_{\mathrm{SHE}}^{0} -\frac{RT}{F} \ln\left(\frac{1}{a_{\mathrm{H^+}}}\right) \]

\[ \mathrm{pH}=-\log(a_{\mathrm{H^+}}) \]

298 Kda bu munosabat taxminan:

\[ E_{\mathrm{RHE}} = E_{\mathrm{SHE}}^{0} -0.059\,\mathrm{pH} \]

va potensial konversiyasi:

\[ U_{\mathrm{RHE}} = U_{\mathrm{SHE}}+0.059\,\mathrm{pH} \]

shaklida yoziladi.

\(R\) universal gaz doimiysi, \(T\) mutlaq harorat, \(F\) Faraday doimiysi va \(a_{\mathrm{H^+}}\) proton faolligidir. Tadqiqotda bu munosabat pH va elektrod potensiali sirt qoplanishi barqarorligiga birgalikda qanday taʼsir qilishini tushuntirish uchun ishlatiladi.

CHE ostidagi adsorbsiya Gibbs energiyalari

n ta adsorbsiyalangan vodorod uchun manba:

\[ \Delta G_{*nH,\mathrm{RHE}} = G_{*nH} - \left( G_*+\frac{n}{2}G_{\mathrm{H_2}} \right) + neU_{\mathrm{RHE}} \]

ifodasidan foydalanadi.

n ta gidroksil uchun:

\[ \Delta G_{*nOH,\mathrm{RHE}} = G_{*nOH} + \frac{n}{2}G_{\mathrm{H_2}} - \left( G_*+nG_{\mathrm{H_2O}} \right) - neU_{\mathrm{RHE}} \]

va sirt kislorodi uchun:

\[ \Delta G_{*nO,\mathrm{RHE}} = G_{*nO} + nG_{\mathrm{H_2}} - \left( G_*+nG_{\mathrm{H_2O}} \right) - 2neU_{\mathrm{RHE}} \]

ishlatiladi.

Bu yerda \(G_{*nX}\) — n ta X turi adsorbsiyalangan yuzaning Gibbs energiyasi; \(G_*\) — toza yuzaning Gibbs energiyasi; \(G_{\mathrm{H_2}}\) va \(G_{\mathrm{H_2O}}\) esa referens vodorod va suv turlarining Gibbs energiyalari.

Bu tenglamalar *H qaytaruvchi sharoitlarda, *OH va *O esa oksidlovchiroq potensiallarda yoki yuqoriroq pH sharoitlarida barqarorlashishga moyilligini hisoblash imkonini beradi.

CCvac va CCsolv oʻrtasidagi farq

CCvac doimiy zaryad hisobini vakuum sirt modeli bilan bajaradi. Erituvchi qutblanishi va elektrolit ekranlanishi hisobga olinmaydi.

CCsolv esa doimiy zaryad doirasini saqlab qoladi, lekin elektrolitni qutblanishi mumkin boʻlgan uzluksiz muhit sifatida ifodalovchi yashirin erituvchi modelini qoʻshadi. Ushbu usul haqiqiy Stern qatlamidagi alohida suv molekulalarini, maxsus ion adsorbsiyasini yoki yoʻnalgan vodorod bogʻlarini hal qilmaydi; biroq oʻrtacha elektrostatik erituvchi taʼsirlarini CCvacga qaraganda realroq ifodalaydi.

Tadqiqot natijalaridan biri shuki, CP hisobini taxminiy baholash kerak boʻlsa, boshlangʻich nuqta CCvac emas, kamida CCsolv boʻlishi kerak.

Grand-kanonik doimiy potensial yondashuvi

Grand-kanonik ansamblda elektrodning elektron elektrokimyoviy potensiali doimiy saqlanadi, elektronlar soni oʻzgarishi mumkin. Shu tariqa yuzadagi adsorbsiya, zaryad oʻtkazilishi yoki dipol hosil boʻlishi elektrod potensialini sunʼiy oʻzgartirish oʻrniga elektron rezervuari bilan zaryad almashinuvi orqali muvozanatlanadi.

Manbada yashirin erituvchining dielektrik funksiyasi:

\[ \varepsilon[\rho] = 1+ (\varepsilon_{\mathrm{bulk}}-1)S[\rho] \]

shaklida aniqlanadi. Bu yerda \(\rho\) — Kohn–Sham elektron zichligi; \(\varepsilon_{\mathrm{bulk}}\) — hajmiy erituvchining dielektrik doimiysi; \(S[\rho]\) esa elektron zichligiga bogʻliq erituvchi boʻshligʻining shakl funksiyasi.

Shakl funksiyasi manbada komplementar xato funksiyasi yordamida aniqlanadi:

\[ S[\rho] = \frac{1}{2} \operatorname{erfc} \left( \frac{\log(\rho/\rho_{\mathrm{cut}})} {\sigma\sqrt{2}} \right) \]

\(\rho_{\mathrm{cut}}\) erituvchi boʻshligʻi hosil boʻla boshlaydigan elektron zichligi chegarasini, \(\sigma\) oʻtish sohasining kengligini boshqaradi.

Erituvchi boʻshligʻini hosil qilish energiyasining sirt hissasi:

\[ G_{\mathrm{cavity}} = \tau\int d^3r\,|\nabla S[\rho]| \]

shaklida yoziladi. \(\tau\) samarali sirt tarangligi parametridir.

Poisson–Boltzmann ekranlanishi

Elektrolitdagi harakatchan ionlarning elektr maydonini ekranlash taʼsiri tadqiqotda chiziqlashtirilgan Poisson–Boltzmann yondashuvi bilan ifodalanadi:

\[ \nabla\cdot \left( \varepsilon_0\varepsilon[\rho]\nabla\phi \right) - \frac{S[\rho]\varepsilon_0\varepsilon_{\mathrm{bulk}}} {L_D^2}\phi = -\rho_s \]

Bu yerda:

• \(\varepsilon_0\): vakuumning dielektrik singdiruvchanligi,
• \(\varepsilon[\rho]\): joylashuvga bogʻliq dielektrik funksiya,
• \(\phi\): elektrostatik potensial,
• \(L_D\): Debye ekranlanish uzunligi,
• \(\rho_s\): erigan tizimning umumiy zaryad zichligi.

Manba \(L_D=3\) Å ishlatgan. Bu qiymat 1 M 1:1 elektrolitni va taxminan kuchli kislotali pH 1 yoki kuchli ishqoriy pH 13 sharoitlarini ifodalash uchun tanlangan.

Mutlaq elektrod potensiali va SHEga aylantirish

Hajmiy elektrolit potensiali \(U_{\mathrm{bulk}}\) nol referens sifatida olinib, mutlaq elektrod potensiali:

\[ \Phi_{\mathrm{abs}} = eU_{\mathrm{bulk}}-E_{\mathrm{Fermi}} \]

shaklida hisoblanadi.

SHE shkalasiga aylantirish:

\[ U_{\mathrm{SHE}} = \frac{\Phi_{\mathrm{abs}}-\Phi_{\mathrm{SHE}}}{e} \]

orqali bajariladi. Tadqiqotda mutlaq SHE potensiali uchun \(\Phi_{\mathrm{SHE}}=4{,}44\) eV ishlatilgan; manba adabiyotda taxminan 4,40–4,80 eV oraligʻidagi qiymatlar ham xabar qilinganini taʼkidlaydi.

PZC nima?

Nol zaryad potensiali (PZC) — elektrod yuzasidagi sof ortiqcha yoki yetishmaydigan elektron zaryadi nol boʻladigan elektrod potensiali. Manbada \(q=0\) shartidagi potensial \(U_{\mathrm{PZC}}\) deb aniqlanadi.

Adsorbat yuzaga kelganda elektron zichligini va sirt dipolini oʻzgartirishi mumkin. Natijada toza yuzaning PZCsi bilan adsorbat qoplagan yuzaning PZCsi farq qilishi mumkin. Tadqiqotning asosiy deskriptori:

\[ \Delta U_{\mathrm{PZC}} = U_{\mathrm{PZC},*nX} - U_{\mathrm{PZC},*} \]

mantiqida aniqlanadigan adsorbsiya tufayli PZC oʻzgarishidir.

Potensial ostida Gibbs energiyasi nega parabolik?

Tadqiqotning CP yondashuvida interfeys sigʻimi oʻrganilgan potensial oraligʻida doimiy deb faraz qilinadi:

\[ \frac{\partial^2G} {\partial\Phi_{\mathrm{abs}}^2} = -2AC \]

Bu yerda \(A\) sirt maydoni va \(C\) interfeys sigʻimidir. Ushbu faraz ostida energiya potensialga nisbatan ikkinchi darajali oʻzgaradi:

\[ G(U)=aU^2+bU+c \]

va PZC markazli koʻrinish:

\[ G(U) = -\frac{1}{2}C(U-U_{\mathrm{PZC}})^2 + G_{\mathrm{PZC}} \]

shaklida.

Ushbu parabolaning choʻqqi nuqtasi \(U_{\mathrm{PZC}}\) va \(G_{\mathrm{PZC}}\) qiymatlarini, egriligi esa sigʻimni beradi. CC yondashuvida ishlatiladigan chiziqli potensial tuzatishidan farqli ravishda CP yondashuvidagi ushbu ikkinchi darajali xatti-harakat ikki usul ayniqsa katta PZC oʻzgarishlarida nega farqlanishini tushuntiradi.

Grand-kanonik energiya

DFT-LPCM energiyasida elektron rezervuari bilan zaryad almashinuvi hisobga olinganda grand-kanonik energiya:

\[ \Omega = E_{\mathrm{DFT-LPCM}} - qE_{\mathrm{Fermi}} \]

shaklida hisoblanadi.

Termokimyoviy tuzatishlar qoʻshilganda:

\[ G = \Omega + \mathrm{ZPE} + U_{\mathrm{therm}} - TS \]

hosil boʻladi.

Bu yerda ZPE nol-nuqta energiyasi, \(U_{\mathrm{therm}}\) issiqlik energiyasi hissasi, \(T\) harorat va \(S\) entropiya.

Potensialga bogʻliq adsorbsiya Gibbs energiyasi

Toza va adsorbatli yuzalarning parabolik energiya egri chiziqlari taqqoslanganda manba Tenglama (21)ni beradi:

\[ \Delta G(U) = -\frac{1}{2}C_{*nX} \left( U-U_{\mathrm{PZC},*nX} \right)^2 + G_{\mathrm{PZC},*nX} + \frac{1}{2}C_* \left( U-U_{\mathrm{PZC},*} \right)^2 - G_{\mathrm{PZC},*} \]

Bu tenglama tadqiqotning asosiy matematik bogʻlanishlaridan biridir. Adsorbatli va toza yuzalarning sigʻimlari va PZClari bir-biriga juda yaqin boʻlsa, ikki parabolaning ikkinchi darajali hissalari katta darajada bir-birini bekor qiladi va natija CCga oʻxshash chiziqli xatti-harakatga yaqinlashadi. Aksincha, adsorbsiya PZCni kuchli siljitsa, CP energiyasi bilan CC yondashuvi oʻrtasidagi farq ortadi.

PZC Siljishi Nega CC bilan CP Oʻrtasidagi Moslikni Belgilaydi?

Adsorbsiya PZCni siljitganda toza va adsorbatli yuzalar bir xil qoʻllangan potensialda turli sirt zaryadlariga va turli elektr maydon–dipol oʻzaro taʼsirlariga ega boʻladi; CC hisobi bu potensialga bogʻliq sigʻim javobini bevosita ifodalay olmagani sababli \(|\Delta U_{\mathrm{PZC}}|\) ortgan sayin CC va CP adsorbsiya energiyalari osonroq farqlanadi.

*OH nega ayniqsa sezgir?

Manbada 0,25 ML *H va *O uchun afzal adsorbsiya joylari oʻrganilgan potensial oraligʻida katta darajada potensialdan mustaqil. *OH uchun vaziyat boshqacha. Ayniqsa Pd, Pt va Au ustida CP hisoblari potensial oshgani sari top joydan fcc kabi hollow joylarga oʻtishni bashorat qiladi.

Bu xatti-harakatning fizik izohi adsorbat dipoli bilan interfeys elektr maydoni oʻrtasidagi oʻzaro taʼsirdir. Birinchi tartibli Stark yondashuvi bilan ushbu oʻzaro taʼsir taxminan:

\[ E_{\mathrm{dipol-alan}} \approx -\boldsymbol{\mu}\cdot\mathbf{F} \]

shaklida qaralishi mumkin.

Top joydagi *OH tashqariga koʻproq yoʻnalgan sirtga-normal dipolga ega boʻlishi mumkin. Potensial PZCning bir tomonida boʻlganda elektr maydoni ushbu joyni barqarorlashtirishi, PZCdan oʻtilganda esa maydon yoʻnalishi oʻzgarib, boshqa dipol yoʻnalishiga ega hollow joy barqarorroq boʻlishi mumkin.

Adsorbsiya joyi PZC farqi

Bir xil adsorbatning ikki xil adsorbsiya joyi orasidagi PZC oʻzgarishi manbada:

\[ \Delta\Delta U_{\mathrm{PZC}} = \Delta U_{\mathrm{PZC}}^{(i)} - \Delta U_{\mathrm{PZC}}^{(j)} \]

shaklida aniqlanadi.

\(\Delta\Delta U_{\mathrm{PZC}}\) ortgani sari ikki joyning elektr maydoni bilan oʻzaro taʼsiri koʻproq farqlanadi va potensialga bogʻliq joy almashishi ehtimoli ortadi. Tadqiqot *OH uchun \(\Delta\Delta U_{\mathrm{PZC}}\geq0{,}50\ \mathrm{V_{SHE}}\) qiymatini CP hisobiga oʻtish uchun amaliy ogohlantirish chegarasi sifatida taklif qiladi.

Sirtga-normal dipol qanday hisoblandi?

Adsorbsiya hosil qilgan elektron zichligi oʻzgarishi:

\[ \Delta\rho = \rho_{*X} - \left( \rho_*+\rho_X \right) \]

shaklida hisoblangan.

\(\rho_{*X}\) relaksatsiya qilingan sirt–adsorbat tizimining, \(\rho_*\) ayni geometrik holatdagi toza slabning va \(\rho_X\) izolyatsiyalangan adsorbatning elektron zichligidir.

Uch oʻlchamli zichlik farqi x va y yoʻnalishlarida integrallanib, z yoʻnalishida bir oʻlchamli profil olinadi:

\[ \Delta\lambda_z = \int_0^{a_1} \int_0^{a_2} \Delta\rho\,dy\,dx \]

Soʻng sirt maydoniga normallashtirilgan dipol zichligi:

\[ \Delta P = \frac{1}{A} \int_{z_0}^{z} (z'-z_0)\Delta\lambda_{z'}\,dz' = \frac{\mu_z}{A} \]

shaklida hisoblanadi.

\(\mu_z\) sirtga-normal dipol momentidir. Manbadagi ishora taʼrifiga koʻra musbat qiymat tashqariga, elektrolit tomon; manfiy qiymat esa slab ichiga tomon yoʻnalgan dipolni koʻrsatadi.

Sirt Dipoli Qimmat CP Hisobi Kerakligini Qanday Bashorat Qiladi?

Helmholtz munosabatiga koʻra adsorbsiya hosil qilgan sirtga-normal dipol zichligi PZC siljishi bilan taxminan chiziqli bogʻlangan; shu sababli arzonroq vakuum hisobidan olingan dipol toʻliq doimiy-potensial hisobiga oʻtishdan oldin \(\Delta U_{\mathrm{PZC}}\) kattaligi va demak CC xatosi xavfini tekshirish uchun ishlatilishi mumkin.

Helmholtz munosabati

Manba bu munosabatni:

\[ \Delta U_{\mathrm{PZC}} = \frac{e}{\varepsilon_{\mathrm{eff}}\varepsilon_0} \Delta P_{\mathrm{norm}} \]

shaklida beradi.

Bu yerda:

• \(\Delta U_{\mathrm{PZC}}\): adsorbsiya tufayli PZC oʻzgarishi,
• \(\Delta P_{\mathrm{norm}}\): sirtga-normal induksiyalangan dipol zichligi,
• \(\varepsilon_0\): vakuumning oʻtkazuvchanligi,
• \(\varepsilon_{\mathrm{eff}}\): interfeysning samarali dielektrik doimiysi,
• \(e\): elementar elektr zaryadidir.

Yashirin erituvchi ichidagi barcha tizimlarni qamrab olgan regressiyada manba taxminan:

• qiyalik: 7,91 V Å e⁻¹,
• R² ≈ 0,97,
• MAE ≈ 0,05 V,
• RMSE ≈ 0,07 V

deb xabar qiladi.

Ushbu qiyalikdan \(\varepsilon_{\mathrm{eff}}=20{,}69\) olingan. Alohida metall regressiyalarida samarali dielektrik doimiy taxminan 17,81–26,57 oraligʻida. Bu qiymat metall yuzasining \(\varepsilon\approx1\) muhiti bilan suvning \(\varepsilon=78{,}40\) qiymati oʻrtasida joylashgan samarali interfeys muhitiga mos keladi.

Vakuum dipol hisoblari ham kuchli bogʻlanish bergan. Manba 5B-rasmda global regressiya taxminan R² ≈ 0,91 darajasida ekanini bildiradi. Bu natija muhim, chunki vakuum hisobi yashirin erituvchi yoki CP hisobidan arzonroq.

Bader zaryadi nega har bir metallda bir xil darajada yaxshi deskriptor emas?

Adsorbsiya vaqtida adsorbat bilan metall oʻrtasidagi sof elektron oʻtkazilishi sirt dipolining manbalaridan biridir. Shu sababli adsorbatning Bader zaryadi, ayniqsa sof zaryad oʻtkazilishi ustun boʻlgan tizimlarda PZC oʻzgarishi bilan korrelyatsiya koʻrsatishi mumkin.

Manba Group 11 metallari Cu, Ag va Au uchun Bader zaryadi bilan \(\Delta U_{\mathrm{PZC}}\) oʻrtasida taxminan chiziqli bogʻlanish topgan. Bunga qarama-qarshi ravishda Group 10 metallari Ni, Pd va Pt uchun bogʻlanish kuchsizroq.

Sababi shundaki, Group 10 metallarida nafaqat sof zaryad oʻtkazilishi, balki metall slabining oʻz ichidagi elektron zichligi qayta taqsimlanishi ham muhim. Adsorbat elektron buluti bilan metall sirt elektronlari ustma-ust tushganda Pauli itarilishi elektron zichligini slab ichiga qaytaradi. Ushbu taʼsir adabiyotda pillow effect yoki cushion effect deb ataladi.

Overlap index qanday aniqlandi?

Tadqiqot slab ichidagi zaryad qayta taqsimlanishini adsorbatga sof zaryad oʻtkazilishidan ajratish uchun yangi taqqoslama koʻrsatkichni aniqlaydi.

Slab ichidagi qayta taqsimlangan zaryad:

\[ Q_{\mathrm{slab}} = \int_{z_0}^{z_{\mathrm{int}}} |\Delta\lambda_z|\,dz - \left| \int_{z_0}^{z_{\mathrm{int}}} \Delta\lambda_z\,dz \right| \]

shaklida.

Adsorbat sohasiga oʻtkazilgan sof zaryad:

\[ Q_{\mathrm{ads}} = \int_{z_{\mathrm{int}}}^{z_{\mathrm{plat}}} \Delta\lambda_z\,dz \]

deb hisoblanadi.

Overlap index:

\[ I_O = \frac{Q_{\mathrm{slab}}} {|Q_{\mathrm{ads}}|} \]

shaklida aniqlanadi.

\(I_O\rightarrow0\) boʻlganda adsorbsiya koʻproq sof zaryad oʻtkazilishi bilan tavsiflanadi. \(I_O\) oshgani sari slab ichidagi orbital ustma-ust tushish/Pauli qayta taqsimlanishining nisbiy ahamiyati ortadi.

*O va *OH uchun \(I_O\leq0{,}55\) boʻlgan zaryad-oʻtkazilishi ustun rejim faqat Group 11 metallari orasida kuzatilgan. *H uchun esa Ag va Audagi juda kichik \(|Q_{\mathrm{ads}}|<0{,}1e\) qiymatlari nisbatni son jihatdan sezgir qilgani sababli, ayni chegara umumlashtirilmagan.

Mualliflar \(I_O=0{,}55\) qiymatini barcha tizimlarga qoʻllanadigan universal chegara sifatida taklif qilmaydilar. Overlap index ushbu tadqiqotda taqqoslama va fizik jihatdan talqin qilinadigan deskriptordir; uning kengroq umumlashtirilishi kelajakda oʻrganilishi kerakligi koʻrsatiladi.

Tadqiqot Usuli va Natijalari

Hisoblash dasturi va elektron tuzilish usuli

Spin-qutblangan periodik Kohn–Sham DFT hisoblari VASP 6.4.2 bilan bajarilgan.

Oltita metall va panjara doimiylari

Boshlangʻich hajmiy kristall tuzilmalari Materials Projectdan olingan va optimallashtirilgan.

Sirt slab modellari

Har bir metallning optimallashtirilgan hajmiy kristali (111) tekisligidan kesilib, p(2×2) sirt yacheykasi hosil qilingan. Slab toʻrtta atom qatlamidan iborat. Sirtga-normal z yoʻnalishda periodik tasvirlar oʻrtasidagi oʻzaro taʼsirni kamaytirish uchun taxminan 15 Å vakuum qoʻshilgan.

Geometriyani optimallashtirishda eng pastdagi ikki metall qatlami hajmiy holatda qotirilgan, yuqori qatlamlar va adsorbatlarning relaksatsiyasiga ruxsat berilgan.

Yashirin erituvchi parametrlari

Erituvchi parametrlari uchun sezgirlik nazorati

Bare va OH-qoplangan Ag(111) namunaviy tizimida \(\varepsilon_{\mathrm{bulk}}\), \(\sigma\), \(\rho_{\mathrm{cut}}\) va \(\tau\) qiymatlari standart qiymatlarining ±%50 miqdorida oʻzgartirilgan. Natijada \(\Delta U_{\mathrm{PZC}}\) eng koʻpi bilan ±0,03 V oʻzgargan.

Manba bu natijani ishlatilgan erituvchi parametrlarining asosiy PZC farqi tendensiyalariga taʼsiri kichik va tizimlar oʻrtasida katta darajada sistematik ekanini koʻrsatuvchi sezgirlik nazorati sifatida talqin qiladi. \(L_D\) oʻzgarishining ham \(\Delta U_{\mathrm{PZC}}\)ga taʼsiri eʼtiborga olinmaydigan darajada ekani; kattaroq Debye uzunliklari esa interfeys sigʻimini kamaytirishi qayd etiladi.

Tebranish va termokimyo hisobi

Optimallashtirilgan geometriyalar minimal energiya konfiguratsiyasi ekanini tekshirish uchun chekli farq usuli bilan tebranish chastotalari hisoblangan. Simmetriya cheklovi qoʻllanmagan va siljish amplitudasi 0,02 Å tanlangan.

Izolyatsiyalangan molekulalarda barcha atomlar, slab tizimlarida esa faqat adsorbat atomlari tebratilgan. Taxminan 100 cm⁻¹dan past sunʼiy mavhum chastotalar manbada qoʻllangan yondashuvga koʻra 12 cm⁻¹ haqiqiy chastotalar bilan almashtirilgan.

Izolyatsiyalangan molekulalar ideal gaz yondashuvi bilan, adsorbatlar esa garmonik yondashuv ostida koʻrib chiqilgan.

Charge-interpolation usuli

Barcha tuzilmalar avval vakuumda optimallashtirilgan, keyin PZC sharoitidagi energiya va elektrod potensialini olish uchun yashirin erituvchi muhitida qayta optimallashtirilgan.

Potensialga bogʻliqlikni namunalash uchun neytral slabning elektronlar soni:

• −6×10⁻³ elektron bilan +6×10⁻³ elektron oraligʻida,
• 2×10⁻³ elektron qadamlar bilan

oʻzgartirilgan. Shu tariqa jami yettita zaryad holati hisoblanib, Gibbs energiyasi–potensial munosabati parabolik tarzda moslashtirilgan.

Mualliflar har bir maqsad potensialini alohida yaqinlashtiradigan toʻgʻridan-toʻgʻri hisoblash potentiostati oʻrniga charge-interpolation usulini xarajat, barqarorlik va potensialga bogʻliq geometriya oʻzgarishlarini ushlash oʻrtasidagi amaliy murosa sifatida tanlaganlar.

Adsorbsiya joylari qanday skanerlandi?

*H, *O va *OH adsorbatlari uchun top, bridge, fcc va hcp joylar oʻrganilgan. Boshlangʻich konfiguratsiyalarning bir qismi optimallashtirish vaqtida bir xil yakuniy fcc geometriyasiga yaqinlashgan yoki yuqori qoplanishlarda molekulyar turlar hosil boʻlishini koʻrsatgan; natijalarga faqat barqaror geometriyalar kiritilgan.

0,25 ML darajasida eng barqaror joylar aniqlangandan keyin *H va *O uchun ayni turdagi joylar yuqoriroq qoplanishlarga toʻldirilgan. 1,00 ML darajasiga yetilganda keyingi qatlamlar uchun joy skanerlashi qayta bajarilgan. *OH potensialga kuchliroq sezgirlik koʻrsatgani sababli, yuqori qoplanishlarda bir nechta adsorbsiya joyi alohida oʻrganilgan.

Potensialga bogʻliq *OH joy almashishlari

Manbada potensial oynasi ichida:

• 0,50 ML *OH / Pt(111),
• 0,50 ML *OH / Au(111),
• 1,00 ML *OH / Pd(111),
• 1,00 ML *OH / Au(111)

tizimlarida potensialga bogʻliq adsorbsiya joyi oʻzgarishi kuzatilgan.

Bu natija past qoplanishda katta \(\Delta\Delta U_{\mathrm{PZC}}\) va joy barqarorligi oʻzgarishi koʻrinsa, yuqori qoplanishlarda turli adsorbsiya joylarini CP usuli bilan ochiq baholash kerakligini qoʻllab-quvvatlaydi.

CCvac, CCsolv va CP qaysi adsorbatda qanchalik mos keldi?

CCsolv Qaysi Chegaralarda CP Oʻrniga Ishlatilishi Mumkin?

Ushbu tadqiqotda CCsolv uchun tavsiya qilingan asosiy amaliy mezon \(|\Delta U_{\mathrm{PZC}}|\leq0{,}40\ \mathrm{V_{SHE}}\) yoki unga mos \(|\Delta P_{\mathrm{norm,vac}}|\leq10{,}09\ \mathrm{me\,Å^{-1}}\) shartidir; qiymat kattaroq boʻlsa CP afzal koʻriladi, biroq ish potensiali toza va adsorbatli yuzalar PZClarining energiya-kesishish oʻrta nuqtasiga 0,40 V_SHE ichida boʻlsa, CCsolv ayrim hollarda baribir yetarli boʻlishi mumkin.

Energiyalar kesishadigan xarakterli potensial:

\[ U_{\Delta G=\Delta G(U)} = \frac{ U_{\mathrm{PZC},*nX} + U_{\mathrm{PZC},*} }{2} \]

shaklida toza va adsorbat qoplagan yuzalar PZClarining oʻrta nuqtasi bilan bogʻlangan.

Eksperimental faza oʻtishlari bilan taqqoslash

Tadqiqot oʻz elektrokimyoviy tajribasini oʻtkazmagan; CP va CCsolv yordamida hisoblangan faza-oʻtish potensiallarini avval nashr qilingan eksperimental oʻlchovlar bilan taqqoslagan.

Pt misolida CCsolvning muvaffaqiyati tadqiqotning PZC mezoniga mos: adsorbsiya tufayli \(\Delta U_{\mathrm{PZC}}\) atigi 0,04 V_SHE. Cu misolida qiymat 0,30 V_SHE boʻlib, yana 0,40 V mezonidan past va CCsolv tajribaga yaqin natija beradi.

Ag misolini ehtiyotkorroq talqin qilish kerak. CPning 0,33 V bashorati 0,51 V eksperimental qiymatdan taxminan 0,18 V farq qiladi; manba buni “reasonable agreement” deb baholaydi. Shu sababli Verianla matnida “mukammal moslik” yoki “toʻliq tasdiq” sifatida kengaytirilmaydi.

Pourbaix diagrammalarining roli

Mualliflar pHni 0–14 oraligʻida oʻzgartirib, har bir pH va potensial nuqtasidagi eng barqaror *H, *O va *OH sirt qoplanishini aniqlagan va sirt-qoplanish Pourbaix diagrammalarini tuzgan.

CCvac, CCsolv va CP diagrammalari taqqoslanganda, ayniqsa Ag kabi ayrim metallarda ayni pH/potensial sohasiga turli adsorbat fazalari biriktirilishi koʻrinadi. Bu farq faqat sonli energiya xatosi emas; katalizator qaysi sirt holatida ekani haqidagi fizik talqinni oʻzgartirishi mumkin.

Manba shuningdek Pourbaix diagrammasi oʻzi barcha energiya yaqinliklarini koʻrsatmasligini taʼkidlaydi. Bir nechta qoplanish fazasi bir-biriga yaqin energiyada boʻlsa, faqat rangli barqarorlik sohasiga qarash yetarli emas; bevosita adsorbsiya Gibbs energiyasi egri chiziqlari bilan birga baholash kerak.

Amaliy qaror qoidalari

Tadqiqot qoʻllab-quvvatlaydigan natijalar

• Oltita metallning (111) yuzalarida *H adsorbsiyasi 1,00 MLgacha CCvac, CCsolv va CP oʻrtasida juda yaqin energetik xatti-harakat koʻrsatgan.
• *OH adsorbsiya joyi Pd, Pt va Au kabi ayrim yuzalarda elektrod potensialiga sezilarli darajada bogʻliq boʻlishi mumkin.
• Adsorbsiya tufayli katta PZC oʻzgarishlari CC–CP energiya farqlarining muhim fizik sababidir.
• \(\Delta U_{\mathrm{PZC}}\) sirtga-normal adsorbsiya dipoli bilan kuchli chiziqli bogʻlanish koʻrsatgan.
• Vakuum dipoli qimmatroq yashirin erituvchi yoki CP hisobi kerakmi-yoʻqmi oldindan tekshirish uchun ishlatilishi mumkin boʻlgan arzonroq deskriptordir.
• Sof Bader zaryadi zaryad-oʻtkazilishi ustun Group 11 metallarida PZC siljishi bilan yaxshiroq bogʻlanish koʻrsatadi.
• Group 10 metallari uchun slab ichidagi Pauli/pillow zaryad qayta taqsimlanishi muhim boʻlgani sababli faqat Bader zaryadi dipolning hammasini tushuntirmaydi.
• CCvac erituvchi va doimiy-potensial taʼsirlarini birga hisobga olmagani sababli ayrim sirt-faza oʻtishlarini sifat jihatdan notoʻgʻri bashorat qilishi mumkin.
• CCsolv \(\Delta U_{\mathrm{PZC}}\) kichik boʻlganda CPga sezilarli darajada yaxshiroq yaqinlashadi.

Tadqiqot qoʻllab-quvvatlamaydigan natijalar

• “Barcha elektrokimyoviy hisoblarda CP majburiy” degan xulosa qoʻllab-quvvatlanmaydi.
• “CC har qanday sharoitda CP kabi aniq” degan xulosa qoʻllab-quvvatlanmaydi.
• 0,40 V va 0,50 V mezonlari barcha metallar, adsorbatlar, elektrolitlar va yuzalar uchun universal fizik konstantalar sifatida tasdiqlanmagan.
• Tadqiqot haqiqiy elektrolitdagi atomistik suv tuzilishini yoki maxsus ion adsorbsiyasini ochiq erituvchi molekulalari bilan modellashtirmaydi.
• CCsolv haqiqiy Stern qatlamining molekulyar tuzilishini toʻliq ifodalamaydi.
• Hisoblangan sirt qoplanishlarini barcha haqiqiy elektrokatalizator morfologiyalariga avtomatik umumlashtirib boʻlmaydi; ular bevosita faqat oʻrganilgan ideal (111) slab modellariga tegishli.
• Mualliflar yangi elektrokimyoviy tajriba oʻtkazmagan; eksperimental taqqoslashlar ilgari xabar qilingan adabiyot qiymatlaridir.
• Overlap index uchun universal chegara tasdiqlanmagan.

Tadqiqotning asosiy cheklovlari

• Oltita metallning faqat (111) yuzalari oʻrganilgan.
• Adsorbat toʻplami *H, *O va *OH bilan cheklangan.
• Yashirin erituvchi yondashuvi birinchi suv qatlamining yoʻnalishini, maxsus ion adsorbsiyasini va yoʻnalgan vodorod bogʻlarini atomistik ravishda hal qilmaydi.
• Interfeys sigʻimi oʻrganilgan potensial oraligʻida potensialdan mustaqil deb faraz qilingan.
• Asosiy hisoblarda Debye uzunligi uchun 3 Å tanlangan va oraliq pHlarda ion kuchini doimiy saqlaydigan tayanch elektrolit faraz qilingan.
• BEEF-vdW, GLSSA13 va ishlatilgan slab geometriyasining natijalarga xos model farazlari mavjud.
• Amaliy 0,40 V va 0,50 V chegaralarining murakkabroq elektrokimyoviy interfeyslar va reaksiya tarmoqlarida umumlashtirilishi hali sinovdan oʻtkazilmagan.
• Overlap indexning DFT funksionalidan mustaqil umumiyroq deskriptor ekanligi kelajakdagi tadqiqot mavzusi sifatida qoldirilgan.

Ilmiy jihatdan eng muhim amaliy xulosa

Tadqiqot tavsiya qilgan ish oqimi “har doim eng qimmat CP hisobini bajar” yondashuvi oʻrniga bosqichli qaror jarayonini qoʻyadi. Avval arzonroq hisoblar bilan adsorbsiya sirt dipoli va PZCni qanchalik oʻzgartirishi baholanadi. Oʻzgarish kichik boʻlsa CCsolv ishlatilishi mumkin; oʻzgarish katta boʻlsa yoki adsorbsiya joyi potensialga bogʻliq boʻlib qolsa CPga oʻtiladi.

Bu yondashuv ayniqsa keng material skanlarida muhim. Chunki koʻp sonli metall, qotishma yoki adsorbat nomzodlarining barchasi uchun CP DFT bajarish katta hisoblash xarajatini tugʻdirishi mumkin. Tadqiqotning maqsadi aniqlikdan voz kechish emas, qimmat usulni faqat fizik jihatdan zarur boʻlgan holatlarda ishlatish uchun oldindan tanlash mezonini yaratishdir.

Manba va Usul Haqida Izoh

Original sarlavha: When do constant-charge and constant-potential computational simulations agree? A dipole-based framework for predicting electrochemical surface coverages

Mualliflar: Conor Brennan-Pollak; Dídac A. Fenoll; Max García-Melchor.

Masʼul muallif: Max García-Melchor.

Teng hissa: Manbada teng birinchi mualliflik/teng hissa belgisi mavjud emas.

Muassasalar: School of Chemistry, Trinity College Dublin, Ireland; Center for Cooperative Research on Alternative Energy (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Spain; IKERBASQUE, Basque Foundation for Science, Bilbao, Spain.

Jurnal: Electrochimica Acta.

Nashriyot: Elsevier Ltd.

Jild / maqola raqami: 566 / 148795.

DOI: 10.1016/j.electacta.2026.148795.

Qabul qilib olingan sana: 16 yanvar 2026.

Qayta koʻrib chiqish sanasi: 26 mart 2026.

Qabul qilingan sana: 29 mart 2026.

Onlayn nashr sanasi: PDFga koʻra 2 aprel 2026.

Jild/son nashr sanasi: ScienceDirect rasmiy qaydiga koʻra 1 avgust 2026.

Manba turi: Taqrizdan oʻtgan hisoblash elektrokatalizi / DFT tadqiqot maqolasi.

Taqriz holati: Electrochimica Actada nashr qilingan version-of-record tadqiqot maqolasi.

Litsenziya: Creative Commons Attribution–NonCommercial 4.0 International (CC BY-NC 4.0).

Rasmiy havola: https://doi.org/10.1016/j.electacta.2026.148795

Maʼlumotlar toʻplami: Tadqiqotning optimallashtirilgan tuzilmalar va energiyalarni oʻz ichiga olgan qoʻllab-quvvatlovchi hisoblash maʼlumotlari ioChem-BDda DOI 10.19061/iochem-bd-6-614 bilan ochiq koʻrsatilgan.

Moliyalashtirish: Taighde Éireann – Research Ireland 21/FFP-A/9161 granti Conor Brennan-Pollak va Max García-Melchorga; 20/FFP-P/8740 granti Max García-Melchorga yordam bergan. European Commission Horizon Europe doirasidagi 101191948 raqamli grant Dídac A. Fenoll va Max García-Melchorga, 101126600 raqamli grant Max García-Melchorga yordam bergan.

Hisoblash resurslari: Mualliflar Trinity College Dublin Research IT Unitga hisoblash resurslari uchun minnatdorchilik bildiradilar.

Manfaatlar toʻqnashuvi / competing interest: Max García-Melchor Horizon Europe va Research Ireland moliyaviy yordamlarini potensial competing-interest aloqasi sifatida bildirgan. Boshqa mualliflar tadqiqotga taʼsir qilishi mumkin boʻlgan maʼlum competing financial interest yoki personal relationship mavjud emasligini bayon qiladilar.

Muallif hissalari: Conor Brennan-Pollak: Writing – review & editing, Writing – original draft, Visualization, Methodology, Investigation, Formal analysis. Dídac A. Fenoll: Writing – review & editing, Visualization, Investigation, Formal analysis. Max García-Melchor: Writing – review & editing, Visualization, Validation, Supervision, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Conceptualization.

Asosiy hisoblash usuli: Spin-qutblangan periodik Kohn–Sham DFT; doimiy zaryad CHE yondashuvlari (CCvac va CCsolv) bilan grand-kanonik doimiy-potensial DFTni taqqoslash.

Model tizimlari: Ni(111), Pd(111), Pt(111), Cu(111), Ag(111) va Au(111) yuzalarida *H, *O va *OH adsorbsiyalari.

Asosiy dasturlar: VASP 6.4.2; ASE; VASPsol/GLSSA13; VASPkit. Hajmiy boshlangʻich tuzilmalar Materials Projectdan olingan.

Asosiy metodologik farazlar: Toʻrt qatlamli p(2×2) ideal slab; pastki ikki qatlamni qotirish; yashirin elektrolit; doimiy sigʻim yondashuvi; 3 Å Debye uzunligi; standart termokimyoviy yondashuvlar.

Eksperimental tasdiqlash holati: Tadqiqot yangi tajriba oʻtkazmaydi. Pt, Ag va Cu sirt faza-oʻtish potensiallari ilgari nashr qilingan elektrokimyoviy tajribalar bilan taqqoslangan.

Asosiy talqin chegarasi: \(|\Delta U_{\mathrm{PZC}}|\leq0,40\) V_SHE, \(\Delta\Delta U_{\mathrm{PZC}}\geq0,50\) V_SHE, 10,09 va 12,81 me Å⁻¹ dipol mezonlari ushbu tadqiqotdagi tizimlardan olingan amaliy tanlash mezonlaridir; murakkabroq interfeyslarda universal amal qilishi hali koʻrsatilmagan.