🔬 materials science

Anisotropic Phonon Heat Flow and Thermoelectric Response in Tetragonal GeS $_2$ and GeSe $_2$

تستخدم هذه الدراسة حسابات المبادئ الأولية لتوضيح أن مركبات GeS $_2$ و GeSe $_2$ ذات النظام البلوري الرباعي تُظهر خصائص كهروحرارية متباينة بوضوح، تتميز بانخفاض كبير في الموصلية الحرارية الشبكية عبر المستويات، وقيم معتدلة لمعامل الاستحقاق (zT)، لا سيما لمركب GeS $_2$ من النوع السالب (n-type) عند درجات الحرارة المرتفعة.

المؤلفون الأصليون: Neeraj Kulhari, Krishna Swaroop Sharma, K. C. Bhamu

نُشر 2026-08-05

📖 5 دقيقة قراءة🧠 قراءة متعمّقة

CC BY 4.0

المؤلفون الأصليون: Neeraj Kulhari, Krishna Swaroop Sharma, K. C. Bhamu

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

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

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

متاهة البلورات: حكاية هندستين

في هذه الدراسة، عمل الباحثون كمهندسين معماريين رقميين، حيث قاموا ببناء واختبار مدينتين بلوريتين خياليتين: إحداهما مكونة من الجرمانيوم والكبريت (GeS2)، والأخرى مكون من الجرمانيوم والسيلينيوم (GeSe2). لم يكتفوا ببنائهما فحسب؛ بل قاموا بمحاكاة كيفية اهتزاز الذرات وكيفية انطلاق الإلكترونات داخلها باستخدام نماذج حاسوبية قوية.

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

ازدحام مروري حراري

كانت النتائج دراماتيكية بشكل مفاجئ. أظهرت عمليات المحاكاة الحاسوبية أن هذه المواد انتقائية للغاية بشأن الاتجاه الذي يمكن للحرارة أن تسلكه.

بالنسبة للمادة القائمة على الكبريت (GeS2)، تتدفق الحرارة مثل الطريق السريع عند التحرك جانبياً عبر الطبقات. عند درجة حرارة الغرفة (300 كلفن)، تبلغ الموصلية الحرارية رقماً ضخماً قدره 26.86 واط م⁻¹ ك⁻¹. ولكن إذا حاولت دفع تلك الحرارة نفسها للأعلى عبر الكومة (اتجاه "عبر المستوى")، فإنها تصطدم بازدحام مروري هائل. تنخفض الموصلية الحرارية إلى 1.19 واط م⁻¹ ك⁻¹ فقط. هذا فرق يزيد عن 20 ضعفاً!

أما المادة القائمة على السيلينيوم (GeSe2) فتتصرف بشكل مشابه، وإن كان "الطريق السريع" أبطأ قليلاً. جانبياً، توصل الحرارة بمعدل 18.74 واط م⁻¹ ك⁻¹، ولكن صعوداً وهبوطاً، تنخفض إلى 1.52 واط م⁻¹ ك⁻¹.

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

رقصة الإلكترونات

بينما كانت الحرارة عالقة، تحقق الباحثون من كيفية سلوك الكهرباء. وجدوا أن هذه المواد هي أشباه موصلات، مما يعني أنها يمكن أن توصل الكهرباء إذا أضفت إليها بعض الإلكترونات الإضافية (نوع n) أو أزلت بعضها (نوع p).

ومع ذلك، لم تحصل الكهرباء على نفس معامل "الازدحام المروري" الذي تعرضت له الحرارة. كان بإمكان الإلكترونات التحرك، لكن المادة لم تكن موصلة فائقة. حسب الباحثون أنه بالنسبة لنسخة الكبريت (GeS2)، إذا ضبطت عدد الإلكترونات بدقة ورفعت درجة حرارتها إلى 800 كلفن، يمكنك الحصول على درجة أداء (zT) تبلغ 0.257 عند النظر في اتجاه "عبر المستوى".

بالنسبة لنسخة السيلينيوم (GeSe2)، كانت الدرجة أقل، حوالي 0.066 للنقل من النوع p عند نفس درجة الحرارة.

ماذا يعني هذا (وما لا يعنيه)

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

