Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

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

تنشأ الكهرباء الحديدية (ferroelectricity) في الهياكل المتغايرة من نوع فانديرفالز (vdW)، مثل نيتريد البورون السداسي (hBN) ثنائي الطبقات، نتيجة لإعادة البناء الهيكلي للمجالات المثلثة (تراصفات AB و BA) التي تفصلها انخلاعات الواجهة. وبينما تم توصيف "الكهرباء الحديدية الانزلاقية" جيداً في الطبقات ثنائية التراكب ذات الالتواء الصغير (زاوية الالتواء $\theta < 2^\circ$) وحالات الإجهاد الصغيرة، لا تزال هناك فجوات كبيرة:

• النطاق غير المستكشف: من غير المعروف ما إذا كانت الكهرباء الحديدية تستمر تحت تأثير التشوهات المتغايرة الكبيرة (مثل زوايا الالتواء الكبيرة أو الإجهاد الكبير) والتشوه المتغاير الصغير بالنسبة لتكوينات التماثل العالي بخلاف تراصف AA القياسي.
• القيود الحسابية: تُعد المحاكاة الذرية (باستخدام الجهود بين الذرية) مكلفة حسابياً للهياكل متغايرة التراكب ذات وحدة الخلية الكبيرة. علاوة على ذلك، تفشل الجهود بين الذرية الموجودة في التنبؤ بدقة بإعادة البناء الهيكلي والاستقطاب في تكوينات الالتواء الكبيرة (تحديداً تكوين $\Sigma7$ عند $\approx 21.79^\circ$) تحت ضغط عمودي شديد.
• عدم الدقة المحتملة: غالباً ما تعطي الجهود القياسية قيم استقطاب ومقادقات خاطئة مقارنة بنظرية الدالة الكثافية (DFT)، مما يؤدي إلى تباينات في التنبؤ بتبديل المجالات تحت تأثير المجالات الكهربائية.

2. المنهجية

استخدم المؤلفون نهجاً متعدد الأوجه يجمع بين المحاكاة الذرية، ونظرية الدالة الكثافية (DFT)، وإطار عمل مستمر مبتكر:

• المحاكاة الذرية (LAMPS):

  • استُخدمت جهود "تيرسوف" المعدلة للروابط داخل الطبقة (B-N) وجهود الطبقة البينية المعتمدة على التسجيل (ILP) لتفاعلات فانديرفالز.
  • تمت محاكاة نظام hBN ثنائي الطبقات ذي الالتواء الصغير ($0.299^\circ$) والإجهاد الصغير ($0.42\%$ إجهاد ثنائي المحور متساوي) لتوصيف انخلاعات الواجهة والتبديل الكهربائي الحديدي.
  • تم تحديد أن الجهود الموجودة تفشل في إعادة إنتاج نتائج DFT لتكوين $\Sigma7$ ($21.79^\circ$) الملتوي، خاصة تحت الضغط العمودي.

• نظرية الدالة الكثافية (DFT):

  • حُسبت طاقة التراصف العامة (GSFE) ومشاهد الاستقطاب (PL) لكل من تكوينات $0^\circ$ (AA) و $\Sigma7$ ($21.786789^\circ$).
  • أثبتت أن تكوين $\Sigma7$، تحت ضغط عمودي بنسبة 28%، يُظهر نهايات صغرى طاقية متحللة واستقطاباً متعاقباً في الاتجاه العمودي، مما يستوفي شروط الكهرباء الحديدية.

• نموذج متعدد المقاييس ثابت الإطار والمعتمد على البلورات الثنائية (BFIM):

  • تم تطوير إطار عمل مستمر لسد الفجوة بين دقة DFT والكفاءة الحسابية.
  • الحركية (Kinematics): يعامل الطبقة الثنائية كصفائح مستمرة. يستخدم بلورات ثنائية تعتمد على الصيغة الطبيعية (Smith Normal Form) لتحديد التكوينات الطبيعية والمرجعية والمشوهة.
  • القانون التكويني: تشمل الطاقة الكلية طاقة الإجهاد المرن، وطاقة فانديرفالز البينية (المشتقة من GSFE)، وشغل المجال الكهربائي (المشتق من PL).
  • المعادلات الحاكمة: يحل النموذج الاسترخاء الهيكلي عبر تقليل تدفق التدرج لتابع الطاقة الكلية، مما يلتقط حركة انخلاعات الواجهة تحت المجالات الكهربائية.

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

