🔬 materials science

Mapping the influence of symmetry breaking in structure-property relationships of ABO $_3$ perovskites

تتناول هذه الورقة التحديات الحسابية في نمذجة هياكل البيروفسكايت المشوهة من خلال تطوير إطار عمل فعال باستخدام قوالب الخلايا الفائقة المقيدة وغير المقيدة بالتناظر لاستكشاف التشوهات المعتمدة على التركيب وتأثيرها على استقرار وخصائص المواد الفيزيائية لمركبات ABO $_3$ بيروفسكايت بشكل منهجي.

المؤلفون الأصليون: Panupol Untarabut, Sylvian Cadars, Fabien Pascale, Sébastien Lebègue, Olivier Masson, Samuel Bernard, Assil Bouzid, Santanu Saha

نُشر 2026-07-31

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

CC BY 4.0

المؤلفون الأصليون: Panupol Untarabut, Sylvian Cadars, Fabien Pascale, Sébastien Lebègue, Olivier Masson, Samuel Bernard, Assil Bouzid, Santanu Saha

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

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

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

المبارزة الكبرى للميل الذري

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

في هذه الدراسة، عمل فريق من العلماء كمدربي رقص لهذه الرقصة الذرية. أرادوا معرفة أي حركات الرقص (التشوهات) تجعل المادة أكثر استقراراً وأيها يمنحها أفضل خصائص توفير الطاقة. وللقيام بذلك، لم يكتفوا بالتخمين؛ بل بنوا ساحة لعب رقمية ضخمة. لقد قاموا بمحاكاة 15 "وصفة" مختلفة من البيروفسكايت، عبر خلط ومطابقة ذرات مختلفة في المواقع (A) و(B) في الصيغة الكيميائية $ABO_3$ .

قواعد ساحة الرقص

اكتشف الباحثون أن "أسلوب الرقص" يعتمد بشكل كبير على رقم بسيط يسمى عامل التسامح غولدشميت ( $t_G$ ). يمكنك التفكير في $t_G$ كمقياس لمدى ملاءمة الراقصين لأماكنهم المخصصة.

الضغط الضيق ( $t_G < 1.00$ ): عندما يكون الراقصون أصغر قليلاً من المساحة، يضطرون للتجمع معاً. هذا يجبر الهيكل بأكمله على الميل والدوران بجنون لملء الفجوات. وجدت الدراسة أنه بالنسبة لهذه المواد ذات "الضغط الضيق"، فإن حركة الرقص الأكثر استقراراً وطاقة هي نوع معين من الميل يسمى نمط $a^+b^-b^-$ (والذي يترجم إلى شكل معيني الوجوه/أورثورومبي). الأمر يشبه اتفاق الراقصين جميعاً على الميل في نمط متناسق ومتعرج للحفاظ على التوازن.
الملاءمة الواسعة ( $t_G > 1.00$ ): عندما يكون الراقصون أكبر من المساحة، لا يحتاجون للميل كثيراً. بدلاً من ذلك، يفضلون التذبذب بشكل فردي. في هذه الحالات، تظل المادة مكعبة (مربعة) في الغالب، لكن الذرة المركزية (الموقع B) تنزاح قليلاً عن المركز، مثل راقص يخطو خطوة خارج الخط ليحضر غرضاً ما.

لعبة المحاكاة

لإيجاد هذه الإجابات، لعب الفريق لعبة عالية المخاطر من "الاستكشاف الهيكلي". استخدموا استراتيجيتين رئيسيتين:

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

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

ما يعنيه هذا للمستقبل

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

انظر إلى المكونات (احسب $t_G$ ).
توقع مقدار ميل الذرات بناءً على ذلك الرقم.
ابنِ صندوقاً بحجم محدد ( $\sqrt{2} \times \sqrt{2} \times 2$ ) وهو مصمم بدقة لاستيعاب هذا الميل.
اترك الكمبيوتر يريح الهيكل.

هذه "الوصفة الذكية" وجدت باستمرار الهياكل ذات الطاقة الأدنى والأكثر استقراراً لجميع المواد التي اختبروها، متفوقة على طرق التخمين العشوائي.

الصورة الأكبر

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

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

الخلاصة

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

ملخص تقني: رسم خرائط تأثير كسر التماثل في علاقات البنية والخواص لمركبات البيروفسكايت من نوع ABO₃

