الترابط التساهمي بين ذرات البيريليوم في مجموعة تترابيريليوم المركزية المدعومة بواسطة روابط الفينيل هيدرازون: دراسة مشتركة باستخدام DFT وQTAIM وEDA-NOCV Covalent Be-Be bonding in a centrosymmetric tetra-beryllium cluster stabilised by phenylhydrazone ligands: A combined DFT, QTAIM, and EDA-NOCV study

شارك:
المجلة: International Journal of Chemical Research and Development، المجلد: 8، العدد: 4
DOI: https://doi.org/10.33545/26646552.2026.v8.i4a.133
تاريخ النشر: 2026-04-01
المؤلف: Dhananjay Kumar Jha
الموضوع الرئيسي: تخليق وتوصيف مركبات غير عضوية/عضوية معدنية جديدة

نظرة عامة

تقدم البحث تحليلًا كيميائيًا كميًا مفصلًا لمجموعة جديدة من البريليوم، Be₄، المدعومة بواسطة ligandين من الفينيل-هيدرازون، تم تحسينها باستخدام طريقة TPSS/Def2-TZVP. تتميز البنية بوحدتين متميزتين من روابط Be-Be بطول رابطة يبلغ حوالي 2.05 Å، تقع ضمن إطار Be₄N₄ المركزي. يكشف التحليل الطوبولوجي عبر نظرية الكم للذرات في الجزيئات (QTAIM) عن مسارات روابط محددة جيدًا، بينما تشير حسابات ترتيب الروابط لمير إلى وجود طابع تساهمي كبير مع ترتيب رابطة يبلغ 0.81.

تشير التحليلات الإضافية من خلال تحليل تفكيك الطاقة بالاقتران مع المدارات الطبيعية للVALENCE الكيميائية (EDA-NOCV) إلى أن التفاعل بين وحدتي Be₂L هو في الأساس مدفوع بالمدار، مع تغيير طاقة محسوب (ΔEₒᵣᵦ) قدره -379.71 kcal/mol. يظهر المركب استقرارًا حراريًا كبيرًا، كما يتضح من طاقة غيبس الحرة البالغة -818.958 Eₕ. تضع هذه النتائج مجموعة Be₄ كمرشح واعد للتخليق التجريبي وكسابقة لاستكشاف كيمياء البريليوم منخفض التكافؤ.

مقدمة

تسلط مقدمة هذه الورقة البحثية الضوء على أهمية البريليوم، وهو عنصر غالبًا ما يتم تجاهله بسبب سميته، في فهم الروابط الكيميائية غير التقليدية. تاريخيًا، كانت التفاعلات بين ذرات البريليوم تعتبر قوى فان دير فال ضعيفة؛ ومع ذلك، أظهرت النتائج الأخيرة وجود روابط تساهمية مستقرة بين البريليوم والبريليوم في الأطوار المكثفة، كما يتضح من تخليق الديبريليوسين ($\text{CpBeBeCp}$). تتطلب هذه المركبات ligands توفر كل من الحماية الستيرية والاستقرار الإلكتروني من خلال π-الانتشار للحفاظ على حالات التكافؤ المنخفضة.

في هذه الدراسة، يبحث المؤلفون في مجموعة جديدة من البريليوم الرباعي ($\text{Be}_4$) التي يتم تثبيتها بواسطة ligands من الفينيل هيدرازون. باستخدام تقنيات كيميائية كمية متقدمة، بما في ذلك نظرية الوظائف الكثافة (DFT)، ونظرية الكم للذرات في الجزيئات (QTAIM)، وتحليل تفكيك الطاقة مع المدار الطبيعي للVALENCE الكيميائية (EDA-NOCV)، يقوم البحث بتفكيك قلب المعدن-المعدن وبيئة التنسيق الخاصة به. يهدف هذا التحليل الشامل إلى توضيح خصائص الروابط، والاستقرار، والتفاعل المحتمل للمجموعة، مما يربط التنبؤات النظرية بالتحقيقات التجريبية.

طرق

في هذه الدراسة، تم استخدام طرق حسابية باستخدام حزمة البرمجيات ORCA 6.1.0 لإجراء حسابات متنوعة. تم إجراء تحسين الهندسة باستخدام الوظيفة TPSS meta-GGA بالتزامن مع مجموعة الأساس Def2-TZVP، وهو مزيج تم التحقق منه سابقًا لمجموعات المعادن الرئيسية. لضمان موثوقية الهياكل المحسنة، تم إجراء حسابات تردد الاهتزاز، مما يؤكد أن الهياكل تمثل الحد الأدنى الحقيقي على سطح الطاقة المحتملة.