• اكتشاف الكهرباء الحديدية ذات الالتواء الكبير: تثبت هذه الدراسة أن الكهرباء الحديدية لا تقتصر على الالتواءات الصغيرة، بل تستمر في تكوينات hBN ثنائية الطبقات القريبة من تراصف $\Sigma7$ ($21.79^\circ$)، بشرط خضوع النظام لضغط عمودي محدد.
• توصيف الانخلاعات الدوامية:
  • في أنظمة الالتواء الصغير، تكون الانخلاعات بشكل أساسي من النوع اللولبي وتشكل شبكات مثلثية مستقيمة.
  • في أنظمة الإجهاد الصغير (ثنائي المحور المتساوي)، تشكل الانخلاعات شبكة مثلثية دوامية (حلزونية) نتيجة تحول انخلاعات الحافة بالقرب من وصلات AA لتفريغ طاقة الإجهاد.
  • في أنظمة الالتواء الكبير ($\Sigma7$)، تمتلك انخلاعات الواجهة متجه بورغرز أصغر بكثير ($\approx 0.55$ Å) مقارًة بأنظمة الالتواء الصغير، ومع ذلك لا تزال تشكل شبكة مثلثية.
• تطوير نموذج BFIM: أنشأ المؤلفون نموذجاً مستمراً مستنداً إلى DFT ينجح في التنبؤ بإعادة البناء الهيكلي والتبديل الكهربائي الحديدي في الهياكل ذات وحدة الخلية الكبيرة حيث تكون المحاكاة الذرية باهظة التكلفة أو غير دقيقة.
• تحديد القيود المحتملة: تسلط الورقة الضوء على أن الجهود بين الذرية القياسية (ILP) غير موثوقة لـ hBN ذي الالتواء الكبير تحت الضغط، مما يستلهّم الحاجة إلى مناهج متعددة المقاييس مستندة إلى DFT.

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

• التشوهات المتغايرة الصغيرة:
  • أكدت المحاكاة الذرية أن تطبيق مجال كهربائي عمودي يؤدي إلى توسع أو انكماش مجالات AB و BA، مدفوعاً بانحناء انخلاعات الواجهة.
  • في الأنظمة المجهدة، تؤدي الشبكة الانخلاعية الدوامية إلى تغيير أكثر وضوحاً في مساحة المجالات مقارنة بالأنظمة الملتوية.
• التشوهات المتغايرة الكبيرة ($\Sigma7$):
  • كشفت DFT أن تكوين $\Sigma7$ يمتلك نهايتين صغريتين متحللتين بكثافات استقطاب متضادة ($\pm 1.35 \times 10^{-4}$ C/m$^2$) تحت ضغط بنسبة 28%.
  • نجح نموذج BFIM في إعادة إنتاج تشكيل شبكة انخلاع مثلثية في التكوين القريب من $\Sigma7$ ($21.96^\circ$ التواء).
  • تنبأ النموذج بأن هذه الأنظمة ذات الالتواء الكبير تُظهر تبديلاً كهربائياً حديدياً، رغم أن نسبة تغير المساحة ($A_{AB}/A_{BA}$) أصغر بنحو 4 مرات من أنظمة الالتواء الصغير بسبب ضعف مشهد الاستقطاب.
• التحقق من النموذج:
  • تم ضبط معاملات نموذج BFIM باستخدام بيانات تجريبية (Lv et al.) لـ hBN ذي الالتواء الصغير، مما أظهر توافقاً ممتازاً في نسب مساحة المجالات تحت مجالات كهربائية متغيرة.
  • تنبأ النموذج بشكل صحيح بالاعتماد غير الخطي لعزم ثنائي القطب الصافي على المجال الكهربائي المطبق.

5. الأهمية

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