Quantum speed-up for solving the one-dimensional Hubbard model using quantum annealing

ملخص تقني: التسريع الكمي لحل نموذج هوبارد أحادي البعد باستخدام التلدين الكمي

بيان المشكلة
يُعد نموذج هوبارد أحادي البعد إطاراً جوهرياً في فيزياء المادة المكثفة لفهم أنظمة الإلكترونات المترابطة. وبينما يُعتبر هذا النموذج "قابلاً للحل تحليلياً" (integrable) — أي توجد نتائج تحليلية دقيقة لحالته الأرضية في الحد الديناميكي الحراري عبر معادلات "بيث-أنساتز" (Bethe-ansatz) — فإن حل هذه المعادلات للأنظمة ذات الحواف المفتوحة يطرح تحديات حوسبية كبيرة. وتحديداً، بينما يمكن إيجاد طاقة الحالة الأرضية في وقت حدودي عبر حل معادلات "بيث-أنساتز" غير الخطية المترابطة، فإن تحديد دالة الموجة الكاملة للحالة الأرضية (المعاملات) يتطلب حساب عدد هائل من الحدود، مما يستليزم موارد حوسبية أسية. يستقصي المؤلفون ما إذا كانت الخوارزميات الكمية، وتحديداً التلدين الكمي (quantum annealing)، يمكن أن توفر تسريعاً على الطرق الكلاسية لإيجاد الحالة الأرضية لأنظمة هوبارد أحادية البعد عند نصف الامتلاء.

المنهجية
يقوم المؤلفون بتنفيذ محاكاة قائمة على البوابات لعملية التلدين الكمي على محاكي حاسوب كمي عام (JUQCS) لحل هاميلتونيان هوبارد. تتضمن المنهجية ثلاث مراحل رئيسية:

1. رسم خريطة الهاميلتونيان: يتم تحويل هاميلتونيان هوبارد الفيرميوني إلى هاميلتونيان كيوبت باستخدام تحويل "جوردان-ويجنر". يقوم هذا التحويل بتحويل مؤثرات الخلق والفناء الفيرميونية إلى مصفوفات "باولي" تعمل على $2L$ من الكيوبتات (حيث $L$ هو عدد مواقع الشبكة)، مع الحفاظ على علاقات التبادل المضاد.
2. تحضير الحالة الأولية: تبدأ عملية التلدين بالحالة الأرضية لهاميلتونيان القفز غير المتفاعل ($H_I$). يستخدم المؤلفون خوارزمية تعتمد على دورات "جيفنز" (Givens rotations) لتحضير هذه الحالة الأولية، والتي تقابل ملء أدنى حالات الزخم الطاقية لعدد معين من الإلكترونات ذات اللف المغزلي العلوي والسفلي.
3. ديناميكيات التلدين: يتطور النظام تحت تأثير هاميلتونيان يعتمد على الزمن $H(s) = (1-s)H_I + sH_H$، حيث تتغير $s$ من 0 إلى 1. يتم محاكاة التطور الزمني باستخدام تفكيك "سوزوكي-تروتر" (تحديداً تقريب من الدرجة الثانية) لتقسيم التطور المستمر إلى خطوات زمنية منفصلة. يتم تفكيك المؤثرات الوحدوية لكل خطوة إلى بوابات عالمية أحادية وثنائية الكيوبت (بوابات Hadamard، و $R_Z$، و $R_{ZZ}$، و $X$).
4. تحليل القياس (Scaling Analysis): أُجريت المحاكاة لأحجام أنظمة تصل إلى $L=20$ (40 كيوبت). مقياس الأداء هو الطاقة المتبقية $\Delta E = \langle \psi(1) | H_H | \psi(1) \rangle - E_0$، حيث $E_0$ هي طاقة الحالة الأرضية الدقيقة المحسوبة عبر "بيث-أنساتز". يحلل المؤلفون كيف يتدرج زمن التلدين الإجمالي المطلوب ($T_A$) مع حجم النظام ($L$) لتحقيق دقة معينة.

المساهمات والنتائج الرئيسية

• التنفيذ القائم على البوابات: تقدم الورقة بناءً ملموساً لمحاكاة التلدين الكمي لنموذج هوبارد على حاسوب كمي قائم على البوابات، مع تفصيل تفكيك الدائرة والمتطلبات الموردية (أعداد البوابات).
• تدرج زمن التلدين: بالنسبة للأنظمة عند نصف الامتلاء ($2N_\downarrow = N = L$)، يجد المؤلفون أن الطاقة المتبقية تتدرج كـ $\Delta E \propto L^{1.234} / T_A^2$ لجداول التلدين الخطية. ولتحقيق دقة ثابتة $\Delta E^*$، يتدرج زمن التلدين المطلوب كـ $T_A \propto L^{0.62}$ لأحجام الأنظمة التي تقل عن عتبة معينة $L^*$.
• نظام الاقتران القوي: تؤكد عمليات المحاكاة لقوى التفاعل القوية ($U/t_H = 8, 16$) أن بداية نظام التدرج تظل حميدة، حيث يظل زمن التلدين محدوداً بدالة خطية لحجم النظام ($T_A \propto L^{0.995}$) لأي حجم $L$ كبير جداً.
• المقارنة مع الطرق الكلاسيكية: يقارن المؤلفون نتائجهم بخوارزميات البحث عن الجذور الكلاسيكية لمعادلات "بيث-أنساتز". وبينما إيجاد الطاقة عبر "بيث-أنساتز" يتم في وقت حدودي ($T_{BA} \propto L^{3.3}$)، فإن الحصول على معاملات الحالة الأرضية الكاملة مكلف أسياً. في المقابل، يوفر نهج التلدين الكمي الحالة الأرضية بتدرج زمني وموردِيّ تحت خطي إلى خطي بالنسبة لحجم النظام.

الأهمية والادعاءات
يدعي البحث أن التلدين الكمي يوفر تسريعاً كمياً جوهرياً على الخوارزميات الكلاسيكية لإيجاد الحالة الأرضية لنموذج هوبارد أحادي البعد.

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

يظل المؤلفون متواضعين فيما يتعلق بالتطبيق العملي المباشر، مدركين أن الحسابات الحالية التي تتضمن آلاف البوابات تتجاوز قدرات الأجهزة المتاحة حالياً، لكن تدرج الموارد يبشر بالخير للهياكل المستقبلية القادرة على تصحيح الأخطاء.