Quantum-echo Markov process for combinatorial optimization
تقدم هذه الورقة عملية ماركوف ذات صدى كمي للتحسين التوافقي تستفيد من الديناميكيات الكمية لهندسة نوى انتقال مهيكلة، مما يثبت أن الجمع بين الاستكشاف المدفوع كمياً والاستغلال الجشع يوازن بفعالية بين عدم التمركز في فضاء هامينج والتموضع في فضاء الطاقة لتعزيز أداء التحسين.
يعد حل الألغاز المعقدة جزءاً أساسياً من كيفية تنقلنا في العالم، بدءاً من تنظيم مسار تسليم وصولاً إلى جدولة غرف العمليات في المستشفيات. هذه هي المسائل التوافقية (combinatorial problems)، حيث يتمثل الهدف في إيجاد أفضل ترتيب واحد من بين عدد هائل من الاحتمالات. لعقود من الزمن، بحث العلماء في ميكانيكا الكم طلباً للمساعدة، آملين أن يتمكن السلوك الغريب للجسيمات من استكشاف مساحات البحث الضخمة هذه بشكل أسرع من أي حاسوب كلاسيكي. يستخدم نهجان بارزان، يُعرفان بالتقليب الكمي (quantum annealing) وخوارزمية التحسين التقريبي الكمي (quantum approximate optimization algorithm)، حركات كمية محكومة لتوجيه النظام نحو الحل. ومع ذلك، أظهرت الأبحاث الحديثة أنه عندما تُستخدم هذه الأدوات الكمية مع موارد محدودة — أي أنها تعمل لفترة زمنية قصيرة أو بعدد ثابت من الخطوات — فإنها غالباً ما تتعثر؛ إذ تميل إلى النظر فقط في الخيارات القريبة، مما يؤدي إلى تفويت الحلول الأفضل التي تقع بعيداً، أو تقفز بجنون بحيث تغير تكلفة الحل بشكل جذري يجعلها غير مفيدة.
اقترح باحث في جامعة واسيدا طريقة جديدة لتسخير هذه الموارد الكمية المحدودة، ليس لإيجاد الإجابة النهائية مباشرة، بل لتعمل كدليل متطور لعملية البحث. لقد طور طريقة تسمى "عملية ماركوف للصدى الكمي" (quantum-echo Markov process). تخيل مسافراً يحاول العثور على أدنى نقطة في سلسلة جبلية ضبابية شاسعة. قد يتفحص المتسكع البسيط الأرض المحيطة بقدميه مباشرة فقط، مما يعرضه لخطر الوقوع في وادٍ صغير. أما القافز المتهور فقد يقفز عبر السلسلة بأكملها، لكنه معرض بقدر ما هو معرض للوقوع على قمة عالية كما هو معرض للوقوع في وادٍ منخفض. أراد الباحث طريقة يمكنها أخذ المسافر بعيداً عن مكانه الحالي دون إرساله ليطير إلى ارتفاع أعلى وأسوأ بكثير. ولتحقيق ذلك، استخدم تسلسلاً كمياً محدداً: التحرك للأمام في الزمن، تطبيق دفعة محلية صغيرة، ثم التحرك للخلف في الزمن. تسمح تقنية "الصدى" هذه للنظام باستكشاف تكوينات بعيدة في مساحة البحث مع الحفاظ على التغييرات في التكلفة الإجمالية صغيرة ويمكن التحكم فيها.
اختبر الباحث هذا النهج على نوعين مختلفين من المشاهد الرياضية. الأول كان نموذج "إيسينج" عشوائي (random Ising model)، والذي يحاكي نظاماً معقداً تتفاعل فيه الأجزاء مع بعضها البعض بطرق محددة، مما يخلق تضارساً وعراً من التلال والوديان. والثاني كان نموذج الطاقة العشوائي (random energy model)، وهو مشهد أكثر فوضوية حيث لا توجد علاقة بين ارتفاع التضاريس والموقع، مما يمثل اختباراً صارماً لقدرة الطريقة على إيجاد بنية حيث لا توجد بنية طبيعية. ومن خلال تشغيل عمليات محاكاة لأنظمة تصل إلى أربعة عشر متغيراً، لاحظوا أنه مع زيادة مدة الحركة الكمية أو عدد الخطوات في خوارزميتهم، تصبح العملية فعالة بشكل ملحوظ. بدأت في الوصول إلى تكوينات مختلفة تماماً عن نقطة البداية، ومع ذلك ظلت تكلفة هذه التكوينات الجديدة قريبة من التكلفة الأصلية. وهذا مزيج نادر: القدرة على السفر بعيداً دون دفع ثمن باهظ.
اكتشف الباحث أن هذا النجاح يأتي من آليتين متمايزتين تعملان معاً. فالقدرة على الوصول إلى أماكن بعيدة تنبع من الطريقة التي تنتشر بها المعلومات الكمية، مما يربط بفعالية بين أجزاء متباعدة في مساحة البحث. أما القدرة على البقاء قريباً في التكلفة فتنشأ من ارتباط دقيق يولده المسار الكمي بين موقع النظام وطاقته. في نموذج "إيسজ" العشوائي، يعد هذا الارتباط نتيجة طبيعية لتطور النظام ببطء بما يكفي لاحترام بنيته الأساسية. وفي نموذج الطاقة العشوائي الأكثر فوضوية، يتم إنشاء هذا الارتباط من خلال ضبط معايير الدائرة الكمية بدقة. ووجد الباحث أن هذا التوازن دقيق؛ فإذا أصبحت العملية شديدة التركيز على إبقاء التكلفة منخفضة، فإنها تفقد قدرتها على الاستكشاف، وتتوقف عملية البحث.
ولوضع هذا الدليل الكمي في العمل، طبق الباحث استراتيجية تحسين تكرارية. لقد ترك العملية الكمية تقترح تكويناً جديداً، ولكن لم يقبل الحركة إلا إذا كانت تحسن جودة الحل أو تحافظ عليها. وعندما اختبر ذلك على سلسلة مغناطيسية بسيطة ونموذج "إيسلج" عشوائي معقد، وجدوا أن طريقة الصدى الكمي تفوقت على عمليات البحث العشوائية القياسية، خاصة عند البحث عن حلول عالية الجودة. ومع ذلك، لاحظوا أيضاً حداً: إذا أصبحت العملية الكمية مقيدة للغاية، فإنها ستفشل في الهروب من الفخاخ المحلية. ولحل ذلك، دمجوا خطوات الصدى الكمي مع تقنية كلاسيكية تُعرف باسم "الانحدار الجشع" (greedy descent). فبعد أن تقترح العملية الكمية مكاناً جديداً، يقوم حاسوب كلاسيكي فوراً باتخاذ سلسلة من الخطوات الصغيرة نحو الأسفل لإيجاد أفضل حد أدنى محلي من تلك النقطة الجديدة.
أثبت هذا النهج الهجين أنه الأكثر قوة. فقد وفرت الديناميكيات الكمية الاستكشاف اللازم للقفز خارج الوديان المحلية، بينما ضمن الانحدار الجشع استغلال كل فرصة للتحسين بمجرد الهبوط في منطقة جديدة. وفي عمليات المحاكاة، أدت إضافة خطوة الانحدار الجشع إلى تحسين معدل النجاح وسرعة إيجاد أفضل الحلول بشكل كبير، حتى في الحالات التي كانت فيها العملية الكمية وحدها تعاني. تشير النتائج إلى أن الموارد الكمية المحدودة، عندما يتم هندستها بشكل صحيح، يمكن أن تعمل كأداة أولية قوية للتحسين التكراري. فبدلاً من محاولة حل المشكلة بأكملها في قفزة كمية واحدة، تستخدم هذه الطريقة الديناميكيات الكمية لتوليد حركات مهيكلة وذكية يمكن للحاسوب الكلاسيكي بعد ذلك صقلها. وتوضح الدراسة أن هذا التوازن بين الاستكشاف البعيد والبقاء قريباً هو المفتاح لإطلاق إمكانات الحواسيب الكمية لحل مشكلات التحسين في العالم الحقيقي، مما يوفر مساراً واعداً لاستخدام الأجهزة الكمية المحدودة اليوم لمعالجة أصعب ألغاز الغد.
بيان المشكلة تتناول الورقة البحثية تحدي استخدام الديناميكيات الكمية ذات الموارد المحدودة (تحديداً التلدين الكمي (QA) وخوارزمية التحسين التقريبي الكمي (QAOA)) في التحسين التوليفي. وبينما يعد التلدين الكمي وQAOA واعدين، تشير الدراسات الحديثة إلى أنهما في ظل ظروف الزمن المحدود أو العمق الثابت، يظهران هياكل محلية فعالة. وتحديداً، غالباً ما تقتصر ديناميكياتهما على جيران محددين في الرسم البياني الأساسي، مما يمنعهما من الهروب من النهايات الصغرى المحلية بفعالية، أو على العكس من ذلك، يعملان كعملية أخذ عينات عشوائية إذا استكشفا مساحات واسعة جداً. وتتمثل المشكلة الجوهرية في هندسة نواة انتقال (transition kernels) لعمليات البحث التكرارية تكون غير محلية في فضاء هامينج (Hamming space) (لاستكشاف التكوينات البعيدة) ولكنها محلية في فضاء الطاقة (energy space) (لتجنب التغيرات الكبيرة والضارة في الطاقة)، مما يسمح بالهروب الفعال من النهايات الصغرى المحلية دون عشوائية البحث.
المنهجية يقترح المؤلف إنشاء عملية ماركوف للصدى الكمي لبناء نوى الانتقال هذه. تُعرف العملية بواسطة مؤثر وحدوي U^ (مشتق من QA أو QAOA) ومجموعة من المؤثرات الوحدوية الهيرميتية المحلية {V^i} (المختارة كقلبات سبين فردية σ^ix).
ديناميكيات الصدى الكمي: يتم توليد الانتقال عبر المتسلسلة U^V^iU^†. بدءاً من حالة القاعدة الحسابية ∣z⟩، يتطور النظام للأمام عبر U^†، ثم يخضع لاضطراب محلي V^i، ويتطور للخلف عبر U^.
مصفوفة الانتقال: يتم تحديد احتمال الانتقال من التكوين z إلى z′ من خلال قياس الحالة ∣ψi(z)⟩=U^V^iU^†∣z⟩ في القاعدة الحسابية.
المحلية في فضاء هامينج (dH): متوسط مسافة هامينج بين التكوينات الحالية والمرشحة.
المحلية في فضاء الطاقة (dE): متوسط مربع التغير في الطاقة (ΔC2) الناتج عن الانتقال.
النماذج: تم اختبار الإطار العملي على نموذج آيزينج العشوائي (Random Ising Model) (حيث ترتبط بنية الطاقة وهامينج ببعضهما) ونموذج الطاقة العشوائي (REM) (حيث لا توجد علاقة بينهما).
استراتيجيات التحسين: تم اقتراح خوارزميتين للتحسين التكراري:
التحسين المحلي للصدى الكمي: يقبل التكوين المرشح فقط إذا كان ΔC≤0.
التحسين المحلي للصدى الكمي مع النزول الجشع (Greedy Descent): بعد انتقال الصدى الكمي، يتم إجراء عملية نزول جشع بقلب سبين فردي للوصول إلى نهاية صغرى محلية قبل التكرار التالي.
المساهمات والآليات الرئيسية
هندسة الانتقالات المهيكلة: تثبت الورقة أن زيادة زمن التلدين τ في QA أو عدد الطبقات p في QAOA تزيد في آن واحد من عدم المحلية في فضاء هامينج (عدم التمركز). ومع ذلك، فإن سلوك المحلية في فضاء الطاقة يختلف باختلاف النموذج: في نموذج الطاقة العشوائي (REM)، تقل dE (مما يحسن المحلية) مع زيادة τ أو p، بينما في نموذج آيزينج العشوائي، تكون dE مستقلة تقريباً من p وτ. ومن المثير للاهتمام أنه في نموذج REM، يمكن أن تكون هذه الانتقالات أكثر محلية في فضاء الطاقة من مجرد قلبات السبين الفردية البسيطة، رغم كونها أكثر عدم محلية في فضاء هامينج.
آلية عدم التمركز في هامينج: ينشأ هذا من انتشار المؤثر (operator spreading). مع زيادة τ أو p، تتطور سلاسل باولي من الرتب المنخفضة في صورة هيزنبرغ إلى سلاسل عالية الرتبة تشمل سبينات عديدة، مما يسمح للانتقال بالربط بين تكوينات متباعدة في فضاء هامينج.
آلية التمركز في الطاقة: ينشأ هذا من ارتباط ديناميكي مُولد بين دالة التكلفة (Cz) والوزن الهامينجي (wz) للقيم الذاتية للهاملتوني الفعال U^†C^U^.
في QA، ينشأ هذا الارتباط من الديناميكيات شبه الأديباتية (quasiadiabatic) الخشنة قبل حل الفروقات الأسية الضئيلة في مستويات الطاقة.
في QAOA، يتم توليده عبر التحسين المتغير لمعلمات الدائرة، مما يعزز الحدود منخفضة الرتبة في الهاملتوني الفعال، ويخلق ارتباطاً خطياً بين Cz وwz.
دور النزول الجشع: وجدت الدراسة أنه بينما توفر الديناميكيات الكمية الاستكشاف الضروري (عدم المحلية)، فإن التمركز المفرط في فضاء الطاقة قد يؤدي إلى تدهور الأداء عبر كبح الانتقالات الضرورية. إن دمج خطوة النزول الجشع (الاستغلال/Exploitation) يحسن الأداء بشكل كبير، مما يسلط الض الضوء على الأدوار المتكاملة للديناميكيات الكمية للاستكشاف (Exploration) والخطوات الجشعة الكلاسيكية للاستغلال.
النتائج
ضبط المحلية: تظهر المحاكاة العددية على نموذج آيزينج العشوائي وREM أن dH تزديد مع زيادة τ أو p. بالنسبة لـ REM، تقل dE مع زيادة τ أو p، ويستطيع تقريب الطور العشوائي (random-phase approximation) إعادة إنتاج النتيجة الدقيقة بدقة. أما بالنسبة لـ نموذج آيزينج العشوائي، تظل dE مستقلة تقريباً من τ وp، ويظهر انحراف بين النتائج الدقيقة وتقريبات الطور العشوائي عند قيم τ وp الصغيرة.
سلسلة فيرمو-ماغنيتيكية (Ferromagnetic Chain): في سلسلة آيزينج الفيرومغناطيسية، يتناسب عدد التكرارات المطلوبة للوصول إلى فجوة المثالية المستهدفة خطياً مع حجم النظام N. وتتفوق عملية الصدى الكمي على قلبات السبين الفردية العشوائية للفجوات المثالية الصغيرة. كما تبين أن الانجراف الديناميكي هو تقريباً حاصل ضرب الانجراف الحراري في "نسبة العيب المحلي"، والتي تقترب من الواحد الصحيح مع زيادة عمق الدائرة.
نموذج آيزينج العشوائي:
بدون النزول الجشع: تزداد احتمالية النجاح (psuc) في البداية مع زيادة τ أو p ولكنها قد تنخفض عند τ الكبيرة (في QA) بسبب تحول مصفوفة الانتقال إلى قطرية كتلية (block-diagonal) في الحد الأديباتي، مما يؤدي لحبس النظام. يشير هذا إلى أن التمركز المفرط في فضاء الطاقة قد يكون ضاراً.
مع النزول الجشع: تقترب احتمالية النجاح بسرعة من الواحد حتى لقيم τ أو p أصغر. والأهم من ذلك، في QA، تظل احتمالية النجاح هي الواحد حتى عند τ الكبيرة حيث يفشل النسخ غير الجشع. كما يتناقص وقت الوصول (hit time) مع زيادة p أو τ.
القابلية للتوسع (Scalability): يظل المعامل ξ (مؤشر المحلية في فضاء الطاقة) مستقلاً تقريباً عن حجم النظام في نموذج آيزينج العشوائي، مما يشير إلى أن الآلية تستمر في الأنظمة الأكبر.
الأهمية والادعاءات تؤسس هذه الورقة عملية ماركوف للصدو الكمي كإطار عمل لهندسة نوى انتقال ذات خصائص محلية قابلة للضبط. وتدعي الورقة أن الديناميكيات الكمية متعددة الأجسام ذات الموارد المحدودة يمكن أن تعمل كأداة حسابية أولية للتحسين التكراري، بشرط استخدامها لتوليد مقترحات ضمن حلقة هجينة كمية-كلاسيكية بدلاً من استخدامها كمحسنات مستقلة.
يشير المؤلف بتواضع إلى أنه بينما تشير نتائجه في نموذج آيزينج العشوائي إلى القابلية للتوسع، فإن القابلية للتوسع التقاربية في نموذج REM وإمكانية تدريب معلمات QAOA عند أحجام الأنظمة الكبيرة (مشكلة الهضاب المسطحة/Barren Plateaus المحتملة) تظل أسئلة مفتوحة. كما يذكر أن زمن التلدين الأمثل أو العمق للمشاكل الأكثر صعوبة، وتأثير قواعد القبول الأكثر تطوراً (مثل السماح بالحركات الصاعدة في الطاقة)، تتطلب مزيداً من الاستقصاء. يوفر هذا العمل مساراً لاستخدام الديناميكيات الكمية محدودة الزمن ليس فقط لإيجاد الحل النهائي، بل لتحسين التحسين التقريبي بشكل منهجي من خلال الاستكشاف المهيكل.