Interplay of interlayer distance and in-plane lattice relaxations in encapsulated twisted bilayers

ملخص تقني: التداخل بين مسافة الطبقة البينية والاسترخاء الشبكي داخل المستوي في الطبقات الثنائية الملتوية المغلفة

بيان المشكلة
تُظهر الهياكل غير المتجانسة لـ "فان دير فالس" الملتوية، ولا سيما ثنائيات الطبقات من ثنائيات كالكوجينيدات المعادن الانتقالية (TMDs)، شبكات "مواريه" (moiré) ميزوسكوبية تغير الخصائص الإلكترونية والإلكترونية البصرية بشكل كبير. يؤدي تشكل هذه الشبكات إلى دفع الاسترخاء الذري لتقليل الطاقة الكلية، مما يحول أنماط "مواريه" الصلبة إلى مصفوفات من النطاقات المتوافقة التي تفصل بينها جدران نطاق ممتصة للإجهاد. وبينما تم توصيف العبور من أنماط "مواريه" الصلبة إلى هياكل النطاقات المسترخية للطبقات الثنائية المعلقة (مع زوايا عبور $\theta^* \approx 2.5^\circ$ للحالات المتوازية و $\approx 1^\circ$ للحالات غير المتوافقة في المحاذاة)، تشير التجارب الحديثة إلى أن تغليف طبقات الـ TMD في نتريد البورون السداسي (hBN) يعزز الاسترخاء الشبكي ويزيح زوايا العبور هذه إلى قيم أعلى. إن النماذج النظرية الحالية للاسترخاء الشبكي لا تأخذ في الحسبان القيود الميكانيكية التي تفرضها طبقات التغليف، وتحديداً تعديل المسافات بين الطبقات عند الواجهة.

المنهجية
يطور المؤلفون نموذجاً نظرياً يصف الاسترخاء الشبكي في طبقات الـ TMD الملتوية بحيث يدمج صراحةً تأثيرات التغليف. يتضمن النهج ما يلي:

1. صياغة طاقة الالتصاق: استخدام كثافة طاقة التصاق تعتمد على التراص (stacking-dependent)، $W_{ad}(r_0, d)$، والتي تعتمد على الإزاحة داخل المستوي $r_0$ والمسافة بين الطبقات $d$.
2. نمذجة التغليف: معاملة الطبقة الثنائية الملتوية كأنها مغروسة بين صفائح متعددة الطبقات (التغليف). يتم نمذجة التفاعل بين طبقة الـ TMD الثنائية وصفائح التغليف عبر جهد توافقي يتميز بمعامل فعال $k$، والذي يمثل صلابة الواجهة بين الطبقة الثنائية والصفائح.
3. تعديل المسافة بين الطبقات: اشتقاق تعبير جديد لتعديل المسافة بين الطبقات، $d_{enc}(r_0)$، من خلال موازنة الضغط المعتمد على التراص من الطبقات الملتوية ($P_{in}$) مقابل ضغط الاستعادة من صفائح التغليف ($P_{out}$). ينتج عن ذلك كبح لتغير المسافة بين الطبقات مقارنة بالطبقات الثنائية المعلقة.
4. حساب الاسترخاء داخل المستوي: تقليل دالة الطاقة الكلية التي تضم الطاقة المرنة داخل المستوي وطاقة الالتصاق المحلية (باستخدام $d_{enc}$ المعدلة) لتحديد مجالات الإزاحة داخل المستوي. تستخدم الحسابات توسيعات متسلسلة فوريه عبر ناقلات التردد لشبكة "مواريه" الفائقة.
5. القياس الكمي: تحليل قوة الاسترخاء الشبكي باستخدام النسبة $R_\nu(\theta)$، وهي مساحة المناطق ذات التراص من نوع XX بالنسبة لمساحة خلية "مواريه" الفائقة، عبر زوايا التواء مختلفة ومعاملات صلابة الواجهة.

المساهمات والنتائج الرئيسية

• تطوير النموذج: يقدم البحث معاملًا ظاهريًا واحدًا، $k$، يميز صلابة الواجهة بين الطبقة الثنائية الملتوية ومادة التغليف. يتحكم هذا المعامل في كبح تعديل المسافة بين الطبقات.
• تأثير التغليف: يوضح النموذج أن التغليف يكبح نطاق المسافات المتاحة بين الطبقات. هذا الكبح يزيد من فرق الطاقة بين التراصات المفضلة طاقيًا (مثل التراص المروحي 3R) والتراصات غير المفضلة (XX).
• إزاحة زاوية العبور: النتيجة الأساسية هي أن زيادة صلابة الواجهة ($k$) ترفع زاوية العبور ($\theta^*$) بين نظامي الاسترخاء الضعيف والقوي. بالنسبة للطبقات الثنائية المعلقة ($k=0$)، تُحسب $\theta^*$ بنحو $\approx 3.8^\circ$. أما بالنسبة للواجهات الصلبة تماماً ($k \gg 4\varepsilon$)، فتزدัง $\theta^*$ لتصل إلى $\approx 4.5^\circ$.
• التوافق الكمي: من خلال ضبط معامل الصلابة لنمذجة واجهة صلبة تماماً (والتي يُرجح أنها ناتجة عن عدم التوافق الشبكي بين الـ TMD والـ hBN)، يتنبأ النموذج بإزاحة في زاوية العبور قدرها $\Delta\theta^* \approx 0.75^\circ$. وهذا يتوافق جيداً مع الملاحظات التجريبية لـ $\Delta\theta^* \approx 1^\circ$ لطبقات WSe$_2$ المتوازية المغلفة بـ hBN.
• سلوك النطاق: تؤكد الدراسة أنه في نظام الاسترخاء القوي (زوايا الالتواء الصغيرة)، تصبح مساحة مناطق التراص XX مستقلة عن زاوية الالتواء، متبعةً مقياس $\theta^2$ لمعامل قوة الاسترخاء، وهو سلوك يظل عالمياً عبر مختلف مستويات قوة التغليف.

الأهمية
يدعي البحث أن أهميته تكمن في توفير إطار نظري يسد الفجوة بين نماذج الطبقات الثنائية المعلقة والنتائج التجريبية للأنظمة المغلفة. ومن خلال تحديد أن صلابة واجهة التغليف تكبح تعديل المسافة بين الطبقات، يشرح المؤلفون سبب إظهار الأنظمة المغلفة استرخاءً شبكيًا معززًا وزوايا عبور أعلى. ينجح النموذج في تحقيق توافق كمي مع البيانات التجريبية الموجودة من خلال معاملة الواجهة كواجهة صلبة فعلياً، مما يشير إلى أن عدم التوافق الشبكي عند واجهة الـ TMD/hBN هو عامل حاسم في تحديد الاستجابة الميكانيكية لهذه الهياكل غير المتجانسة. ويشير المؤلفون إلى إمكانية توسيع هذا الإطار ليشمل أنظمة "توستروينك" (twistronic) أخرى، مثل أنظمة WX$_2$/MoX$_2$، للتنبؤ بالتغيرات في مقدار الإجهاد وسلوك شبكة جدران النطاق.