Efficient Operator Selection and Warm-Start Strategy for Excitations in Variational Quantum Eigensolvers

ملخص تقني لورقة بحثية بعنوان: "الاختيار الفعال للمؤثرات واستراتيجية البدء الدافئ للإثارات في خوارزميات القيم الذاتية الكمومية التباينية" للباحث ماكس هاس وآخرون.

بيان المشكلة

تعد خوارزمية القيم الذاتية الكمومية التباينية (VQE) خوارزمية هجينة رائدة لحل مشكلات الحالة الأرضية في الأنظمة الكمومية متعددة الأجسام، ومع ذلك فهي تواجه عقبات كبيرة في تحقيق التفوق الكمومي. وتتمثل التحديات الرئيسية في:

• الهضاب القاحلة (Barren Plateaus): تلاشي التدرجات بشكل أسي في دالة الخسارة، مما يعيق التقارب.
• عمق الدائرة (Circuit Depth): تتوسع الـ ansätze التقليدية مثل (UCCSD) بشكل متعدد الحدود مع عدد الإلكترونات، مما يؤدي إلى أعماق دوائر كبيرة جداً بالنسبة للأجهزة الكمومية ذات الضجيج والقياس المتوسط (NISQ) الحالية.
• الأعباء الإضافية لاختيار المؤثرات (Operator Selection Overhead): تخفف الاستراتيجيات التكيفية مثل ADAPT-VQE من مشكلات العمق عبر إضافة مؤثرات بشكل تكراري، لكنها تعاني من تكاليف حسابية عالية بسبب التقييم المتسلسل لأعداد هائلة من الدوائر الكمومية أثناء عملية الاختيار.
• الحساسية للتهيئة (Initialization Sensitivity): غالباً ما توفر التهيئة القياسية (مثل المعاملات الصفرية) أو الطرق الاضطرابية مثل نظرية مۆلر-بليسيت من الدرجة الثانية (MP2) نقاط بداية دون المستوى الأمثل، خاصة في الأنظمة ذات الارتباط القوي.

المنهجية

يقترح المؤلفون بروتوكولاً جديداً يجمع بين مُحسّن ExcitationSolve وترتيب الطاقة (ES) لبناء ansätze لـ VQE فعالة وذات بدء دافئ.

1. ExcitationSolve وترتيب الطاقة (ES):

• ExcitationSolve: هو مُحسّن يعيد بناء دالة التكلفة لكل معامل، مما يسمح بالتحسين المباشر دون الحاجة إلى التدرج النازل (gradient descent). وهو يحدد قيمة المعامل المثلى لمؤثر ما عن طريق تقييم مشهد الطاقة عند نقاط قليلة.
• استراتيجية ترتيب الطاقة (ES): بدلاً من اختيار المؤثرات بشكل متسلسل، يقوم ES بتقييم تأثير الطاقة ($\Delta E_j = E_{ref} - E_j$) لجميع المؤثرات في مجموعة (pool) واحدة في وقت واحد مقابل حالة مرجعية. ويتم اختيار المؤثرات التي تتجاوز عتبة $\epsilon_A$ في مسحة واحدة.
• البدء الدافئ (Warm-Start): يوفر ExcitationSolve المعامل الأمثل $\theta_j$ لكل مؤثر مختار على الفور، مما يسمح بتهيئة الـ ansatz بقيم حقيقية غير صفرية بدلاً من الصفر.

2. المعالجة المسبقة الكلاسيكية والاشتقاق التحليلي:

• بالنسبة لخطوة الاختيار الأولية التي تبدأ من حالة مرجعية لـ "هارتري-فوك" (HF)، يوضح المؤلفون أنه يمكن إجراء اختيار الإثارات الثنائية (double excitations) كلاسيكياً بالكامل.
• من خلال اختزال هاميلتوني البنية الإلكترونية إلى فضاء فرعي ثنائي الحالة (حالة HF + إثارة ثنائية واحدة)، اشتق المؤلفون صيغاً تحليلية (المعادلتان 4 و5) لحساب الحد الأقصى لخفض الطاقة ($\Delta E_{min}$) والمعاملات الأولية المثلى ($\theta_{ES}$) باستخدام تكاملات أحادية وثنائية الإلكترون فقط. وهذا يلغي الحاجة إلى الموارد الكمومية أثناء عملية اختيار المؤثر الأولية.

3. التكامل مع OVP-CEOs:

• تم تكييف الطريقة مع مؤثرات تبادل الأزواج ذات المعامل الواحد (OVP-CEOs)، والتي تقلل من عدد بوابات CNOT (من 13 إلى 9) وعمق الدائرة (من 11 إلى 7) مقارنة بالإثارات الكيوبت القياسية.
• لمعالجة التكرار (حيث تعمل كل من النسختين $+$ و $-$من OVP-CEOs بشكل متشابه على حالة HF) ومشكلات خلط المدارات، تم تنفيذ خطوة اختيار إضافية. تحدد هذه الخطوة المتغير الأمثل ($+$أو$-$) لكل زوج من المؤثرات بناءً على تأثيرهم على الحالة المتطورة، مما يضمن التقارب دون مضاعفة عمق الدائرة.

4. الارتباط النظري بالنظرية الاضطرابية:

• ربط المؤلفون نهجهم نظرياً بنظرية "إبستين-نيسبيت" (EN2) من الدرجة الثانية. وأظهروا أن سعات EN2 تتوافق مع توسع تايلور من الدرجة الأولى لسعات ES بدلالة مقياس الارتباط ($b/a$).
• يثبت هذا أن ES و EN2 متطابقان تقريباً في الأنظمة ضعيفة الارتباط، لكن ES يظل قوياً في الأنظمة ذات الارتباط القوي حيث قد تفشل EN2 (و MP2).

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

• اختيار المؤثرات بمسحة واحدة: بروتوكول يبني ansatz تقريبية للحالة الأرضية عبر مسحة واحدة على مجموعة المؤثرات، مما يقلل بشكل كبير من عدد تقييمات الدائرة الكمومية المطلوبة للاختيار.
• البدء الدافئ الكلاسيكي: الاشتقاق الرياضي للصيغ التي تسمح باختيار الإثارات الثنائية ذات الصلة وتهيئة معاملاتها باستخدام الحوسبة الكلاسيكية فقط، مما يوفر نقطة بداية عالية الجودة للمُحسّن الكمومي.
• التسريع التربيعي: الإثبات التجريبي بأن الطريقة المقترحة (ExcitationSolve + ES) تحقق تسارعاً في التقارب تربيعياً فيما يتعلق بحجم المجموعة (pool size) مقارنة بكل من ADAPT-VQE القياسي (التدرج النازل + اختيار التدرج) و ExcitationSolve التكيفي بدون ترتيب.
• التكامل مع OVP-CEO: التكيف الناجح مع OVP-CEOs، مما يوازن بين تقليل عمق الدائرة وبين التكاليف الحسابية الممكنة لاختيار المؤثرات.

النتائج

تم اختبار الطريقة على 12 جزيئاً (تتراوح من 4 إلى 20 كيوبت، مثل H2, LiH, CH4, C2) باستخدام محاكيات الحالة المتجهة (state-vector simulators) الخالية من الضجيج.

• الكفاءة: بالنسبة للجزيئات ذات مجموعات المؤثرات الكبيرة (مثل CH4 التي تحتوي على 804 مؤثراً)، قللت طريقة (ExcitationSolve + ES) الموارد الحسابية بعامل قدره $10^{-2}$ مقارنة بـ Pure ExcitationSolve.
• التقارب: تبني الطريقة ansätze مدمجة تحتوي فقط على المؤثرات ذات الصلة. بالنسبة لـ LiH، استخدم الـ ansatz الذي بناه ES عدد 34 مؤثراً مقارنة بـ 96 في UCCSD ثابت، مما أدى إلى تقارب أسرع وتجنب هضاب التحسين.
• جودة التهيئة: في الأنظمة ذات الارتباط القوي (المحاكية عبر تمديد أطوال الروابط)، قدم ES تخميناً أولياً متفوقاً على MP2 و EN2. وبينما كان MP2 هو الأسوأ باستمرار، حافظ ES على قوته في الحالات التي انحرفت فيها MP2 و EN2 عن الحل الدقيق.
• الإثارات من الرتب العليا: نجحت الطريقة في دمج الإثارات الثلاثية لـ LiH، محققة دقة أفضل برتبتين من الإثارات الثنائية وحدها، مع زيادة طفيفة فقط في تكلفة التقييم.
• أداء OVP-CEO: أدى استخدام OVP-CEOs مع معيار الاختيار المقترح إلى التقارب نحو نفس طاقة الحالة الأرضية مثل الإثارات القياسية، ولكن مع تقليل عمق الدائرة (عدد أقل من بوابات CNOT)، وإن كان ذلك مع زيادة طفيفة في عدد تقييمات الطاقة بسبب حجم المجموعة الأكبر.

الأهمية والادعاءات

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

• القابلية للتوسع: من خلال تقليل التعقيد الحسابي لاختيار المؤثرات وتمكين "البدء الدافئ" عبر المعالجة المسبقة الكلاسيكية، تجعل هذه الطريقة محاكاة الجزيئات الأكبر (تتجاوز حد الـ 12 كيوبت المعتاد في الأدبيات الحديثة) أكثر جدوى.
• الصلة بأجهزة NISQ: يوفر التكامل مع OVP-CEOs مقايضة مواتية لأجهزة NISQ، حيث يمثل عمق الدائرة العامل المحدد الأساسي.
• الرؤية النظرية: ترسي الورقة رابطاً جوهرياً بين اختيار المؤثرات الاستدلالي ونظريات الاضطراب الراسخة، موضحة أن ES ليس مجرد استدلال فعال بل هو تعميم قوي لنظرية EN2.
• قابلية التطبيق الواسعة: بينما صُممت الطريقة خصيصاً لـ VQE، يشير المؤلفون إلى أن استراتيجية المعالجة المسبقة الكلاسيكية والبدء الدافئ يمكن أن تستفيد منها طرق أخرى مكلفة حسابياً (مثل مونت كارلو الكمي، و QSE) وطرق الهاميلتوني المتغيرة التي تعتمد على حالات أولية عالية الجودة.

تخلص الورقة إلى أن نهجهم يجسّر الفجوة بفعالية بين الدوائر الضحلة للطرق التكيفية والخطوات المحدودة للـ ansätze الثابتة، مما يوفر مساراً لحسابات بنية إلكترونية أكثر كفاءة في عصر NISQ.