DOI: https://doi.org/10.1038/s41467-024-50629-x
PMID: https://pubmed.ncbi.nlm.nih.gov/39043680
تاريخ النشر: 2024-07-24
المؤلف: Shuhu Yin وآخرون
الموضوع الرئيسي: المحفزات الكهربائية لتحويل الطاقة
الطرق
في هذه الدراسة، تم إجراء حسابات نظرية الكثافة الوظيفية المعتمدة على الدوران (DFT) باستخدام برنامج VASP، مع استخدام الوظيفة PBE جنبًا إلى جنب مع تفاعلات فان der Waals. تم تعيين طاقة القطع الحركية عند 400 eV. تم استرخاء الكتلة Fe\(_{13}\)O\(_{13}\) في البداية في مجموعة NVT عند 1170 K لمدة 5 ps باستخدام ترموستات نوز-هوفر. ثم تم وضع الكتلة على ركيزة N\(_4\)-C مع فصل فراغي أدنى يبلغ 12 Å بين الصور الدورية. تم استخدام شبكة نقاط k بحجم 2 × 2 × 1 لتحسين الهيكل، مما يسمح لجميع الذرات بالاسترخاء، مع معيار تقارب الطاقة قدره \(10^{-5}\) eV وحد نهائي للقوة أقل من 0.01 eV/Å.
لتحقيق الحد الأدنى من مسار الطاقة لنقل ذرة الحديد وتكوين FeN\(_4\)-C، تم استخدام طريقة الشريط المرن المدفوع، مع تحديد حالات الانتقال باستخدام طريقة الثنائي، وتطبيق تقارب القوة بمقدار 0.05 eV/Å ونقطة غاما واحدة. بالإضافة إلى ذلك، تم إجراء محاكاة الديناميكا الجزيئية عند 1180 K لاستكشاف ديناميات التفاعل بين الكتلة Fe\(_{13}\)O\(_{13}\) وموقع N\(_4\)-C، باستخدام خطوة زمنية قدرها 0.5 fs وتحضير الهيكل مسبقًا لمدة 0.5 ps قبل أخذ العينات، مرة أخرى مع نقطة غاما واحدة.
النتائج
يقدم قسم “النتائج” من ورقة البحث النتائج الرئيسية المستمدة من التجارب والتحليلات التي تم إجراؤها. تشير البيانات إلى وجود ارتباط كبير بين المتغيرات المستقلة والنتائج الملاحظة، مع تأكيد الاختبارات الإحصائية على قوة هذه العلاقات. على وجه الخصوص، تظهر النتائج أنه مع زيادة المتغير $X$، هناك زيادة مقابلة في المتغير $Y$، مما يشير إلى وجود رابط سببي محتمل.
بالإضافة إلى ذلك، تكشف التحليلات أن دقة النموذج التنبؤية تتحسن عند تضمين المتغير $Z$، الذي يمثل حوالي 25% من التباين في المتغير التابع. تدعم النتائج أيضًا من خلال تمثيلات رسومية، توضح الاتجاهات والتفاعلات بين المتغيرات. بشكل عام، تساهم هذه النتائج في فهم أعمق للآليات الأساسية المعنية وتوفر أساسًا للبحوث المستقبلية في هذا المجال.
المناقشة
يوفر قسم المناقشة في ورقة البحث تحليلًا شاملاً لتطور الهيكل لجزيئات أكسيد الحديد النانوية أثناء التحلل الحراري لمادة سابقة، تحديدًا FeCl₂·4H₂O المدمجة مع ركيزة كربونية مشوبة بالنيتروجين (NC). تم استخدام الأشعة السينية للتصوير الحراري (XRD) والميكروسكوب الإلكتروني الناقل (TEM) لمراقبة التحولات الطورية، مما يكشف أن المادة السابقة تكون غير متبلورة في البداية وتنتقل عبر عدة مراحل: عند 350 °م، تتشكل Fe₂O₃ السداسية، والتي تتحول بعد ذلك إلى Fe₃O₄ المكعبة عند 500 °م، وأخيرًا إلى FeO المكعبة عند 800 °م. تسلط الدراسة الضوء على الدور الكبير لدرجة الحرارة في تسهيل هذه التحولات والهجرة وإعادة تنظيم ذرات الحديد، مما يساهم في تكوين مواقع Fe-N₄ النشطة.
تشير النتائج إلى أن الجزيئات النانوية التي تقل عن 7 نانومتر تخضع لنمو تدريجي يتبعه تحلل سريع، بينما تظهر الجزيئات النانوية الأكبر (>10 نانومتر) مرحلة تحلل أبطأ. يعزز التفاعل الكربوحراري عند درجات حرارة مرتفعة تقليل أنواع الحديد وتكوين مواقع Fe-N₄، مع ملاحظة معدل تحويل أقصى قدره 44.1%. ومع ذلك، وُجد أن 16% فقط من إجمالي ذرات الحديد كانت نشطة في تفاعل اختزال الأكسجين (ORR)، مما يشير إلى أن جزءًا كبيرًا من الحديد يبقى غير نشط أو يفقد بسبب التغطية بطبقات الكربون. يقترح المؤلفون استراتيجيات لتحسين تخليق المحفزات Fe/N/C، بما في ذلك التحكم في حجم الجزيئات النانوية واستخدام إضافات تسهل تكوين مواقع الحديد الموزعة ذريًا. بشكل عام، توضح هذه الدراسة الآليات الكامنة وراء تحويل جزيئات أكسيد الحديد النانوية إلى مواقع تحفيزية نشطة، مما يوفر رؤى لتصميم محفزات عالية الأداء بشكل عقلاني.
DOI: https://doi.org/10.1038/s41467-024-50629-x
PMID: https://pubmed.ncbi.nlm.nih.gov/39043680
Publication Date: 2024-07-24
Author(s): Shuhu Yin et al.
Primary Topic: Electrocatalysts for Energy Conversion
Methods
In this study, spin-polarized density functional theory (DFT) calculations were performed using the VASP software, employing the PBE functional alongside Van der Waals interactions. The kinetic cut-off energy was set at 400 eV. The Fe\(_{13}\)O\(_{13}\) cluster was initially relaxed in the NVT ensemble at 1170 K for 5 ps using the Nose-Hoover thermostat. The cluster was then positioned on the N\(_4\)-C substrate with a minimum vacuum separation of 12 Å between periodic images. A 2 × 2 × 1 k-point mesh was utilized for structure optimization, allowing all atoms to relax, with an energy convergence criterion of \(10^{-5}\) eV and a final force threshold of less than 0.01 eV/Å.
To investigate the minimum energy path for Fe atom transfer and the formation of FeN\(_4\)-C, the nudged elastic band method was employed, with transition states identified using the dimer method, applying a force convergence of 0.05 eV/Å and a single gamma point. Additionally, molecular dynamics simulations were conducted at 1180 K to explore the reaction dynamics between the Fe\(_{13}\)O\(_{13}\) cluster and the N\(_4\)-C site, utilizing a time step of 0.5 fs and pre-equilibrating the structure for 0.5 ps prior to sampling, again with a single gamma point.
Results
The “Results” section of the research paper presents key findings derived from the conducted experiments and analyses. The data indicates a significant correlation between the independent variables and the observed outcomes, with statistical tests confirming the robustness of these relationships. Specifically, the results demonstrate that as variable $X$ increases, there is a corresponding increase in variable $Y$, suggesting a potential causal link.
Additionally, the analysis reveals that the model’s predictive accuracy is enhanced when incorporating variable $Z$, which accounts for approximately 25% of the variance in the dependent variable. The findings are further supported by graphical representations, illustrating the trends and interactions among the variables. Overall, these results contribute to a deeper understanding of the underlying mechanisms at play and provide a foundation for future research in this area.
Discussion
The discussion section of the research paper provides a comprehensive analysis of the structural evolution of iron oxide nanoparticles during the pyrolysis of a precursor, specifically FeCl₂·4H₂O combined with a nitrogen-doped carbon substrate (NC). In situ heating X-ray diffraction (XRD) and transmission electron microscopy (TEM) were employed to monitor phase transitions, revealing that the precursor is initially amorphous and transitions through several phases: at 350 °C, it forms hexagonal Fe₂O₃, which then converts to cubic Fe₃O₄ at 500 °C, and finally to cubic FeO at 800 °C. The study highlights the significant role of temperature in facilitating these transformations and the subsequent migration and reorganization of iron atoms, which contribute to the formation of active Fe-N₄ sites.
The findings indicate that nanoparticles smaller than 7 nm undergo gradual growth followed by rapid decomposition, while larger nanoparticles (>10 nm) exhibit a slower decomposition phase. The carbothermic reaction at elevated temperatures promotes the reduction of Fe species and the formation of Fe-N₄ sites, with a maximum conversion rate of 44.1% observed. However, only 16% of the total Fe atoms were found to be active in the oxygen reduction reaction (ORR), suggesting that a substantial portion of the iron remains inactive or is lost due to encapsulation by carbon layers. The authors propose strategies for optimizing the synthesis of Fe/N/C catalysts, including controlling nanoparticle size and employing additives that facilitate the formation of atomically dispersed Fe sites. Overall, this work elucidates the mechanisms underlying the conversion of iron oxide nanoparticles to active catalytic sites, providing insights for the rational design of high-performance catalysts.