Structure–Property Relationships and Force-Field Transferability in 1-Propanol/Water Mixtures from 243 to 298 K

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

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

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

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

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

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

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

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