Revisiting Phase Transitions of Yttrium: Insights from Density Functional Theory

ملخص تقني: إعادة النظر في التحولات الطورية لليتريوم: رؤى من نظرية وظيفية الكثافة

بيان المشكلة
يظل فهم التحولات الطورية الهيكلية للعناصر الأرضية النادرة تحت الضغط العالي تحديًا جوهريًا في فيزياء المادة المكثفة. ويعد اليتريوم (Y)، وهو معدن انتقالي ثلاثي التكافؤ لا يحتوي على إلكترونات 4f، نموذجًا مرجعيًا حاسمًا لاستقصاء هذه الآليات نظرًا لتشابهه الهيكلي مع سلسلة اللانثانيدات. وبينما تشير البيانات التجريبية إلى تسلسل محدد من التحولات الطورية تحت الضغط — من السداسي مغلق الرص (hcp) إلى نوع ساماريوم (Sm-type) عند ~10 جيجا باسكال، ثم لاحقًا إلى السداسي مزدوج مغلق الرص (dhcp) عند ~25 جيجا باسكال — فإن النماذج النظرية القائمة على دوال تقريب التدرج العام (GGA) القياسية (مثل PBE) قد فشلت في إعادة إنتاج ضغوط التحول بدقة. وغالبًا ما كانت الحسابات السابقة تقلل بشكل كبير من الضغوط المطلوبة لهذه التحولات، مما خلق فجوة بين النظرية والتجربة. علاوة على ذلك، فإن آليات الاستقرار الديناميكي والديناميكي الحراري التي تقود هذه التحولات، وخاصة التفاعل بين انتقال الشحنة الإلكترونية واهتزازات الشبكة، تتطلب مزيدًا من التوضيح.

المنهجية
استخدم المؤلفون نظرية وظيفية الكثافة (DFT) للتحقيق بشكل منهجي في سلوك اليتريوم تحت ظروف الضغط المنخفض (<30 جيجا باسكال). واستخدمت الدراسة حزمة محاكاة فيينا من البداية (VASP) مع جهد كاذب من نوع (PAW). وتضمنت المكونات المنهجية الرئيسية ما يلي:

• دوال تبادل الارتباط: تم إجراء تحليل مقارن باستخدام GGA القياسي (PBE)، وmeta-GGA (مثل SCAN و r2SCAN)، والمتغيرات المصححة بالتشتت (DFT-D3 مع تخميد Becke-Johnson).
• النماذج الهيكلية: أُجريت الحسابات لكل من أطوار hcp ($P6_3/mmc$)، وSm-type ($R\bar{3}m$)، وdhcp ($P6_3/mmc$). ولضمان الدقة في حسابات طاقة الحالة الأرضية، تم استخدام خلايا فائقة (2×2×1) تحتوي على 8، 36، و16 ذرة على التوالي، متجاوزة بذلك تقريبات الخلية الأولية التي أظهرت عدم دقة عددية.
• تحليل الخصائص: حسبت الدراسة المحتويات الإنثالبية النسبية لتحديد استقرار الطور، وعلاقات الحجم والضغط، وتجمعات شحنة موليكن لتتبع انتقال الإلكترونات من مدار s إلى مدار d.
• الاستقرار الديناميكي والميكانيكي: تم حساب منحنيات تشتت الفونون باستخدام طريقة الفرق المحدود مع نهج الخلية الفائقة لتحديد عدم الاستقرار الاهتزازي. كما حُسبت الثوابت المرنة ($C_{ij}$) والمعاملات الميكانيكية لتقييم التلين الميكانيكي عند حدود الطور.

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

1. دقة دالة r2SCAN: تثبت الدراسة أن دالة meta-GGA من نوع r2SCAN توفر تنبؤات أكثر دقة بضغوط التحول الطوري مقارنة بـ PBE-GGA.

   • PBE-GGA: توقع التحول من hcp إلى Sm-type عند 0.5 جيجا باسكال (خلية أولية) أو 2.8 جيجا باسكال (خلية فائقة)، والتحول من Sm-type إلى dhcp عند 7.5 جيجا باسكال، وكلاهما ينحرف بشكل كبير عن القيم التجريبية (~10 جيجا باسكال و ~25 جيجا باسكال، على التوالي).
   • r2SCAN: توقع التحول من hcp إلى Sm-type عند 9.2 جيجا باسكال، والتحول من Sm-type إلى dhcp عند 18.6 جيجا باسكال. وتظهر هذه القيم توافقًا ممتازًا مع البيانات التجريبية، مما يؤكد صلاحية r2SCAN كأداة متفوقة لهذا النظام.
   • تصحيحات فان دير فال: لم يؤدِ إدراج تصحيحات DFT-D3 إلى تحسين النتائج بشكل ملحوظ، مما يشير إلى أن قوى التشتت بعيدة المدى تلعب دورًا ثانويًا في التحولات الطورية لليتريوم.

2. آلية التحول الطوري:

   • انتقال الشحنة: أكد تحليل شحنة موليكن حدوث انتقال للشحنة الناجم عن الضغط من مدار 5s إلى مدارات 4d. ومع زيادة الضغط، يفقد مدار 5s الشحنة بينما يكتسب مدار 4d الشحنة، مما يعدل التفاعلات بين الذرات.
   • عدم الاستقرار الاهتزازي: تم تحديد عدم الاستقرار الاهتزازي كالمحرك الرئيسي للتحولات الطورية. فقد أظهرت منحنيات تشتت الفونون لطور hcp ظهور أنماط صوتية مستعرضة لينة على طول الاتجاه [001] ($\Gamma \to A$)، والتي أصبحت تخيلية عند 10 جيجا باسكال، مما يشير إلى التحول من hcp إلى Sm-type. وبالمثل، أظهر طور Sm-type تليناً على طول المسار $S_2 \to F$، مما أدى إلى التحول إلى طور dhcp.
   • التلين الميكانيكي: كشفت حسابات الخصائص المرنة عن تلين ميكانيكي، لا سيما في الثابت المرن $C_{44}$ لطور hcp بعد 8 جيجا باسكال. ويرتبط هذا التلين باللين الملحوظ في الفونونات، مما يشير إلى وجود صلة قوية بين عدم الاستقرار المرن والتحولات الهيكلية.

3. ملفات الحجم والإنثالبيا: أعادت دالة r2SCAN إنتاج ملفات الحجم والضغط بدقة، مطابقة انخفاضات الحجم التجريبية عند التحول الأول (~1.26%) والطبيعة المستمرة للتحول الثاني.

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

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