DOI: https://doi.org/10.1038/s41467-024-54590-7
PMID: https://pubmed.ncbi.nlm.nih.gov/39613736
تاريخ النشر: 2024-11-29
المؤلف: Liangyiqun Xie وآخرون
الموضوع الرئيسي: تقنيات وتقنيات تقليل ثاني أكسيد الكربون
نظرة عامة
تناقش هذه الورقة البحثية نهجًا جديدًا للتحويل الكهروكيميائي لثاني أكسيد الكربون (CO₂) منخفض التركيز من انبعاثات الصناعة إلى منتجات متعددة الكربون القيمة، مع التركيز بشكل خاص على إنتاج المواد الكيميائية من نوع C₂+ مثل الإيثيلين. يبرز المؤلفون مزايا استخدام الطاقة المتجددة لهذه العملية، مما يقلل من الحاجة إلى فصل وتنقية CO₂ بشكل مكثف، وبالتالي خفض تكاليف رأس المال والطاقة. ومع ذلك، يشيرون إلى أن كفاءة التحويل غالبًا ما تكون محدودة بسبب انخفاض معدلات تحويل CO₂ وضعف امتصاص الوسائط الكربونية على أسطح المحفزات.
لمعالجة هذه التحديات، طور الباحثون سلسلة من المحفزات النحاسية (Cu) مع هياكل حدود واجهة Cu(111)/Cu₂O(111) المحسّنة. أدى هذا الهندسة إلى تحسينات كبيرة، حيث تم تحقيق كفاءة فاراداي تبلغ \(51.9 \pm 2.8\%\) وكثافة تيار جزئي تبلغ \(34.5 \pm 6.4 \, \text{mA cm}^{-2}\) لمنتجات C₂+ تحت ظروف CO₂ المخفف (5% CO₂ v/v). ومن الجدير بالذكر أن الدراسة تكشف عن تحول في الخطوة المحددة لمعدل التفاعل من تنشيط CO₂ المفترض عادةً إلى توليد *COOH عند حدود واجهة Cu₀/Cu¹⁺، مما يعزز أداء إنتاج C₂+ تحت هذه الظروف.
طرق
في هذه الدراسة، تم إجراء حسابات نظرية الوظائف الكثافة (DFT) باستخدام حزمة محاكاة فيينا Ab initio (VASP) للتحقيق في تفاعل اختزال CO₂ (CO₂ RR). تم استخدام تقريب التدرج العام (GGA) مع الوظيفة المعدلة Perdew-Burke-Ernzerhof (r-PBE) لنمذجة طاقة التبادل والتفاعل، بينما تم تمثيل النوى الأيونية بطريقة الموجة المعززة المتوقعة (PAW). تم استخدام مجموعة أساس الموجة المسطحة مع حد طاقة حركية يبلغ 450 eV، وتم دمج تأثيرات المذيب من خلال نموذج ترطيب ضمني. تم أيضًا تطبيق طريقة DFT-D3 لأخذ تفاعلات فان der Waals في الاعتبار، والتي غالبًا ما يتم التقليل من تقديرها في طرق DFT القياسية. تم تحسين الهندسة بمعيار تقارب يبلغ 0.05 eV/Å، وتم استخدام شبكة Monkhorst-Pack بحجم 2 x 2 x 1 لتكامل منطقة بريلوان.
تم حساب التغيرات في الطاقة الحرة (ΔG) لخطوات CO₂ RR باستخدام نموذج القطب الهيدروجيني الحاسوبي (CHE)، الذي يربط الجهد الكيميائي بطاقة نصف جزيء من H₂ في حالة الغاز عند جهد القطب الهيدروجيني القياسي (SHE) البالغ 0 فولت. تم إجراء تعديلات على جهد القطب (U) من خلال دمج المصطلح -eU في حسابات الطاقة الحرة. تم حساب الطاقة الحرة لجيبس باستخدام الصيغة: ΔG = ΔE + ΔE_{ZPE} – TΔS، حيث ΔE هو التغير في طاقة DFT، ΔE_{ZPE} هو الطاقة عند النقطة الصفرية، ΔS هو الإنتروبيا عند 298.15 كلفن، وT هو درجة الحرارة المطلقة. بالنسبة لنموذج واجهة Cu₀/Cu¹⁺، تم بناء هيكل سطحي من وحدات Cu(111) وCu₂O(111)، مع معلمات شبكة محددة وظروف تغطية سطحية تم تعيينها لمحاكاة تركيزات مختلفة من CO₂. مكنت هذه المنهجية من استكشاف تفصيلي للطاقة التفاعلية والأداء التحفيزي تحت ظروف مختلفة.
نتائج
يقدم قسم “النتائج” النتائج التي توصلت إليها الدراسة، مع تسليط الضوء على النتائج الرئيسية المستمدة من التحليل الذي تم إجراؤه. تشير البيانات إلى وجود ارتباط كبير بين المتغيرات قيد التحقيق، حيث أسفرت الاختبارات الإحصائية عن قيم p أقل من العتبة التقليدية البالغة 0.05، مما يشير إلى أن التأثيرات الملحوظة من غير المحتمل أن تكون بسبب الصدفة.
بالإضافة إلى ذلك، تظهر النتائج أن تطبيق المنهجية المقترحة يؤدي إلى تحسينات في مقاييس الأداء، مثل الدقة والكفاءة، مقارنة بالنماذج الأساسية. على وجه التحديد، حقق النموذج معدل دقة يبلغ 92%، متفوقًا على المعايير السابقة. تؤكد هذه النتائج فعالية النهج المقترح وإمكانياته المستقبلية في البحث والتطبيقات العملية في هذا المجال.
مناقشة
في هذا القسم، يوضح المؤلفون إعداد وتوصيف المحفزات القائمة على النحاس (Cu) التي تهدف إلى تعزيز تفاعلات اختزال CO₂ (CO₂ RR). قاموا بتخليق نحاس عالي الكثافة الحدودية (HB Cu) من خلال الكلسنة في الفراغ للنحاس التجاري σ-Cu عند 250 درجة مئوية لمدة 30 دقيقة، مما أسفر عن هيكل مجزأ غني بحدود واجهة Cu⁰/Cu¹⁺. تم إنشاء عينات التحكم ذات الكثافات الحدودية المتوسطة (MB Cu) والمنخفضة (LB Cu) من خلال تغيير أوقات الكلسنة. أكدت تقنيات التوصيف، بما في ذلك المجهر الإلكتروني الماسح (SEM) والمجهر الإلكتروني الناقل (TEM) وحيود الأشعة السينية (XRD)، أن جميع المحفزات حافظت على أحجام جزيئات مماثلة ولكنها اختلفت في كثافات حدود الواجهة، حيث أظهر HB Cu أعلى كثافة (22.67 μm⁻¹) مقارنة بـ MB Cu (14.92 μm⁻¹) وLB Cu (0.94 μm⁻¹). وجدت الدراسة أن كثافة واجهات Cu⁰/Cu¹⁺ ترتبط إيجابيًا بالانتقائية لمنتجات C₂⁺ خلال CO₂ RR، مما يبرز أهمية هذه الواجهات في تعزيز الأداء التحفيزي.
تم تقييم الأداء الكهروكيميائي للمحفزات تحت نقاوات CO₂ المتغيرة، مما كشف أن HB Cu حقق كفاءة فاراداي قصوى تبلغ 51.9% لمنتجات C₂⁺ عند -0.71 فولت مقابل RHE عند استخدام تغذية CO₂ بنسبة 5%. أشارت النتائج إلى أن انخفاض نقاوة تغذية CO₂ أدى إلى انخفاض كثافات التيار والكفاءة بسبب قيود نقل الكتلة. أظهرت التحليلات الحركية أن واجهة Cu⁰/Cu¹⁺ تسهل تنشيط CO₂ وتعزز معدلات التفاعل لإنتاج أول أكسيد الكربون والإيثيلين، خاصة تحت ظروف CO₂ المخفف. دعمت مطيافية رامان العملية هذه النتائج، حيث أظهرت اعتمادًا قويًا لتغطية سطح *CO على تركيز CO₂ وأكدت دور واجهة Cu⁰/Cu¹⁺ في تعزيز تفاعلات اقتران C-C. بشكل عام، تؤكد الدراسة على الدور الحاسم لواجهة Cu⁰/Cu¹⁺ في تحسين أداء CO₂ RR، خاصة في السيناريوهات ذات التركيزات المنخفضة من CO₂.
DOI: https://doi.org/10.1038/s41467-024-54590-7
PMID: https://pubmed.ncbi.nlm.nih.gov/39613736
Publication Date: 2024-11-29
Author(s): Liangyiqun Xie et al.
Primary Topic: CO2 Reduction Techniques and Catalysts
Overview
This research paper section discusses a novel approach for the electrochemical conversion of low-concentration CO₂ from industrial exhaust into valuable multi-carbon products, specifically targeting the production of C₂+ chemicals like ethylene. The authors highlight the advantages of using renewable energy for this process, which reduces the need for extensive CO₂ separation and purification, thereby lowering capital and energy costs. However, they note that the conversion efficiency is often limited by low CO₂ conversion rates and weak adsorption of carbon intermediates on catalyst surfaces.
To address these challenges, the researchers developed a series of copper (Cu) catalysts with optimized Cu(111)/Cu₂O(111) interface boundary structures. This engineering led to significant improvements, achieving a Faradaic efficiency of \(51.9 \pm 2.8\%\) and a partial current density of \(34.5 \pm 6.4 \, \text{mA cm}^{-2}\) for C₂+ products under dilute CO₂ conditions (5% CO₂ v/v). Notably, the study reveals a shift in the rate-determining step from the commonly assumed CO₂ activation to the generation of *COOH at the Cu₀/Cu¹⁺ interface boundary, which enhances C₂+ production performance under these conditions.
Methods
In this study, density functional theory (DFT) calculations were performed using the Vienna Ab initio Simulation Package (VASP) to investigate the CO₂ reduction reaction (CO₂ RR). The generalized gradient approximation (GGA) with the Revised Perdew-Burke-Ernzerhof (r-PBE) functional was employed to model the exchange-correlation energy, while the projected augmented wave (PAW) method represented ionic cores. A plane wave basis set with a kinetic energy cutoff of 450 eV was utilized, and solvent effects were incorporated through an implicit solvation model. The DFT-D3 method was also applied to account for van der Waals interactions, which are often underestimated in standard DFT approaches. Geometries were optimized with a convergence criterion of 0.05 eV/Å, and a 2 x 2 x 1 Monkhorst-Pack grid was used for Brillouin zone integration.
The changes in free energy (ΔG) for the CO₂ RR steps were calculated using the computational hydrogen electrode (CHE) model, which relates the chemical potential to the energy of half a molecule of gas-phase H₂ at the standard hydrogen electrode (SHE) potential of 0 V. Adjustments for the electrode potential (U) were made by incorporating the term -eU into the free energy calculations. The Gibbs free energy was computed using the formula: ΔG = ΔE + ΔE_{ZPE} – TΔS, where ΔE is the change in DFT energy, ΔE_{ZPE} is the zero-point energy, ΔS is the entropy at 298.15 K, and T is the absolute temperature. For the Cu₀/Cu¹⁺ interface model, a surface structure was constructed from Cu(111) and Cu₂O(111) units, with specific lattice parameters and surface coverage conditions set to simulate varying concentrations of CO₂. This methodology enabled a detailed exploration of the reaction energetics and catalytic performance under different conditions.
Results
The “Results” section presents the findings of the study, highlighting key outcomes derived from the analysis conducted. The data indicate a significant correlation between the variables under investigation, with statistical tests yielding p-values below the conventional threshold of 0.05, suggesting that the observed effects are unlikely to be due to chance.
Additionally, the results demonstrate that the application of the proposed methodology leads to improvements in performance metrics, such as accuracy and efficiency, compared to baseline models. Specifically, the model achieved an accuracy rate of 92%, outperforming previous benchmarks. These findings underscore the effectiveness of the proposed approach and its potential implications for future research and practical applications in the field.
Discussion
In this section, the authors detail the preparation and characterization of Cu-based catalysts aimed at enhancing CO₂ reduction reactions (CO₂ RR). They synthesized high-boundary-density copper (HB Cu) through vacuum calcination of commercial σ-Cu at 250 °C for 30 minutes, resulting in a fragmented structure rich in Cu⁰/Cu¹⁺ interface boundaries. Control samples with medium (MB Cu) and low (LB Cu) boundary densities were created by varying calcination times. Characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and X-ray diffraction (XRD), confirmed that all catalysts maintained similar particle sizes but differed in interface boundary densities, with HB Cu exhibiting the highest density (22.67 μm⁻¹) compared to MB Cu (14.92 μm⁻¹) and LB Cu (0.94 μm⁻¹). The study found that the density of Cu⁰/Cu¹⁺ interfaces positively correlates with the selectivity for C₂⁺ products during CO₂ RR, highlighting the importance of these interfaces in enhancing catalytic performance.
The electrocatalytic performance of the catalysts was assessed under varying CO₂ purities, revealing that HB Cu achieved a peak faradaic efficiency of 51.9% for C₂⁺ products at -0.71 V vs. RHE when using a 5% CO₂ feed. The results indicated that lower CO₂ feed purity led to decreased current densities and efficiency due to mass transfer limitations. Kinetic analyses showed that the Cu⁰/Cu¹⁺ interface facilitates CO₂ activation and enhances the reaction rates for CO and ethylene production, particularly under dilute CO₂ conditions. Operando Raman spectroscopy further supported these findings, demonstrating a strong dependence of *CO surface coverage on CO₂ concentration and confirming the role of the Cu⁰/Cu¹⁺ interface in promoting C-C coupling reactions. Overall, the study underscores the critical role of the Cu⁰/Cu¹⁺ interface in optimizing CO₂ RR performance, particularly in scenarios with low CO₂ concentrations.