Spin-Orbit Coupling Effects on the Structural and Electronic Properties of Planar Pentagonal p-MS$_{2}$ (M = Si, Ge, and Pb)

ملخص تقني: تأثيرات اقتران المدار-اللف المغزلي على بنية p-MS₂ خماسية السداسية المستوية (M = Si, Ge, Pb)

بيان المشكلة
تقدم المواد ثنائية الأبعاد (2D) مساراً لهندسة البنى الإلكترونية من خلال الحصر الكمي وكسر التماثل. وبينما تمت دراسة عائلات مثل ثنائي كالكوجينيدات المعادن الانتقالية (MX₂) والبنى الخماسية المستوية (مثل penta-NiN₂) بشكل مكثف، فإن الفئة المحددة من كالكوجينيدات المجموعة IVA ذات البنية الخماسية المستوية (p-MX₂) لا تزال غير مستكشفة بشكل كافٍ فيما يتعلق بالتداخل بين الاستقرار الهيكلي واقتران المدار-اللف المغزلي (SOC). ركزت الدراسات السابقة على الخصائص الطوبولوجية أو التغيرات الإلكترونية الناجمة عن SOC بمعزل عن غيرها. ومع ذلك، لم يتم تقديم تحقيق منهجي حول كيفية تأثير SOC على الاستقرار الهيكلي، والمعاملات الهندسية، وإعادة التشكيل الإلكتروني التفصيلية (تحديداً عبر التفكيك المداري الموزع حسب $j$) لأنظمة p-MS₂ (حيث M = Si, Ge, Pb).

المنهجية
تستخدم الدراسة حسابات نظرية الدوال الوظيفية للكثافة (DFT) باستخدام حزمة Quantum ESPRESSO. وتشمل الميزات المنهجية الرئيسية ما يلي:

• الجهود الكاذبة (Pseudopotentials): استُخدمت جهود كاذبة من نوع SSSP Efficiency PBEsol v1.3 لعناصر Si، Ge، Pb، وS.
• الاسترخاء والتقارب: أُجريت عمليات تحسين البنية تحت ضغط قدره 0.5 كيلوبار، مع تقارب القوة عند $5 \times 10^{-4}$ Ry·Bohr$^{-1}$، وتقارب الطاقة عند $3.0 \times 10^{-5}$ Ry. تم تضمين SOC في كل من استرخاء البنية وحسابات الخصائص الإلكترونية للحالة الأرضية.
• تحليل الاستقرار: تم تقييم الاستقرار الحراري عبر محاكاة الديناميكا الجزيئية (MD) عند 300 كلفن باستخدام حاسبة FALCON القائمة على التعلم النشط (انحدار العملية الغاوسية) ضمن مجموعة NVT لمدة 10 بيكوثانية. كما حُسبت طاقات التماسك ($E_c$) لمقارنة الاستقرار الطاقي مقابل المواد ثنائية الأبعاد النظرية والتجريبية المعروفة.
• التحليل الإلكتروني: استُخدمت هياكل النطاق (BS)، وكثافة الحالات (DOS)، وكثافة الحالات المسقطة (PDOS)، والتفكيكات المدارية الموزعة حسب $j$ (الزخم الزاوي الكلي) لتحليل الحالات الإلكترونية بالقرب من مستوى فيرمي.

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

1. الاستقرار الهيكلي والهندسة:

   • عدم استقرار p-SiS₂: تخلص الدراسة إلى أن بنية p-SiS₂ غير مستقرة على الأرجى. يكشف تحليل فرق كثافة الإلكترونات (EDD) عن إزاحة غير متماثلة في كثافة الإلكترونات المشتركة على طول روابط S–S، وتظهر محاكاة الديناميكا الجزيئية (MD) إعادة بناء هيكلية سريعة وتقلبات في الطاقة خلال أول 3 بيكوثانية عند 300 كلفن.
   • استقرار p-GeS₂ و p-PbS₂: يحافظ كل من p-GeS₂ و p-PbS₂ على سلامتهما الخماسية المستوية تحت محاكاة الديناميكا الجزيئية. وتقع طاقات تماسكها (4.14 إلكترون فولت/ذرة و3.94 إلكترون فولت/ذرة بدون SOC، على التوالي) ضمن نطاق المواد المستقرة المعروفة في النماذج النظرية والمواد التجريبية ثنائية الأبعاد.
   • الانكماش الناجم عن SOC: يتسبب تضمين SOC في انكماش هيكلي طفيف؛ حيث يقل ثابت الشبكة ($a$) وأطوال روابط M–S ($d_1$) (بنسبة تتراوح بين 0.1% إلى 0.26%)، بينما تظل أطوال روابط S–S ($d_2$) دون تغيير تقريباً. هذا الانكماش يكون أكثر وضوحاً مع الذرات الأثقل (Pb > Ge).

