Simulation of additive binding energies in asphalt using quantum-selected configuration interaction (QSCI)

ملخص تقني: محاكاة طاقات الترابط المضافة في الأسفلت باستخدام التفاعل التكويني المختار كموميًا (QSCI)

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

المنهجية
تقدم الدراسة QuantumPave، وهو سير عمل هجين (كمي-كلاسيكي) مصمم لحساب طاقات ترابط المواد المضافة. تدمج المنهجية ثلاثة نهج متكاملة:

1. اختيار النظام: اختار المؤلفون معقداً من 24 ذرة (بيريدين-فينول) ($C_{11}H_{11}NO$) كنموذج تمثيلي. يتميز هذا النظام برابطة هيدروجينية قوية (O-H$\cdots$N) وأنظمة إلكترونية مترافقة من نوع $\pi$، مما يحاكي المجموعات الوظيفية التي تم تحديدها في دراسات الديناميكا الجزيئية لتقادم الأسفلت.
2. تحسين الهندسة (Geometry Optimization): استُخدمت جهود التفاعل بين الذرات القائمة على التعلم الآلي (ORB v3) لتحسين الهندسة بسرعة، مما وفر توازناً بين الكفاءة الحسابية والدقة الكيميائية.
3. حساب طاقة الترابط: تم استخدام ثلاث طرق لحساب طاقات الترابط:
   • التعلم الآلي (ORB v3): استُخدم للفحص السريع والتقديرات الأولية.
   • نظرية الدالة الكثافية (DFT): عملت حسابات B3LYP-D3/def2-SVP كمرجع كلاسيكي، مع دمج تصحيحات التشتت لتفاعلات التراص $\pi$-$\pi$ والروابط الهيدروجينية.
   • التفاعل التكويني المختار كمومياً (QSCI): يُعرف أيضاً بـ "التشخيص المعتمد على العينات" (SQD)، وهو النهج المستوحى من الكم والأساسي في هذه الدراسة. يستخدم تقريب الفضاء النشط (10 إلكترونات في 10 مدارات، مُسقطة على 20 كيوبت) لتقليل التعقيد مع استهداف التفاعلات التي تهيمن عليها التكافؤات.
4. التنفيذ الكمي: نُفذت خوارزمية QSCI على معالج فائق التوصيل IQM Emerald بسعة 54 كيوبت. يتضمن سير العمل ما يلي:
   • أخذ العينات (Sampling): يقوم المعالج الكمي بأخذ عينات من التكوينات الإلكترونية المهيمنة (سلاسل البتات/bitstrings) من الحالة المُعدة.
     extrapolation.
   • تخفيف الخطأ: بدلاً من استخدام "الاستقراء الخالي من الضجيج" (zero-noise extrapolation)، تستخدم الطريقة استعادة التكوين المتسق ذاتياً. هذه العملية التكرارية تقوم بإصلاح التكوينات غير الصالحة (التي تنتهك تناظر عدد الجسيمات) باستخدام متوسط إشغالات المدارات من العينات المحتفظ بها.
   • التشخيص (Diagonalisation): تقوم موارد الحوسبة عالية الأداء (HPC) الكلاسيكية بعملية التشخيص داخل الفضاء الفرعي المحدد.
   • تصحيح التشتت: تمت إضافة تصحيحات D3BJ التجريبية إلى نتائج QSCI لمراعاة تفاعلات "فان دير فالس" بعيدة المدى التي لا تغطيها البنية الإلكترونية للفضاء النشط وحدها.

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

• دقة الأجهزة: على معالج IQM Emerald، أعادت طريقة SQD مرجع (CASCI) الخاص بالفضاء النشط الكلاسيكي بدقة تامة. وقد أعطى كلاهما طاقة ترابط بلغت -3.52 كيلو سعرة/مول (-0.153 إلكترون فولت).
• المتانة تجاه الضجيج: بدلاً من أن يؤدي ضجيج الجهاز إلى تدهور النتيجة، فقد ساهم في توسيع توزيع أخذ العينات ليشمل الفضاء النشط بالكامل. وبناءً على ذلك، استعاد التشخيص الكلاسيكي طاقة الفضاء النشط بدقة دون الحاجة إلى تقنية الاستقراء الخالي من الضجيج.
• المقارنة بالتجربة: قيمة الفضاء النشط المحسوبة (-3.52 كيلو سعرة/مول) تُظهر ترابطاً أقل من المحتوى الحر الحراري التجريبي الذي يبلغ حوالي -6.25 كيلو سعرة/مول. ويعزو المؤلفون هذا التفاوت إلى الطبيعة المدمجة للفضاء النشط، الذي يلتقط الارتباطات التكافلية الداخلية ولكنه يغفل المساهمات من الفضاء المداري الأوسع المشارك في الرابطة الهيدروجينية الكهروستاتيكية.
• المقارنة بـ DFT: بالغت طريقة B3LYP-D3 في تقدير طاقة الترابط (-9.08 كيلو سعرة/مول)، وهو ما يتوافق مع الميول المعروفة للدوال الهجينة للمبالغة في تقدير الروابط الهيدروجينية المتوسطة.

المساهمات الرئيسية

• إثبات المفهوم للحوسبة الفائقة المرتكزة على الكم: تُثبت الورقة نموذجاً قابلاً للتكرار بالكامل وصغير النطاق للحوسبة الفائقة المرتكزة على الكم، حيث يقوم المعالج الكمي (QPU) بأخذ عينات من التكوينات بينما يقوم الحاسوب الكلاسيكي (HPC) بعملية التشخيص.
• التحقق من صحة QSCI/SQD على الأجهزة: توفر الدراسة دليلاً تجريبياً على قدرة QSCI على تحقيق دقة عددية مقابل حد الفضاء النشط الكلاسيكي على الأجهزة الكمية الحالية ذات الضجيج (NISQ).
• سير عمل مقاوم للضجيج: تسلط الدراسة الضوء على سير عمل يكون فيه ضجيج الجهاز مفيداً للتقارب عبر توسيع التوزيع المأخوذ، مما يلغي الحاجة إلى تقنيات معقدة مثل الاستقراء الخالي من الضجيج في هذا النطاق المحدد.
• التطبيق في علم المواد: نجحت الدراسة في تطبيق هذه الخوارزميات الكمية على مشكلة ذات صلة بصناعة الأسفلت، مما يسد الفجوة بين تجارب الحوسبة الكمية ومشكلات المواد الصناعية.

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

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

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