On Performance and Limitations of NISQ Hardware for Simulations of Quantum Wave Packet Dynamics

ملخص تقني: حول أداء وقيود أجهزة NISQ في محاكاة ديناميكيات الحزمة الموجية الكمومية

بيان المشكلة
تخضع المحاكاة العددية الدقيقة للعمليات الديناميكية الكمومية، مثل الاهتزازات الجزيئية، والتشتت، والنفق الكمومي، لمعادلة شرودنجر المعتمدة على الزمن. تواجه النهج الكلاسيكية "لعنة الأبعاد"، حيث تزداد التكلفة الحسابية بشكل أسي مع عدد درجات الحرية، مما يجعل الانتشار طويل الأمد للأنظمة التي تتجاوز بضعة أبعاد أمراً مستعصياً. وبينما توفر الحواسيب الكمومية منصة طبيعية لمحاكاة الأنظمة الكمومية نظراً لتوسع فضاء الحالة فيها بشكل أسي، إلا أن الانتشار الصريح للحزم الموجية في الزمن الحقيقي باستخدام الدوائر الكمومية الرقمية لا يزال غير مستكشف بعمق مقارنة بمشكلات البنية الإلكترونية. علاوة على ذلك، تعتمد طرق المحاكاة الكمومية الحالية غالباً على دالات طاقة وضع تحليلية أو تعاني من توسع أسي في تفكيك المؤثرات (عادةً $\mathcal{O}(4^n)$ لـ $n$ من الكيوبتات)، مما يحد من قابليتها للتطبيق على الأجهزة الحالية ذات المقياس المتوسط والضجيج (NISQ).

المنهجية
يقترح المؤلفون ترميزاً يعتمد على الشبكة للدالة الموجية على سجل مكون من $n$ من الكيوبتات، حيث تمثل $M = 2^n$ من نقاط الشبكة الإحداثي المكاني المنفصل. يتم تنفيذ التطور الزمني باستخدام نهج "المشغل المنقسم" (split-operator approach)، والذي يفصل مؤثري الطاقة الحركية ($\hat{T}$) وطاقة الوضع ($\hat{V}$).

1. الطاقة الحركية: يتم تطبيق مؤثر الطاقة الحركية باستخدام تحويل فوريه الكمومي (QFT). يقوم هذا التحويل بنقل الدالة الموجية إلى أساس الزخم حيث يكون $\hat{T}$ قطرياً، مما يسمح بالتطبيق بكفاءة عبر عمق دائرة يتوسع حدودياً بمقدار $\mathcal{O}(n^2)$.
2. طاقة الوضع: على عكس العمل السابق الذي اقتصر على طاقات الوضع التحليلية، يستخدم هذا البحث تمثيلاً عاماً لمؤثر طاقة الوضع كتركيبة خطية من سلاسل "باولي-Z" المتوافقة. ولأن طاقة الوضع تكون قطرية في أساس الموضع، فإن تفكيكها يتطلب $2^n$ من الحدود فقط (تتضمن عمليات $I$ و $Z$ فقط) بدلاً من تفكيك "باولي" الكامل $\mathcal{O}(4^n)$. هذا يقلل من توسع المؤثر إلى $\mathcal{O}(2^n)$ ويسمح بالترميز المباشر لمشاهد طاقة الوضع المنفصلة التعسفية.
3. التنفيذ: يتم تقريب الناشر العالمي عبر تفكيك "تروتر" (Trotter decomposition) من الدرجة الأولى. أُجريت عمليات المحاكاة على محاكيات كلاسيكية خالية من الضجيج (Qiskit SDK) ونُفذت على أجهزة NISQ، وتحديداً معالجات IBM فائقة التوصيل (Torino, Boston, Miami) وجهاز IonQ Forte ذي الأيونات المحاصرة.

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

• الكفاءة الخوارزمية: يوضح البحث طريقة تقلل من توسع تفكيك مؤثر طاقة الوضع من $\mathcal{O}(4^n)$ إلى $\mathcal{O}(2^n)$ من خلال استغلال الطبيعة القطرية للوضع في أساس الموضع.
• القابلية للتعميم: يتيح هذا النهج محاكاة طاقات وضع منفصلة تعسفية، متجاوزاً الأنظمة المقتصرة على الأشكال الدالية التحليلية المعروفة.
• قياس الأداء (Benchmarking): تقدم الدراسة قياساً شاملاً لديناميكيات الحزمة الموجية عبر أعداد كيوبتات تتراوح من 2 إلى 5 (ما يعادل 4 إلى 32 نقطة شبكة) على منصات أجهزة متعددة حالية.

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

• الأداء على النطاقات الصغيرة (2-3 كيوبت): أعادت جميع المنصات المختبرة (IBM و Ion-Q) إنتاج الديناميكيات المرجعية نوعياً. وفي حالات الـ 2-كيوبت، أظهرت أجهزة IBM الأحدث (Boston, Miami) توافقاً يقارب الكمي مع المحاكي.
• القابلية للتوسع والحساسية للضجيج: مع زيادة عدد الكيوبتات إلى 4 و 5، تباين الأداء بشكل كبير بين المنصات.
  • أجهزة IBM فائقة التوصيل: أظهرت النتائج تدهوراً سريعاً مع زيادة عدد الكيوبتات. في محاكاة الـ 5 كيوبت، انهارت الدالة الموجية على أجهزة IBM إلى توزيع منتظم بحلول الخطوة الزمنية الثالثة، مما يشير إلى فقدان بنية الحزمة الموجية بسبب تراكم أخطاء البوابات ومحدودية عمق الدائرة.
  • جهاز IonQ ذي الأيونات المحاصرة: حافظ معالج IonQ Forte على توافق أفضل بكثير مع المرجع الكلاسيكي عبر جميع أعداد الكيوبتات، بما في ذلك حالة الجسيم الحر المكونة من 5 كيوبتات. ويُعزى هذا الأداء إلى دقة أعلى لبوابات الكيوبت الثنائية والاتصال "الكل-إلى-الكل" المتأصل في بنيات الأيونات المحاصرة.
• خصوصية المشكلات: في سيناريوهات النفق والمتذبذب التوافقي، تفوق IonQ باستمرار على الأجهزة فائقة التوصيل، لا سيما في نطاق الـ 5 كيوبت حيث أظهرت أجهزة IBM انحرافات جوهرية.

الأهمية والادعاءات
تخلص الورقة إلى أنه بينما يمكن لأجهزة NISQ الحالية إعادة إنتاج الديناميكيات الكمومية نوعياً للمجموعات الصغيرة من الأنظمة، فإن دقة المحاكاة المعتمدة على الزمن حساسة للغاية لعمق الدائرة وضجيج الأجهزة. وتبرز النتائج أن بنيات الأيونات المحاصرة توفر حالياً أداءً متفوقاً لهذه المحاكاة الديناميكية المحددة مقارنة بالمعالجات فائقة التوصيل. ويؤكد المؤلفون أن التنفيذات القابلة للتوسع لديناميكيات الحزمة الموجية الكمومية ستتطلب تصميمات دوائر تقلل العمق وتكون قوية ضد ضجيج الأجهزة، مشددين على أن ترميز طاقة الوضع المقترح $\mathcal{O}(2^n)$ هو خطوة ضرورية نحو استخدام طاقات الوضع التعسفية على الأجهزة القريبة.