2025-11-21T07:28:22.047250

Thermoelectric transport properties of electron doped pyrite FeS2

Mukherjee, Subedi
Pyrite FeS$_2$ has been investigated for a wide range of applications, including thermoelectrics due to previous observation of large thermopower at room-temperature. However, the values of thermopower reported in the literature is extremely sensitive to the nature of sample -- whether they are natural or lab grown, bulk crystals or thin films -- and an ambiguity in the magnitude and sign of thermopower of pure FeS$_2$ exists. Variation in the magnitude of room-temperature thermopower has also been observed in Co-doped samples. Therefore, it is of interest to clarify the intrinsic thermopower of this system that could be measured in more pure samples. In this paper, we investigate the thermoelectric properties of Co-doped FeS$_2$ using first principles calculations. We apply three different doping schemes to understand the effect of electron doping in FeS$_2$, namely explicit Co-substitution, jellium doping and electron addition within rigid band approximation (RBA) picture. The calculated thermopower is less than $-50$ $μ$V/K for all values of Co doping that we studied, suggesting that this system may not be useful in thermoelectric applications. Interestingly, we find that RBA substantially overestimates the magnitude of calculated thermopower compared to the explicit Co-substitution and jellium doping schemes. The overestimation occurs because the changes in the electronic structure due to doping-induced structural modification and charge screening is not taken into account by the rigid shift of the Fermi level within RBA. RBA is frequently used in first principles investigations of the thermopower of doped semiconductors, and Co-substituted FeS$_2$ illustrates a case where it fails.
academic

خصائص النقل الكهروحراري لـ FeS₂ الم�عم بالإلكترونات

المعلومات الأساسية

  • معرّف الورقة: 2311.06135
  • العنوان: Thermoelectric transport properties of electron doped pyrite FeS2
  • المؤلفون: Anustup Mukherjee, Alaska Subedi (CPHT, CNRS, École Polytechnique, Institut Polytechnique de Paris)
  • التصنيف: cond-mat.mtrl-sci
  • تاريخ النشر: 13 نوفمبر 2023
  • رابط الورقة: https://arxiv.org/abs/2311.06135

الملخص

يتم دراسة البيريت FeS₂ على نطاق واسع للتطبيقات الكهروحرارية نظراً للقوة الكهروحرارية الكبيرة المرصودة في درجة حرارة الغرفة. ومع ذلك، فإن قيم القوة الكهروحرارية المبلغ عنها في الأدبيات حساسة للغاية لخصائص العينة (طبيعية أو معدة في المختبر، بلورات سائبة أو أغشية رقيقة)، مما يؤدي إلى غموض في حجم وإشارة القوة الكهروحرارية للـ FeS₂ النقي. كما توجد تباينات في حجم القوة الكهروحرارية في درجة حرارة الغرفة للعينات المعدة بـ Co. تدرس هذه الورقة الخصائص الكهروحرارية لـ FeS₂ المعدى بـ Co باستخدام حسابات من مبادئ أولى، مستخدمة ثلاث خطط تطعيم مختلفة: الاستبدال الصريح لـ Co، والتطعيم بـ jellium، والتقريب الصارم للنطاق (RBA). تبين الحسابات أن القوة الكهروحرارية عند جميع تركيزات Co المعدى أقل من -50 μV/K، مما يشير إلى أن هذا النظام قد لا يكون عملياً للتطبيقات الكهروحرارية. ومن المثير للاهتمام أن RBA يبالغ بشكل كبير في تقدير القوة الكهروحرارية لأنه لا يأخذ في الاعتبار التعديلات الهيكلية الناجمة عن التطعيم وتأثيرات الحماية الكهربائية على البنية الإلكترونية.

الخلفية البحثية والدافع

تعريف المشكلة

  1. عدم الاتساق في البيانات التجريبية: تتراوح قيم القوة الكهروحرارية المقاسة لـ FeS₂ من -115 إلى 524 μV/K مع عدم اتساق في الإشارة
  2. الاعتماد على خصائص العينة: تعتمد الخصائص الكهروحرارية بشدة على نوع العينة (طبيعية مقابل اصطناعية، سائبة مقابل رقيقة)
  3. دقة الطرق النظرية: تستخدم معظم الدراسات النظرية الموجودة التقريب الصارم للنطاق (RBA)، لكن دقتها موضع شك

