🔬 materials science

Influence of ligand field and correlation on the electronic structure of NiO and CoO from DFT+DMFT calculations

تستخدم هذه الدراسة حسابات نظرية الكثافة الوظيفية مضافة إليها نظرية المجال المتوسط الديناميكي ذات الشحنة ذاتية الاتساق (charge self-consistent DFT+DMFT) لاستقصاء كيفية تأثير البنية البلورية، والمجالات الليغاندية، وقوى الارتباط المتغيرة (بما في ذلك ارتباطات أكسجين 2p) على البنية الإلكترونية والدوال الطيفية لأكسيد النيكل (NiO) وأكسيد الكوبالت (CoO) في الحالة البارامغناطيسية.

المؤلفون الأصليون: Daniel Mutter, Frank Lechermann, Daniel F. Urban, Christian Elsässer

نُشر 2026-05-06

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: Daniel Mutter, Frank Lechermann, Daniel F. Urban, Christian Elsässer

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

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

هذه الورقة البحثية هي غوص عميق لفهم كيفية سلوك هذه الإلكترونات بدقة في أكسيد النيكل (NiO) وأكسيد الكوبالت (CoO). استخدم المؤلفون طريقة محاكاة حاسوبية قوية تسمى DFT+DMFT (وهي مزيج من نظريتين فيزيائيتين متقدمتين) لإنشاء "فيلم" لمستويات طاقة الإلكترونات، والتي يسمونها الدوال الطيفية (spectral functions). فكر في هذا كخريطة توضح الأماكن التي تحب الإلكترونات التواجد فيها، وكمية الطاقة التي تحتاجها للقفز إلى مكان جديد.

إليك تفصيل لنتائجهم باستخدام تشبيهات بسيطة:

1. المكونان الرئيسيان: "الزحام" و"شكل الغرفة"

نظر الباحثون في شيئين رئيسيين يغيران كيفية سلوك الإلكترونات:

الزحام (قوة الارتباط): مدى قوة دفع الإلكترونات لبعضها البعض. في المحاكاة، قام المؤلفون بضبط قرص تحكم يسمى $U$ لمعرفة ما يحدث عندما تصبح الإلكترونات أكثر "تذمراً" وترغب في البقاء متباعدة.
شكل الغرفة (المجال الرابط - Ligand Field): تترتب الذرات في هذه المواد في أشكال هندسية محددة. قارن المؤلفون بين شكلين:
- ملح الصخر (Rock-Salt - RS): مثل المكعب حيث تحاط ذرة المعدن بذرات الأكسجين من جميع الجوانب (ثماني أوجه). هذا هو الشكل الطبيعي والمستقر لهذه المواد.
- الزنك بلايند (Zincblende - ZB): مثل شكل الهرم (رباعي الأوجه). وهو شكل غير مستقر، "افتراضي"، اختبره المؤلفون لمعرفة كيف يغير الشكل رقصة الإلكترونات.

2. تصحيح "التفاعل الذاتي" (SIC)

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

التشبيه: تخيل حفلة صاخبة. النموذج القياسي يعتقد أن ذرات الأكسجين تصرخ على نفسها، مما يخلق الكثير من الضجيج في الخلفية الذي يفسد الموسيقى. طريقة الـ SIC تقوم بخفض هذا الضجيج المحدد.
النتيجة: عندما قاموا بتفعيل ميزة "إلغاء الضجيج" هذه للأكسجين، تطابقت المحاكاة مع التجارب الواقعية بشكل أفضل. فقد تنبأت بشكل صحيح بحجم "الفجوة" (الطاقة اللازمة لجعل المادة موصلة للكهرباء)، لتصل مباشرة إلى نطاق 5-6 إلكترون فولت الموجود في المختبرات الحقيقية.

3. النيكل مقابل الكوبالت: فرق "الخطوة الواحدة"

النيكل والكوبالت جيران في الجدول الدوري. الفرق الوحيد هو أن الكوبالت لديه إلكترون واحد أقل في "فرقة الرقص" الخاصة به من النيكل.

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

4. الشكل يهم: المكعبات مقابل الأهرامات

عندما قارنوا بين هيكل ملح الصخر (المكعب) المستقر وهيكل الزنك بلايند (الهرمي) غير المستقر:

مستويات الطاقة: تغير "الغرف" (الأشكال الهندسية) مستويات طاقة مدارات الإلكترونات. في المكعب، تكون بعض المدارات ذات طاقة أقل؛ أما في الهرم، فينقلب الترتيب.
الاستقرار: يوضح المؤلفون أن الشكل الهرمي أقل استقراراً لهذه المواد لأنه يجبر المزيد من الإلكترونات على الدخول في حالات "الترابط العكسي" (antibonding).
- التشبيه: فكر في حالات "الترابط" ككراسي مريحة حيث تكون الإلكترونات سعيدة. أما حالات "الترابط العكسي" فهي مثل الجلوس على مقعد مهتز. الشكل الهرمي يجبر المزيد من الإلكترونات على الجلوس على مقاعد مهتزة، مما يجعل الهيكل بأكمله مهتزاً. وهذا يفسر لماذا يشكل أكسيد النيكل مكعبات بشكل طبيعي، بينما الشكل الهرمي نادر وغير مستقر.

5. "القمر الصناعي" و"الفجوة"

تبحث الورقة في سمات محددة في خريطة طاقة الإلكترونات:

قمة القمر الصناعي (The Satellite Peak): "صدى" عالي الطاقة ناتج عن تفاعلات إلكترونية قوية. وجد المؤلفون أن طريقتهم يمكنها رؤية هذا الصدى، لكنه كان باهتاً في بعض الحسابات.
فجوة النطاق (The Band Gap): الرقم الأكثر أهمية. إنها جدار الطاقة الذي يمنع تدفق الكهرباء (مما يجعل المادة عازلة).
- النتيجة: بدون تصحيح الأكسجين (SIC)، توقعت نماذجنا الحاسوبية جداراً منخفضاً جداً (حوالي 2.5-3 إلكترون فولت). ومع تصحيح الأكسجين، قفز الجدار إلى 5.1 إلكترون فولت، وهو ما يطابق ما يقيسه العلماء في التجارب الحقيقية.

ملخص

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

ملخص تقني: تأثير مجال الليجاند والارتباط على البنية الإلكترونية لمركبي NiO وCoO

بيان المشكلة
تُظهر أكاسيد المعادن الانتقالية (TMOs) مثل NiO وCoO سلوكيات إلكترونية معقدة مدفوعة بالتفاعل بين ارتباط الإلكترونات، ومجالات الليجاند، والتركيب البلوري، والتهجين. ورغم أن هذه المواد واعدة في تطبيقات الكهروحرارية، والإلكترونيات الضوئية، والحفز، إلا أن بنيتها الإلكترونية يصعب نمذجتها بدقة. غالبًا ما تفشل نظرية الدالة الوظيفية للكثافة (DFT) القياسية باستخدام تقريب الكثافة الموضعية (LDA) أو تقريب التدرج المعمم (GGA) في إعادة إنتاج الأطياف التجريبية، حيث تتنبأ بشكل خاطئ بحالات معدنية للأطوار العازلة لهذه المركبات. وعلى الرغم من أن تقنية DFT+U تقدم تحسينات، إلا أنها تفتقر إلى المعالجة الديناميكية لتقلبات الشحنة واللف المغزلي الضرورية للأطوار البارامغناطيسية عالية الحرارة. ويكمن التحدي في التقاط الارتباطات القوية في مدارات $3d$ للمعدن الانتقالي، والتهجين مع مدارات الأكسجين $2p$ ، وتأثير بيئات المجال البلوري المختلفة (ثماني الأوجه مقابل رباعي الأوجه) على الدوال الطيفية في آن واحد.

المنهجية
يستخدم المؤلفون مزيجًا من نظرية الدالة الوظيفية للكثافة ونظرية المجال المتوسط الديناميكي (DFT+DMFT) المتسق شحنيًا لاستقصاء مركبي NiO وCoO البارامغناطيسيين. وتستخدم الدراسة الإطار المنهجي المحدد التالي:

الإعداد الحسابي: تُجرى الحسابات باستخدام كود الجهد الكاذب ذو القاعدة المختلطة (MBPP) لجزء الـ DFT، ومحلل (CTHYB) للزمن المستمر ضمن حزمة TRIQS لمشكلة الشوائب في DMFT.
الفضاء المرتبط: يتم تحديد الفضاء المرتبط باستخدام مدارات $3d$ الشبيهة بالذرية المتمركزة على أيونات المعدن الانتقالي، والمُسقطة على قاعدة كوهن-شام مبتورة.
معاملات التفاعل: يتم حل مشكلة هوبارد متعددة المدارات المعممة باستخدام هاميلتوني سلايتر المتباين دورانياً. وتغير الدراسة معامل هوبارد $U$ (6 إلكترون فولت و9 إلكترون فولت لـ CoO؛ 10 إلكترون فولت لـ NiO) وتفاعل هوند $J$ (حوالي 1 إلكترون فولت).
تصحيح التفاعل الذاتي (SIC): لمعالجة آثار الارتباط في مدارات الأكسجين $2p$ ، والتي تهجن بقوة مع حالات المعدن الانتقالي $3d$ ، يطبق المؤلفون مخطط جهد كاذب لتصحيح التفاعل الذاتي (DFT+sicDMFT). ويتضمن ذلك تعديل الجهود الكاذبة للأكسجين بمعاملات وزن محددة ( $w_l=1.0$ لـ $2s$ و $w_l=0.8$ لـ $2p$ ).
الهياكل: تقارن الدراسة بين هيكلين بلوريين: هيكل ملح الصخر (RS) المستقر ذي مجالات الليجاند ثمانية الأوجه، وهيكل الزنك بليند (ZB) شبه المستقر ذي مجالات الليجاند رباعية الأوجه.
التحليل: يتم اشتقاق الدوال الطيفية للتردد الحقيقي عبر الاستمرار التحليلي (طريقة الإنتروبيا القصوى) من ترددات ماتسوبارا التخيلية. ويركز التحليل على الدوال الطيفية الكلية، والمُسقطة على العناصر (المعدن الانتقالي والأكسجين)، والمُحللة على مستوى المدار المحلي.

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

