Quantum Metropolis-Hastings via Penalised Qubitized Walks: Spectral Filtering and Circuit Implementation

1. بيان المشكلة

تعد خوارزمية "متروبوليس-هستينغز" (Metropolis-Hastings) حجر الزاوية في طرق "مونت كارلو" بسلاسل ماركوف (MCMC)، وهي مستخدمة على نطاق واسع في الفيزياء، والاستدلال البايزي، والتعلم الآلي. ومع ذلك، تعاني طرق MCMC الكلاسيكية من بطء في الاختلاط (slow mixing) في الأسطح الطاقية عالية الأبعاد، أو متعددة الأنماط، أو الوعرة، حيث تتعثر السلسلة في النهايات الصغرى المحلية. ويتناسب زمن الاختلاط عكسياً مع الفجوة الطيفية لمصفوفة الانتقال.

تعد المشيات الكمومية (Quantum walks)، وتحديداً تلك القائمة على إطار "سيجيدي" للكمومة (Szegedy's quantization)، بتقديم تسريع تربيعي في زمن الاختلاط عن طريق تضخيم الفجوة الطيفية. ومع ذلك، يوجد عائق كبير: تتضمن خوارزمية متروبوليس-هستينغز القياسية حركات مرفوضة، مما يؤدي إلى عدم قابلية انعكاس جوهرية. إن تنفيذ منطق الرفض هذا مباشرة في دائرة كمومية وحدوية (unitary) أمر غير بسيط.

بينما اقترح كلاودون وآخرون (2025) إطاراً نظرياً للتغلب على ذلك من خلال العمل على فضاء حالة مزدوج (edge state space) واستخدام المشيات المكمتة (qubitized walks)، حدد مؤلفو هذه الورقة خللاً حرجاً في تنفيذ الدائرة المباشر لها:

• عدم التبادلية (Non-commutativity): يتطلب البناء النظري الإسقاط على فضاء القيم الذاتية 1 لمؤثر المشية ($W$) و تقييد الحالة إلى صورة شبه عملية جزئية معينة ($\text{Im}(\boxtimes)$). هذان الفضاءان الجزئيان لا يتبادلان، مما يعني أن عمليات الترشيح المتتالية تتداخل مع بعضها البعض، مما يمنع استخراج الحالة المستقرة الصحيحة.
• التعددية (Degeneracy): تكون القيمة الذاتية 1 عادةً متعددة في فضاء هيلبرت الكامل، مما يجعل من الصعب عزل الحالة المستهدفة الفريدة دون آليات إضافية.

2. المنهجية

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

أ. المشية المكمتة المعاقبة (The Penalised Qubitized Walk)

لحل مشكلات عدم التبادلية والتعددية، قدم المؤلفون مؤثر المشية المكمتة المعاقب ($V$):
$$V = \left( \Pi_{\boxtimes} + e^{i\phi}(\text{Id} - \Pi_{\boxtimes}) \right) W$$

• الآلية: يقوم هذا المؤثر بتطبيق دفعة طور (phase kick) قدرها $e^{i\phi}$ لأي حالة خارج صورة شبه العملية الجزئية $\boxtimes$، بينما يترك الحالات داخل $\text{Im}(\boxtimes)$ دون تغيير.
• الأثر: هذا يكسر التعددية العالمية للقيمة الذاتية 1. حيث تظل الحالة المستقرة المستهدفة (التي تقع في $\text{Im}(\boxtimes)$) عند طور 0، بينما يتم إزاحة جميع القيم الذاتية الأخرى ذات القيمة 1 إلى أطوار غير صفرية. وهذا يسمح لمرشح طيفي بعزل الحالة المستهدفة دون الحاجة إلى عمليات إسقاط متزامنة وغير متبادلة.

ب. الترشيح الطيفي عبر تقدير الطور الكمومي (QPE)

يستخدم سير العمل تقدير الطور الكمومي (QPE) لترشيح القيم الذاتية لـ $V$.

• استراتيجية الترشيح: بما أن الحالة المستهدفة تقابل طوراً ذاتياً قدره 0 (القيمة الذاتية 1)، فإن الخوارزمية تقوم باختيار لاحق (post-selection) بناءً على قياس سجل الطور للقيمة $|0\rangle$.
• قيود الدقة: تحلل الورقة المقايضة بين عدد الكيوبتات المساعدة للطور ($m$) والفجوة الطيفية ($\delta$). لتمييز الحالة المستهدفة عن النمط الأول المثار، يجب أن تكون دقة الطور $2^{-m}$ أصغر من الفجوة الزاوية $\Delta \approx \sqrt{\delta}$.
• التنفيذ شبه الكلاسيكي: لتقليل العبء على الكيوبتات في المحاكاة، استخدم المؤلفون تقنية تقدير الطور التكراري شبه الكلاسيكي (باستخدام كيوبت مساعد واحد معاد تدويره) لمحاكاة نموذج إيسينج، مما يضحي بعرض الدائرة مقابل عمقها.

ج. سير العمل الخوارذجي الكامل

يتكون سير العمل المقترح من الخطوات التالية:

1. التهيئة: إعداد حالة بذرة في $\text{Im}(\boxtimes)$ (على سبيل المثال، توزيع متساوٍ فوق الحواف الموجهة).
2. بناء المشية المعاقبة: بناء المؤثر $V$ باستخدام أوراكل (oracles) لنواة الاقتراح ($T$) ونواة القبول ($A$).
3. ترشيح QPE: تطبيق QPE على $V$ والاختيار اللاحق لسجل الطور ليكون $|0\rangle$. هذا يعزل الحالة الوحيدة ذات القيمة الذاتية 1.
4. فك الترميز (Decoding): تطبيق تحويل فك ترميز من "المزدوج إلى الرأس" (dual-to-vertex) لرسم الخريطة من تمثيل الحافة إلى تمثيل الرأس.
5. أخذ العينات: قياس سجل الموقع للحصول على عينات من التوزيع المستقر المستهدف $\pi$.

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

1. تحديد عقبات التنفيذ: تثبت الورقة بصرامة أن البناء النظري لـ كلاودون وآخرين لا يمكن تنفيذه مباشرة بسبب عدم تبادل مؤثرات الإسقاط المطلوبة.
2. حل المشية المعاقبة: يعد تقديم المؤثر المعاقب $V$ مساهمة مبتكرة تعمل على رفع تعددية القيم الذاتية وتضمن أن الفضاء الجزئي المستهدف ثابت تحت عملية الترشيح، مما يتيح تنفيذاً قوياً للدائرة.
3. بناء صريح للدائرة: يقدم المؤلفون تفاصيل الدائرة الكاملة، بما في ذلك بناء أوراكل الاقتراح والقبول لطوبولوجيا محددة (الشبكات الدورية ومكعبات هامينج)، ودمج آلية العقوبة.
4. إطار محاكاة شامل: قام المؤلفون بمحاكاة الخوارزمية الكاملة باستخدام موارد حوسبة عالية الأداء، مع تحليل تأثير دقة الطور المحدودة على جودة الحالة المُعدة.

4. النتائج

تحقق المؤلفون من صحة الخوارزمية على نظامين: جهد البئر المزدوج ثنائي الأبعاد (شبكة 4x4) وسلسلة إيسينج لـ 4 سبينات.

• جهد البئر المزدوج:

  • دور العقوبة: بمقارنة المشية المعاقبة ($m=4$) مع مشية غير معدلة ($m=5$)، فشلت المشية غير المعدلة في الوصول إلى التوزيع الصحيح رغم دقة الطور الأعلى، حيث أظهرت سوء تخصيص هائل للكتلة الاحتمالية بين البئرين. أما المشية المعاقبة فقد نجحت في إعادة إنتاج التوزيع ثنائي الأنماط المستهدف.
  • الاعتماد على الدقة: مع زيادة عدد الكيوبتات المساعدة للطور ($m$)، تحسنت دقة الحالة (fidelity) تجاه التوزيع المستهدف بشكل رتيب، وانخفضت مسافة التباين الكلي بشكل ملحوظ.

• نموذج إيسينج:

  • تحليل الفجوة الطيفية: تم اختبار الخوارزمية عبر درجات حرارة عكسية ($\beta$) متفاوتة. مع زيادة $\beta$، تضيق الفجوة الطيفية، مما يتطلب دقة طور أعلى ($m$) للحفاظ على الدقة.
  • نظام نقص الدقة: عندما تكون دقة الطور غير كافية، لا تكون الحالة الناتجة عشوائية بل تحتفظ بـ معلومات مهيكلة. فهي تلتقط بشكل صحيح معاملات النظام العالمية (مثل المغناطيسية) ولكنها تظهر أخطاء في الملاحظات الحساسة (مثل عد جدران النطاقات) وتقدير الطاقة.
  • ترشيح الطاقة: يعمل المرشح الطيفي بشكل هرمي، حيث يقمع قطاعات الطاقة العالية أولاً. وحتى قبل الوصول للتقارب الكامل، تكون الحالة منحازة نحو التكوينات منخفضة الطاقة، مما يشير إلى فائدة محتملة كبداية عالية الجودة لعمليات التحسين الكلاسيكية.

5. الأهمية والآفاق المستقبلية

• الجدوى: تثبت هذه العمل أن النسخة الكمومية من متروبوليس-هستينغز سليمة نظرياً وقابلة للتنفيذ عملياً في نظام مقاوم للأخطاء (fault-tolerant)، بشرما تمت معالجة تعددية القيم الذاتية عبر العقوبة.
• متطلبات الموارد: تتطلب الخوارزمية حواسيب كمومية مقاومة للأخطاء نظراً لعمق عمليات QPE والمشيات. هي ليست مناسبة حالياً للأجهزة ذات القدرة المحدودة (NISQ).
• الإمكانات الهجينة: يقترح المؤلفون استراتيجية هجينة كمومية-كلاسيكية: استخدام الخوارزمية الكمومية لإعداد حالة أولية عالية الجودة تلتقط بنية الحوض العالمية (للتغلب على بطء الاختلاط)، يتبعها MCMC كلاسيكي للتحسين المحلي.
• القيود: يعتمد التنفيذ الحالي على المحاكاة الكلاسيكية لمتجهات الحالة، مما يحد من حجم النظام. يظل الاختيار الأمثل لطور العقوبة $\phi$ تجريبياً ويتطلب مزيداً من التحليل النظري.

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