أهمية البحث

  • يتكون FeS₂ من عناصر وفيرة وغير سامة على الأرض، مع إمكانيات في تطبيقات الطاقة مثل الخلايا الكهروضوئية وأقطاب البطاريات والتطبيقات الكهروحرارية
  • من الضروري توضيح خصائصه الكهروحرارية الجوهرية لتقييم قيمة التطبيق العملي
  • توفير إرشادات منهجية لدراسات الخصائص الكهروحرارية النظرية للأنظمة المعدى

حدود الطرق الموجودة

  • يحاكي التقريب الصارم للنطاق (RBA) تأثيرات التطعيم فقط من خلال الإزاحة الصارمة لمستوى فيرمي
  • يتجاهل التغييرات الهيكلية الناجمة عن التطعيم وتأثيرات الحماية الكهربائية
  • يفتقر إلى المقارنة المنهجية بين خطط التطعيم المختلفة

المساهمات الأساسية

  1. المقارنة المنهجية لثلاث خطط تطعيم: أول مقارنة بين الاستبدال الصريح لـ Co والتطعيم بـ jellium و RBA في حسابات الخصائص الكهروحرارية لـ FeS₂
  2. الكشف عن حدود RBA: اكتشاف أن RBA يبالغ بشكل كبير في تقدير القوة الكهروحرارية مع شرح الأسباب من منظور البنية الإلكترونية
  3. توضيح الأداء الكهروحراري لـ FeS₂ المعدى بـ Co: إثبات أن القوة الكهروحرارية عند جميع تركيزات التطعيم أقل من 50 μV/K، مما يشير إلى عدم ملاءمة هذا النظام للتطبيقات الكهروحرارية
  4. توفير إرشادات منهجية: توفير رؤى منهجية مهمة لدراسات الخصائص الكهروحرارية النظرية للأنظمة المعدى

شرح تفصيلي للطرق

الطرق الحسابية

استخدام إطار نظرية الدالة الكثافة (DFT) مع حزمة محاكاة Vienna Ab initio (VASP):

  • دالة التبادل والارتباط: التقريب الكثافة المحلية (LDA)
  • طريقة الكمون الزائف: طريقة الموجات المستوية المعززة بالإسقاط (PAW)
  • قطع الطاقة: ≥420 eV
  • شبكة نقاط k: 8×8×8 للحلقة الذاتية، 24×24×24 لحساب كثافة الحالات

ثلاث خطط تطعيم

1. التقريب الصارم للنطاق (RBA)

محاكاة الإلكترونات الإضافية من خلال إزاحة مستوى فيرمي في البنية الإلكترونية لـ FeS₂، المقابلة لقيمة تطعيم Co المكافئة.

2. التطعيم بـ Jellium

إضافة صريحة للإلكترونات كشحنة خلفية في خلية FeS₂ الأولية، مع إجراء حسابات البنية الإلكترونية.

3. استبدال Co الصريح

استبدال تدريجي لذرات Fe بذرات Co، مع الأخذ في الاعتبار تأثيرات الاستبدال الكيميائي الفعلي.

حساب خصائص النقل

استخدام رمز BoltzTraP2 لحل معادلة النقل البولتزمان الخطية تحت تقريب وقت الاسترخاء الثابت (CRTA):

معامل سيبيك: S(μ,T)=ναjσαiS(\mu, T) = \frac{\nu_{\alpha j}}{\sigma_{\alpha i}}

معامل هول: RH(μ,T)=ναβkσαjσiβR_H(\mu, T) = \frac{\nu_{\alpha \beta k}}{\sigma_{\alpha j}\sigma_{i\beta}}

حيث يتم تعريف موتر النقل كـ: σαβ(μ,T)=1Ω∫σαβ(ϵ)[−∂ϵfμ(ϵ)]dϵ\sigma_{\alpha\beta}(\mu, T) = \frac{1}{\Omega}\int \sigma_{\alpha\beta}(\epsilon)[-\partial_\epsilon f_\mu(\epsilon)]d\epsilon

الإعداد التجريبي

نموذج البنية

  • البنية البلورية: البيريت FeS₂، مجموعة الفضاء Pa3
  • تكوين الخلية الأولية: 4 ذرات Fe (موضع 4a)، 8 ذرات S (موضع 8c)
  • تركيزات التطعيم: تم دراسة تطعيم Co بنسب 1%، 5%، 10%، 15%، 25%، 50%

الاعتبارات المغناطيسية

تم الأخذ في الاعتبار كل من الحالات غير المغناطيسية والحديدية المغناطيسية، حيث تشير التجارب إلى أن FexCo1-xS₂ حديدية مغناطيسية عند x≥0.1.

