Lattice thermal conductivity decomposition: Peierls vs. non-Peierls contributions

ملخص تقني: تفكيك الموصلية الحرارية للشبكة: مساهمات بيريلز وغير بيريلز

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

المنهجية
تستخدم الدراسة محاكاة الديناميكا الجزيئية الكلاسيكية (MD) مدمجة مع صيغة "جرين-كوبو" لحساب موتر الموصلية الحرارية ($\kappa_{\alpha\beta}$). يركز التحليل على ثلاثة أنظمة بلورية متميزة:

1. الأرجون الصلب (SA): شبكة مكعبة مركزية الوجه (fcc) أحادية الذرة تم نمذجتها باستخدام جهد "لينارد-جونز".
2. الأرجون الصلب بكتل متبادلة (SAAM): نموذج يستخدم نفس جهد "لينارد-جونز" ولكن مع كتل ذرية متبادلة (نسبة 10) لإنشاء شبكة رباعية الأضلاع تحتوي على فروع صوتية وبصرية.
3. $\alpha$-الكوارتز: بلورة ثلاثية الأضلاع تحتوي على تسعة ذرات في وحدة الخلية، تم نمذجتها باستخدام جهد "بيست-كرايمر-فان سانتن" (BKS).

لكل نظام، يتم حساب الموصلية الحرارية باستخدام ثلاثة تعريفات لتيار الحرارة:

• تيار الحرارة الكامل ($J_{full}$): التعبير الكلاسيكي الكامل الذي يتضمن جميع الحدود.
• تيار الحرارة التربيعي ($J_{quad}$): مفكوك تايلور من الدرجة الثانية لتيار الحرارة بالنسبة للإزاحات والسرعات الذرية.
• تيار حرارة بيريلز ($J_{Peierls}$): مجموعة فرعية من التيار التربيعي يتم الحصول عليها عن طريق التحويل إلى إحداثيات الأنماط الطبيعية والاحتفاظ فقط بالحدود القطرية في مؤشر فرع الفونون.

بالإضافة إلى ذلك، يتم تقييم تقريب زمن الاسترخاء (RTA) عن طريق تطبيق نهج (GK) على تقلبات طاقات الأنماط الطبيعية الفردية. أُجريت عمليات المحاكاة عند درجات حرارة محددة (20 كلفن لـ SA/SAAM، و300 كلفن لـ $\alpha$-الكوارتز) وعند ضغط صفري باستخدام حزمة (LAMMPS). وتم تحديد عدم الخطية في الأنظمة باستخدام مقياس نسبي ($\eta$) يقارن متوسط طاقة الوضع المعدلة بالحد التوافقي.

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

• مقارنة تعريفات تيار الحرارة: بالنسبة لجميع الأنظمة الثلاثة، تختلف الموصلية الحرارية المحسوبة باستخدام تيارات الحرارة التربيعية وبيريلز اختلافاً طفيفاً فقط عن بعضها البعض. ومع ذلك، فإن الموصلية الحرارية المحسوبة باستخدام تيار الحرارة الكامل تكون أعلى باستمرار من تقريبات التربيعي/بيريلز.
  • في SA وSAAM، تكون قيم التربيعي/بيريلز أقل بنسبة 10-13% تقريباً من قيم التيار الكامل.
  • في $\alpha$-الكوارتز، يكون الفرق بين التيار الكامل والتربيعي/بيريلز أصغر بكثير، ويقع ضمن حدود عدم اليقين العددي.
• المساهمات الصوتية مقابل البصرية:
  • في نموذج SAAM، تهيمن الأنماط الصوتية على الموصلية الحرارية، بينما تكون المساهمات البصرية مهملة (خاصة على طول المحور $c$).
  • في $\alpha$-الكوارتز، تتجاوز مساهمة الفونونات البصرية في موصلية بيريلز الحرارية تلك الخاصة بالأنماط الصوتية على طول المحور $a$ (البصري $\approx$ 6.7 W m$^{-1}$K$^{-1}$ مقابل الصوتي $\approx$ 2.1 W m$^{-1}$K$^{-1}$). وعلى طول المحور $c$، تكون المساهمات البصرية والصوتية متقاربة.
• تقريب زمن الاسترخاء (RTA): يقلل (RTA) بشكل منهجي من تقدير الموصلية الحرارية في جميع الأنظمة الثلاثة مقارنة بحسابات (GK) الكاملة. بالنسبة لـ SA، يعطي (RTA) قيماً أقل بنسبة ~15% من نتيجة بيريلز.
• عدم الخطية: إن قيمة عدم الخطية المقدرة ($\eta$) صغيرة لجميع الأنظمة (1% لـ SA/SAAM و1.6% لـ $\alpha$-الكوارتز)، ومع ذلك فهي كافية لإنتاج موصلية حرارية محدودة وذات حجم ملموس.

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

ومن النتائج الهامة حدوث عكس في آلية نقل الحرارة المهيمنة في $\alpha$-الكوارتز مقارنة بنموذج SAAM: حيث تساهم الأنماط البصرية أكثر من الأنماط الصوتية في $\alpha$-الكوارتز، بينما يحدث العكس في نموذج SAAM. ويعزو المؤلفون ذلك إلى العدد الأكبر من الفروع البصرية في $\alpha$-الكوارتز (24 فرعاً) مقارنة بـ SAAM (3 فروع).

علاوة على ذلك، تسلط الدراسة الضوء على قصور في الأدبيات السابقة (تحديداً المرجع [5] بخصوص $\alpha$-الكوارتز) التي استخدمت ملاءمة تجريبية (heuristic fitting) لـ (HCCF) لتقسيم المساهمات الصوتية والبصرية. ويجادل المؤلفون بأن نهج الملاءمة هذا قد يؤدي إلى تقديرات غير صحيحة لأنه يحاول ملاءمة دوال ذات تكاملات زمنية محدودة مع بيانات تحتوي على مكونات تذبذبية تتلاشى تكاملاتها الزمنية (وبالتالي لا تساهم في الموصلية الحرارية). ومن خلال التقييم المباشر لمساهمات الأنماط دون استخدام ملاءمة ارتجالية، تشير هذه الدراسة إلى أن المساهمات البصرية في $\alpha$-الكوارتز أكبر بكثير مما تم الإبلاغ عنه سابقاً.

أخيراً، يشير المؤلفون إلى أنه بينما تظل التقريبات صالحة لأنظمة خلايا الوحدة الصغيرة المدروسة هنا، فإن نتائجهم الأخيرة على الأنظمة ذات خلايا الوحدة الكبيرة (مثل $\beta$-HMX و $\alpha$-RDX) تظهر أن موصلية بيريلز الحرارية يمكن أن تكون أصغر بكثير من موصلية التيار الكامل، مما يشير إلى أن صلاحية تقريب بيريلز تعتمد على النظام.