Molecular dynamics simulation of high slip flow of water confined between graphene nanochannels at experimentally accessible strain rates

ملخص تقني لورقة بحثية: محاكاة الديناميكا الجزيئية للتدفق عالي الانزلاق للماء المحصور بين قنوات الجرافين النانوية عند معدلات انفعال يمكن الوصول إليها تجريبياً.

1. بيان المشكلة

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

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

2. المنهجية

استخدم المؤلفون طريقة دالة الارتباط للزمن العابر (TTCF) لسد الفجوة بين نطاقات الاستجابة الخطية وغير الخطية، مما سمح بمحاكاة معدلات القص التي يمكن الوصول إليها تجريبياً.

• الإطار النظري (TTCF):

  • على عكس طريقة NEMD القياسية (المتوسط المباشر أو DAV)، التي تحسب متوسط مسار واحد عبر الزمن، تقوم طريقة TTCF بحساب استجابة النظام من خلال ربط كمية فيزيائية $B$ عند الزمن $t$ بدالة التبديد $\Omega$ عند $t=0$ (عند تشغيل المجال).
  • المعادلة الأساسية هي: $\langle B(t) \rangle = \langle B(0) \rangle + \int_0^t ds \langle B(s)\Omega(0) \rangle$.
  • تستفيد هذه الطريقة من طبيعة "الخلط" (mixing) في النظام للوصول إلى الحالة المستقرة حتى عندما تكون القوة الدافعة ضعيفة، مما يتجاك بفعالية قيود نسبة الإشارة إلى الضجيج الناتجة عن المتوسط المباشر.

• إعداد المحاكاة:

  • النظام: 320 جزيء ماء محصورة بين جدران من الجرافين (3 طبقات لكل منها، بتراصف A-B-A).
  • النماذج:
    • الماء: نموذج النقاط الشحنية البسيطة الممتدة (SPC/E) مع قيود SHAKE.
    • الجرافين: جهد Tersoff ثلاثي الأجسام.
    • التفاعلات: جهد لندلارد-جونز (قطع عند 10 أنجستروم) وطريقة (PPPM) للقوى الكولومبية بعيدة المدى.
  • الظروف: الحفاظ على درجة الحرارة عند 300 كلفن باستخدام منظم حراري "لانجفنز" (مطبق على الطبقات الداخلية للجدار)؛ عرض القناة $h \approx 2$ نانومتر.
  • البروتوكول:
    • المسارات الأم (Mother Trajectories): تم توليد 8,000 مسار توازن مستقل.
    • المسارات الابنة (Daughter Trajectories): من هذه المسارات، تم إطلاق 4 ملايين مسار غير متوازن إجمالاً عن طريق تطبيق مجالات قص.
    • معدلات القص: شملت ست درجات أسية، من $\dot{\gamma} = 5 \times 10^5 \, \text{s}^{-1}$ إلى $5 \times 10^{11} \, \text{s}^{-1}$.
    • التحقق: تمت مقارنة النتائج مع طريقة DAV القياسية (عند معدلات قص عالية) والديناميكا الجزيئية المتوازنة (EMD) باستخدام طريقة Green-Kubo/Varghese.

3. المساهمات الرئيسية

1. أول تطبيق لـ TTCF على نظام الماء والجرافين: تعد هذه الدراسة أول دراسة تطبق طريقة TTCF على نظام الماء والجرافين للوصول إلى معدلات الانفعال ذات الصلة التجريبية عبر NEMD.
2. الربط بين المحاكاة والتجربة: نجحت الدراسة في حساب أطوال الانزلاق ومعاملات الاحتكاك عند معدلات انفعال ($10^5 - 10^9 \, \text{s}^{-1}$) يمكن الوصول إليها في التجارب، وهو النطاق الذي تفشل فيه محاكاة NEMS المباشرة بسبب الضجيج.
3. التحقق من المنهجية: أثبت المؤلفون أن نتائج TTCF تتطابق مع نتائج DAV عالية نسبة الإشارة إلى الضبوح (عند القص العالي) ومع الطرق المتوازنة (عند القص المنخفض)، مما يؤكد تفوق TTCF كأداة للأنظمة ذات الانزلاق العالي والانفعال المنخفض.
4. المعايرة الكمية: قدمت دراسة شاملة تقارن النتائج الحسابية مع مجموعة واسعة من الأدبيات التجريبية (التدفق الشعري، مجهر القوة الذرية AFM، والتدفق المدفوع بالتبخر).

4. النتائج الرئيسية

• التحقق من TTCF مقابل DAV:
  • عند معدلات قص عالية ($\dot{\gamma} = 5 \times 10^{10} \, \text{s}^{-1}$)، أنتجت كل من TTCF وDAV ملفات تعريف سرعة وأطوال انزلاق متطابقة (حوالي 45-49 نانومتر)، مما يؤكد دقة الخوارزمية.
  • أظهرت بيانات TTCF ضجيجاً حرارياً أقل بكثير مقارنة بـ DAV.
• الأداء عند معدلات القص المنخفضة:
  • مع انخفاض معدلات القص دون $10^9 \, \text{s}^{-1}$، فشلت طريقة DAV تماماً بسبب ضعف نسبة الإشارة إلى الضجيج (أصبحت الأخطاء غير قابلة للإدارة).
  • وعلى العكس من ذلك، ظل الخطأ الإحصائي في TTCF مستقراً أو تحسن بالنسبة للإشارة، مما سمح بحساب دقيق لأطوال الانزلاق وصولاً إلى $\dot{\gamma} = 5 \times 10^5 \, \text{s}^{-1}$.
• طول الانزلاق ومعامل الاحتكاك:
  • طول الانزلاق ($L_s$): تراوحت القيم المحسوبة بين 49.5 نانومتر (النطاق الخطي) و 102.7 نانومتر (النطاق غير الخطي عند $\dot{\gamma} \sim 10^{11} \, \text{s}^{-1}$).
  • معامل الاحتكاك ($\xi$): تم حساب معامل احتكاك نافييه كـ $\xi = \tau_{zx} / v_s$.
    • النطاق الخطي (EMD/TTCF): $\xi \approx 1.56 \times 10^4 \, \text{kg}/(\text{m}^2\text{s})$.
    • الانتقال غير الخطي: لوحظ انخفاض سريع في $\xi$ بين $5 \times 10^9$ و $5 \times 10^{10} \, \text{s}^{-1}$.
• المقارنة مع الأدبيات:
  • تتوافق أطوال الانزلاق المحسوبة (نطاق 29-80 نانومتر لنماذج SPC/E) بشكل جيد مع الدراسات الحسابية السابقة (مثل Kannam et al. و Celebi et al.).
  • تُظهر النتائج توافقاً ملحوءاً مع البيانات التجريبية من التدفق الشعري (مثل Xie et al. و Ashok Keerthi et al.) ودراسات AFM، رغم الاختلافات في الشحنة السطحية وتأثيرات الركيزة في التجارب.

5. الأهمية

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