Energy-Conserving Contact Dynamics of Nonspherical Rigid-Body Particles

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

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

• الديناميكا الجزيئية الذرية/المجردة (MD): رغم دقتها، إلا أن النماذج الذرية الكاملة مكلفة حاسوبياً للأنظمة ذات المقياس الميكروني. وغالباً ما تفشل طرق التجريد (Coarse-grained) في التقاط تفاعلات التلامس قصيرة المدى بدقة، مما يغير خشونة السطح الفعالة، والقوى المماسية، والاحتكاك.
• نماذج العناصر المنفصلة (DEMs): تعتمد نماذج DEM القياسية غالباً على نقاط تلامس منفصلة، مما يؤدي إلى مشكلتين حرجتين:
  1. انقطاعات القوة: يؤدي الكشف غير السلس لنقاط التلامس إلى قفزات مفاجئة في القوى المحسوبة.
  2. التداخلات غير الفيزيائية: فرض التلامس المنفصل يتم محلياً فقط، مما يسمح لأجزاء بعيدة من الجسيمات بالاختراق أثناء الحركات الانتقالية أو الدورانية الكبيرة.
• حفظ الطاقة: غالباً ما تنتهك هذه المشكلات حفظ الطاقة، مما يجعل عمليات المحاكاة غير موثوقة لدراسة الديناميكيات غير المتوازنة والخصائص الديناميكية الحرارية.

2. المنهجية

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

أ. صياغة قوة التلامس

يستخدم إطار العمل نموذج تلامس معدل من نمط "هيرتز" (Hertzian-style) يتضمن مكونات عمودية ومماسية:

• القوة العمودية ($F_n$): دالة مجزأة تتضمن حداً طارداً خطياً ($-k_n \delta_n$)، حداً مخمداً ($-c_n v_n$)، وحداً جاذباً للانفصالات التي تتجاوز سمك طبقة الجلد ولكنها تقع ضمن نصف قطر القطع ($r_c$).
• القوة المماسية ($F_t$): تم نمذجتها لالتقاط الاحتكاك وتبدد الطاقة، مع تقييدها بمعامل احتكاك حركي ($\mu$).
• طبقة الجلد (Skin Layer): تُعامل الجسيمات على أنها تمتلك "طبقة جلد" (هندسة إزاحة)، وهي مكافئة رياضياً لمجموع مينكوفسكي (Minkowski sum) لشكل الجسيم وكرة أو قرص. يضمن هذا انتقالاً سلساً للقوى.

ب. خوارزميات الكشف عن التلامس

يميز إطار العمل بين الهندسات ثنائية الأبعاد وثلاثية الأبعاد:

• ثنائي الأبعاد (رأس–حدود): بالنسبة للمضلعات، تقوم الخوارزمية بتقييم المسافات بين كل رأس من جسيم ما وكل حافة من الجسيم الآخر. وتحدد أقرب نقطة على الحدود لكل رأس.
• ثلاثي الأبعاد (رأس–سطح و حافة–حافة):
  • رأس–سطح: تحسب أقصر مسافة من رأس أحد الأشكال متعددة السطوح إلى سطح شكل آخر.
  • حافة–حافة: من الضروري أن هذا النوع يكشف عن التلامسات حيث تتقاطع أو تخترق حافتان بعضهما البعض، وهو سيناريو يغفله الكشف (رأس-سطح).
• التعامل مع الأزواج المكررة: تحتفظ الخوارمة بـ أزواج التلامس المكررة (على سبيل المثال، حيث يتطابق تلامس حافة-حافة مع تلامس رأس-سطح). هذا أمر ضروري لضمان استمرارية القوة؛ فمع تحرك الجسيمات، تتباعد هذه النقاط التي كانت متطابقة في البداية، وتجاهلها سيؤدي إلى انقطاعات في القوة.

ج. التنفيذ

• مدمج في حزمة body الخاصة ببرنامج LAMMPS.
• يستخدم وحدات "لينارد-جونز" المختزلة.
• يستفيد من حزمة coxeter بلغة بايثون لتوليد مدخلات الأشكال المعقدة.
• يدعم القابلية للتوسع في التوازي والتصور عبر برنامج OVITO.

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

