Thermal conductivity and tunable thermal anisotropy of magnetic CrSBr monolayer

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

تتناول الورقة البحثية التباينات الكبيرة في القيم المسجلة للموصلية الحرارية الشبكية ($\kappa$) لطبقة أحادية (ML) من CrSBr، وهي مادة مغناطيسية ثنائية الأبعاد من نوع فان دير فالس. قدمت الدراسات النظرية السابقة نتائج متفاوتة للغاية:

• Xuan et al. سجلوا قيمًا منخفضة للغاية ($\kappa_{xx} \approx 0.8$ W m$^{-1}$K$^{-1}$).
• Liu et al. سجلوا قيمًا متوسطة ($\kappa_{xx} \approx 25.5$ W m$^{-1}$K$^{-1}$).
• Han et al. سجلوا قيمًا عالية ($\kappa_{xx} \approx 80.8$ W m$^{-1}$K$^{-1}$).

استخدمت هذه الدراسات أكواد نظرية الكثافة الوظيفية (DFT) متشابهة، لكنها اختلفت في حل معادلة بولتزمان للنقل (BTE) وفي تضمين عمليات التشتت غير المرنة من الرتب العليا (مثل عمليات الفونون الرباعية). علاوة على ذلك، ظل الأصل الفيزيائي لتباين الخواص الحرارية القوي داخل المستوي ($\kappa_{xx} \neq \kappa_{yy}$) والقدرة على ضبط هذا التباين عبر المحفزات الخارجية (مثل الإجهاد أو الحجم) غير واضح.

2. المنهجية

استخدم المؤلفون حسابات المبادئ الأولى لحل هذه التناقضات وتحليل خصائص النقل:

• البنية الإلكترونية: أُجريت الحسابات باستخدام برنامج VASP مع دالة PBE-GGA، وتصحيحات فان دير فالس DFT-D3، و DFT+U ($U_{eff} = 4.0$ eV) لمعالجة إلكترونات Cr $d$ الموضعية.
• ثوابت القوة: تم حساب ثوابت القوة البينية من الدرجة الثانية والثالثة (IFC2 و IFC3) باستخدام طريقة الفروق المحدودة (Phonopy و Thirdorder.py) على خلية فائقة بأبعاد $6 \times 5 \times 1$.
  • خطوة حاسمة: تم فرض قواعد المجموع الدوراني باستخدام hiPhive لضمان التشتت التربيعي الصحيح لأنماط الصوتيات الالتوائية (ZA)، والتي تكون عرضة لعدم الاستقرار العددي في المواد ثنائية الأبعاد.
  • تم إهمال التصحيحات غير التحليلية لإعادة إنتاج كبت انقسام LO-TO بشكل صحيح في الأنظمة ثنائية الأبعاد.
• النقل الحراري: استُخدمت حزمة almaBTE لحل معادلة BTE الخطية بشكل تكراري، متجاوزةً تقريب زمن الاسترخاء (RTA). يسمح هذا بتضمين عمليات التشتت الطبيعي (Normal scattering)، وهي عمليات بالغة الأهمية في المواد ثنائية الأبعاد.
• التحليل المقارن: قارنت الدراسة بين الحلول الكاملة لـ BTE وتقريبات RTA، وبحثت في تأثيرات:
  • درجة الحرارة (100–150 كلفن، تحت درجة حرارة كوري).
  • الإجهاد (أحادي المحور وثنائي المحور حتى $\pm 5\%$).
  • حجم الشريحة ($L = \infty$، و $10$ $\mu$m، و $100$ nm) لمحاكاة تشتت الحدود.

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

أ. حل التناقضات في الموصلية الحرارية

يوضح المؤلفون أن التباين الكبير في الأدبيات السابقة ينبع إلى حد كبير من استخدام تقريب RTA.

• فشل RTA: يقلل RTA من تقدير الموصلية الحرارية بنسبة $\approx 25\%$ بالنسبة لـ $\kappa_{xx}$ ويشوه نسبة التباين بشكل كبير؛ حيث يعامل عمليات التشتت الطبيعي (Normal scattering) بشكل خاطئ كعمليات مقاومة.
• نتائج BTE الكاملة: من خلال حل BTE الكامل، حصل المؤلفون على $\kappa_{xx} \approx 86.3$ W m$^{-1}$K$^{-1}$ و $\kappa_{yy} \approx 43.1$ W m$^{-1}$K$^{-1}$ عند 150 كلفن. وتتوافق هذه القيم بشكل وثيق مع دراسة Han et al. (الذين ضمنوا تشتت الفونون الرباعي)، مما يشير إلى أن عمليات الفونون الرباعية ليست المصدر الأساسي للتشتت في CrSBr، بل المعالجة الدقيقة لعمليات الفونون الطبيعية الثلاثية هي السبب.

