Water cavitation results from the kinetic competition of bulk, surface and surface-defect nucleation events

ملخص تقني: نتائج تكهف المياه من التنافس الحركي بين أحداث النواة في الكتلة، والسطح، وعيوب السطح

بيان المشكلة
توجد المياه عند الضغوط السالبة في حالة شبه مستقرة حتى تصل إلى حالة التوازن عبر التكهف (cavitation)—وهو نشوء فقاعات البخار. تجريبياً، يظهر عتبة ضغط التكهف ($p_{cav}$) تباينًا واسعًا اعتمادًا على نقاء الماء، وزوايا التلامس السطحي، وجودة السطح. وبينما تتنبأ نظرية النواة الكلاسيكية (CNT) بأن المياه النقية في الكتلة يمكن أن تتحمل ضغوطًا سالبة تصل إلى ما دون $-100$ميجا باسكال، فإن الملاحظات التجريبية تظهر غالبًا حدوث التكهف عند ضغوط أعلى بكثير (على سبيل المثال،$-30$ ميجا باسكال أو حتى ضغوط موجبة في الحاويات القياسية). يشير هذا التباين إلى أن مسارات التكهف تتنافس: النواة المتجانسة في الكتلة مقابل النواة غير المتجانسة عند الأسطح وعيوب السطح. لقد كان غياب إطار حركي موحد يصف التنافس بين هذه المسارات الثلاثة (الكتلة، والسطح، والعيب) عائقاً أمام فهم التكهف في كل من الأنظمة الاصطناعية والبيولوجية.

المنهجية
يسد المؤلفون هذه الفجوة من خلال صياغة نموذج حركي يدمج نظرية النواة الكلاسيكية (CNT) مع محاكاة الديناميكا الجزيئية (MD) الذرية.

• إعداد المحاكاة: تستخدم الدراسة نموذج سطح مستوٍ يتكون من طبقتين من الطبقات أحادية الجزيء ذاتية التجميع (SAMs) من الألكانات الهيدروكسيلية. يتم ضبط قطبية السطح، وبالتالي زاوية تلامس الماء ($\theta$)، عن طريق تغيير مقياس الشحنات الجزئية للمجموعات الرأسية. النظام مشبع بالماء (نموذج SPC/E) ويخضع لضغوط سالبة.
• البروتوكول الحركي: لتحديد ترددات محاولة النواة ($\kappa$)، يستخدم المؤلفون بروتوكول محاكاة بزيادة تدريجية للضغط حيث ينخفض الضغط السالب خطياً مع مرور الوقت ($p(t) = \dot{p}t$). يتم تحديد التكهف عندما يزدل طول صندوق المحاكاة بنسبة 50%.
• الإطار النظري: يتم نمذجة معدل التكهف الإجمالي ($k_{tot}$) كمجموع لمعدلات مستقلة لكل من الكتلة ($k_{3D}$)، والسطح الأملس ($k_{2D}$)، وعيب السطح ($k_{def}$). يفترض النموذج أن هذه المسارات لا تتداخل. يستنتج المؤلفون تعبيرات لحواجز الطاقة الحرة ($G^*$) وترددات المحاولة لكل مسار، باستخدام بيانات الديناميكا الجزيئية (MD) لتعيين معاملات ترددات المحاولة ($\kappa_{3D}$ و $\kappa_{2D}$) التي يصعب الحصول عليها تحليلياً.

المساهمات والنتائج الرئيسية

1. تحديد ترددات المحاولة: حدد المؤلفون كثافة تردد محاولة الكتلة ($\kappa_{3D}$) وكثافة تردد محاولة السطح ($\kappa_{2D}$) كدالة لزاوية التلامس $\theta$ من محاكاة الديناميكا الجزيئية. ووجدوا أن $\kappa_{2D}$ يتغير مع $\theta$، حيث يظهر ذروة حول $75^\circ$، رغم أن هذا التباين يعتبر ثانوياً مقارنة بالاعتماد الأسي لحاجز الطاقة الحرة.
2. التنافس الحركي والتحول: يكشف النموذج عن انتقال حركي حاد بين تكهف الكتلة وتكهف السطح.
   • الأسطح الخالية من العيوب: بالنسبة للأسطح المحبة للماء ذات زوايا التلامس التي تقل عن العتبة الحرجة $\theta^* \approx 50^\circ - 60^\circ$، يهيمن تكهف الكتلة، ويكون $p_{cav}$ متوافقاً مع تنبؤات نظرية النواة الكلاسيكية (حوالي $-100$ ميجا باسكال).
   • الأسطح الكارهة للماء: عندما تتجاوز $\theta$ القيمة $\theta^*$، يصبح مسار التكهف السطحي هو المسار المهيمن. ويرتفع ضغط التكهف النموذجي بشكل ملحوظ ليصل إلى حوالي $-30$ ميجا باسكال للأسطح شديدة الكره للماء.
   • الاعتماد على حجم النظام: تعتمد زاوية التحول $\theta^*$ بشكل ضعيف فقط على حجم النظام ووقت المراقبة، حيث تتغير ضمن نطاق ضيق.
3. دور عيوب السطح: توضح الدراسة أن العيوب السطحية الكارهة للماء على المقياس النانوي تعمل كنوى تكهف عالية الكفاءة.
   • يمكن لعيب نانوي واحد حتى أن يهيمن على حركية التكهف في نظام ماكروسكوبي (عياني).
   • العيوب التي تحتوي على فقاعات بخار موجودة مسبقاً يمكن أن ترفع ضغط التكهف الحرج بشكل أكبر بكما لو كانت الأسطح ملساء.
   • يوفر النموذج إطاراً موحداً للتنبؤ بـ $p_{cav}$ بناءً على حجم العيب ($R_{def}$)، وزاوية تلامس العيب ($\theta_{def}$)، وعدد العيوب ($N_{def}$).
4. التحقق من الصحة: أكدت محاكاة الديناميكا الجزيئية الصريحة للأنظمة التي تحتوي على حفرة كارهة للماء واحدة (عيب) التنبؤات النظرية. الأنظمة التي تحتوي على عيوب بنصف قطر 1 نانومتر و2 نانومتر تعرضت للتكهف عند ضغوط أعلى بكثير ($-66$ميجا باسكال و$-45$ميجا باسكال، على التوالي) مقارنة بالأسطح الخالية من العيوب ($-95$ ميجا باسكال)، وهو ما يتوافق نوعياً مع النموذج.

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

• حد "الكتلة النقية" ($p_{cav} \approx -100$ ميجا باسكال) لا يمكن تحقيقه إلا في البيئات الخالية من العيوب وشديدة المحبة للماء (مثل تجاويف الكوارتز المجهرية).
• في معظم السيناريوهات العملية، يحكم التكهف خصائص السطح والعيوب بدلاً من الديناميكا الحرارية للكتلة.
• يوفر النموذج أداة تنبؤية لفهم التكهف في تطبيقات تتراوح من الطباعة بنفث الحبر والتنظيف بالموجات فوق الصوتية إلى الظواهر البيولوجية مثل صعود العصارة في الأشجار وميكانيكا الروبيان القافز (snapping shrimp).
• تشير النتائج إلى أن الامتزاز السطحي (على سبيل المثال، للجزيئات البرمائية) يمكن أن يثبط التكهف عن طريق تغيير زوايا التلامس أو تغليف العيوب، مما يوفر آليات محتملة للتحكم في استقرار التكهف.