بيان المشكلة
تعد أكاسيد البيروفسكايت (ABO₃) بالغة الأهمية لتطبيقات الطاقة نظرًا لمرونتها التركيبية، مما يسمح باستقرار كل من الأطوار عالية التماثل (المكعبة) والأطوار منخفضة التماثل. وبينما تهيمن الأطوار منخفضة التماثل، التي تتميز بتميل ثماني الأوجه وإزاحات موقع B، في الظروف المحيطة، ركزت الدراسات الحسابية تاريخيًا على الطور المكعب المثالي. وتظل نمذجة التشوهات في الأطوار منخفضة التماثل تحديًا بسبب غياب إطار كمي قادر على استيعاب التشوهات المعتمدة على التركيب والناشئة عن أنماط متعددة من الإزاحات الذرية غير المتكافئة. إن النهج الحالي، مثل طريقة الفونون اللين (soft phonon method) أو استرخاء الخلايا الفائقة الكبيرة غير المقيدة، مكلف حاسوبيًا. علاوة على ذلك، ورغم وجود قوالب بنيوية مكسورة التماثل (مثل قوالب جليزر - Glazer templates)، إلا أن تبنيها على نطاق واسع يعوقه غياب التحقق المنهجي عبر تركيبات متنوعة وعدم وضوح كيفية ارتباط أنماط التشوه بالاستقرار الديناميكي الحراري والخواص الفيزيائية.

المنهجية
قام المؤلفون بتطوير إطار عمل حسابي لاستكشاف المشهد البنيوي بشكل منهجي لـ 15 تركيبة ممثلة لبيروفسكايت ABO₃ (حيث A = Ca, Sr, Ba؛ و B = Ti, Zr, Hf, Sn, Nb). استخدمت الدراسة نظرية الدوال الوظيفية للكثافة (DFT) باستخدام دالة PBE وجهود ONCV الكاذبة (ONCV pseudopotentials). تضمنت المنهجية ثلاث فئات رئيسية من القوالب البنيوية:

القوالب المقيدة بالتماثل:
- قوالب جليزر (Glazer Templates): 15 نظام تميل فريد يمثل مجموعات من دورات ثمانية الأوجه في الطور (in-phase - IR) وخارج الطور (out-of-phase - OR) على طول المحاور البلورية.
- القوالب القياسية: خمسة أطوار مرجعية (مكعبة، رباعية الأضلاع I/II، معينة، سداسية) تمثل الأطوار المعروفة والتشوهات المحددة مثل إزاحة موقع B.
القوالب غير المقيدة:
- قوالب الخلايا الفائقة: خلايا فائقة كبيرة (تختلف في الحجم والشكل والزوجية - مثل $2\times2\times2$ ، $3\times2\times2$ ، $\sqrt{2}\times\sqrt{2}\times1$ القائمة على) مع اضطرابات ذرية عشوائية لأخذ عينات من سطح طاقة الوضع دون قيود التماثل.
- القوالب المشوهة: نماذج محددة تعزل أنواع التشوه الفردية (IR، أو OR، أو إزاحة موقع B) لقياس مساهماتها الفردية.
التحليل والتحقق:
- تم إجراء التصنيف الديناميكي الحراري عند درجة حرارة 0 كلفن عبر حسابات الطاقة النسبية ( $\Delta E$ ).
- تم حساب الطاقة الحرة المعتمدة على درجة الحرارة ( $\Delta G(T)$ ) باستخدام التقريب شبه التوافقي (الفونونات) لتركيبات مختارة (CaZrO₃, SrTiO₃, BaTiO₃).
- تم تعريف واصفات بنيوية لقياس تميل ثماني الأوجه بين الأوجه ( $\langle \sin(\theta_{tilt}) \rangle$ ) والتشوهات داخل ثماني الأوجه (إزاحة موقع B $\langle \Delta r_B \rangle$ ومعلم الترتيب $\langle q \rangle$ ).
- تم حساب الخواص الفيزيائية (فجوة HOMO-LUMO $E_g$ ومعامل المرونة الحجمي $B_0$ ) لمختلف القوالب.

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