2. الخصائص الإلكترونية والتحول من معدن إلى شبه موصل:

   • p-GeS₂: يظل معدنياً حتى مع وجود SOC. تأثير SOC على بنية النطاق و DOS ضئيل، حيث تهيمن مدارات Ge-$p$ و S-$p$ على الحالات الإلكترونية بالقرب من مستوى فيرمي دون إعادة تشكيل كبيرة.
   • p-PbS₂: يُظهر تحولاً دراماتيكياً من معدن إلى شبه موصل مدفوعاً بـ SOC. يفتح SOC فجوة نطاق شبه مباشرة تبلغ حوالي 0.475 إلكترون فولت. ويُعزى هذا التحول إلى الـ SOC القوي لذرة Pb الثقيلة، والذي يعدل تشتت النطاق بشكل كبير ويعزز تموضع الحالات الإلكترونية بالقرب من مستوى فيرمي.

3. التفكيك المداري وإعادة بناء الترابط:

   • التحليل الموزع حسب $j$: يكشف تحليل مداري مفصل موزّع حسب $j$ (الزخم الزاوي الكلي) أن SOC يرفع درجة التحلل للمدارات $p$ ($l \ge 1$). في p-PbS₂، يحدث تهجين قوي بين حالات S-$p$ و Pb-$p$. وتحديداً، عند قمة حزمة التكافؤ (VBM)، تساهم حالة Pb-$p$ ($j=0.5$) بشكل كبير، بينما عند قاع حزمة التوصيل (CBM)، تهيمن خليط من حالات S-$p$ و Pb-$p$ ($j=1.5$).
   • طابع الترابط: يشير تعزيز التموضع وخلط المدارات إلى إعادة بناء في خصائص الترابط، مما يؤدي إلى تقليل التداخل المداري وانخفاض طفيف في الاستقرار الطاقي (تنخفض طاقة التماسك إلى حوالي 86.55% من قيمة غير SOC لـ p-PbS₂).

4. تباين الخواص في p-PbS₂:

   • يُظهر التوزيع المكاني للحالات الإلكترونية في p-PbS₂ تباين خواص ملحوظاً على طول روابط S–S. تتوزع حالات قمة حزمة التكافؤ (VBM) بشكل أساسي على طول اتجاه بلوري واحد ($\vec{i}$)، بينما تُظهر حالات قاع حزمة التوصيل (CBM) مساهمة معززة على طول الاتجاه العمودي ($\vec{ii}$) مع طابع ترابط $\pi$ متميز. يشير هذا إلى أن روابط S–S ليست متكافئة إلكترونياً.

الأهمية والادعاءات
يزعم البحث أن هذه النتائج توفر رؤية نقدية لإعادة التشكيل الهيكلي والإلكتروني المدفوع بـ SOC في الكالكوجينيدات خماسية السداسية المستوية. وتؤكد الدراسة أن:

• SOC هو عامل حاسم في استقرار ووجود طبقات p-MS₂ أحادية الطبقة، مما قد يجعل نظائرها الأخف (مثل p-SiS₂) غير مستقرة بينما يعمل على استقرار النظائر الأثقل من خلال التعديلات الهيكلية.
• بالنسبة للعناصر الثقيلة مثل Pb، فإن SOC ليس مجرد اضطراب، بل هو آلية تقود تحولاً جوهرياً من معدن إلى شبه موصل وتغير طبيعة الترابط الكيميائي.
• التباين الإلكتروني الفريد والحالات الإلكترونية المحددة لـ p-PbS₂ تجعل منه مرشحاً واعداً لتطبيقات استشعار الغاز، مع احتمال وجود المواقع النشطة بالقرب من ذرات Pb.

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