الدوال الطيفية لـ RS-CoO:
- تحدد الدراسة قمة منخفضة الطاقة تحت مستوى فيرمي تقابل الحالة النهائية $d^7\underline{L}$ (حيث $\underline{L}$ هي فجوة الليجاند). وتنتقل هذه القمة من $-1.0 $إلكترون فولت إلى$ -1.6 $إلكترون فولت مع زيادة$ U$ من 6 إلى 9 إلكترون فولت.
- يؤدي تطبيق O( $2p$ )-SIC إلى إزاحة هذه القمة إلى $-0.6 $إلكترون فولت، مما يتوافق بشكل أفضل مع بيانات UPS تجريبية محددة، رغم أن بيانات XPS تشير إلى موقع قريب من$ -1.9$ إلكترون فولت.
- يؤدي تضمين O( $2p$ )-SIC إلى تحسين وصف فجوة النطاق بشكل كبير. فبدون SIC، يتم تقدير الفجوة بين القمم المشغولة وغير المشغولة بأقل من قيمتها. ومع استخدام $U=9$ إلكترون فولت وSIC، تصبح المسافة المحسوبة 5.1 إلكترون فولت، وهو ما يقع ضمن النطاق التجريبي لـ 5–6 إلكترون فولت. وقد أبلغت دراسات سابقة لـ DFT+DMFT بدون SIC عن قيم منخفضة تترا-وح 2.5–3.8 إلكترون فولت.
- حزمة هوبرد السفلى (القمة التابعة عند $\approx -9.5$ إلكترون فولت في التجارب) لا تظهر بوضوح في الدالة الطيفية الكلية بدون SIC، ولكنها تظهر كبروز طفيف مع SIC، وإن كانت مزاحة إلى $-7$ إلكترون فولت.
التأثيرات الهيكلية والتركيبية (RS مقابل ZB، NiO مقابل CoO):
- التحقق من نظرية المجال البلوري: تؤكد النتائج تنبؤات نظرية المجال البلوري (CF). في هياكل RS (ثمانية الأوجه)، تكون حالات $t_{2g}$ أقل طاقة من $e_g$ . وفي هياكل ZB (رباعية الأوجه)، ينعكس الترتيب ( $e$ أقل من $t_2$ ).
- الاحتلال: يُظهر RS-NiO تكوين $t_{2g}^6 e_g^2$ ، بينما يُظهر RS-CoO تكوين $t_{2g}^5 e_g^2$ . وتُظهر هياكل ZB تكوين $e^4 t_2^4$ (لـ NiO) و $e^4 t_2^3$ (لـ CoO).
- الاستقرار: وُجد أن هياكل ZB أقل استقرارًا بسبب احتلال حالات مضادة للترابط. في ZB، تكون حالات $t_2$ مضادة للترابط؛ حيث تمتلك CoO وNiO في ZB 3 و4 إلكترونات في هذه الحالات، على التوالي، مقارنة بإلكترونين في حالات $e_g$ المضادة للترابط في طور RS المستقر.
- التشابهات الطيفية: رغم الاختلافات الهيكلية، تُظهر الدوال الطفية المحلية لأطوار RS وZB تشابهات نوعية عند مقارنة المجموعات المتقابلة (على سبيل المثال، RS- $t_{2g}$ تشبه ZB- $e$ ).
- CoO مقابل NiO: تُظهر أطياف CoO قممًا أكثر حدة وزيادة في الوزن الطيفي عند الطاقات المنخفضة مقارنة بـ NiO. ويعزو المؤلفون القمة الحادة منخفضة الطاقة في CoO إلى حالات $d^7\underline{Z}$ (حالات Zhang-Rice)، المماثلة لـ $d^8\underline{Z}$ في NiO، والناتجة عن الفجوة الإضافية في التكوين التكافؤي.
تأثير إعدادات الارتباط:
- تأثير SIC: يؤدي تطبيق O( $2p$ )-SIC باستمرار إلى إزاحة النطاقات غير المشغولة المسقطة على المعدن الانتقالي (حزم هوبرد العليا) والنطاقات المسقطة على الأكسجين إلى طاقات أعلى. كما يزيد من فرق الطاقة ( $\Delta E_{pd}$ ) بين حالات المعدن الانتقاني غير المشغولة وحالات الأكسجين المشغولة، مما يقلل من تقلبات الشحنة.
- اعتماد U: تزيد زيادة $U$ من المسافة بين مراكز الثقل للأطياف المسقطة على المعدن الانتقالي (المشغولة وغير المشغولة)، وهو ما يتوافق مع نموذج هوبارد.
- طاقة انتقال الشحنة ( $\Delta$ ): تميز الدراسة بين طاقة انتقال الشحنة الفعالة $\Delta$ (المشتقة من مستويات DFT أحادية الجسيم) وبين الفجوة الطيفية. وتتبع قيم $\Delta$ المحسوبة لـ RS-NiO (4.3 إلكترون فولت) و RS-CoO (4.9 إلكترون فولت) التسلسل التجريبي (5.4 إلكترون فولت و6.1 إلكترون فولت، على التوالي) لكنها أقل تقديرًا قليلاً.

الأهمية والادعاءات
يزعم البحث أن نهج DFT+sicDMFT يوفر إطارًا قويًا لوصف البنية الإلكترونية لأكاسيد المعادن الانتقالية شديدة الارتباط في أطوارها البارامغناطيسية. وتتمثل المساهمة الرئيسية في إثبات أن مراعاة الارتباطات في مدارات الأكسجين $2p$ صراحةً عبر SIC أمر بالغ الأهمية للحصول على فجوات نطاق وميزات طيفية تتوافق مع أطياف الاستثارة التجريبية (خاصة لـ RS-CoO). نجحت الدراسة في فصل تأثيرات التركيب البلوري (مجال الليجاند) وقوة الارتباط الإلكتروني، مما يدعم تطبيق نظرية المجال البلوري على الميزات الطيفية لهذه الأنظمة المعقدة. ويضع المؤلفون هذه النتائج كواصفات يمكن أن تدعم البحث وتحسين أكاسيد المعادن الانتقالية لأجهزة تحويل الطاقة، مشيرين إلى أن السمات المميزة لأطياف استثارة الإلكترونات ترتبط بالنشاط الحفزي. لا يقترح العمل بروتوكولات تجريبية جديدة، بل يقدم منهجية حسابية منقحة لتفسير البيانات التجريبية الموجودة وتوجيه فهم وظائف المواد.