ب. أصل التباين الحراري

تحدد الدراسة أن التباين القوي ($\kappa_{xx}/\kappa_{yy} \approx 1.8$) ينشأ من تأثير مشترك لـ:

1. سرعات الفونونات: السرعات الجماعية على طول محور $x$ ($v_x$) أعلى منها على طول محور $y$ ($v_y$) عند الترددات المنخفضة ($\omega \lesssim 4$ THz).
2. أعمار الفونونات: تمتلك الفونونات المنتشرة على طول محور $x$ أعمارًا أطول بكثير ($\tau_x \gg \tau_y$) عند الترددات المنخفضة بسبب انخفاض معدلات التشتت.

ج. القابلية للضبط عبر تعديل الحجم

من النتائج الرئيسية أن التباين الحراري ليس ثابتًا جوهريًا، بل يمكن ضبطه عبر التحكم في حجم الشريحة.

• الآلية: يعمل تشتت الحدود على كبح الفونونات ذات المسار الحر المتوسط (MFP) الطويل. وبما أن $\kappa_{xx}$ يهيمن عليها الفونونات ذات الـ MFP الأطول من $\kappa_{yy}$، فإن تقليل حجم الشريحة يقلل من $\kappa_{xx}$ بشكل أكثر حدة.
• النتيجة: مع تناقص حجم الشريحة من اللانهاية إلى 100 نانومتر، تنخفض نسبة التباين من $\sim 1.8$ إلى $\sim 1.6$. وعند المقاييس الصغيرة جدًا ($L \sim 10^{-8}$ m)، يصبح النقل أقل تباينًا، مقتربًا من حد RTA حيث يتم كبح العمليات الطبيعية.

د. الاستقرار المغناطيسي تحت الإجهاد

بحث المؤلفون فيما إذا كان الإجهاد الميكانيكي يمكن أن يحفز انتقال طور من الحالة الأرضية الفيرومغناطيسية (FM) إلى الحالة الفيرومغناطيسية المضادة (AFM).

• النتيجة: تظل الحالة (FM) هي الحالة الأرضية عبر جميع أنواع الإجهاد أحادي وثنائي المحور ($\pm 5\%$). وتظل فرق الطاقة ($\Delta E_{FM-AFM}$) سالبة، مما يشير إلى أن التشوه الميكانيكي وحده لا يمكنه تغيير النظام المغناطيسي في طبقة CrSBr الأحادية ضمن هذا النطاق.

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

المعلم (Parameter)	القيمة / الملاحظة
الموصلية الحرارية (150 K)	$\kappa_{xx} \approx 86.3$ W m$^{-1}$K$^{-1}$؛ $\kappa_{yy} \approx 43.1$ W m$^{-1}$K$^{-1}$
نسبة التباين ($\kappa_{xx}/\kappa_{yy}$)	1.8 (BTE كامل) مقابل 1.5 (RTA)
خطأ RTA	يقلل من تقدير $\kappa_{xx}$ بنسبة $\sim 25\%$؛ ويفشل في التقاط التباين الصحيح.
الحالة الأرضية المغناطيسية	الفيرومغناطيسية (FM) مستقرة؛ و AFM أعلى طاقة بمقدار 49 meV/f.u.
تأثير الإجهاد	لم يُلاحظ انتقال من FM $\to$ AFM حتى إجهاد $\pm 5\%$
تأثير الحجم	تقليل $L$ من $\infty$ إلى 100 nm يقلل $\kappa$ بنحو $\sim 7\times$ ويخفض نسبة التباين.

5. الأهمية

• معيار منهجي: تسلط الورقة الضوء على الضرورة القصوى لحل معادلة BTE الكاملة (بما في ذلك التشتت الطبيعي) بدلاً من الاعتماد على RTA للحصول على تنبؤات دقيقة للنقل الحراري في المواد المغناطيسية ثنائية الأبعاد.
• تصميم المواد: تثبت أن CrSBr موصل حراري عالي التباين، وهو أمر مفيد لتطبيقات الإدارة الحرارية حيث يتطلب الأمر تدفقًا حراريًا اتجاهيًا.
• أداة تحكم: إن اكتشاف أن الحصر الهندسي (حجم الشريحة) يمكنه ضبط التباين الحراري يقدم درجة جديدة من الحرية لهندسة الخصائص الحرارية في الأجهزة ثنائية الأبعاد دون الحاجة إلى تطعيم كيميائي أو مجالات خارجية معقدة.
• المتانة المغناطيسية: إن تأكيد استقرار الحالة (FM) ضد الإجهاد يسهل تصميم المستشعرات القائمة على الإجهاد، حيث لن يتم تغيير النظام المغناطيسي عن غير قصد بسبب الإجهاد الميكانيكي.