Resource-efficient quantum-selected configuration interaction for molecular properties
تقترح هذه الورقة إطار عمل للتفاعل التشكلي المختار كمومياً وفعال في استهلاك الموارد، يقوم ببناء هاميلتوني مدمج لتقليل تعقيد الدائرة بشكل كبير، مما يتيح حسابات دقيقة للخصائص الجزيئية مثل طاقات الحالة الأرضية وعزوم ثنائية القطب لمركبات أحادي فلوريد المجموعة IIIA على أجهزة كمومية ذات مقياس متوسط مشوبة بالضجيج وتصل إلى 20 كيوبت.
المؤلفون الأصليون:Suprava Sahoo, Abdul Kalam, Kenji Sugisaki, V. S. Prasannaa, B. P. Das
تخيل أنك تحاول حل لغز "بازل" كوني عملاق، حيث كل قطعة فيه هي جسيم طاقة صغير غير مرئي، والصورة التي تحاول إكمالها هي كيفية سلوك جزيء ما. هذا هو عالم الكيمياء الكمومية، وهو مجال يستخدم فيه العلماء القواعد الغريبة لعالم الكم للتنبؤ بكيفية تلاحم الذرات مع بعضها البعض. لعقود من الزمن، كان "المعيار الذهبي" لحل هذه الألغاز هو طريقة تسمى "التفاعل التكويني الكامل" (Full Configuration Interaction - FCI). فكر في (FCI) على أنها محاولة للنظر في كل طريقة ممكنة لتوافق قطع اللغز معاً للوصول إلى الصورة المثالية الوحيدة. المشكلة؟ كلما أصبح الجزيء أكبر، انفجر عدد الترتيبات الممكنة مثل حبات الفشار في مقلاة ساخنة، ليصبح ضخماً لدرجة أن أسرع الحواسيب الفائقة في العالم لا تستطيع معالجة هذه الأرقام خلال عمر كامل.
هنا يأتي دور الحاسوب الكمومي، وهو آلة بُنيت لتتحدث اللغة الأصلية لهذه الجسيمات. فبدلاً من محاولة حساب كل الاحتمالات باستخدام حاسوب عادي، يمكن للحاسوب الكمومي محاكاة الجزيء مباشرة. ومع ذلك، فإن هذه الآلات تمر حالياً بـ "سنوات مراهقتها"؛ فهي قوية ولكنها "ضوضائية"، مما يعني أنها ترتكب الأخطاء بسهء إذا طلبت منها القيام بالكثير من المهام في وقت واحد. وللحصول على إجابة واضحة، يحتاج العلماء إلى تبسيط المشكلة دون فقدان التفاصيل المهمة. وهنا يكمن التحدي: كيف تخبر حاسوباً كمومياً ضوضائياً أن يركز فقط على قطع اللغز الأكثر أهمية حتى لا يغمره الأمر ويعطيك نتيجة مشوشة؟
هذا بالضبط ما سعى الباحثون في هذه الورقة البحثية لحله. فقد طوروا استراتيجية جديدة ذكية تسمى "eos-QSCI" (التفاعل التكويني المختار الكمومي القائم على فرز مؤثرات باولي الفعالة). تخيل أنك طاهٍ يحاول صنع حساء مثالي، لكن ليس لديك سوى قدر صغير ووقت محدود. بدلاً من رمي كل التوابل الموجودة في خزانة المؤن (والذي سيكون نهج "الهاملتوني الكامل")، تقوم أولاً بتذوق المرق لترى أي التوابل تُحدث فرقاً فعلياً. ثم تستبعد تلك التي لا تغير النكهة كثيراً وتطهو فقط بالتوابل الأساسية القليلة.
طبق المؤلفون طريقة "تذوق الاختبار" هذه على مجموعة من الجزيئات تسمى "أحادية فلوريد المجموعة IIIA"، والتي تشمل فلوريد البورون (BF)، وفلوريد الألمنيوم (AlF)، وفلوريد الثاليوم (TlF) الثقيل والنسبي. ووجدوا أنه من خلال تحديد والاحتفاظ فقط بـ "مؤثرات الاستثارة" الأكثر هيمنة (وهي المعادل الكمومي لأهم التوابل)، استطاعوا تقليص حجم المشكلة بشكل هائل. فبالنسبة لنظام مكون من 20 كيوبت لجزيء (TlF)، قللت طريقتهم عدد الخطوات المطلوبة في الدائرة الكمومية بأكثر من 98%. إنه يشبه اختزال ماراثون بطول 100 ميل إلى ركض لمسافة ميلين، مع الوصول في النهاية إلى نفس النتيجة.
وعندما اختبروا ذلك على حاسوب كمومي حقيقي (معالج Marrakesh من IBM)، كانت النتائج مبهرة. فحتى مع الضوضاء الطبيعية للآلة، أنتج نهجهم المبسط طاقات الحالة الأرضية وعزوم ثنائية القطب الكهربائية (وهي مقياس لكيفية توزيع شحنة الجزيء) بدقة تتجاوز 99.99% مقارنة بالقيم النظرية المثالية. في الواقع، بالنسبة لجزيء (TlF) الثقيل، احتاجوا فقط إلى عينة تبلغ حوالي 0.88% من إجمالي التكوينات الممكنة للحصول على نتيجة شبه مثالية. وهذا يشير إلى أنه من خلال التعامل بذكاء مع الأجز much التي يجب التركيز عليها في اللغز الكمومي، يمكننا إجراء محاكاة كيميائية معقدة على آلات اليوم غير المثالية، مما يفتح الباب لفهم كل شيء من المواد الجديدة إلى الفيزياء الأساسية دون الحاجة إلى حاسوب مستقبلي خالٍ من العيوب.
ملخص تقني: التفاعل التكويني المختار كمومياً بكفاءة في الموارد لخصائص الجزيئات
بيان المشكلة تواجه المحاكاة الكمومية للبنى الإلكترونية الجزيئية تحديات كبيرة في القابلية للتوسع على الأجهزة الكمومية ذات المقياس المتوسط ذات الضجيج (NISQ). وبينما تقدم طريقة التفاعل التكويني المختار كمومياً (QSCI) نهجاً هجيناً واعداً يجمع بين الحوسبة الكمومية والكلاسيكية من خلال تحديد دالات الأساس المحددات الهامة عبر التطور في الوقت الفعلي، فإن تنفيذها باستخدام الهاملتوني الإلكتروني الكامل يؤدي إلى زيادة سريعة في تعقيد الدارات. وتحديداً، تتطلب طريقة QSCI القائمة على محاكاة الهاملتوني (HSB-QSCI) عدداً من حدود "باولي" يتناسب مع O(N4) وبوابات ثنائية الكيوبت تتناسب مع O(N5) لنظام مكون من N كيوبت تحت تحويل "جوردان-ويجنر". هذا الطلب العالي على الموارد يحد من إمكانية تطبيق الطريقة على الأنظمة الصغيرة ويعيق حساب الخصائص الجزيئية التي تتجاوز طاقات الحالة الأرضية، مثل العزوم الكهربائية ثنائية القطب الدائمة (PDMs)، لا سيما في الأنظمة النسبية حيث تكون حدود الهاملتوني أكثر عدداً.
المنهجية: eos-QSCI يقترح المؤلفون إطار عمل التفاعل التكويني المختار كمومياً القائم على فرز مؤثرات باولي الفعالة (eos-QSCI) المصمم لتقليل تعقيد الهاملتوني قبل المحاكاة الكمومية. تمر الطريقة بثلاث مراحل متميزة:
المرحلة الأولى: الفرز الفيرميوني: يتم فرز الهاملتوني ذو التكميم الثاني من خلال فحص تأثير سلاسل المؤثرات الفيرميونية على المحدد المرجعي (عادةً حالة "ديراك-هارتري-فوك"). يتم استبعاد المؤثرات التي تلغي الحالة المرجعية (مثل تلك التي تتضمن الإلغاء في المدارات الافتراضية). يفرض هذا قواعد اختيار قائمة على الإشغال، مما يقلل عدد الحدود قبل رسم الخرائط إلى الكيوبتات عبر تحويل "جوردان-ويجنر".
المرحلة الثانية: الاختيار القائم على فقدان الوفاء: يتم ترتيب مؤثرات "باولي" المتبقية حسب أهميتها باستخدام مقياس فقدان وفاء الحالة المرجعية. لكل مؤثر Pα، يتم تحليل التطور قصير المدى ∣Ψα(t)⟩=e−ihαPαt∣Φ0⟩. يعمل فقدان الوفاء ωα(t)، المتناسب مع تباين المؤثر بالنسبة للحالة المرجعية، كوزن. يتم ترتيب المؤثرات بترتيب تنازلي وفقاً لـ ωα، ويتم حساب وزن تراكمي معياري. يتم اقتطاع الهاملتوني للاحتفاظ فقط بأعلى K من الحدود بحيث يتجاوز الوزن التراكمي حداً معيناً (مثلاً 0.99).
المرحلة الثالثة: التنفيذ والمعالجة اللاحقة: يتم تقسيم الهاملتوني المقتطع إلى مجموعات من الملاحظات التبادلية (SCOs) لبناء دوائر تطور زمني فعالة. يتم تحسين هذه الدوارات باستخدام روتين خطوط الأنابيب (Qiskit و Pytket) وتنفيذها على الأجهزة الكمومية. تحدد التكوينات التي تم أخذ عينات منها فضاء هاملتوني فرعي مختزل، والذي يتم قطره كلاسيكياً للحصول على طاقات الحالة الأرضية والدوال الموجية. بعد ذلك، يتم تقييم العزوم الكهربائية ثنائية القطب الدائمة (PDMs) من هذه الدوال الموجية.
المساهمات الرئيسية
اختزال الهاملتوني: يبني إطار العمل هاملتونياً فعالاً ومدمجاً من خلال تحديد مؤثرات الإثارة الفيرميونية المهيمنة، مما يقلل بشكل كبير من عدد حدود "باولي"، وعمق الدارة، وعدد بوابات الكيوبت الثنائية مقارنة بـ HSB-QSCI كامل الهاملتوني.
الخصائص الجزيئية النسبية: توسع الطريقة نطاق QSCI لتشمل مشكلات البنية الإلكترونية النسبية، حيث نجحت في حساب كل من طاقات الحالة الأرضية والعزوم الكهربائية ثنائية القطب الدائمة (PDMs) لمونوفلوريدات المجموعة IIIA (BF, AlF, GaF, InF, TlF).
تحليل القابلية للتوسع: يوضح المؤلفون أن متطلبات الموارد لـ eos-QSCI تتوسع بشكل أفضل من نهج الهاملتوني الكامل، محققة تحسناً يقترب من التربيع في مقياس حدود الهاملتوني (O(N2.3) مقابل O(N4)).
النتائج تم التحقق من صحة إطار العمل من خلال عمليات المحاكاة والتنفيذ على الأجهزة باستخدام معالج IBM Marrakesh (بمعمارية Heron r2) لأنظمة AlF و TlF مع مساحات نشطة تصل إلى 20 كيوبت.
كفاءة الموارد: بالنسبة لنظام TlF المكون من 20 كيوبت، حقق الهاملتوني المختزل انخفاضاً أكبر من 98% في كل من عمق الدارة وعدد بوابات الكيوبت الثنائية مقارنة بنهج الهاملتوني الكامل. كما انخفض عدد حدود الهاملتوني بنسبة 97.14% (لـ AlF) و 98.00% (لـ TlF).
الدقة:
المحاكاة: انحرفت طاقات الحالة الأرضية المحسوبة عن قيم التفاعل التكويني للمساحة النشطة الكاملة (CASCI) بأقل من 1.6 mHa (عتبة الدقة الكيميائية) عبر جميع أحجام الأنظمة. تراوحت طاقات الارتباط المستردة بين 90.02% و 98.46%. وكانت انحرافات الـ PDM ضمن نطاق 1.49×10−2 a.u. لـ AlF و 7.0×10−3 a.u. لـ TlF.
التنفيذ على الأجهزة: في عملية حساب TlF المكون من 20 كيوبت على معالج IBM Marrakesh، استردت الطريقة 99.93% من طاقة ارتباط CASCI. وطابقت العزوم الكهربائية ثنائية القطب (PDMs) الناتلة قيم CASCI بنسبة 99.99%. وحتى مع ضجيج الأجهزة، ظلت انحرافات الطاقة أقل من 0.01% بالنسبة لـ CASCI.
اختزال فضاء المحددات: ركزت الطريقة أخذ العينات على المحددات المهيمنة، مما قلل من تكلفة القطر الكلاسيكي. بالنسبة لنظام TlF المكون من 20 كيوبت، تطلب الأمر حوالي 0.88% فقط من إجمالي فضاء محددات CASCI لتحقيق دقة عالية، مما قلل حجم مصفوفة القطر من 181×181 إلى 43×43.
الأهمية يزعم البحث أن إطار عمل eos-QSCI يحسن بشكل كبير من عملية تطبيق HSB-QSCI على أجهزة NISQ. ومن خلال تقليل تعقيد الهاملتوني، وموارد الدائرة الكمومية، وتكاليف المعالجة اللاحقة الكلاسيكية في آن واحد، يتيح هذا النهج محاكاة فعالة للموارد للأنظمة الجزيئية النسبية التي كانت غير قابلة للتطبيق سابقاً بسبب قيود عمق الدارة. إن النجاح في التنفيذ على الأجهزة باستخدام 20 كيوبت يثبت قابلية التوسع لهذا النهج لحساب كل من الطاقات والخصائص الجزيئية مثل الـ PDMs، والتي تعد بالغة الأهمية للتطبيقات في التقنيات الكمومية والفيزياء الأساسية.