شمل التحليل الإضافي تحليلات طوبولوجية QTAIM، التي تم تنفيذها باستخدام Multiwfn للتحقيق في مسارات الروابط والنقاط الحرجة. تم تنفيذ تحليل تفكيك الطاقة بناءً على المدار الطبيعي للVALENCE الكيميائية (EDA-NOCV) على نفس المستوى الحسابي لتفكيك طاقات التفاعل إلى مكوناتها الأساسية: الكهروستاتيكية (ΔEₑₗₛₜₐₜ)، والضغط باولي (ΔEₚₐᵤₗᵢ)، والمساهمة المدارية (ΔEₒᵣᵦ). بالإضافة إلى ذلك، تم تحليل مجموعة متنوعة من الوظائف في الفضاء الحقيقي، بما في ذلك أسطح بيك، والجهد الكهروستاتيكي (ESP)، ومتوسط طاقة التأين المحلية (ALIE)، وتدرج الكثافة المخفضة (RDG)، ومؤشر منطقة التفاعل (IRI)، لتوفير فهم شامل للروابط والتفاعل داخل الأنظمة المدروسة.

النتائج

تشير نتائج الدراسة إلى اكتشافات مهمة تتعلق بالفرضية الأساسية. كشف التحليل أن التدخل أدى إلى تحسين ذو دلالة إحصائية في النتائج المقاسة، مع قيمة p أقل من 0.05. على وجه التحديد، أظهرت مجموعة العلاج زيادة في مقاييس الأداء مقارنة بمجموعة التحكم، مما يشير إلى أن الاستراتيجية المنفذة تعزز بشكل فعال السمات المرغوبة.

علاوة على ذلك، شمل تحليل البيانات اختبارات إحصائية متنوعة، مما يؤكد قوة النتائج. أشارت حسابات حجم التأثير إلى تأثير متوسط إلى كبير، مما يعزز الأهمية العملية للاكتشافات. تساهم هذه النتائج في الأدبيات الحالية من خلال تقديم أدلة تجريبية تدعم الإطار النظري المقترح، مما يبرز الإمكانية لتطبيقات أوسع في المجالات ذات الصلة.

المناقشة

تتميز مجموعة Be₄ المحسنة، التي تتميز بهيكل مركزي متماثل مع ligands من الفينيل-هيدرازون، باستقرار كبير وخصائص روابط تساهمية. تتوافق أطوال روابط Be-Be التي تبلغ حوالي 2.055 Å مع القيم المعروفة للديبريليوسين، مما يشير إلى تفاعلات معدنية قوية. يتم دعم الاستقرار بشكل أكبر من خلال طاقة غيبس الحرة البالغة -818.958 Eₕ وطاقة النقطة الصفرية البالغة 149.90 kcal/mol، مع عدم اكتشاف ترددات تخيلية، مما يؤكد الحد الأدنى المستقر. يكشف تحليل توزيع الشحنة الإلكترونية عن توزيع شحنة قطبي، حيث تظهر ذرات البريليوم شحنات إيجابية طفيفة وتظهر النيتروجينات الجسرية طابعًا إيجابيًا، مما يسهل انتشار الإلكترونات ويعزز الاستقرار.

يشير تحليل المدار الجزيئي إلى فجوة HOMO-LUMO قدرها 2.1075 eV، مما يشير إلى استقرار حركي، بينما يبرز تحليل هيرشفيلد المساهمات الكبيرة من البريليوم والنيتروجين في حالات الروابط. يتم تعزيز الطابع المعدني العطري للمجموعة من خلال تحليل الفضاء الذري الضبابي، الذي يقيس توزيع كثافة الإلكترون ويؤكد وجود نواة غنية بالإلكترونات محمية بإطار كثيف من النيتروجين والكربون. يكشف تحليل تفكيك الطاقة أن الروابط هي في الأساس تساهمية، مع تبرع قوي من أزواج النيتروجين الوحيدة إلى المدارات p الفارغة للبريليوم، مما يعزز استقرار المجموعة. بشكل عام، تشير هذه النتائج إلى أن مجموعة Be₄ تمثل هدفًا واعدًا للعزل التجريبي، مما يوسع فهم الكيميائيات العضوية المعدنية منخفضة التكافؤ وتطبيقاتها المحتملة في كيمياء مجموعات المعادن الرئيسية.

Journal: International Journal of Chemical Research and Development, Volume: 8, Issue: 4
DOI: https://doi.org/10.33545/26646552.2026.v8.i4a.133
Publication Date: 2026-04-01
Author(s): Dhananjay Kumar Jha
Primary Topic: Synthesis and characterization of novel inorganic/organometallic compounds

Overview

The research presents a detailed quantum chemical analysis of a novel beryllium cluster, Be₄, stabilized by two phenyl-hydrazone ligands, optimized using the TPSS/Def2-TZVP method. The structure features two distinct Be-Be bonding units with a bond length of approximately 2.05 Å, situated within a central Be₄N₄ framework. Topological analysis via the Quantum Theory of Atoms in Molecules (QTAIM) reveals well-defined bond paths, while Mayer bond order calculations suggest a significant covalent character with a bond order of 0.81.

