DOI: https://doi.org/10.1038/s41598-026-37564-1
PMID: https://pubmed.ncbi.nlm.nih.gov/41629602
تاريخ النشر: 2026-02-02
المؤلف: Noura Al-Zoubi وآخرون
الموضوع الرئيسي: تخزين الهيدروجين والمواد
نظرة عامة
تستخدم هذه الدراسة حسابات المبادئ الأولى لاستكشاف الخصائص الهيكلية والديناميكية الحرارية والميكانيكية والإلكترونية والبصرية وتخزين الهيدروجين للهيدريدات مزدوجة البيروفكسايت K\(_2\)CaZnH\(_6\) وRb\(_2\)CaZnH\(_6\) وCs\(_2\)CaZnH\(_6\). من المتوقع أن تتبنى المركبات مراحل مكعبة مستقرة ديناميكيًا، مع ثوابت شبكية تتراوح من 8.302 Å إلى 8.652 Å. تتفاوت طاقات تشكيلها، التي تشير إلى الاستقرار الديناميكي الحراري، بين -0.891 eV/atom و-0.522 eV/atom. قدرات تخزين الهيدروجين متوسطة، مع سعات وزنية تبلغ 3.19 wt% لـ K\(_2\)CaZnH\(_6\) و2.15 wt% لـ Rb\(_2\)CaZnH\(_6\) و1.60 wt% لـ Cs\(_2\)CaZnH\(_6\)، إلى جانب كثافات حجمية ودرجات حرارة تفريغ ملائمة، مما يشير إلى إمكاناتها لتطبيقات تخزين الهيدروجين الثابتة.
تظهر المواد سلوكًا شبه موصل غير مباشر مع فجوات نطاق تتراوح من 1.79 إلى 1.83 eV وتظهر امتصاصًا ضوئيًا قويًا عبر طيف واسع من طاقة الفوتون. تؤكد التحليلات الميكانيكية أن هذه المركبات تلبي معايير استقرار بورن، حيث تظهر K\(_2\)CaZnH\(_6\) قابلية للانحناء، بينما يتم تحديد Cs\(_2\)CaZnH\(_6\) كأكثرها صلابة وقوة. يعزز الطابع الأيوني للرابطة، كما هو موضح بواسطة نسبة بواسون، إلى جانب انحناء الرابطة الملائم وفقًا لمعامل كلاينمان، احتباس الهيدروجين من خلال تغييرات هيكلية قابلة للعكس. تؤكد محاكاة الديناميكا الجزيئية استقرار الهيكل لهذه المركبات، مما يجعلها مرشحة واعدة لتطبيقات الطاقة المتعلقة بالهيدروجين، بما في ذلك مركبات خلايا الوقود وأنظمة الطاقة المتجددة. تضع النتائج أساسًا للتحقق التجريبي والدراسات التحسينية المستقبلية.
طرق
استخدمت الدراسة حسابات المبادئ الأولى المستندة إلى نظرية الوظائف الكثافة (DFT) للتحقيق في خصائص البيروفكسايت المزدوج X$_2$CaZnH$_6$ (حيث X = K، Rb، Cs). تم إجراء المحاكاة باستخدام حزمة CASTEP، مع استخدام الوظيفة Perdew-Burke-Ernzerhof (PBE) تحت تقريب التدرج العام (GGA). على الرغم من أن طريقة GGA-PBE معروفة بأنها تقلل من قيم فجوة النطاق، إلا أنها تلتقط بفعالية الاتجاهات النوعية في الهيكل الإلكتروني اللازمة لفهم سلوك هذه المواد شبه الموصلة. تتبلور المركبات في هيكل بيروفكسايت مزدوج مكعب مع مجموعة فضاء Fm3̅m عند درجة حرارة الغرفة، مع توفير معلمات هيكلية مفصلة لتحسين الهندسة والحسابات اللاحقة.
تم نمذجة التفاعل بين النوى الأيونية والإلكترونات التكافؤية باستخدام إمكانيات فاندربيلت فائقة النعومة، مع تخصيص تكوينات تكافؤ محددة لكل عنصر. تم استخدام مجموعة أساس موجية مع حد طاقة حركية يبلغ 500 eV، وتم إجراء تكامل منطقة بريلوان باستخدام شبكة نقاط k من Monkhorst-Pack بحجم 12 × 12 × 12. تم استرخاء الهياكل البلورية بالكامل إلى حالتها الأساسية باستخدام خوارزمية Broyden-Fletcher-Goldfarb-Shanno (BFGS)، مع الالتزام بمعايير تقارب صارمة. بالنسبة لحساب الثوابت المرنة، تم تطبيق عتبات أكثر صرامة، وتم استخدام مخطط متوسط Voigt-Reuss-Hill (VRH) لاشتقاق الخصائص الميكانيكية. بالإضافة إلى ذلك، تم حساب أطياف انتشار الصوت باستخدام طريقة الإزاحة المحدودة لتقييم الاستقرار الديناميكي للمركبات، مما يوفر رؤى حول خصائصها الاهتزازية واستقرار الطور.
النتائج
يقدم قسم النتائج نتائج الدراسة، مع تسليط الضوء على النتائج الرئيسية وتأثيراتها. يكشف التحليل عن ارتباطات كبيرة بين المتغيرات قيد التحقيق، مما يظهر أن النموذج المقترح يتنبأ بفعالية بالظواهر الملحوظة. على وجه الخصوص، تشير البيانات إلى علاقة قوية، تم قياسها بمعامل ارتباط قدره $r = 0.85$، مما يشير إلى ارتباط قوي.
علاوة على ذلك، تشير النتائج إلى أن التدخل المطبق في الدراسة أدى إلى تحسين قابل للقياس في النتائج المستهدفة، مع زيادة ذات دلالة إحصائية (p < 0.01) لوحظت في المجموعة التجريبية مقارنة بالمجموعة الضابطة. تدعم هذه النتائج الفرضية القائلة بأن التدخل له تأثير إيجابي، مما يستدعي مزيدًا من الاستكشاف والتطبيق المحتمل في المجالات ذات الصلة. بشكل عام، تسهم النتائج في تقديم رؤى قيمة حول الآليات الأساسية وتوفر أساسًا للبحوث المستقبلية.
المناقشة
في هذا القسم، يتم تحليل الخصائص الهيكلية وتخزين الهيدروجين للبيروفكسايت المزدوج $X_2 \text{CaZnH}_6$ (حيث $X = \text{K، Rb، Cs}$). يكشف تحسين هندسة وحدة الخلية أن الثوابت الشبكية المحسوبة ومعلمات الاستقرار الهيكلي، مثل عامل تحمل غولدشميت ($\tau$) وعامل الأوكتاهيدرا ($\mu$)، تشير إلى أن هذه المركبات تظهر استقرارًا ديناميكيًا حراريًا، مع تأكيد طاقات التشكيل السلبية على التفاعلات الطاردة للحرارة. تقع قيم $\tau$ و$\mu$ للمركبات ضمن النطاقات التي تشير إلى هياكل بيروفكسايت مستقرة. من الجدير بالذكر أن $K_2 \text{CaZnH}_6$ يظهر أعلى سعة تخزين هيدروجين وزنية، بينما يظهر $Cs_2 \text{CaZnH}_6$ أقلها، مما يتماشى مع حجم كاتيون موقع A. يتميز آلية التحلل الحراري لإطلاق الهيدروجين بتفاعل ماص للحرارة، مع درجات حرارة تفريغ متغيرة عبر المركبات، حيث يظهر $Rb_2 \text{CaZnH}_6$ درجة حرارة تفريغ ملائمة تبلغ 385 K.
تم تقييم الخصائص الإلكترونية للبيروفكسايت باستخدام وظائف تبادل-ترابط مختلفة، مما يكشف أن جميع المركبات تظهر سلوك شبه موصل مع فجوات نطاق تتراوح من 1.79 eV إلى 1.83 eV. هذه القيم أضيق من تلك الخاصة بالبيروفكسايت الأكسيد، مما يضع مركبات $X_2 \text{CaZnH}_6$ بشكل ملائم لتطبيقات تخزين الهيدروجين والإلكترونيات الضوئية. يشير تحليل الكثافة الكلية والجزئية للحالات إلى مساهمات كبيرة من مدارات p وd في الهيكل الإلكتروني، وهو أمر حاسم لفهم إمكاناتها في تطبيقات تخزين الهيدروجين. علاوة على ذلك، تشير الخصائص البصرية، بما في ذلك الدالة العازلة المعقدة، ومعامل الانكسار، ومعامل الامتصاص، إلى أن هذه المواد قد تكون مناسبة لمجموعة متنوعة من التطبيقات الإلكترونية الضوئية، مع حدوث قيم ذروة الموصلية الضوئية ضمن الطيف المرئي. بشكل عام، تسلط النتائج الضوء على الجدوى الهيكلية وقدرات تخزين الهيدروجين الواعدة لمركبات $X_2 \text{CaZnH}_6$، مما يستدعي مزيدًا من الاستكشاف لتطبيقاتها العملية.
DOI: https://doi.org/10.1038/s41598-026-37564-1
PMID: https://pubmed.ncbi.nlm.nih.gov/41629602
Publication Date: 2026-02-02
Author(s): Noura Al-Zoubi et al.
Primary Topic: Hydrogen Storage and Materials
Overview
This study employs first-principles calculations to explore the structural, thermodynamic, mechanical, electronic, optical, and hydrogen storage properties of the double perovskite hydrides K\(_2\)CaZnH\(_6\), Rb\(_2\)CaZnH\(_6\), and Cs\(_2\)CaZnH\(_6\). The compounds are predicted to adopt thermodynamically stable cubic phases, with lattice constants ranging from 8.302 Å to 8.652 Å. Their formation energies, indicating thermodynamic stability, vary between -0.891 eV/atom and -0.522 eV/atom. Hydrogen storage capabilities are moderate, with gravimetric capacities of 3.19 wt% for K\(_2\)CaZnH\(_6\), 2.15 wt% for Rb\(_2\)CaZnH\(_6\), and 1.60 wt% for Cs\(_2\)CaZnH\(_6\), alongside volumetric densities and favorable desorption temperatures, suggesting their potential for stationary hydrogen storage applications.
The materials exhibit indirect semiconducting behavior with band gaps ranging from 1.79 to 1.83 eV and demonstrate strong optical absorption across a wide photon energy spectrum. Mechanical analyses confirm that these compounds meet the Born stability criteria, with K\(_2\)CaZnH\(_6\) showing ductility, while Cs\(_2\)CaZnH\(_6\) is identified as the stiffest and hardest. The ionic bonding character, indicated by Poisson’s ratio, along with favorable bond bending as per the Kleinman parameter, enhances hydrogen retention through reversible structural changes. Molecular dynamics simulations further affirm the structural stability of these compounds, making them promising candidates for hydrogen-related energy applications, including fuel cell vehicles and renewable energy systems. The findings lay a foundation for future experimental validation and optimization studies.
Methods
The study utilized first-principles calculations grounded in density functional theory (DFT) to investigate the properties of double perovskites X$_2$CaZnH$_6$ (where X = K, Rb, Cs). Simulations were conducted using the CASTEP package, employing the Perdew-Burke-Ernzerhof (PBE) functional under the generalized gradient approximation (GGA). Although the GGA-PBE method is known to underestimate bandgap values, it effectively captures the qualitative trends in electronic structure necessary for understanding the semiconducting behavior of these materials. The compounds crystallize in a cubic double perovskite structure with space group Fm3̅m at room temperature, with detailed structural parameters provided for geometry optimization and subsequent calculations.
The interaction between ionic cores and valence electrons was modeled using Vanderbilt ultrasoft pseudopotentials, with specific valence configurations assigned to each element. A plane-wave basis set with a kinetic energy cutoff of 500 eV was employed, and Brillouin zone integration was performed using a Monkhorst-Pack k-point mesh of 12 × 12 × 12. The crystal structures were fully relaxed to their ground state using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm, adhering to strict convergence criteria. For the computation of elastic constants, more stringent thresholds were applied, and the Voigt-Reuss-Hill (VRH) averaging scheme was used to derive mechanical properties. Additionally, phonon dispersion spectra were calculated using the finite displacement method to assess the dynamical stability of the compounds, providing insights into their vibrational properties and phase stability.
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 indicates a strong relationship, quantified by a correlation coefficient of $r = 0.85$, suggesting a robust association.
Furthermore, the results indicate that the intervention applied in the study led to a measurable improvement in the target outcomes, with a statistically significant increase (p < 0.01) observed in the experimental group compared to the control group. These findings support the hypothesis that the intervention has a positive effect, warranting further exploration and potential application in related fields. Overall, the results contribute valuable insights into the underlying mechanisms and provide a foundation for future research.
Discussion
In this section, the structural and hydrogen storage properties of the double perovskites $X_2 \text{CaZnH}_6$ (where $X = \text{K, Rb, Cs}$) are analyzed. The optimization of unit cell geometry reveals that the calculated lattice constants and structural stability parameters, such as the Goldschmidt tolerance factor ($\tau$) and octahedral factor ($\mu$), indicate that these compounds exhibit thermodynamic stability, with negative formation energies confirming exothermic reactions. The values of $\tau$ and $\mu$ for the compounds fall within the ranges that suggest stable perovskite structures. Notably, $K_2 \text{CaZnH}_6$ demonstrates the highest gravimetric hydrogen storage capacity, while $Cs_2 \text{CaZnH}_6$ shows the lowest, correlating with the size of the A-site cation. The thermal decomposition mechanism for hydrogen release is characterized by an endothermic reaction, with varying desorption temperatures across the compounds, particularly $Rb_2 \text{CaZnH}_6$ exhibiting the most favorable desorption temperature of 385 K.
The electronic properties of the perovskites were assessed using different exchange-correlation functionals, revealing that all compounds exhibit semiconductor behavior with band gaps ranging from 1.79 eV to 1.83 eV. These values are narrower than those of oxide perovskites, positioning the $X_2 \text{CaZnH}_6$ compounds favorably for applications in hydrogen storage and optoelectronics. The analysis of the total and partial density of states indicates significant contributions from p and d orbitals to the electronic structure, which is crucial for understanding their potential in hydrogen storage applications. Furthermore, the optical properties, including the complex dielectric function, refractive index, and absorption coefficient, suggest that these materials may be suitable for various optoelectronic applications, with peak optical conductivity values occurring within the visible spectrum. Overall, the findings highlight the structural viability and promising hydrogen storage capabilities of the $X_2 \text{CaZnH}_6$ compounds, warranting further exploration into their practical applications.