Universal Dynamic Scaling for Sub-Yield Aging and Over-Yield Rejuvenation in A Cyclically Driven Metallic Glass
تحل هذه الدراسة الجدل القائم حول الأصول المتباينة للتقادم دون حد الخضوع وتجدد ما فوق حد الخضوع في الزجاج المعدني المدفوع دورياً، وذلك من خلال الكشف عن ساعة ديناميكية موحدة حيث يتم دفع تجدد ما فوق حد الخضوع بواسطة عمل القص البلاستيكي العياني، بينما يحكم التقدم دون حد الخضوع تقلبات الإجهاد المرن الموضعية من الدرجة الثانية، مما يؤسس لهوية قياس عالمية تربط بين الريولوجيا غير الخطية وديناميكيات تقلب الإجهاد.
تخيل مادة تبدو وتُلمس كالمعدن الصلب، ومع ذلك تترتب ذراتها بفوضى عارمة تشبه السائل المتجمد. هذا هو الزجاج المعدني. وخلافًا للمعادن العادية، التي تمتلك هياكل بلورية منتظمة ومتكررة تسمح لها بالانثناء والتمدد دون انكسار، يفتقر هذا الزجاج إلى هذا النظام الداخلي. بدلاً من ذلك، تكون ذراته محشورة معًا في مشهد طاقي وعر وغير مستوٍ. وعندما تضغط عليها أو تسحبها، يجب على الذرات أن تقفز من أحد الوديان العميقة في هذا المشهد إلى وادٍ آخر لتتحرك. لقد عرف العلماء منذ فترة طويلة أن كيفية حركتها تعتمد على مدى قوة الضغط الممارس عليها. فإذا ضغطت بلطف، تميل المادة إلى الاستقرار بشكل أعمق في وديانها، لتصبح أكثر صلابة واستقرارًا، وهي عملية تُسمى "التعتيق" (aging). أما إذا ضغطت بقوة كافية، فإنك تهز الذرات خارج وديانها العميقة وتدفعها إلى نقاط أعلى طاقة وأقل عمقًا، مما يجعل المادة أكثر ليونة وقابلية للتشكيل، وهي عملية تُسمى "التجديد" (rejuvenation). ولسنوات، تنازع الباحثون حول ما إذا كان هذان السلوكان المتضادان — الاستقرار مقابل الاهتزاز — ينبعان من نفس الفيزياء الأساسية أم أنهما مدفوعان بقواعد مجهرية مختلفة تمامًا.
لقد استخدم فريق من الباحثين الآن عمليات محاكاة حاسوبية قوية لحسم هذا الجدل، وكشفوا أن كلاً من التعتيق والتجديد محكومان بساعة واحدة موحدة تملي على المادة كيفية التنقل في مشهد الطاقة الخاص بها. ومن خلال دراسة نموذج لزجاج معدني من النحاس والزركونيوم، اكتشفوا أن الفرق بين السلوكين ليس تغيراً في القواعد الأساسية، بل هو تحول في كيفية حقن الطاقة في النظام. فعندما يتم الضغط على المادة بما يتجاوز نقطة الخضوع (yield point) — وهو الحد الذي تبدأ عنده بالتدفق بشكل دائم — يكون التجديد مدفوعًا مباشرة بالعمل واسع النطاق للتشوه اللدن. ومع ذلك، عندما يتم الضغط على المادة بلطف تحت هذا الحد، فإن عملية التعتيق لا تكون مدفوعة بالضغط الإجمالي، بل بتقلبات الإجهاد الموضعية الصغيرة التي تعمل كمحفز ميكانيكي، حيث تدفع الذرات بالكاد بما يكفي لمساعدتها على الاستقرار.
قام الباحثون بمحاكاة سلوك 16,000 ذرة في كتلة من هذا الزجاج، وإخضاعها لدورات متكررة من التمدد والضغط. وقد قاموا بتغيير شدة التمدد وسرعة تطبيق القوة. وفي عمليات المحاكاة، عندما كان التمدد لطيفًا، كانت المادة تفقد الطاقة باستمرار، حيث تهبط إلى حالات أكثر استقرارًا وعمقًا بمرور الوقت، وهذا أكد تأثير التعتيق. ولكن عندما زادوا شدة التمدد لتتجاوز عتبة معينة، انقلب السلوك؛ بدأت المادة في اكتساب الطاقة، وتتسلق خارج وديانها العميقة لتستقر في حالة ذات طاقة أعلى وأكثر سيولة. وظهرت النتائج الأكثر إثارة للاهتمام عند الحد الفاصل بين هذين العالمين؛ ففي مستويات الإجهاد المتوسطة، كانت المادة تكتسب الطاقة في البداية، كما لو كانت تخضع لعملية تجديد، ثم تعود لتفقد الطاقة وتبدأ في التعتيق مجددًا بعد عدد معين من الدورات. هذا السلوك غير الرتيب، حيث تستيقظ المادة أولاً ثم تعود للنوم، قد لوحظ من قبل ولكن كان يفتقر إلى تفسير واضح.
ومن خلال تحليل البيانات، وجد الفريق علاقة رياضية دقيقة تربط بين النظامين. فقد حددوا أن المعدل الذي تتجدد به المادة تحت الإجهاد الثقيل يتبع قانون قوة (power law) محددًا، بينما يتبع المعدل الذي تعتق به المادة تحت الإجهاد الخفيف قانون قوة مختلفًا. والأهم من ذلك، وجدوا أن الأس (exponent) الذي يحكم عملية التعتيق هو ضعف الأس الذي يحكم عملية التجديد تمامًا. هذات النسبة البسيطة كشفت عن الآلية الخفية: فالاهتزاز الشديد للمادة يخلق تموجات صغيرة موضعية في الإجهاد. وتحت الأحمال الثقيلة، يتم طغيان هذه التموجات بفعل التدفق الهائل للمادة، أما تحت الأحمال الخفيفة، حيث يتوقف التدفق العام، تصبح هذه التموجات الصغيرة هي القوة المهيمنة. فهي تعمل كمصدر للضجيج الميكانيكي، مما يمنع الذرات من العلوق في أعمق الحفر ويبطئ عملية التعتيق. وكلما زادت سرعة دورات المادة، أصبحت هذه التموجات أقوى، مما يجعل من الصعب على الذرات الاستقرار.
يوفر هذا الاكتشاف إطارًا واحدًا للتنبؤ بكيفية تطور المادة. وقد أظهر الباحثون أن نقطة التحول بين التجديد الأولي والتعتيق اللاحق يمكن التنبؤ بها من خلال موازنة الطاقة المحقونة بواسطة القوة الخارجية مقابل الطاقة المبددة بواسطة تقلبات الإجهاد الداخلية هذه. الأمر كما لو أن المادة تمتلك مؤقتًا داخليًا يحصي عدد الدورات اللازمة قبل أن تسيطر التموجات الموضعية الصغيرة وتدفع النظام نحو حالة مستقرة. هذا الاكتشاف يجسّر الفجوة بين التدفق واسع النطاق للمادة وبين الارتجافات المجهرية لذراتها. ويشير إلى أنه من خلال التحكم الدقيق في سرة وكثافة التدريب الميكانيكي، يمكن للمهندسين تغيير خصائص الزجاج المعدني حسب الحاجة، سواء بجعله أكثر صلابة أو أكثر ليونة. إن هذا العمل يثبت أنه حتى في مادة تتسم بالفوضى مثل الزجاج، فإن المسار من الفوضى إلى النظام محكوم بمجموعة ثابتة وعالمية من القواعد.
ملخص تقني: القياس الديناميكي الشامل للشيخوخة دون حد الخضوع وتجديد البنية فوق حد الخضوع في زجاج معدني مدفوع دورياً
بيان المشكلة يمثل التطور الهيكلي للمواد الصلبة غير المتبلورة، وتحديداً الزجاج المعدني (MGs)، تحت تأثير التشوه الميكانيكي الدوري، ازدواجية جوهرية. فاعتماداً على سعة الانفعال (γA)، يؤدي التحفيز الدوري إما إلى شيخوخة ميكانيكية (نطاق ما دون الخضوع)، حيث يسترخي النظام نحو مستويات طاقة أدنى وأكثر استقراراً، أو تجديد هيكلي (نطاق فوق الخضوع)، حيث يُدفع النظام نحو حالات طاقة أعلى وأقل استقراراً. وثمة جدل طويل لم يُحل حول ما إذا كانت هذه المسارات المتعارضة تنشأ من آليات مجهرية متميزة أم أنها تشترك في أصل فيزيائي واحد. علاوة على ذلك، كان هناك غياب ملحوظ لإطار كمي موحد يربط حركية الاسترخاء في نطاق ما دون الخضوع بالريولوجيا غير الخطية لتدفق ما فوق الخضوع، لا سيما فيما يتعلق بديناميات الانتقال عند سعات انفعال متوسطة حيث يحدث تطور غير رتيب للطاقة.
المنهجية استخدم المؤلفون محاكاة ديناميكا جزيئية (MD) كلاسيكية مكثفة باستخدام حزمة LAMMPS على نظام زجاج معدني ثنائي نموذجي، Cu64Zr36، يتكون من 16,000 ذرة. تم نمذجة التفاعلات بين الذرات باستخدام جهد "طريقة الذرة المضمنة" (EAM) الذي تم التحقق من صحته جيداً.
بروتوكول المحاكاة: تم إعداد التكوينات الأولية من الطور السائل عبر التبريد السريع من 2000 كلفن إلى 300 كلفن بمعدل 1 كلفن/بيكو ثانية، متبوعاً بالاسترخاء الحراري.
التشوه: تم تطبيق شد وانضغاط أحادي المحور دوري على طول اتجاه x باستخدام بروتوكول انفعال جيبي، γ(t)=γAsin(2πt/tp)، يغطي نطاقاً واسعاً من سعات الانفعال (γA) وفترات التحميل (tp).
تحليل البيانات: تم تتبع التغير في الطاقة الكامنة (ΔU) كدالة لعدد دورة التحميل (Ncycle)، مع تسجيلها حصرياً عند حالات الانفعال الصفري لاستبعاد المساهمات المرنة العابرة. كما تم حساب المتوسطات الإحصائية عبر 10 تكوينات مستقلة تم إنشاؤها عبر طريقة "الإيزوكونفيغوراشن" (isoconfigurational method).
المساهمات والنتائج الرئيسية
تحديد نطاقات القياس المتميزة: يحدد البحث الاستجابة الديناميكية عبر فضاء المعلمات (γA, tp)، ويحدد أربعة مجالات تشغيلية:
المنطقة 1 (ما دون الخضوع): شيخوخة متسارعة ميكانيكياً رتيبة.
المنطقتان 2 و3 (العبور): منطقة انتقال ديناميكي تتميز بمسارات غير رتيبة (تجديد أولي يليه شيخوخة).
المنطقة 4 (فوق الخضوع): تجديد مستدام فوق حد الخضوع يصل إلى هضبة طاقة في الحالة المستقرة.
قوانين القياس العالمية: كشف المؤلفون عن أسّين عالميين متميزين يحكمان هذه النطاقات:
تجديد ما فوق الخضوع (المنطقة 4): تتناسب هضبة طاقة الحالة المستقرة (ΔUp) مع العمل البلاستيكي الكلي المستمد من نموذج "هيرشل-بولكي" (H-B) الريولوجي. علاقة القياس هي ΔUp∝γAγ˙α، مع أس تخلخل القص α≈0.17.
شيخوخة ما دون الخضوع (المنطقة 1): يتبع تطور الطاقة منحنى رئيسي (master curve) عند قياسه بـ γ˙β/γA، مما يعطي أس شيخوخة β≈0.34. يتم مطابقة التطور باستخدام دالة اضمحلال أسي ممدد.
هوية القياس الموحدة (β≈2α): تتمثل إحدى النتائج المركزية في العلاقة الكمية β≈2α (0.34 ≈ 2 × 0.17). يقترح البحث آلية فيزيائية موحدة لتفسير ذلك:
الاستجابة من الدرجة الأولى (فوق الخضوع): التجديد مدفوع مباشرة بالعمل البلاستيكي القصي الكلي (∝γ˙α).
الاستجابة من الدرجة الثانية (ما دون الخضوع): في غياب التدفق البلاستيكي العالمي، تتم عملية الشيخوخة الهيكلية عبر تقلبات الإجهاد المرن الموضعية من الدرجة الثانية (⟨δσlocal2⟩). تعمل هذه التقلبات، التي ترث مقياس الإجهاد من المصفوفة المدفوعة (∝γ˙α)، كمصدر ضجيج ميكانيكي. وتتدرج طاقة الاضطراب وفق ⟨δσlocal2⟩∝(γ˙α)2=γ˙2α.
الآلية: يعمل الاضطراب المكثف عالي المعدل على زيادة احتمالية الهروب من مستويات طاقة دنيا عميقة، مما يقلل من كفاءة الاحتفاظ بالاسترخاء الهيكلي. وقد أثبت المؤلفون وجود مبرر ديناميكي يشير إلى أن β≈2α، مما يوحي بأصل فيزيائي محتمل يربط بين الاسترخاء في ما دون الخضوع والريولوجيا في ما فوق الخضوع.
ديناميات العبور (Nc): بالنسب, لسابعات الانفعال المتوسطة (المنطقتان 2 و3)، يظهر النظام عبوراً عابراً من التجديد إلى الشيخوخة عند عدد دورات مميز Nc. صاغ المؤلفون Nc كتوازن ديناميكي بين حقن العمل البلاستيكي من الدرجة الأولى وتبديد تقلبات الإجهاد المرن من الدرجة الثانية. كما اشتقوا علاقة قياس عالمية: Ncγ˙2α=G(γ˙α/γA)، حيث G هي دالة قياس عالمية. تعمل هذه "الساعة الديناميكية" على تحديد الانتقال من الإثارة قصيرة المدى (التجديد) إلى الاسترخاء الجماعي طويل المدى (الشيخوخة).
الأهمية والادعاءات يزعم البحث حل الازدواجية بين شيخوخة ما دون الخضوع وتجديد ما فوق الخضوع من خلال إنشاء ساعة ديناميكية موحدة تتحكم في التنقل عبر مشهد الطاقة عبر كلا النطاقين.
الجسر النظري: يربط البحث مباشرة بين الريولوجيا غير الخطية (إطار هيرشل-بولكي) وديناميات تقلبات الإجهاد، مما يوفر إطاراً واحداً للتنبؤ بالتطور الهيكلي.
الرؤية الميكانيكية: يعيد تعريف شيخوخة ما دون الخضوع ليس مجرد استرخاء سلبي، بل كعملية محكومة بتقلبات الإجهاد المرن الموضعية من الدرجة الثانية التي تعمل كـ "محفز ميكانيكي" ومصدر ضجيج.
القدرة التنبؤية: من خلال تحديد هوية القياس β≈2α والانهيار الموحد لمنحنيات العبور، يقدم هذا العمل بروتوكولاً تنبؤياً لتطويع الحالة الهيكلية للمواد غير المتبلورة من خلال التدريب الميكانيكي الدوري.
أشار المؤلفون إلى أنه بينما تعد قيم الأسس المحددة (α≈0.17,β≈0.34) خاصة بظروف NPT أحادية المحور لنموذج Cu64Zr36، فإن إطار العمل القائم على التقلبات وعلاقة القياس β≈2α يمثلان اتصالاً جوهرياً بين تقلبات الإجهاد وديناميات الاسترخاء يستحق المزيد من التحقق عبر أوضاع التحميل المختلفة.