الاستقرار الديناميكي الحراري و $t_G$ : يوجد ارتباط كمي قوي بين عامل تولرانس جولدشميت ( $t_G$ $t_{G}$ ) وطاقة الاستقرار للأطوار المشوهة.
- $t_G$ المنخفض ( $<1.00$ ): الأنظمة مثل CaZrO₃ تظهر طاقات استقرار كبيرة (تصل إلى $\sim$ 200 ملي إلكترون فولت/ذرة) بالنسبة للطور المكعب. الحالة الأرضية هي باستمرار طور المعين $Pnma $(جليزر$ a^+b^-b^-$)، مدفوعًا بشكل أساسي بدورات ثمانية الأوجه الجماعية.
- $t_G$ المرتفع ( $\ge 1.00$ ): الأنظمة مثل SrTiO₃ و BaTiO₃ تظهر فروق طاقة ضئيلة ( $<15$ ملي إلكترون فولت/ذرة) بين الأطوار. بالنسبة لـ BaTiO₃، تصبح إزاحة موقع B هي آلية الاستقرار المهيمنة، بينما يظل SrTiO₃ قريبًا من الحد المكعب مع عدم استقرار دوراني ضعيف.
دور زوجية الخلية الفائقة: يؤثر اختيار حجم وزوجية الخلية الفائقة بشكل كبير على القدرة على التقاط حالات الطاقة المنخفضة.
- الخلايا الفائقة ذات الزوجية الزوجية (مثل $2\times2\times2$ ، $\sqrt{2}\times\sqrt{2}\times2$ ) تلتقط بفعالية أنماط الدوران المناهضة للطور المتماسكة (مثل $Pnma$) وتصل إلى الحد الأدنى العالمي.
- الخلايا الفائقة ذات الزوجية الفردية (مثل $3\times2\times2$ ) تعيق الدورات المثالية خارج الطور، مما يؤدي إلى أنماط تميل مختلطة ( $\pm$ ) وطاقات أعلى، مما يفشل في الوصول إلى الحالة الأرضية الحقيقية في الأنظمة ذات $t_G$ المنخفض.
تأثيرات درجة الحرارة:
- CaZrO₃: يظل طور المعين $Pnma$ هو الحالة الأرضية حتى 2000 كلفن؛ لا يتم التنبؤ بأي انتقال طوري ضمن التقريب التوافقي.
- SrTiO₃: مشهد الطاقة الحرة مسطح للغاية. الطور المكعب مستقر عند درجات الحرارة المنخفضة في هذا الحساب (على عكس بعض التقارير التجريبية التي تشير إلى انتقال عند 105 كلفن)، حيث تظل الأطوار منخفضة التماثل متقاربة جدًا في الطاقة.
- BaTiO₃: يُتوقع حدوث انتقال من الرباعي الأضلاع إلى المكعب حول 880 كلفن (التقريب التوافقي)، مدفوعًا بالإنتروبيا الاهتزازية، رغم أن هذه الدرجة أعلى من القيمة التجريبية (398 كلفن)، مما يشير إلى الحاجة إلى تصحيحات لا توافقية.
علاقات البنية والخواص:
- فجوة النطاق ( $E_g$ ): تزداد عمومًا مع انخفاض التماثل. تضيق الفجوة بين المكعب والمعين مع زيادة $t_G$ . إزاحة موقع B (مثل $P4mm $) تغير بشكل كبير$ E_g$ في مركبات Ti/Sn.
- معامل المرونة الحجمي ( $B_0$ ): يظهر الطور المعين عادةً أعلى $B_0$ ، بينما يظهر الطور المكعب أدنى قيمة. تميل التشوهات التي تتضمن إزاحة موقع B إلى تقليل $B_0$ .

بروتوكول التحسين المقترح
بناءً على النتائج، يقترح المؤلفون بروتوكولًا قويًا ومنخفض التكلفة حاسوبيًا للاستكشاف البنيوي السريع:

حساب $t_G$ لتقدير ثابت الشبكة الأولي.
إنشاء خلية فائقة $2\times2\times2$ من الطور المكعب المثالي.
تطبيق نمط دوران مستهدف (تحديدًا $+--$ أو $a^+b^-b^-$ ) مع زاوية ميل أولية مشتقة من نموذج خطي بناءً على $t_G$ .
إدخال إزاحات عشوائية صغيرة لمواقع A و B (0.2 أنجستروم) لبذر إزاحة موقع B.
إجراء استرخاء بنيوي كامل.
هذا سير العمل يحدد باستمرار الأطوار الأدنى طاقة عبر الـ 15 تركيبة، متفوقًا على أخذ العينات العشوائية وقوالب جليزر القياسية من حيث الكفاءة والموثوقية، خاصة في الأنظمة ذات $t_G$ المنخفض.

الأهمية والادعاءات
يدعي البحث تقديم إطار حسابي موحد يسد الفجوة بين النماذج عالية التماثل المثالية والواقع المعقد لمركبات البيروفسكايت المشوهة منخفضة التماثل. ويؤكد المؤلفون ما يلي:

الفهم المنهجي: توضح الدراسة الأدوار المحددة لكل من التميل بين الأوجه (inter-octahedral tilting) مقابل الإزاحات داخل الأوجه (intra-octahedral displacements)، وتظهر التحول من الاستقرار المدفوع بالدوران (عند $t_G$ منخفض) إلى الاستقرار المدفوع بالإزاحة (عند $t_G$ مرتفع).
التقدم المنهجي: يقدم البروتوكول المقترح حلاً عمليًا وعالي الإنتاجية لتوليد هياكل منخفضة التماثل ذات صلة فيزيائية، متجاوزًا الاختناقات الحسابية لأخذ العينات العشوائية للخلايا الفائقة والقيود المفروضة على قوالب التماثل الموجودة.
تصحيح قواعد البيانات: من خلال تمكين التوليد الفعال للنماذج منخفضة التماثل، يهدف هذا العمل إلى تخفيف التحيز الحالي في قواعد بيانات المواد مفتوحة المصدر نحو الأطوار المكعبة المثالية، مما يثري مجموعة النماذج ذات الصلة طاقيًا لعمليات الفحص والتصميم المستقبلية.
التنبؤ بالخواص: يوضح العمل أن التشوهات التي تحكم الاستقرار الديناميكي الحراري تتحكم أيضًا بشكل مباشر في الخواص الإلكترونية والميكانيكية، مما يشير إلى أن النمذجة الدقيقة لكسر التماثل ضرورية للتنبؤ بأداء المواد في التطبيقات الواقعية.

أبحاث مشابهة