Application and Performance Assessment of Annealing Methods for Electrostatic-Energy-Based Configuration Search in Mixed Crystals

ملخص تقني: تقييم أداء وتطبيق طرق التلدين للبحث عن التكوينات القائمة على الطاقة الكهروستاتيكية في البلورات المختلطة

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

المنهجية
يقترح المؤلفون إطار عمل يصيغ البحث عن تكوينات الطاقة الكهروستاتيكية المنخفضة كمسألة تحسين توليفية.

1. الصياغة: يتم تمثيل حالات الإشغال للمواقع الاستبدالية بواسطة متغيرات ثنائية ($\sigma_i = 0$ أو $1$). تم إسقاط طاقة إيوالد على هاميلتوني من نوع إيزينج (Ising-type Hamiltonian). ولفرض نسب استبدال محددة، تمت إضافة حد قيد إلى الهاميلتوني، مما نتج عنه مسألة تحسين مقيدة.
2. الخوارزميات: تم تنفيذ طريقتين للتلدين لحل هذه المسألة:
   • التلدين المحاكي (SA): باستخدام مكتبة dwave-samplers ، استخدمت الدراسة كلاً من جداول التبريد الهندسي والخطي. يستكشف النظام مشهد الطاقة عبر قبول حالات ذات طاقة أعلى باحتمالات تحكمها توزيعات بولتزمان، مع التبريد تدريجياً للتقارب نحو حالات الطاقة المنخفضة.
   • التلدين الكمي (QA): تم تنفيذه على نظام D-Wave Advantage 4.1 QPU. تم دمج المسألة في المعالج الكمي باستخدام AutoEmbeddingComposite. يتطور البحث من هاميلتوني المجال المستعرض إلى هاميلتوني المسألة. تبحث الدراسة في تأثير وقت التلدين وحدوث "كسر السلاسل" (حيث تفشل الكيوبتات الفيزيائية التي تمثل متغيراً منطقياً واحداً في الاصطفاف).
3. التحقق: تم تقييم أداء التلدين المحاكي (SA) والتلدين الكمي (QA) كمياً مقابل عمليات البحث الشاملة لثلاثة أنظمة مستهدفة متزايدة التعقيد:
   • نطاق صغير: $\text{CaYAlO}_4$ (70 تكويناً ممكناً).
   • نطاق متوسط: $\beta\text{-KSbF}_4$ (794 تكويناً ممكناً).
   • نطاق كبير: $\text{Ba-doped SiAlON}$ (127,400 تكويناً ممكناً).

المساهمات والنتائج الرئيسية
تقدم الدراسة مقارنة كمية منهجية لطرق التلدين للفحص المسبق لتكوينات البلورات المختلطة:

• أداء التلدين المحاكي (SA):

  • النطاق الصغير ($\text{CaYAlO}_4$): حقق SA تسريعاً بنحو 26.6 إلى 43.8 مرة مقارنة بالبحث الشامل (اعتماداً على جدول التبريد) ونجح في تحديد جميع الهياكل ذات الطاقة الأدنى دون إغفال.
  • النطاق المتوسط ($\beta\text{-KSbF}_4$): حقق SA تسريعات تتراوح بين 18.8 و248 مرة تقريباً. أظهر التبريد الخطي سرعة فائقة (248 ضعفاً) مقارنة بالتبريد الهندسي (18.8 ضعفاً) مع الاستمرار في التقاط جميع الهياكل ذات الطاقة الأدنى.
  • النطاق الكبير ($\text{Ba-doped SiAlON}$): حقق SA تسريعاً بنحو 286 مرة (التبريد الهندسي) ونجح في التقاط جميع الهياكل ذات الطاقة الأدنى. وعند زيادة عدد العينات إلى 40,000، وصل التسريع إلى 668 مرة، رغم أن ذلك تطلب موارد حاسوبية أكبر بكثير.
  • جداول التبريد: قدم التبريد الخطي عموماً أوقات حساب أسرع ولكنه تطلب ضبطاً دقيقاً لأحجام العينات أو عدد المسحات لتجنب إغفال الهياكل منخفضة الطاقة في الأنظمة الكبيرة. كان التبريد الهندسي أكثر متانة ضد الإغفال ولكنه كان أبطأ حسابياً بسبب العدد الأكبر من عمليات قلب السبين (spin flips) المقبولة لزيادة الطاقة.

• أداء التلدين الكمي (QA):

  • النطاق الصغير ($\text{CaYAlO}_4$): أظهر QA كفاءة عالية، محققاً تسريعاً قدره 116 مرة والتقاط جميع الهياكل ذات الطاقة الأدنى دون إغفال. لم يتم ملاحظة أي كسر في السلاسل.
  • النطاق المتوسط ($\beta\text{-KSbF}_4$): كان أداء QA محدوداً. بينما سمح زيادة وقت التلدين (حتى 2,000 ميكروثانية) بالتقاط جميع الهياكل ذات الطاقة الأدنى، انخفض التسريع إلى عامل قدره 1.29 فقط مقارنة بالبحث الشامل. والأهم من ذلك، حدثت كسور في السلاسل في 5.55% إلى 33% من الحلول. أظهرت الدراسة وجود علاقة طردية بين كسر السلاسل ومتوسط طاقة الحلول، مما يشير إلى أن كسور السلاسل تعيق البحث عن حالات الطاقة المثلى.
  • النطاق الكبير ($\text{Ba-doped SiAlON}$): لم يكن من الممكن إجراء QA لهذا النظام بسبب قيود التضمين.

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

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