Hydrogen Bonding Governs Vapor–Liquid Equilibria in Alcohol-Based Fuels: A Thermodynamic Equation-of-State Approach

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

قام فريق من الباحثين من جامعة شنغهاي للمحيطات ومركز هلمهولتز-تسن دريسدن-روسندورف بتطوير طريقة جديدة لرسم خرائط هذه التفاعلات لثلاثة كحوليات محددة تُستخدم في الصناعة وهي: الأيزوبوتانول، و1-بوتانول، و2-إيثيل-1-هكسانول. وتعد هذه المواد مكونات حيوية في إنتاج البلاستيك والمذيبات والوقود المتجدد. ركز الباحثون على المخاليط الناتجة عن دمج هذه الكحوليات معاً، وهو سيناريو يتكرر كثيراً في المصانع الصناعية ولكنه كان من الصعب التنبؤ به بدقة. ومن خلال إنشاء إطار رياضي موحد يأخذ في الاعتبار الطريقة التي تتشبث بها هذه الجزيئات ببعضها البعض، تمكن الفريق من وصف سلوكها بدقة متناهية. ويوفر عملهم دليلاً موثوقاً للمهندسين عند تصميم المعدات اللازمة لفصل وتنقية هذه المواد الكيميائية القيمة.

يكمن التحدي في دراسة هذه المخاليط في الاختلافات الدقيقة بين الجزيئات. فبينما يتشابه الأيزوبوتانول و1-بوتانول بشكل كبير، حيث يتشاركان في الوزن الكيميائي نفسه، فإن 2-إيثيل-1-هكسانول أكبر بكثير وأكثر ضخامة. وهذا الحجم الزائد يخلق عائقاً فيزيائياً، أو ما يعرف بالإعاقة الفراغية، حول الجزء من الجزيء الذي يقوم بعملية الالتصاق. تخيل أنك تحاول مصافحة شخص ما بينما ترتدي قفازاً سميكاً وضخماً؛ الاتصال لا يزال ممكناً، لكنه أضعف قليلاً ويصعب تشكيله مقارنة بالقبضة باليد المجردة. في حالة 2-إيثيل-1-هكسانول، يجعل هذا "القفاز" من الصعب على الجزيئات تكوين روابط قوية مع جيرانها مقارنة بجزيئات البوتانول الأكثر نحافة. وكان الباحثون بحاجة إلى طريقة يمكنها استيعاب هذا الفارق الدقيق دون الحاجة إلى مجموعة جديدة تماماً من القواعد لكل خليط على حدة.

ولحل هذه المشكلة، استخدم الفريق نموذجاً ديناميكياً حرارياً يسمى معادلة الحالة "المكعبة-المضافة-الارتباط" (Cubic-Plus-Association). يجمع هذا النهج بين طريقة قياسية لحساب سلوك الغازات والسوائل ونظرية محددة تأخذ في الاعتبار "الارتباط" أو التلاصق بين الجزيئات. وقد عامل الباحثون كل جزيء كحول على أنه يمتلك موقعين محددين: موقع يقدم ذرة هيدروجين للارتباط، وموقع آخر يستقبل ذرة هيدروجين من جار له. نجح هذا النموذج البسيط ذو الموقعين مع جميع الكحوليات الثلاثة، بما في ذلك 2-إيثيل-1-هكسانول الضخم. وكان مفتاح نجاحهم هو تعديل قوة الرابطة للجزيء الضخم؛ حيث وجدوا أن طاقة الارتباط، التي تقيس مدى تمسك الجزيئات ببعضها البعض، كانت أقل بالنسبة لـ 2-إيثيل-1-هكسانول مقارنة بالنوعين الآخرين. وقد أكد هذا أن الضخامة الفيزيائية للجزيء تضعف بالفعل الترابط الهيدروجيني، وهو تفصيل غالباً ما تغفله النماذج الأبسط.

بعد تحديد قواعد الجزيئات الفردية، انتقل الفريق إلى دراسة المخاليط. واختبروا مدى قدرة نموذجهم على التنبؤ بسلوك ثلاثة أزواج مختلفة: الأيزوبوتانول مع 1-بوتانول، والأيزوبوتانول مع 2-إيثيل-1-هكسانول، و1-بوتانول مع 2-إيثيل-1-هكسانول. وقارنوا حساباتهم بمجموعة كبيرة من البيانات التجريبية الجديدة التي جمعها علماء آخرون، والتي غطت نطاقاً واسعاً من درجات الحرارة والضغوط. كانت النتائج دقيقة بشكل مذهل. فبالنسبة لخليط البوتانولين، اللذين يتشابهان جداً، لم يحتج النموذج إلا إلى تعديل طفيف ليتطابق مع البيانات الواقعية. أما بالنسبة للخليط الذي يتضمن 2-إيثيل-1-هكسانول الضخم، فقد ظل النموذج يعمل بشكل جيد، رغم أنه تطلب تعديلاً أكثر تعقيداً قليلاً لمراعاة الطبيعة المتغيرة للتفاعلات عند درجات حرارة مختلفة.

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

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