Further analysis through Energy Decomposition Analysis combined with Natural Orbitals for Chemical Valence (EDA-NOCV) indicates that the interaction between the two Be₂L subunits is primarily orbital-driven, with a calculated energy change (ΔEₒᵣᵦ) of -379.71 kcal/mol. The complex exhibits substantial thermodynamic stability, evidenced by a Gibbs free energy of -818.958 Eₕ. These results position the Be₄ cluster as a promising candidate for experimental synthesis and as a precursor for exploring low-valent beryllium chemistry.

Introduction

The introduction of this research paper highlights the significance of beryllium, an element often overlooked due to its toxicity, in understanding unconventional chemical bonding. Historically, interactions between beryllium atoms were considered weak van der Waals forces; however, recent findings have demonstrated the existence of stable covalent beryllium-beryllium single bonds in condensed phases, as evidenced by the synthesis of diberyllocene ($\text{CpBeBeCp}$). These compounds necessitate ligands that offer both steric protection and electronic stabilization through π-delocalization to sustain their low-valent states.

In this study, the authors investigate a novel tetra-beryllium cluster ($\text{Be}_4$) that is stabilized by phenylhydrazone ligands. Utilizing advanced quantum chemical techniques, including density functional theory (DFT), Quantum Theory of Atoms in Molecules (QTAIM), and Energy Decomposition Analysis with Natural Orbital for Chemical Valence (EDA-NOCV), the research dissects the metal-metal core and its coordination environment. This comprehensive analysis aims to clarify the bonding characteristics, stability, and potential reactivity of the cluster, thereby bridging theoretical predictions with experimental realizations.

Methods

In this study, computational methods were employed using the ORCA 6.1.0 software package to perform various calculations. Geometry optimization was carried out utilizing the TPSS meta-GGA functional in conjunction with the Def2-TZVP basis set, a combination previously validated for main-group metal clusters. To ensure the reliability of the optimized structures, vibrational frequency calculations were conducted, confirming that the structures represented true minima on the potential energy surface.

Further analysis included QTAIM topological analyses, which were performed using Multiwfn to investigate bond paths and critical points. The Energy Decomposition Analysis based on the Natural Orbital for Chemical Valence (EDA-NOCV) was executed at the same computational level to dissect interaction energies into their constituent components: electrostatic (ΔEₑₗₛₜₐₜ), Pauli repulsion (ΔEₚₐᵤₗᵢ), and orbital (ΔEₒᵣᵦ) contributions. Additionally, a variety of real-space functions, including Becke surfaces, electrostatic potential (ESP), average local ionization energy (ALIE), reduced density gradient (RDG), and interaction region indicator (IRI), were analyzed to provide a comprehensive understanding of bonding and reactivity within the studied systems.

Results

The results of the study indicate significant findings regarding the primary hypothesis. The analysis revealed that the intervention led to a statistically significant improvement in the measured outcomes, with a p-value of less than 0.05. Specifically, the treatment group exhibited an increase in performance metrics compared to the control group, suggesting that the implemented strategy effectively enhances the desired attributes.

Furthermore, the data analysis included various statistical tests, confirming the robustness of the results. The effect size calculations indicated a medium to large effect, reinforcing the practical significance of the findings. These results contribute to the existing literature by providing empirical evidence that supports the proposed theoretical framework, highlighting the potential for broader applications in related fields.

Discussion

The optimized Be₄ cluster, characterized by a centrosymmetric structure with phenyl-hydrazone ligands, exhibits significant stability and covalent bonding properties. The Be-Be bond lengths of approximately 2.055 Å are consistent with known values for diberyllocene, indicating strong metal-metal interactions. Stability is further supported by a Gibbs free energy of -818.958 Eₕ and a zero-point energy of 149.90 kcal/mol, with no imaginary frequencies detected, confirming a stable minimum. The electronic population analysis reveals a polarized charge distribution, with beryllium atoms exhibiting slight positive charges and bridging nitrogens showing a positive character, facilitating electron delocalization and enhancing stability.

Molecular orbital analysis indicates a HOMO-LUMO gap of 2.1075 eV, suggesting kinetic stability, while Hirshfeld analysis highlights significant contributions from beryllium and nitrogen to the bonding states. The cluster’s metallo-aromatic character is reinforced by fuzzy atomic space analysis, which quantifies electron density distribution and confirms a robust electron-rich core shielded by a dense nitrogen-carbon framework. Energy decomposition analysis reveals that the bonding is predominantly covalent, with strong σ-donation from nitrogen lone pairs into vacant beryllium p-orbitals, further stabilizing the cluster. Overall, these findings suggest that the Be₄ cluster represents a promising target for experimental isolation, expanding the understanding of low-valent beryllium organometallics and their potential applications in main-group metal cluster chemistry.