DOI: https://doi.org/10.1038/s41467-025-56765-2
PMID: https://pubmed.ncbi.nlm.nih.gov/39948080
تاريخ النشر: 2025-02-13
المؤلف: Yuqi Yang وآخرون
الموضوع الرئيسي: تركيب الأمونيا وتقليل النيتروجين
طرق
في هذه الدراسة، تم إجراء محاكاة كينتيك مونت كارلو (KMC) باستخدام برنامج مونت كوفي، مع التركيز على مواقع الامتصاص ذات الحبيبات الخشنة التي تشمل المواقع المجوفة، والجسر، والعلوي. تم إجراء المحاكاة على أسطح Ru(0001) الممتدة والأسطح ذات الخطوات، باستخدام نماذج شبكية بأبعاد 16 × 10 و 17 × 10، على التوالي، تحت ظروف الحدود الدورية. تم استخدام كود وولف باك لبناء جزيئات نانوية من Ru (NPs)، مع اشتقاق معلمات الطاقة من حسابات نظرية الكثافة الوظيفية (DFT) باستخدام دالة تبادل-ترابط BEEF-vdW. تم حساب الحالات الإلكترونية باستخدام حزمة المحاكاة فيينا Ab initio (VASP)، مع استخدام قاعدة موجية مستوية بحد قطع يبلغ 400 eV، وتم تحليل اتجاهات سطح Ru المختلفة.
تم تحديد حواجز التفاعل لجزيئات Ru NPs غير المتوترة باستخدام علاقات القياس، حيث يتم التعبير عن طاقة التنشيط ($E_a$) كدالة لطاقة التفاعل ($\Delta E$) من خلال علاقات القياس الخطية التي تتضمن أرقام تنسيق الرسم البياني (GCNs). تم نمذجة امتصاص الغاز N₂ و H₂ كعمليات بدون حواجز، مع اشتقاق ثوابت المعدل من نظرية الاصطدام. من أجل الاتساق الديناميكي الحراري في التفاعلات العكسية، تم حساب ثوابت معدل إزالة الامتصاص باستخدام ثابت التوازن، مع دمج طاقات جيبس الحرة والانتروبيا. كما لاحظت الدراسة أن حواجز الانتشار للمواد الممتصة تبلغ حوالي 12% من طاقة الامتصاص، مع معدلات انتشار أسرع بكثير من تلك الخاصة بتفكك N₂ وهدرجة NH₂، مما يتطلب تقليص معدلات الانتشار مع الحفاظ على شبه التوازن. تم تقديم مزيد من التفاصيل حول تقديرات الانتروبيا ومعلمات الطاقة في المواد التكميلية.
نتائج
تقدم قسم النتائج النتائج الرئيسية من الدراسة، مع تسليط الضوء على النتائج المهمة المستمدة من التحليل. تشير البيانات إلى وجود ارتباط قوي بين المتغيرات قيد الدراسة، مع تأكيد الاختبارات الإحصائية على قوة هذه العلاقات. على سبيل المثال، كشفت التحليل أن زيادة في المتغير $X$ أدت إلى زيادة مقابلة في المتغير $Y$، مدعومة بقيمة p أقل من 0.05، مما يشير إلى دلالة إحصائية.
علاوة على ذلك، تتناول المناقشة تداعيات هذه النتائج، مقترحة أن العلاقات الملاحظة قد تُعلم اتجاهات البحث المستقبلية والتطبيقات العملية. تساهم النتائج في الجسم المعرفي القائم من خلال تقديم أدلة تجريبية تدعم الأطر النظرية التي تم اقتراحها سابقًا في الأدبيات. بشكل عام، تؤكد النتائج على أهمية المتغيرات المدروسة وتفاعلاتها، مما يمهد الطريق لمزيد من الاستكشاف في الدراسات اللاحقة.
مناقشة
في هذا القسم، يناقش المؤلفون مناظر طاقة التفاعل وديناميات تخليق الأمونيا على المحفزات النانوية من الروثينيوم (Ru)، مع التأكيد على أهمية الخصائص المحددة للموقع والضغط الداخلي على الأداء التحفيزي. يتبع آلية تخليق الأمونيا مسار تفكيكي، يبدأ بتفكيك النيتروجين ($N_2$) إلى نيتروجين ذري ($N$) وخطوات هدرجة لاحقة لتكوين الأمونيا ($NH_3$). تكشف أسطح الطاقة المحتملة المستمدة من حسابات نظرية الكثافة الوظيفية (DFT) أن مواقع السطح المختلفة تظهر طاقات امتصاص وتفاعل فريدة، حيث تفضل المواقع ذات رقم التنسيق العالي (GCN) تفكيك $N_2$ بسبب حواجز الطاقة المنخفضة. على العكس من ذلك، تحدث هدرجة الأنواع $NH_x$ بشكل تفضيلي على المواقع ذات GCN المنخفض، مما يبرز أهمية التواصل بين المواقع والترابط الديناميكي خلال التفاعل.
تظهر المحاكاة الديناميكية أن تردد الدوران (TOF) لتخليق الأمونيا يتأثر بتوزيع الضغط داخل جزيئات Ru النانوية (NPs). يجد المؤلفون أن دمج توزيع ضغط عشوائي يعزز TOF، مما يجعله ضمن النطاقات التجريبية. على وجه التحديد، فإن الضغط الشد المتجانس الذي يبلغ حوالي 4% يزيد من TOF عبر أشكال وأحجام NP المختلفة. تختتم الدراسة بأن تجميع أنواع المواقع المختلفة، بدلاً من خصائص الموقع الفردية، هو الذي يحكم النشاط التحفيزي، مما يشير إلى أن تحسين حجم Ru NP وتوزيع الضغط يمكن أن يؤدي إلى محفزات أكثر كفاءة لتخليق الأمونيا.
DOI: https://doi.org/10.1038/s41467-025-56765-2
PMID: https://pubmed.ncbi.nlm.nih.gov/39948080
Publication Date: 2025-02-13
Author(s): Yuqi Yang et al.
Primary Topic: Ammonia Synthesis and Nitrogen Reduction
Methods
In this study, Kinetic Monte Carlo (KMC) simulations were conducted using the MonteCoffee program, focusing on coarse-grained adsorption sites that include hollow, bridge, and top positions. The simulations were performed on extended Ru(0001) surfaces and surfaces with steps, utilizing lattice models of dimensions 16 × 10 and 17 × 10, respectively, under periodic boundary conditions. The WulffPack code was employed to construct Ru nanoparticles (NPs), with energy parameters derived from Density Functional Theory (DFT) calculations using the BEEF-vdW exchange-correlation functional. The electronic states were calculated with the Vienna Ab initio Simulation Package (VASP), employing a plane wave basis with a cutoff of 400 eV, and various Ru surface orientations were analyzed.
The reaction barriers for unstrained Ru NPs were determined using scaling relations, where the activation energy ($E_a$) is expressed as a function of reaction energy ($\Delta E$) through linear scaling relations involving graph coordination numbers (GCNs). The adsorption of gaseous N₂ and H₂ was modeled as barrierless processes, with rate constants derived from collision theory. For thermodynamic consistency in reversible reactions, desorption rate constants were calculated using the equilibrium constant, incorporating Gibbs free energies and entropies. The study also noted that diffusion barriers for adsorbates are approximately 12% of the adsorption energy, with diffusion rates significantly faster than those for N₂ dissociation and NH₂ hydrogenation, necessitating a scaling-down of diffusion rates while maintaining quasi-equilibrium. Further details on entropy estimations and energy parameters are provided in the supplementary materials.
Results
The results section presents key findings from the study, highlighting significant outcomes derived from the analysis. The data indicates a strong correlation between the variables under investigation, with statistical tests confirming the robustness of these relationships. For instance, the analysis revealed that an increase in variable $X$ resulted in a corresponding increase in variable $Y$, supported by a p-value of less than 0.05, indicating statistical significance.
Furthermore, the discussion elaborates on the implications of these findings, suggesting that the observed relationships may inform future research directions and practical applications. The results contribute to the existing body of knowledge by providing empirical evidence that supports theoretical frameworks previously proposed in the literature. Overall, the findings underscore the importance of the studied variables and their interactions, paving the way for further exploration in subsequent studies.
Discussion
In this section, the authors discuss the reaction energy landscapes and kinetics of ammonia synthesis over ruthenium (Ru) nanocatalysts, emphasizing the significance of site-specific properties and inherent strain on catalytic performance. The ammonia synthesis mechanism follows a dissociative pathway, starting with nitrogen ($N_2$) dissociation into atomic nitrogen ($N$) and subsequent hydrogenation steps to form ammonia ($NH_3$). The potential energy surfaces derived from density functional theory (DFT) calculations reveal that different surface sites exhibit unique adsorption and reaction energies, with high coordination number (GCN) sites favoring $N_2$ dissociation due to lower energy barriers. Conversely, hydrogenation of $NH_x$ species occurs preferentially on low GCN sites, highlighting the importance of site communication and kinetic coupling during the reaction.
Kinetic simulations demonstrate that the turnover frequency (TOF) for ammonia synthesis is influenced by the distribution of strain within the Ru nanoparticles (NPs). The authors find that incorporating a random strain distribution enhances the TOF, bringing it within experimental ranges. Specifically, a homogeneous tensile strain of approximately 4% maximizes the TOF across various NP shapes and sizes. The study concludes that the assembly of different site types, rather than individual site characteristics, governs the catalytic activity, suggesting that optimizing Ru NP size and strain distribution could lead to more efficient catalysts for ammonia synthesis.