ومع ذلك، تشير الورقة البحثية بحذر إلى أنه على الرغم من أن هذه بداية رائعة، إلا أن هذه المواد ليست بعد "الكأس المقدسة" للطاقة. تُعت-بر قيمة zT البالغة 0.257 قيمة "متوسطة". ولكي تكون مفيدة حقاً في المحركات أو محطات الطاقة في العالم الحقيقي، يبحث العلماء عادة عن أرقام أقرب إلى 1.0 أو أعلى. كما لاحظ الباحثون أن نتائجهم تعتمد على محاكاة حاسوبية، فهم لم يبنوا آلة مادية بهذه البلورات بعد لإثبات نجاحها في العالم الحقيقي.

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

الخلاصة

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

ملخص تقني: تدفق الحرارة الفونوني متباين الخواص والاستجابة الكهروحرارية في GeS₂ و GeSe₂ رباعي الأضلاع

بيان المشكلة
تتطلب المواد الكهروحرارية، التي تحول الحرارة مباشرة إلى كهرباء، معامل جدارة (zT) مرتفعاً، وهو المعرف بنسبة عامل القدرة ( $\sigma S^2$ ) إلى الموصلية الحرارية الكلية ( $\kappa_e + \kappa_l$ ). ويتمثل التحدي المركزي في تصميم المواد الكهروحرارية في الاقتران المتبادل بين هذه المعلمات: فزيادة الموصلية الكهربائية ( $\sigma$ ) غالباً ما تؤدي إلى زيادة الموصلية الحرارية الإلكترونية ( $\kappa_e$ ). وبينما يعد تقليل الموصلية الحرارية للشبكة ( $\kappa_l$ ) من خلال الهيكلة النانوية أو السبائك استراتيجية شائعة، إلا أن هناك حاجة لاستكشاف مركبات غير سامة ومتوفرة بكثرة في الأرض تمتلك بطبيعتها $\kappa_l$ منخفضة وخصائص إلكترونية مواتية. وتحديداً، لا تزال الخصائص الإلكترونية والبصرية ونقل الفونونات في ثنائي كالكوجينيد الجرمانيوم رباعي الأضلاع (GeS₂ و GeSe₂) في الحالة الصلبة غير مستكشفة بشكل كافٍ من منظور كهروحراري، رغم تباينها الهيكلي وإمكانية فصل النقل الإلكتروني عن نقل الفونونات.

المنهجية
استخدم المؤلفون إطار عمل شاملاً من المبادئ الأولى (first-principles) لاستقصاء الخصائص الهيكلية والإلكترونية والحرارية لمركبات GeS₂ و GeSe₂ رباعية الأضلاع (المجموعة الفراغية $P42/nmc$).

البنية الإلكترونية: أُجريت الحسابات باستخدام نظرية وظيفية الكثافة (DFT) مع دالة PBE وتصحيحات تشتت Grimme D2 للاسترخاء الهيكلي والخصائص المرنة. وللحصول على فجوات نطاق دقيقة، أُجريت حسابات باستخدام دالة الهجين المختلطة HSE03. كما حُسبت معاملات النقل الإلكتروني باستخدام AMSET (نقل بولتزمان من المبادئ الأولى)، مع دمج آليات تشتت تعتمد على الحالة بما في ذلك جهود تشوه الصوت (acoustic deformation potentials)، والفونونات البصرية القطبية، والشوائب المتأينة.
ديناميكا الشبكة والنقل الحراري: تم استقصاء الديناميكا الهارمونية وغير الهارمونية للشبكة باستخدام نظرية الاضطراب الوظيفي للكثافة (DFDT) وطرق الإزاحة المحدودة. وحُسبت الموصلية الحرارية للشبكة ( $\kappa_l$ ) باستخدام كود ShengBTE، عبر حل معادلة بولتزمان لنقل الفونونات ضمن تقريب زمن الاسترخاء (RTA).
تحليل الترابط: تم تحليل الأصل المجهري لسلوك النقل باستخدام LOBSTER لتقييم واصفات "تعداد هاملتون لمدارات البلورات المتكامل" (ICOHP) و"تعداد تداخل مدارات البلورات" (ICOOP)، مما يميز الطبيعة التساهمية-الأيونية لروابط Ge–X.
التقييم الكهروحراري: تم حساب معامل الجدارة ( $zT$ ) من خلال دمج تينسورات النقل الإلكتروني لـ AMSET مع تينسورات الموصلية الحرارية للشبكة لـ ShengBTE، مما يضمن الاتساق بين المكونات (باقتران المعاملات الإلكترونية داخل المستوي مع $\kappa_l$ داخل المستوي، والمعاملات عبر المستوي مع $\kappa_l$ عبر المستوي).

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

