Pressure-Driven Evolution of the Physical Properties of Layered Y2GaB and Y2InB MAB Phases: A Comparative First-Principles Study

ملخص تقني: التطور المدفوع بالضغط للخصائص الفيزيائية لأطوار MAB من نوع Y2GaB وY2InB الطبقية

بيان المشكلة
تُظهر البوريدات الثلاثية الطبقية المعروفة باسم أطوار MAB مزيجًا فريدًا من السلوك الإلكتروني المعدني والصلابة الشبيهة بالسيراميك نظرًا لبنيتها ذات الطبقات النانوية، حيث تتواجد روابط معدنية-بورون قوية داخل المستويات مع تفاعلات بين الطبقات أضعف. وبينما شخصت الدراسات السابقة مركبات MAB محددة عند الضغط الجوي، إلا أن هناك نقصًا في الفهم الموحد حول كيفية تأثير الضغط الهيدروستاتيكي واستبدال الموقع A (تحديدًا Ga مقابل In) على الاستقرار البنيوي، والتباين الميكانيكي، والطوبولوجيا الإلكترونية، والخصائص الحرارية الفيزيائية. علاوة على ذلك، أفادت الأعمال الحسابية السابقة حول Y2GaB وY2InB بوجود أنماط فونون ذات ترددات سالبة، مما يثير تساؤلات حول استقرارها الديناميكي ومدى حساسية هذه التقييمات للمعاملات الحسابية (معالجة وحدة الخلية مقابل الخلية الفائقة). تعالج هذه الدراسة هذه الفجوات من خلال المقارنة المنهجية للتطور المعتمد على الضغط لمركبي Y2GaB وY2InB من 0 إلى 20 جيجا باسكال لتحديد نوافذ ضغط صالحة ميكانيكيًا وتوضيح العلاقات بين البنية والخصائص.

المنهجية
تستخدم الدراسة حسابات المبادئ الأولية القائمة على نظرية دالة الكثافة (DFT) باستخدام كود CASTEP مع تقريب التدرج المعمم (GGA) من نوع Perdew-Burke-Ernzerhof (PBE). وتشمل السمات المنهجية الرئيسية ما يلي:

• تحسين البنية: تم تحسين البنى السداسية من نوع $P6_3/mmc$ تحت ضغوط هيدروستاتيكية (0، 5، 10، 15، و20 جيجا باسكال) باستخدام خوارزمية BFGS.
• الاستقرار المرن: تم اشتقاق الثوابت المرنة عبر طريقة الإجهاد والانفعال المحدودة. وتم تقييم الاستقرار باستخدام معايير مصححة بالضغط للتناظر السداسي، مع التمييز بين استقرار الانفعال المتجانس واستقرار الفونون.
• اشتقاق الخصائص: تم حساب المعاملات متعددة البلورات (متوسطات Voigt-Reuss-Hill)، والسرعات الصوتية، ودرجات حرارة ديباي، وتقديرات الموصلية الحرارية (نماذج Clarke وSlack) للحالات المستقرة مرنًا.
• التحليل الإلكتروني والبصري: تم حساب بنى النطاقات، وكثافات الحالات (DOS)، وأسطح فيرمي مع وبدون اقتران المدار المغزلي (SOC) باستخدام Quantum ESPRESSO. كما تم تحليل الخصائص البصرية بالنسبة لاستقطابات المجال الكهربائي على طول [100] و[001].
• الترابط والديناميكا: استُخدمت خرائط فرق كثافة الإلكترونات، وتحليل شحنة Mulliken/Hirshfeld، وتركيزات تداخل الروابط لتوصيف الترابط. كما تم حساب تشتت الفونون باستخدام طريقة الإزاحة المحدودة لكل من حسابات وحدة الخلية والخلية الفائقة الموسعة لتقييم الاستقرار الديناميكي والحساسية العددية.

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

• التطور البنيوي: يحتفظ كلا المركبين ببنية $P6_3/mmc$ تحت الضغط. يعد الانضغاط متباين الخواص، حيث ينكمش معامل الشبكة القاعدي ($a$) بشكل أكبر من المعامل المحوري ($c$). وبناءً على ذلك، تزدل نسبة $c/a$ لكليهما، وينخفض حجم وحدة الخلية بشكل رتيب (حوالي 21.8% لـ Y2GaB و20.3% لـ Y2InB خلال 0–20 جيجا باسكال).
• الاستقرار المرن والخصائص الميكانيكية:
  • نوافذ الاستقرار: يستوفي Y2GaB معايير الاستقرار المرني المصححة بالضغط عند 0 و5 و10 جيجا باسكال، ولكنه يصبح غير مستقر عند 15 و20 جيجا باسكال. أما Y2InB فهو مستقر مرنًا فقط عند الضغط الجوي (0 جيجا باسكال).
  • الاستجابة الميكانيكية: ضمن نطاقه المستقر، يُظهر Y2GaB زيادة في معاملات الحجم، والقص، ويونغ مع الضغط. ويُتوقع أن يكون كلا المركبين مطيليين (نسبة بوج $G/B < 0.57$). يُظهر Y2GaB تباينًا مرنًا عاليًا عند 0 جيجا باسكال ($A^U = 5.25$)، والذي يتناقص مع الضغط، بينما يعد Y2InB أقل تباينًا عند الضغط الجوي.
  • الخصائص الحرارية الفيزيائية: بالنسبة لـ Y2GaB المستقر، يؤدي الضغط إلى زيادة سرعات الصوت، ودرجة حرارة ديباي، ودرجات حرارة الانصهار المقدرة. ترتفع الموصلية الحرارية الشبكية المقدرة بنموذج Slack عند 300 كلفن من 0.605 واط م⁻¹ كلفن⁻¹ (0 جيجا باسكال) إلى 1.543 واط م⁻¹ كلفن⁻¹ (10 جيجا باسكال). يُظهر كلا المركبين موصلية حرارية دنيا منخفضة جدًا، مما يشير إلى إمكانات لاستخدامها في تطبيقات الحواجز الحرارية عند درجات الحرارة المنخفضة، وهو ما يحده انخفاض درجات انصهارهما نسبيًا.
