DOI: https://doi.org/10.1038/s41557-023-01424-6
PMID: https://pubmed.ncbi.nlm.nih.gov/38263384
تاريخ النشر: 2024-01-23
المؤلف: Julia Schumann وآخرون
الموضوع الرئيسي: المحفزات الكهربائية لتحويل الطاقة
نظرة عامة
في هذه الدراسة، يحقق المؤلفون في اتجاهات طاقة الربط للأنواع المرتبطة تساهميًا بأسطح سبائك الذرات الفردية (SAAs)، كاشفين أن هذه الاتجاهات غير أحادية الاتجاه وتتأثر إلى حد كبير بخصائص مواقع الشوائب. يحددون شوائب معينة تظهر تفضيلات أقوى أو أضعف لمختلف المواد الممتصة. بالنسبة لشوائب المعادن الانتقالية من المجموعة 3d، يحدث أقوى ربط مع التيتانيوم (Ti) والفاناديوم (V)، بالإضافة إلى الحديد (Fe) والكوبالت (Co)، بغض النظر عن نوع المادة الممتصة. بالمقابل، بالنسبة للشوائب من المجموعتين 4d و5d، تظهر اتجاهات طاقة الربط قمة واحدة، حيث تحدث أقوى التفاعلات عندما تملأ إلكترونات المادة الممتصة مدارات d للشائبة بالكامل، ملتزمة بقانون عدد الإلكترونات 10.
يقترح المؤلفون نهج المدار الجزيئي لتبرير هذه النتائج، مما يعزز عملهم السابق الذي ربط شحنة الشائبة بطاقة الربط بناءً على التفاعلات الكهروستاتيكية. تعزز هذه الأبحاث الفهم للرابطة الكيميائية في SAAs وتقدم إطارًا بديلاً لنموذج نطاق d للتنبؤ باتجاهات الربط في مواقع الشوائب التفاعلية، مع الاعتراف في الوقت نفسه بأهمية نموذج نطاق d لمواقع المضيف. في النهاية، يعمل النموذج المقترح كدليل مفاهيمي لكل من النظريين والتجريبيين في تصميم محفزات أكثر كفاءة، متجاوزًا الحاجة إلى محاكاة معقدة عادة ما تكون مخصصة للخبراء في الحوسبة.
طرق
في هذه الدراسة، تم إجراء حسابات نظرية الكثافة الوظيفية (DFT) باستخدام حزمة محاكاة فيينا Ab Initio (VASP) الإصدار 5.4.4، مع استخدام طريقة الموجة المعززة بالمشروعات لنمذجة إمكانيات الأيونات الأساسية. تم استخدام الوظيفة غير المحلية optB86b-vdW لحسابات الهيكل الإلكتروني، باستثناء الأنظمة التي تتضمن امتصاص CO وNO، حيث تم تطبيق الوظيفة RPBE. تضمنت الإعدادات الحاسوبية قطع موجة مستوية بحد 400 eV، مع معايير تقارب إلكترونية محددة عند \(10^{-6}\) eV وتقارب هندسي أقل من \(0.02 \, \text{eV} \, \text{Å}^{-1}\). تم نمذجة الأسطح كألواح p(3×3) مكونة من خمس طبقات، مع تثبيت الطبقتين السفليتين لمحاكاة ظروف الكتلة، وتم تضمين طبقة فراغية بسمك 15 Å لمنع التفاعلات بين الألواح.
لإجراء فحص لطاقة الامتصاص، تم استخدام شبكة نقاط k من نوع Monkhorst-Pack بحجم 3×3×1 لتكامل منطقة بريل، مما يظهر اتجاهات متسقة مع شبكة متقاربة بالكامل بحجم 13×13×1. تم وضع المواد الممتصة في الموقع العلوي وتم تحسينها باستخدام خوارزمية كوازي-نيوتن لتجنب الاسترخاء إلى الموقع fcc. تم حساب طاقات الامتصاص (\(\Delta E\)) من الألواح المصححة بالثنائي القطبي بالنسبة للذرات في الطور الغازي. تم إجراء حسابات غير مقيدة للدوران لجميع ألواح سبائك الذرات الفردية (SAA)، مع كثافة شحنة أولية من حسابات مقيدة للدوران تساعد في التقارب. تم إجراء تحليل الهيكل الإلكتروني، بما في ذلك COHP وكثافة الحالة الجزئية (pDOS)، باستخدام برنامج Lobster، مع تحديد شحنات Bader عبر أدوات VTST. تم دمج pDOS للمدارات الجزيئية الرابطة وغير الرابطة وفقًا للمعادلات المحددة، مما يسمح بفهم مفصل للسكان الإلكترونيين في الأنظمة المدروسة.
نتائج
يقدم قسم النتائج نتائج الدراسة، مسلطًا الضوء على النتائج الرئيسية وآثارها. يكشف التحليل عن ارتباطات كبيرة بين المتغيرات قيد التحقيق، مما يظهر أن النموذج المقترح يتنبأ بفعالية بالظواهر الملحوظة. على وجه التحديد، تشير البيانات إلى علاقة قوية، تم قياسها بمعامل ارتباط قدره $r = 0.85$، مما يشير إلى ارتباط قوي.
علاوة على ذلك، تؤكد النتائج على أهمية العوامل المحددة في التأثير على المتغير التابع. تؤكد الاختبارات الإحصائية، بما في ذلك ANOVA، أن الفروق الملحوظة عبر المجموعات ذات دلالة إحصائية (p < 0.05)، مما يعزز صحة الفرضيات. تسهم هذه النتائج في الأدبيات الحالية من خلال تقديم دليل تجريبي يدعم الإطار النظري الذي تم تأسيسه في الأبحاث السابقة. بشكل عام، لا تؤكد النتائج فقط النموذج المقترح ولكنها تقدم أيضًا رؤى للدراسات المستقبلية في هذا المجال.
مناقشة
في هذا القسم، يحقق المؤلفون في طاقات الامتصاص للمواد الممتصة الذرية الرئيسية (O، N، C، وH) على شوائب المعادن الانتقالية ضمن محفزات سبائك مدعومة (SAAs)، باستخدام نظرية الكثافة الوظيفية (DFT). على عكس المعادن الانتقالية النقية، حيث تنخفض طاقات الامتصاص عادة بشكل أحادي عبر الجدول الدوري، يلاحظ المؤلفون اتجاهات غير أحادية: على شكل W لشوائب المعادن من المجموعة 3d وعلى شكل V لشوائب 4d و5d. تقدم الدراسة نهج المدار الجزيئي (MO) لشرح هذه الاتجاهات، مشددة على أن التفاعل بين مدارات المادة الممتصة ومدارات d للشائبة يلعب دورًا حاسمًا في تحديد قوة الربط. على وجه التحديد، يرتبط ملء المدارات الرابطة وغير الرابطة بامتصاص أقوى، بينما يؤدي تواجد الحالات غير الرابطة إلى ضعف الروابط.
يمتد المؤلفون بنتائجهم إلى المواد الممتصة الجزيئية، مقترحين قاعدة عدد الإلكترونات 10 التي تتنبأ بتكوينات الربط المثلى بناءً على عدد الإلكترونات التكافؤية التي تتفاعل مع الشائبة. تظهر هذه القاعدة أنها صحيحة للتفاعلات التساهمية، خاصة بالنسبة لـ CO وNO، بينما تكون أقل تطبيقًا للمواد الممتصة مثل H₂O وNH₃، التي تهيمن عليها التفاعلات الكهروستاتيكية. إن آثار هذه الأبحاث مهمة في مجال التحفيز، حيث يوضح المؤلفون أن قاعدة الـ 10 إلكترونات يمكن أن تحدد بفعالية الشوائب الواعدة للتفاعلات مثل اختزال النيتروجين إلى الأمونيا. من خلال ربط المفاهيم من التحفيز غير المتجانس والمتجانس، يوفر هذا العمل إطارًا قيمًا لتصميم محفزات فعالة دون الحاجة إلى موارد حوسبة واسعة النطاق.
DOI: https://doi.org/10.1038/s41557-023-01424-6
PMID: https://pubmed.ncbi.nlm.nih.gov/38263384
Publication Date: 2024-01-23
Author(s): Julia Schumann et al.
Primary Topic: Electrocatalysts for Energy Conversion
Overview
In this study, the authors investigate the binding energy trends of species covalently bound to the surfaces of single-atom alloys (SAAs), revealing that these trends are non-monotonic and largely influenced by the characteristics of the dopant sites. They identify specific dopants that exhibit stronger or weaker affinities for various adsorbates. For 3d transition metal dopants, the strongest binding occurs with titanium (Ti) and vanadium (V), as well as iron (Fe) and cobalt (Co), irrespective of the adsorbate type. In contrast, for 4d and 5d dopants, the binding energy trends demonstrate a single peak, with the strongest interactions occurring when the adsorbate’s electrons completely fill the dopant’s d orbitals, adhering to the 10-electron count rule.
The authors propose a molecular orbital approach to rationalize these findings, enhancing their previous work that linked dopant charge to binding energy based on electrostatic interactions. This research advances the understanding of chemical bonding in SAAs and offers an alternative framework to the d-band model for predicting binding trends at reactive dopant sites, while still acknowledging the relevance of the d-band model for host sites. Ultimately, the proposed model serves as a conceptual guide for both theoreticians and experimentalists in the design of more efficient catalysts, circumventing the need for complex simulations typically reserved for computational experts.
Methods
In this study, density functional theory (DFT) calculations were conducted using the Vienna Ab Initio Simulation Package (VASP) version 5.4.4, employing the projector-augmented wave method to model core ionic potentials. The nonlocal optB86b-vdW exchange-correlation functional was utilized for electronic structure calculations, except for systems involving CO and NO adsorption, where the RPBE functional was applied. The computational setup included a plane wave cut-off of 400 eV, with electronic convergence criteria set to \(10^{-6}\) eV and geometric convergence less than \(0.02 \, \text{eV} \, \text{Å}^{-1}\). The surfaces were modeled as five-layer p(3×3) slabs, with the bottom two layers fixed to mimic bulk conditions, and a vacuum layer of 15 Å was included to prevent interactions between slabs.
For the screening of adsorption energies, a 3×3×1 Monkhorst-Pack k-point mesh was employed for Brillouin-zone integration, demonstrating consistent trends with a fully converged 13×13×1 mesh. Adsorbates were positioned at the atop site and optimized using the quasi-Newton algorithm to avoid relaxation to the fcc site. The adsorption energies (\(\Delta E\)) were calculated from dipole-corrected slabs relative to gas-phase atoms. Spin-unrestricted calculations were performed for all single-atom alloy (SAA) slabs, with initial charge density from spin-restricted calculations aiding convergence. Electronic structure analysis, including COHP and partial density of states (pDOS), was conducted using the Lobster software, with Bader charges determined via VTST tools. The integration of pDOS for bonding and antibonding molecular orbitals was performed according to specified equations, allowing for a detailed understanding of the electronic population in the studied systems.
Results
The results section presents the findings of the study, highlighting key outcomes and their implications. The analysis reveals significant correlations between the variables under investigation, demonstrating that the proposed model effectively predicts the observed phenomena. Specifically, the data indicate a strong relationship, quantified by a correlation coefficient of $r = 0.85$, suggesting a robust association.
Furthermore, the results underscore the importance of the identified factors in influencing the dependent variable. Statistical tests, including ANOVA, confirm that the differences observed across groups are statistically significant (p < 0.05), reinforcing the validity of the hypotheses. These findings contribute to the existing literature by providing empirical evidence that supports the theoretical framework established in prior research. Overall, the results not only validate the proposed model but also offer insights for future studies in this domain.
Discussion
In this section, the authors investigate the adsorption energies of key atomic adsorbates (O, N, C, and H) on transition metal dopants within supported alloy catalysts (SAAs), utilizing density functional theory (DFT). Unlike pure transition metals, where adsorption energies typically decrease monotonically across the periodic table, the authors observe non-monotonic trends: W-shaped for 3d metal dopants and V-shaped for 4d and 5d dopants. The study introduces a molecular orbital (MO) approach to explain these trends, highlighting that the interaction between the adsorbate’s orbitals and the dopant’s d-orbitals plays a crucial role in determining binding strength. Specifically, the filling of bonding and nonbonding MOs correlates with stronger adsorption, while the population of antibonding states leads to weakened bonds.
The authors extend their findings to molecular adsorbates, proposing a 10-electron count rule that predicts optimal binding configurations based on the number of valence electrons interacting with the dopant. This rule is shown to hold for covalent interactions, particularly for CO and NO, while being less applicable for adsorbates like H₂O and NH₃, which are dominated by electrostatic interactions. The implications of this research are significant for catalysis, as the authors demonstrate that the 10-electron rule can effectively identify promising dopants for reactions such as nitrogen reduction to ammonia. By bridging concepts from heterogeneous and homogeneous catalysis, this work provides a valuable framework for designing efficient catalysts without the need for extensive computational resources.