تحسين البنية

  • تتبع ثوابت الشبكة للمركبات المختلطة قانون Vegard
  • تم إرخاء المواضع الذرية لتقليل القوى الذرية
  • المعامل الداخلي لـ FeS₂ هو u=0.3821

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

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

القوة الكهروحرارية لـ Fe₀.₇₅Co₀.₂₅S₂ (درجة حرارة الغرفة)

خطة التطعيمالحالة غير المغناطيسية (μV/K)الحالة الحديدية المغناطيسية (μV/K)
RBA-46.71-46.71
Jellium-14.68-14.69
استبدال Co-19.67-20.49

القوة الكهروحرارية لـ Fe₀.₅Co₀.₅S₂ (درجة حرارة الغرفة)

خطة التطعيمالحالة غير المغناطيسية (μV/K)الحالة الحديدية المغناطيسية (μV/K)
RBA-39.60-39.60
Jellium-17.52-17.52
استبدال Co-25.58-25.58

الاكتشافات الرئيسية

  1. المبالغة المنهجية في RBA: يبالغ RBA بشكل كبير في تقدير القيمة المطلقة للقوة الكهروحرارية عند جميع تركيزات التطعيم
  2. الاتساق في الاعتماد على درجة الحرارة: تظهر جميع الطرق الثلاث أن القوة الكهروحرارية تزداد مع زيادة درجة الحرارة
  3. نوع الناقلات: معامل سيبيك سالب في جميع الحالات، مما يشير إلى أن الإلكترونات هي الناقلات الرئيسية
  4. الأداء الكهروحراري المحدودة: القوة الكهروحرارية في درجة حرارة الغرفة أقل من 50 μV/K عند جميع تركيزات التطعيم

تحليل البنية الإلكترونية

أسباب المبالغة في RBA

من خلال تحليل كثافة الحالات (DOS) وبنية النطاق:

  • في RBA: يقع مستوى فيرمي في منطقة شديدة الانحدار من كثافة الحالات
  • في التطعيم بـ Jellium: ينتقل مستوى فيرمي إلى منطقة كتف أوسع
  • في استبدال Co: يقع مستوى فيرمي في منطقة نسبياً مسطحة

نظراً لأن معامل سيبيك S ∝ d ln N(EF)/dE، فإن كثافة الحالات الأكثر انحداراً بالقرب من مستوى فيرمي في RBA تؤدي إلى قيم S أكبر.

النتائج عند تركيزات تطعيم مختلفة

التركيبRBA (μV/K)Jellium (μV/K)
Fe₀.₉₉Co₀.₀₁S₂-32.45-28.52
Fe₀.₉₅Co₀.₀₅S₂-10.47-10.24
Fe₀.₉Co₀.₁S₂-17.11-10.13
Fe₀.₈₅Co₀.₁₅S₂-25.41-11.75

الأعمال ذات الصلة

الحالة الحالية للدراسات التجريبية

  • البيانات التاريخية: تتراوح معاملات سيبيك المبلغ عنها في النصف الأول من القرن العشرين من -115 إلى 524 μV/K
  • الدراسات الحديثة:
    • Kato وآخرون: عينات n من النوع الطبيعي S = -300 μV/K
    • Willeke وآخرون: بلورات أحادية عالية الحركة تصل إلى S = -320 μV/K
    • دراسات الأغشية الرقيقة: تقارير عامة عن قيم S موجبة (60-75 μV/K)

الدراسات النظرية

  • الحسابات المبكرة: استخدم Gudelli وآخرون و Harran وآخرون RBA للتنبؤ بـ S محسّن يتجاوز -400 μV/K
  • مساهمة هذا العمل: أول مقارنة منهجية بين خطط التطعيم المختلفة، مع الكشف عن حدود RBA

الخلاصات والمناقشة

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

  1. الأداء الكهروحراري المحدودة: FeS₂ المعدى بـ Co عند جميع تركيزات التطعيم المدروسة، القوة الكهروحرارية في درجة حرارة الغرفة أقل من 50 μV/K، وهي أقل بكثير من معايير المواد الكهروحرارية الممتازة (>200 μV/K)
  2. حدود طريقة RBA: يبالغ التقريب الصارم للنطاق بشكل كبير في تقدير القوة الكهروحرارية لأنه يتجاهل:
    • التعديلات الهيكلية الناجمة عن التطعيم
    • تأثيرات الحماية الكهربائية
    • التغييرات في اتساع النطاق والانقسام النسبي
  3. الأهمية المنهجية: FeS₂ المعدى بـ Co هو مثال نموذجي لفشل RBA، مما يوفر تحذيراً مهماً لدراسات الخصائص الكهروحرارية للأنظمة المعدى