• البنية الإلكترونية: كلا المادتين عبارة عن معادن متعددة النطاقات. يمتلك Y2GaB تسعة أوراق لسطح فيرمي، ويمتلك Y2InB ثمانية عند الضغط الجوي. يعدل اقتران المدار المغزلي (SOC) تشتت النطاقات المحلي ولكنه لا يفتح فجوة أو يغير عدد تقاطعات سطح فيرمي. لا يُلاحظ أي انتقال ليفشيتز (إعادة بناء طوبولوجية لسطح فيرمي) بين 0 و10 جيجا باسكال لـ Y2GaB، رغم أن الجيب المركزي يتقلص بشكل كبير. تهيمن مدارات Y-4d على الحالات الإلكترونية بالقرب من مستوى فيرمي، مع مساهمات ثانوية من مدارات B-2p والمدارات p للموقع A.
• طبيعة الترابط: يكشف التحليل عن ترابط مختلط أيوني-تساهمي بين Y وB، يتميز بتراكم الشحنة حول B ونقصانها حول Y. وتُظهر عناصر الموقع A (Ga/In) ترابطًا معدنيًا أكثر انتشارًا.
• الخصائص البصرية: يُظهر كلا المركبين استجابات درود معدنية وتفاوتًا بصريًا قويًا بين الاستقطابات [100] و[001]. ويعملان كعاكسات فعالة في المنطقة المرئية وامتصاصات في الأشعة فوق البنفسجية. يحفز الضغط إزاحات غير منتظمة في الميزات البصرية، مما يرفع عمومًا الطاقات المميزة لـ Y2GaB.
• ديناميكا الفونون: تُظهر حسابات فونون وحدة الخلية ترددات موجبة لكلا المركبين عند 0 و10 جيجا باسكال، مما يشير إلى استقرار ديناميكي ضمن هذا المعالجة المقيدة. ومع ذلك، تكشف حسابات الخلية الفائقة الموسعة عن أنماط تخيلية، مما يسلط الض الضوء على حساسية استقرار الفونون منخفض التردد لحجم الخلية ومعالجة ثابت القوة. ومن الجدير بالذكر أن Y2InB يظل غير مستقر مرنًا عند 10 جيجا باسكال رغم إظهاره ترددات موجبة في طيف فونون وحدة الخلية، مما يؤكد التمييز بين معايير الاستقرار المرني والاستقرار الديناميكي.

الأهمية والادعاءات
تضع الورقة إطارًا مقارنًا لفهم كيفية تعديل استبدال الموقع A (Ga مقابل In) والضغط الهيدروستاتيكي لخصائص أطوار MAB. وتتمثل مساهماتها الرئيسية في:

1. تحديد نوافذ الاستقرار: تحدد النطاقات الضغطية التي يكون فيها Y2GaB صالحًا ميكانيكيًا (0–10 جيجا باسكال) وتوضح أن Y2InB مستقر فقط عند الضغط الجوي وفقًا للمعايير المعتمدة.
2. فك الارتباط بين أنواع الاستقرار: توضح الدراسة صراحةً أن الاستقرار المرني (مقاومة الانفعال المتجانس) والاستقرار الديناميكي (مقاومة الإزاحات الذرية) هما أمران متمايزان ويجب تقييمهما بشكل منفصل، مع ملاحظة حساسية حسابات الفونون لحجم الخلية الفائقة.
3. العلاقات بين البنية والخصائص: تفصل الدراسة كيف يعزز الضغط بشكل منهجي الصلابة، والنقل الحراري، والاستجابة البصرية في Y2GaB مع تقليل التباين المرني، مما يوفر خط أساسًا لأي تحقق تجريبي مستقبلي.
4. النطاق المتواضع: يذكر المؤلفون صراحةً عدم وجود مشتقات MABene مثبتة تجريبيًا لهذه المركبات؛ حيث تركز الدراسة على الأطوار الأم كضرورة تأسيسية لأي تقييم مستقبلي للمشتقات ذات الأبعاد المختزلة. وتُقدم قيم الموصلية الحرارية كتقديرات نماذج (Slack/Clarke) وليس كتنبؤات نقل دقيقة، مع الإقرار بغياب حسابات صريحة لتشتت الإلكترون-فونون واللانمطية.