الاستقرار الهيكلي والميكانيكي: تم التأكد من أن كلاً من GeS₂ و GeSe₂ رباعي الأضلاع مستقران ميكانيكياً وديناميكياً. وتظهر الهياكل تبايناً قوياً، يتميز بثوابت مرونة كبيرة داخل المستوي ( $C_{11}$ ) وثوابت أصغر بكثير خارج المستوي ( $C_{33}$ ). ومن الجدير بالذكر أن كلا المادتين تظهران قنوات تشوه قص لينة (منخفضة $C_{66}$ )، وهو أمر ملائم لتقليل الموصلية الحرارية للشبكة.
الخصائص الإلكترونية: أسفرت حسابات HSE03/Wannier interpolation عن فجوات نطاق أشباه موصلات تبلغ 2.48 إلكترون فولت لـ GeS₂ و 1.23 إلكترون فولت لـ GeSe₂. وكشف تحليل LOBSTER عن ترابط مختلط (تساهمي-أيوني)، حيث أظهرت روابط Ge–S مقدار ICOHP أكثر استقراراً (−5.27 eV/bond) مقارنة بروابط Ge–Se (−4.74 eV/bond).
الموصلية الحرارية للشبكة متباينة الخواص: كشفت الدراسة عن تباين واضح في نقل الحرارة. عند 300 كلفن، تبلغ الموصلية الحرارية للشبكة داخل المستوي ( $\kappa_{ab}$ ) نحو 26.86 W m⁻¹ K⁻¹ لـ GeS₂ و 18.74 W m⁻¹ K⁻¹ لـ GeSe₂. وفي المقابل، فإن الموصلية عبر المستوي ( $\kappa_c$ ) أقل بشكل جذري، حيث بلغت 1.19 W m⁻¹ K⁻¹ لـ GeS₂ و 1.52 W m⁻¹ K⁻¹ لـ GeSe₂. ويُعزى هذا التثبيط في النقل عبر المستوي إلى الإطار رباعي الأوجه الطبقي، وانخفاض سرعات المجموعة الفونونية خارج المستوي، واللانمطية (anharmonicity) المتوسطة. وقد أدى استبدال S بـ Se إلى خفض تردد الفونون العلوي من ~13.6 THz إلى ~10.3 THz.
الأداء الكهروحراري: من خلال الجمع بين بيانات النقل الإلكتروني والشبكي، كانت قيم $zT$ المحسوبة متوسطة. لوحظ أعلى أداء في النقل من النوع n عبر المستوي لـ GeS₂ عند 800 كلفن، حيث حقق $z_c T$ قدره 0.257 عند تركيز حاملات قدره $1 \times 10^{19}$ cm⁻³. وبالنسبة لـ GeSe₂، يعطي تقدير PBE-AMSET قيمة $z_c T$ تبلغ 0.066 للنقل من النوع p عبر المستوي عند 800 كلفن. ويشير المؤلفون إلى أن نتائج GeSe₂ يجب أن تُعامل كقيم نوعية نظراً لكون فجوة النطاق تقترب من الصفر في وصف PBE، مما قد يتطلب تصحيح الفجوة للحصول على دقة كمية.

الأهمية والادعاءات
يزعم البحث أن مركبات GeX₂ رباعية الأضلاع تمثل فئة من المواد الكهروحرارية ذات التباين القوي، حيث تستفيد الاستجابة عبر المستوي بشكل كبير من تثبيط نقل الحرارة للشبكة. وتوضح الدراسة أنه بينما يعد انخفاض $\kappa_l$ شرطاً ضرورياً، إلا أنه ليس كافياً لتحقيق $zT$ مرتفع؛ إذ يتم التحكم في الأداء في النهاية من خلال التوازن المتزامن بين عامل القدرة، وتدهور الحركية عند تركيزات الحاملات العالية، والنقل الحراري الإلكتروني. ويخلص المؤلفون إلى أن مركبات GeS₂ و GeSe₂ رباعية الأضلاع النقية تظهر قيم $zT$ محسوبة متوسطة، مع كون قناة النقل الأكثر واعدة هي عبر المحور $c$ . ويقترحون أنه يمكن تحقيق تحسينات مستقبلية من خلال حسابات التشتت المصححة بالفجوة، وهندسة العيوب، والسبائك، والتحكم في الإجهاد للاستفادة بشكل أكبر من الاستجابة الحرارية المواتية عبر المستوي.

متاهة البلورات: حكاية هندستين

ازدحام مروري حراري

رقصة الإلكترونات

ماذا يعني هذا (وما لا يعنيه)

الخلاصة

أبحاث مشابهة