القيود

  1. قيود الطرق الحسابية: قد يؤدي استخدام LDA إلى التقليل من تقدير فجوة النطاق، مما يؤثر على دقة الخصائص الكهروحرارية
  2. نطاق تركيز التطعيم: تقتصر خطة استبدال Co الصريحة على حجم الخلية الأولية، مما يمنع دراسة تركيزات تطعيم منخفضة جداً
  3. تقريب وقت الاسترخاء: قد يؤثر اعتماد تقريب وقت الاسترخاء الثابت على دقة خصائص النقل

الاتجاهات المستقبلية

  1. تحسين الطرق الحسابية: استخدام دوال تبادل وارتباط أكثر دقة (مثل HSE06)
  2. توسيع الأنظمة المدروسة: دراسة تأثيرات التطعيم بمعادن انتقالية أخرى
  3. التحقق التجريبي: تحضير عينات عالية الجودة للتحقق من التنبؤات النظرية

التقييم المتعمق

المميزات

  1. الابتكار المنهجي: أول مقارنة منهجية بين ثلاث خطط تطعيم، توفير معيار منهجي مهم للمجال
  2. الرؤى الفيزيائية العميقة: شرح عميق للآلية الفيزيائية للمبالغة في RBA من منظور البنية الإلكترونية
  3. الصرامة الحسابية: الأخذ في الاعتبار التأثيرات المغناطيسية، استخدام معاملات حسابية مناسبة ومعايير تقارب
  4. القيمة العملية الواضحة: توضيح أن FeS₂ المعدى بـ Co غير مناسب للتطبيقات الكهروحرارية، مما يتجنب الاستثمار التجريبي غير الفعال

أوجه القصور

  1. قيود الطرق النظرية: قد لا تكون LDA دقيقة بما يكفي، خاصة بالنسبة للأنظمة الإلكترونية المترابطة
  2. نطاق درجة الحرارة المحدود: التركيز الأساسي على خصائص درجة حرارة الغرفة، مع نقص الدراسات المنهجية في درجات الحرارة العالية
  3. آلية التطعيم الموحدة: دراسة التطعيم الإلكتروني فقط، دون الأخذ في الاعتبار احتمالية التطعيم بالفجوات

التأثير

  1. المساهمة الأكاديمية: توفير إرشادات منهجية مهمة لدراسات الخصائص الكهروحرارية النظرية للأنظمة المعدى
  2. القيمة العملية: مساعدة الباحثين على تجنب إهدار الموارد على أنظمة غير واعدة
  3. قابلية إعادة الإنتاج القوية: وصف تفاصيل حسابية كافية، سهلة إعادة الإنتاج والتوسع

السيناريوهات المعمول بها

  1. البحث النظري: توفير مرجع منهجي لدراسات الخصائص الكهروحرارية للأنظمة المعدى الأخرى
  2. فحص المواد: استبعاد المواد المرشحة غير المناسبة في الفحص عالي الإنتاجية لمواد كهروحرارية
  3. التوجيه التجريبي: توفير أساس نظري لاختيار استراتيجيات التطعيم المناسبة للباحثين التجريبيين

المراجع

تستشهد الورقة بـ 38 مرجعاً مهماً، تغطي الخصائص الأساسية لـ FeS₂ والدراسات التجريبية وطرق الحساب النظري وغيرها، مما يوفر أساساً أدبياً متيناً للبحث. من الجدير بالملاحظة بشكل خاص الاستشهاد بأحدث الأعمال التجريبية لـ Xi وآخرين، والتي يبدو أنها حلت الخلاف السابق حول نوع الناقلات في FeS₂.

التقييم الشامل: هذه ورقة بحثية نظرية عالية الجودة، لا توضح فقط الخصائص الكهروحرارية لـ FeS₂ المعدى بـ Co، بل الأهم من ذلك أنها تكشف عن حدود طريقة RBA المستخدمة على نطاق واسع، مما يوفر مساهمة منهجية مهمة للمجال. تتمتع نتائج البحث بقيمة مهمة في توجيه الفحص النظري والدراسة التجريبية لمواد كهروحرارية.