1. إطار عمل موحد: طريقة قوية لمحاكاة الجسيمات الصلبة المحدبة التعسفية في كل من الأبعاد الثنائية والثلاثية مع الحفاظ الصارم على الطاقة الإجمالية.
2. كشف شامل للتلامس: يضمن تضمين تفاعلات (حافة-حافة) في الأبعاد الثلاثية والاحتفاظ بأزواج التلامس المكررة منع اختراق الجسيمات وضمان بقاء القوى مستمرة.
3. التحقق من حفظ الطاقة: تثبت الاختبارات الصارمة أن إطار العمل يحقق انحرافات في الطاقة في حدود $10^{-6}$ إلى $10^{-5}$، وهي مقاربة للديناميكا الجزيئية الذرية القياسية، حتى عند كسور التعبئة العالية.
4. توسيع القدرات: على عكس طرق مونت كارلو (MC) التي تكتفي بأخذ عينات من حالات التوازن، يتيح هذا الإطار دراسة الديناميكيات غير المتوازنة، والحركية، وخصائص النقل.

4. النتائج

تم التحقق من صحة إطار العمل من خلال عمليات محاكاة مكثفة لمختلف أشكال الجسيمات:

• محاكاة ثنائية الأبعاد (مضلعات):

  • تم اختبار المثلثات، المربعات، الخماسي، والسداسي عند كسور تعبئة ($\eta$) مختلفة.
  • الطاقة: لوحظ حفظ مثالي للطاقة في مجموعات NVE.
  • التعبئة: نجح في رصد الانتقال إلى التراصات المنتظمة (مثل تشكيل السداسيات لتبليطات مثالية) والإحباط الهندسي في الخماسي.
  • المقاييس: تم حساب دالات التوزيع الشعاعي $g(r)$ ودالات التوزيع الزاوي $g(\theta)$.

• محاكاة ثلاثية الأبعاد (متعددة السطوح وقضبان):

  • تم اختبار خليط من المكعبات، المكعبات الكبيرة، رباعيات الأوجه، المنشورات السداسية، والقضبان.
  • الطاقة: ظلت فروق الطاقة ضئيلة ($\sim 10^{-5}$) حتى عند كسور تعبئة عالية ($\eta \approx 0.72$).
  • الانتشار:
    • قياس الانتشار في الإطار المختبري ($D^{Lab}_T$) وإطار الجسم ($D^{Body}_T, D^{Body}_R$).
    • وُجد أن الأشكال متباينة الخواص تؤدي إلى انتشار اتجاهي (على سبيل المثال، تنتشر القضبان بشكل أسرع طولياً؛ وتدور المنشورات السداسية بشكل أسرع حول محورها).
  • التبلور: محاكاة التجميع الذاتي إلى هياكل مكعبة بسيطة (SC)، ومكعبة مركزية الوجه (FCC)، ومكعبة مركزية الجسم (BCC). تتبع الإطار بنجاح التطور الزمني لنسب التبلور.
  • معادلة الحالة (EOS): حساب الضغط المختزل ($p^*$) مقابل كسر التعبئة ($\eta$). طابقت النتائج عند الصلابة العالية ($k_n$) بيانات المرجع للجسيمات الصلبة من محاكاة مونت كارلو، مما يؤكد صحة تعريفات التلامس.

5. الأهمية والآفاق المستقبلية

• الأثر العلمي: يسد هذا العمل الفجوة بين كفاءة النماذج المجردة والدقة الفيزيائية المطلوبة لديناميكيات التلامس. إنه يوفر أداة موثوقة لاستقصاء الظواهر غير المتوازنة (مثل تدفق القص، والتحشر/الانسداد، وحركية التجميع الذاتي) التي لا يمكن الوصول إليها عبر تقنيات مونت كارلو القياسية.
• التطبيقات: إطار العمل قابل للتطبيق على:
  • التجميع الذاتي الغروي: فهم كيف يحدد تباين الشكل هيكل البلورات.
    টি الجريان الحبيبي: نمذجة الاحتكاك، وسلاسل القوة، وتبدد الطاقة في المواد السائبة.
  • الهيدروديناميكا: يمكن للتوسعات المستقبلية دمج تفاعلات الهيدروديناميكا وقوى التزييت، وهي ضرورية لتجميع الجسيمات النانوية.
• القابلية للتوسع: يسمح التصميم النمطي بدمج عدم التجانس السطحي (مواقع تفاعل محددة) والأشكال البيولوجية المعقدة في المستقبل، مما يسهل تصميم المواد المستجيبة والآلات الجزيئية.

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