Chemical equilibrium constants of diatomic molecules from a q-deformed hyperbolic interaction potential
تستخدم هذه الدراسة جهد تفاعل زائدي مُعدل بـ q لحساب والتحقق من ثوابت الاتزان الكيميائي لتكوين BrF وClF وIF وSO مقابل بيانات NIST-JANAF، وذلك عبر طريقة نيكيفوروف-أوفاروف البارامترية، مما أظهر دقة عالية ضمن نطاقات درجات حرارة محددة.
المؤلفون الأصليون:E. S. Eyube, D. Festus, A. Zakariya, Saïdou Alioum, P. A. Njidda, I. Joel, T. G. Nyagong, J. Abe
في عالم الذرات غير المرئي، لا تُعد الجزيئات مجرد خرزات ساكنة، بل هي أنظمة ديناميكية تترابط عبر قوى غير مرئية تتمدد وتتضغط مثل النوابض. ولفهم كيفية سلوك هذه الهياكل المتناهية الصغر، يعتمد العلماء على خرائط رياضية تسمى دالات طاقة الوضع. تصف هذه الخرائط كيف تتغير الطاقة بين ذرتين كلما اقتربتا من بعضهما البعض أو ابتعدتا، مما يحدد كل شيء بدءاً من مدى قوة تمسكهما ببعضهما وصولاً إلى كيفية اهتزازهما عند التسخين. ومن خلال معرفة مستويات الطاقة هذه، يمكن للباحثين التنبؤ بكيفية تصرف الجزيئات في مجموعات كبيرة، وحساب الخصائص المعتمدة على درجة الحرارة التي تحرك التفاعلات الكيميائية. إن هذا الارتباط بين الرقصة المجهرية للذرات الفردية والسلوك العياني للغازات أمر ضروري لمجالات تتراوح من علوم الغلاف الجوي إلى هندسة الاحتراق، مما يسمح لنا بالتنبؤ بما إذا كان التفاعل سيحدث وإلى أي مدى سيمضي في ظل ظروف مختلفة.
لقد سعى فريق من الباحثين مؤخراً إلى تحسين هذه الخرائط باستخدام نموذج رياضي محدد يُعرف باسم "جهد التفاعل الزائدي المشوه بـ q" (q-deformed hyperbolic interaction potential). يقدم هذا النموذج طريقة مرنة لوصف القوى المعقدة وغير الخطية التي تربط الذرات معاً، لا سيما بالنسبة للجزيئات ثنائية الذرة، والتي تتكون من ذرتين فقط. ركز الفريق على أربعة جزيئات محددة: بروم فلوريد، وكلور فلوريد، ويود فلوريد، وأكسيد الكبريت. هذه ليست اختيارات عشوائية؛ فالثلاثة الأولى هي مركبات هالوجينية مكونة من ذرات هالوجين مختلفة، بينما يعد أكسيد الكبريت لاعباً رئيسياً في عمليات الاحتراق وكيمياء الغلاف الجوي. أراد الباحثون معرفة ما إذا كان نموذجهم الرياضي المحدد يمكنه التنبؤ بدقة بثوابت الاتزان الكيميائي لهذه الجزيئات. وثابت الاتزان هو في الأساس رقم يخبر العلماء مقدار ما سيكون موجوداً من ناتج مقارنة بمكوناته الأولية عند درجة حرارة معينة، وهو معلومة حاسمة لفهم الاستقرار الكيميائي وإنتاجية التفاعلات.
للقيام بذلك، قام الفريق أولاً بحل المعادلات الأساسية التي تحكم حركة الذرات داخل هذه الجزيئات. وقد استخدموا تقنية رياضية متخصصة لإيجاد مستويات الطاقة الدقيقة المرتبطة باهتزاز الذرات، مع معاملة الرابطة بينهما كنظام مهتز. وبمجرد حصولهم على مستويات الطاقة هذه، قاموا ببناء وصف إحصائي للجزيئات، يُعرف باسم "دالة التجزئة" (partition function)، والتي تجمع كل الطرق الممكنة لتخزين الجزيئات للطاقة. تضمن هذا الحساب اهتزازات الذرات، ودورانها في الفضاء، وحركتها عبر الهواء. ومن خلال دالة التجزئة هذه، قام الباحثون بحساب طاقة جيبس الحرة، وهي كمية ثرموديناميكية تحدد اتجاه ومدى التفاعل الكيميائي. ولجعل النموذج يتوافق مع الواقع بأكبر قدر ممكن، طبق الباحثون عامل تصحيح صغيراً لكل جزيء على قيم طاقة جيبس المحسوبة، مما أدى إلى ضبط النتائج بدقة بحيث تتطابق مع البيانات المرجعية المعروفة. ثم استخدموا قيم طاقة جيبس هذه لاستخلاص ثوابت الاتزان لتكوين كل جزيء من متفاعلاته متماثلة النواة.
وقد تمت مقارنة نتائج هذا العمل بقاعدة بيانات موثوقة للقيم التجريبية تُعرف باسم "NIST-JANAF"، والتي تعد المعيار الذهبي للبيانات الثرموكيميائية. كشفت المقارنة أن النموذج الجديد دقيق بشكل ملحوظ، لكن أداءه يعتمد بشدة على الجزيء المحدد ونطاق درجة الحرارة المدروس. فبالنسبة لبروم فلوريد، أنتج النموذج نتائج قريبة جداً من البيانات المرجعية عبر نطاق واسع من درجات الحرارة، من 300 كلفن حتى 6000 كلفن، بمتوسط فرق يقل عن واحد بالمائة. وبالنسبة لكلور فلوريد، عمل النموذج بشكل جيد جداً حتى 4600 كلفن، لكن الدقة بدأت في التراجع عند درجات الحرارة الأعلى. وأظهر يود فلوريد نمطاً مشابهاً، حيث ظهرت أكبر التباينات عند درجات الحرارة الأدنى، بينما قدم أكسيد الكبريت حالة فريدة حيث عانى النموذج بشكل كبير تحت 900 كلفن ولكنه توافق بشكل شبه مثالي مع البيانات المرجعية عند درجات الحرارة المرتفعة.
تثبت هذه النتائج أن "جهد التفاعل الزائدي المشوه بـ q" هو أداة قوية للتنبؤ بالسلوك الكيميائي، ولكنه ليس حلاً واحداً يناسب الجميع. توضح الدراسة أنه بينما يمكن للنموذج إعادة إنتاج بيانات الاتزان المرجعية بدقة عالية عبر فترات درجات حرارة واسعة لبعض الجزيئات، إلا أن له حدوداً محددة تصبح فيها الانحرافات ملحوظة. فبالنسبة لبروم فلوريد، يعد النموذج موثوقاً عبر النطاق المختبر بالكامل. أما بالنسبة للآخرين، فقد حدد الباحثون نوافذ درجات الحرارة المحددة التي تظل فيها التنبؤات جديرة بالثقة والتي تبدأ عندها في التباعد عن القيم الراسخة. يوفر هذا العمل تقييماً كمياً واضحاً للمواضع التي ينجح فيها هذا النهج الرياضي والمواضع التي يخفق فيها، مما يقدم طريقة منقحة لحساب ثوابت الاتزان للجزيئات ثنائية الذرة دون الحاجة إلى الاعتماد حصرياً على القياسات التجريبية لكل حالة جديدة.
ملخص تقني: ثوابت الاتزان الكيميائي للجزيئات ثنائية الذرة من جهد تفاعل هيبيربولي (زائدي) مشوه بـ q
بيان المشكلة يعتمد التنبؤ بثوابت الاتزان الكيميائي للأنظمة الغازية بشكل أساسي على مستويات الطاقة الجزيئية الدقيقة، والتي تُشتق من جهود التفاعل. وبينما توجد نماذج عديدة لجهد طاقة التفاعل تحليلية وشبه تحليلية (مثل نماذج مورس، كراتزر، ودينج-فان)، إلا أن هناك حاجة لتوسيع نطاق هذه النماذج لتتجاوز الثوابت الطيفية إلى الكميات الديناميكية الحرارية العيانية، وتحديداً طاقات غيبس الحرة وثوابت الاتزان. وعلى الرغم من دراسة الجهود المشوهة بـ q في الأطياف الجزيئية والوظائف الديناميكية الحرارية، فإن تطبيقها في التحديد الميكانيكي الإحصائي لثوابت الاتزان الكيميائي عبر فترات درجات حرارة ممتدة لا يزال غير مستكشف بشكل كافٍ. تعالج هذه الدراسة الفجوة في تطبيق جهد التفاعل الهيبربولي المشوه بـ q (qHIP) للتنبؤ بثوابت الاتزان لتكوين جزيئات محددة ثنائية الذرة (BrF، ClF، IF، وSO) من متفاعلاتها متماثلة النواة.
المنهجية يستخدم البحث إطاراً نظرياً يعتمد على جهد التفاعل الهيبربولي المشوه بـ q (qHIP)، والمعرف بقوة التفاعل (U0)، ومعلمة التغطية (α)، ومعلمة التشويه (q). تتبع المنهجية الخطوات التالية:
حل معادلة شرودنجر: يتم حل معادلة شرودنجر الشعاعية التي تحتوي على qHIP باستخدام طريقة نيكيفوروف-أوفاروف البارامترية (PNUS) للحصول على القيم الذاتية للطاقة الاهتزازية (Eν). وقد تم إهمال حد الطرد المركزي للتركيز على الحركة الاهتزازية.
تحديد المعلمات: لا يتم تعيين معلمات الجهد بشكل عشوائي، بل تُشتق من الثوابت الجزيئية القابلة للقياس: طول الرابطة عند التوازن (re)، وطاقة التفكك (De)، والعدد الموجي الاهتزازي (ωe). ويتم ذلك من خلال فرض شروط التوازن، والتفكك، والانحناء على الجهد، باستخدام دالة لامبرت W لإيجاء معلمة التشويه.
بناء دالة التجزئة:
الاهتزازية: يتم بناء دالة التجزئة الاهتزازية (Qa) من القيم الذاتية للطاقة المنفصلة. ولتعويض عدم الانتظام (اللانمطية) من الرتب العليا والاقتران الاهتزازي-الدوراني الذي لا يستوعبه الجهد بشكل صريح، تم إدخال معلمة تصحيح عديمة الأبعاد (η) تعتمد على الجزيء. يتم تحسين هذه المعلمة عن طريق تقليل الانحراف بين طاقات غيبس الحرة المحسوبة والمرجعية. ويتم تقييم المجموع المحدود باستخدام طريقة الجمع بواسون تحت تقريب الحد الكلاسيكي (مع الاحتفاظ بالنمط fourier j=0 فقط).
الحساب الديناميكي الحراري: تُستخدم دالة التجزئة الجزيئية الكلية لحساب طاقة غيبس الحرة (G). ويتم حساب التغير في طاقة غيبس الحرة (ΔG) لتفاعل التكوين (X+Y→XY) باستخدام آنثالبي التكوين القياسية والأنثالبي المولية الحساسة المشتقة من دوال التجزئة.
تحديد ثابت الاتزان: يتم اشتقاق ثابت الاتزان الكيميائي (Kp) من ΔG باستخدام العلاقة lnKp=−ΔG/RT.
التحقق من الصحة: تتم مقارنة قيم log10Kp المحسوبة مع بيانات NIST-JANAF الكيميائية الحرارية. ويتم تقييم الدقة باستخدام الانحراف المطلق النسبي (RAD) ومتوسط الانحراف المطلق النسبي (ARAD) عبر فترات درجات حرارة محددة.
المساهمات الرئيسية
توسيع نطاق qHIP: يعمل هذا العمل على توسيع تطبيق جهد التفاعل الهبربولي المشوه بـ q من أطياف الطاقة الجزيئية إلى التنبؤ بثوابت الاتزان الكيميائي، مما يؤسس رابطاً مباشراً بين نموذج التشويه بـ q وديناميكا التفاعل الحرارية.
الإطار التحليلي: يوفر صياغة تحليلية كاملة لدالة التجزئة الاهتزازية للجهود المشوهة بـ q، مع دمج معلمة تصحيح لمواءمة النتائج النظرية مع البيانات الديناميكية الحرارية المرجعية.
التقييم المنهجي: تقيم الدراسة بشكل منهجي أداء qHIP لأربعة أنظمة ثنائية الذرة مختلفة (BrF، ClF، IF، وSO)، مع تحديد نطاقات درجات الحرارة التي يتوافق فيها النموذج بشكل وثيق مع البيانات المرجعية.
النتائج حسبت الدراسة ثوابت الاتزان لـ BrF، وClF، وIF، وSO عبر نطاقات درجات حرارة واسعة وقارنتها ببيانات NIST-JANAF. وتشير النتائج إلى أداء يعتمد على نوع الجزيء ودرجة الحرارة:
BrF: أظهر أعلى مستوى من الاتساق، مع متوسط انحراف مطلق نسبي (ARAD) قدره 0.345% على كامل النطاق من 300 إلى 6000 كلفن. وتظل الانحرافات أقل من 1% في معظم هذه الفترة.
ClF: أظهر توافقاً جيداً بمتوسط انحراف مطلق نسبي قدره 0.751% عبر 300–4600 كلفن. ومع ذلك، تزدادة الانحرافات بشكل كبير فوق 4600 كلفن، لتصل إلى 19.063% عند 6000 كلفن.
IF: أظهر انحرافاً كبيراً عند أدنى درجة حرارة (8.075% عند 300 كلفن) ولكنه يتحسن بسرعة. وفوق الفترة 400–6000 كلفن، بلغ متوسط الانحراف المطلق النسبي 0.871%، مع بقاء الانحرافات صغيرة بشكل عام.
SO: أظهر انحرافات جوهرية عند درجات الحرارة المنخفضة (187.144% عند 300 كلفن)، والتي تتناقص مع ارتفاع درجة الحرارة. ويحقق النموذج توافقاً وثيقاً من 900 كلفن فصاعداً، بمتوسط انحراف مطلق نسبي قدره 0.905% عبر الفترة 900–6000 كلفن.
الأهمية والادعاءات يدعي البحث أن جهد التفاعل الهيبربولي المشوه بـ q يوفر إطاراً ميكانيكياً إحصائياً قابلاً للتطبيق للتنبؤ بثوابت الاتزان الكيميائي. وتكمن أهمية العمل في قدرته على إعادة إنتاج بيانات اتزان NIST-JANAF المرجعية بانحرافات متوسطة صغيرة ضمن فترات درجات حرارة محددة للجزيئات المختارة.
ويخلص المؤلفون بتواضع إلى أنه بينما تعتبر الصياغة فعالة ضمن النطاقات المحددة (على سبيل المثال، 300–6000 كلفن لـ BrF، و900–6000 كلفن لـ SO)، فإن حالات الابتعاد الملحوظة عند أطراف هذه النطاقات (درجات الحرارة العالية لـ ClF، ودرجات الحرارة المنخفضة لـ SO) تحدد الحدود الحرارية للصياغة الحالية. ولا يدعي البحث قابلية التطبيق الشاملة عبر جميع درجات الحرارة أو جميع الجزيئات، بل يحدد بدقة الفترات التي تتوافق فيها ثوابت qHIP مع البيانات الكيميائية الحرارية الراسخة، مما يؤكد جدوى نموذج التشويه بـ q في الديناميكا الحرارية لتوازن التفاعل ضمن تلك الحدود.