Nonlinear dynamic elastic moduli from equilibrium stress fluctuations

ملخص تقني: معاملات المرونة الديناميكية غير الخطية من تقلبات الإجهاد عند التوازن

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

المنهجية
يشتق المؤلفون التعبيرات المطلوبة بدءاً من المعادلات المجهرية للحركة لنظام من الجسيمات المتفاعلة.

1. معادلات الحركة: يستخدم الاشتقاق معادلات DOLLS/SLLOD للحركة. يقيد المؤلفون تحليلهم بالتشوهات غير الدورانية (حيث يكون موتر تدرج السرعة $L$ متماثلاً)، وهو شرط تتطابق عنده صياغات DOLLS وSLLOD. يسمح هذا باستخدام إطار هاملتوني معزز بعنصر اقتران بين زخم الجسيمات ومجال سرعة التدفق.
2. المتغيرات المسحوبة (Pulled-Back Variables): للتعامل مع اعتماد الانفعال على الزمن، يستخدم المؤلفون تمثيلاً "مسحوباً" لمتغيرات فضاء الطور (الإحداثيات $\vec{R}_a$ والزخم $\vec{P}_a$). يقوم هذا التحويل برسم خريطة للحالة المشوهة الحالية عائدة إلى تكوين مرجعي، مما يضمن بقاء المتغيرات مستمرة حتى عند تطبيق انفعال مفاجئ.
3. التوزيع غير المتوازن: يبني المؤلفون دالة التوزيع الاحتمالي غير المتوازنة عبر تطوير التوزيع الكلاسيكي (Canonical) المتوازن الأولي. ويعبرون عن المتوسط غير المتوازن لموتر الإجهاد كمتوسط فوق المجموعة المتوازنة، بوزن يعتمد على الأس الخاص بتغير الطاقة الناتج عن تطور نقطة الطور ($\Delta H_\Gamma$).
4. توسيع الاضطراب: يتم توسيع موتر الإجهاد وتغير الطاقة في متسلسلة قوى بالنسبة لموتر الانفعال لاغرانج النهائي $E$. يفصل المؤلفون المساهمات إلى حدود خطية وغير خطية، محددين موترات "بورن-كينيتك" ($C^{BK}$ و $N^{BK}$) والتي تمثل ثوابت المرونة اللحظية المستمدة من طاقات الوضع والحركة.
5. اشتقاق دالة الارتباط: من خلال تعويض هذه التوسيعات في تعبير متوسط الإجهاد غير المتوازن والاحتفاظ بالحدود حتى الدرجة الثانية في الانفعال، يشتق المؤلفون صيغاً صريحة للمعاملات الديناميكية كدوال ارتباط تعتمد على الزمن تتضمن موتر الإجهاد ومصطلحات "بورن-كينيتك".

المساهمات والنتائج الرئيسية
النتيجة الأساسية للورقة هي اشتقاق تعبيرات تقلب مغلقة لكل من المعاملات المرونية الديناميكية الخطية وغير الخطية:

• المعاملات الديناميكية الخطية ($C_{ijkl}$): يستعيد المؤلفون صيغة التقلب اللزوجة المرنة الخطية المعروفة (المعادلة 45)، والتي تعبر عن المعامل كمزيج من المتوسط المتوازن لموتر "بورن-كينيتك" ودوال الارتباط الزمني لموتر الإجهاد.
• المعاملات الديناميكية غير الخطية ($N_{ijklmn}$): المساهمة المركزية هي اشتقاق صيغة التقلب لموتر المعامل الديناميكي غير الخطي من الدرجة السادسة (المعادلة 46). يتضمن هذا التعبير:
  • المتوسط المتوازن لموتر "بورن-كينيتك" غير الخطي ($N^{BK}$).
  • الارتباطات المتبادلة بين موتر الإجهاد وموتر "بورن-كينيتك" الخطي.
  • ارتباطات من رتب أعلى تتضمن حاصل ضرب موترات الإجهاد والفرق بين موتر "بورن-كينيتك" في أوقات مختلفة.
  • حدود تتناسب مع $\beta^2$ (مربع مقلوب درجة الحرارة) تتضمن ارتباطات ثلاثية للإجهاد.

تستوفي الموترات المشتقة خصائص تناظر محددة؛ حيث تمتلك موترات "بورن-كينيتك" تناظرات صغرى وكبرى، بينما تحقق المعاملات الديناميكية الزمنية الكاملة تناظرات صغرى فقط (وتناظراً تبادلياً محدداً للموتر غير الخطي).

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

يؤكد المؤلفون على الفائدة العملية لهذه النتائج:

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

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