Quantum computer-based simulation of Stark many-body localization in a 1D Fermi-Hubbard model
تُظهر هذه الورقة محاكاة قائمة على الحاسوب الكمي لظاهرة التمركز الجماعي لستارك في نموذج "فيرمي-هوبارد" أحادي البعد مكون من 12 كيوبت على أجهزة IBM، وذلك باستخدام تقنيات تجميع وتحسين متقدمة لتقليل تعقيد الدارات بشكل كبير ورصد الانتقال من الديناميكيات الحرارية إلى الديناميكيات المتمركزة بنجاح.
المؤلفون الأصليون:Abdul Kalam, Prasenjit Deb, Akitada Sakurai, Tapan Mishra, V. S. Prasannaa, B. P. Das
تخيل الكون كأنه ساحة رقص عملاقة وفوضوية. عادةً، عندما تضع بضعة راقصين على هذه الساحة، فإنهم يصطدمون ببعضهم البعض، ويتبادلون الشركاء، وفي النهاية يمتزجون تماماً لدرجة أنك لا تستطيع معرفة من أين بدأ كل منهم. في الفيزياء، يُسمى هذا الامتزاج الفوضوي "التحرّر الحراري" (thermalization)، وهو الكيفية التي تصل بها معظم الأشياء في الطبيعة إلى حالة التوازن، مثل فنجان قهوة ساخن يبرد ليصل إلى درجة حرارة الغرفة. ولكن ماذا لو كانت ساحة الرقص تحتوي على جدران غير مرئية أو منحدر مائل غريب يجبر الراقصين على البقاء في زواياهم الصغيرة الخاصة، رافضين الامتزاج؟ هذا هو لغز "توطين الأجسام المتعددة" (Many-Body Localization - MBL). إنها ظاهرة ينسى فيها النظام الكمي كيفية التحرر حرارياً، محتفظاً بذاكرة طويلة المدى لموقعه الأصلي. وبينما عرف العلماء منذ زمن طويل أن الفوضى العشوائية (الاضطراب) يمكن أن تسبب ذلك، تشير فكرة أحدث إلى أن حتى منحدرًا مستوياً ومنتظماً تماماً (مجالاً كهربائياً موحداً) يمكنه حبس الجسيمات بنفس الفعالية. يُسمى هذا "Stark MBL". إن فهم هذا الأمر يعد أمراً بالغ الأهمية لأنه يتحدى قواعدنا الأساسية حول كيفية انتقال الطاقة والمعلومات، وقد يكون المفتاح لبناء حواسيب كمومية أفضل لا تفقد بياناتها.
في هذه الورقة البحثية، قرر فريق من الباحثين اختبار فكرة "Stark MBL" هذه باستخدام حاسوب كمومي حقيقي ومادي، بدلاً من مجرد محاكاة على حاسوب محمول عادي. لقد أعدوا نموذجاً يُسمى "نموذج فيرمي-هوبارد" (Fermi-Hubbard model)، وهو يشبه ملعباً رقمياً للإلكترونات التي تقفز بين مواقع الشبكة، لكنهم أضافوا مجالاً كهربائياً موحداً قوياً لإمالة الملعب بأكمله. كان الهدف هو معرفة ما إذا كانت الجسيمات ستظل عالقة في مكانها (موطنة) أم أنها ستتحرك بحرية وتمتزج في النهاية (تتحرر حرارياً). كان التحدي هو أن الحواسيب الكمومية الحالية "ضوضائية" قليلاً—مثل محاولة حل لغز معقد بينما يقوم شخص ما بهز الطاولة. وللحصول على إجابة واضحة، كان على الفريق أن يكون ذكياً للغاية في برمجته. لقد استخدموا تقنية رسم خرائط خاصة لترجمة مشكلة الإلكترونات إلى "كيوبتات" (qubits) -وهي بتات الحاسوب الكمومي- ثم طبقوا روتيناً لـ "ضغط الدائرة" (circuit compression). فكر في الأمر كأنك تأخذ كرة ضخمة ومتشابكة من الخيوط وتصغرها لتصبح عقدة صغيرة ومرتبة دون فقدان النمط. كان هذا التحسين نجاحاً باهراً، حيث قلل عدد العمليات المعقدة ثنائية الكيوبت بنسبة 88% تقريباً وقلل عمق الدائرة بنسبة 87%، مما جعل التجربة ممكنة على هذه الأجهزة غير المثالية اليوم.
كانت النتائج انتصاراً واضحاً لنظرية "المنحدر المائل". فعندما طبق الباحثون مجالاً كهربائياً ضعيفاً، تصرفت الجسيمات بشكل طبيعي: كانت تقفز وتختلط، وفقد النظام ذاكرته عن مكان بدايته، تماماً مثل ساحة رقص عادية. ومع ذلك، عندما رفعوا قوة المجال الكهربائي إلى مستوى عالٍ (تحديداً قوة ميل قدرها 15)، تغير السلوك بشكل دراماتيكي. توقفت الجسيمات عن الحركة؛ وبقيت تماماً حيث وُضعت، رافضة القفز إلى المواقع المجاورة. احتفظ النظام بـ "ذاكرة" قوية لحالته الأولية، ونما التشابك بين الجسيمات ببطء شديد، وهو العلامة المميزة للتوطين. تحقق الفريق من هذه النتائج عبر مقارنة نتائج حاسوبهم الكمومي الضوضائي مع عمليات محاكاة مثالية وخالية من الضوضاء وحسابات رياضية دقيقة، وكانت النتائج متطابقة بشكل ملحوظ. لقد رصدوا تأثير "التجميد" هذا باستخدام عدة قياسات مختلفة، مثل تتبع عدد الجسيمات التي بقيت في جانب واحد من السلسلة وقياس "عدم توازن الشحنة" (مدى عدم انتظام توزيع الجسيمات).
تستبعد الورقة البحثية صراحةً فكرة أن هذا التوطين ناتج عن اضطراب عشوائي أو فوضى في النظام؛ بل تؤكد أن مجالاً كهربائياً حتمياً وموحداً كافٍ لحبس الجسيمات. الباحثون واثقون من هذه النتائج لأن نتائج أجهزتهم كانت تحاكي بدقة عمليات المحاكاة الدقيقة، رغم أنه كان عليهم محاربة الضوضاء في المعالج الكمومي. لم يكتفوا بمجرد الاقتراح بأن هذا يحدث، بل قاموا بقياسه في الديناميكيات الزمنية الحقيقية على نظام مكون من 12 كيوبت. تُظهر الدراسة أنه باستخدام حيل التحسين الصحيحة، يمكن للحواسيب الكمومية الحالية محاكة حالات كمومية معقدة غير متحررة حرارياً بنجاح. وهذا يثبت أنه يمكننا استخدام هذه الآلات الضوضائية لاستكشاف أسئلة فيزيائية عميقة، مثل كيف تظل الجسيمات عالقة في عالم منظم تماماً، مما يمهد الطريق لتجارب مستقبلية قد تساعد يوماً ما في بناء تقنيات كمومية أكثر استقراراً.
ملخص تقني: محاكاة حاسوبية كمومية لظاهرة التموضع في الأجسام المتعددة تحت تأثير "ستارك" (Stark Many-Body Localization)
بيان المشكلة يمثل التموضع في الأجسام المتعددة (MBL) انهياراً في عملية التوازن الحراري في الأنظمة الكمومية المعزولة، ويتميز بالاحتفاظ بذاكرة الحالة الأولية ونمو بطيء للتشابك. وبينما يرتبط هذا النوع تقليدياً بالاضطراب الناجم عن التغيرات العشوائية (Quenched Disorder)، تشير أعمال نظرية حديثة إلى أن مجالاً كهربائياً خطياً منتظماً (جهد ستارك) يمكن أن يحفز ظاهرة الـ MBL في الأنظمة الخالية من الاضطراب، وهو ما يُعرف بـ "Stark MBL". ويتمثل التحدي الرئيسي في دراسة هذه الظواهر في محدودية قدرات الحواسيب الكلاسيكية في محاكاة ديناميكيات الأنظمة متعددة الأجسام المتفاعلة، لا سيما مع زيادة حجم النظام ووقت التطور. علاوة على ذلك، فإن محاكاة الـ MBL على أجهزة الكموم الحالية ذات النطاق المتوسط المشوب بالضجيج (NISQ) تعوقها قيود عمق الدارة وعدد البوابات، مما يؤدي إلى فقدان الترابط السريع وتراكم الأخطاء. يتناول هذا العمل جدوى محاكاة Stark MBL في نظام فرميوني مترابط مكون من 12 كيوبت باستخدام أجهزة كمومية من شركة IBM، مع التركيز بشكل خاص على التغلب على قيود الموارد في المعالجات المشوبة بالضجيج.
المنهجية تستخدم الدراسة نموذج "فيرمي-هوبارد" المائل أحادي البعد (t-FHM)، حيث يتم تطبيق جهد "ستارك" خطي على شبكة "فيرمي-هوبارد" القياسية. تتضمن الهاملتوني (Hamiltonian) القفز بين الجيران الأقرب (J)، والتفاعل الكولومي في الموقع (U)، ومصطلح ميل "ستارك" (Δ0).
لمحاكاة هذا النظام على حاسوب كمومي، استخدم المؤلفون المسار التقني التالي:
الرسم (Mapping): تم استخدام تحويل "جوردان-ويغنر" (Jordan-Wigner) الموزع حسب اللف المغزلي لربط الهاملتوني الفرميوني بالكيوبتات. وبخلاف عمليات الربط التقليدية المتداخلة، يقوم هذا النهج بتعيين الفرميونات ذات اللف المغزلي للأعلى والأسفل إلى سلاسل تماثل منفصلة (الكيوبتات الزوجية والفردية)، مما يقلل من طول سلاسل "باولي" المرتبطة بعمليات القفز.
بناء الدارة: تم تقريب التطور الزمني باستخدام تفكيك "تروتر-سوزوكي" من الدرجة الأولى. تم تنظيم الدارة في ثلاث كتل: جهد "ستارك"، التفاعل في الموقع، والقفز. وللتعامل مع مصطلحات القفز غير الموضعية الناتجة عن رسم "جوردان-ويغنر" الموزع على أجهزة ذات اتصال بين الجيران الأقرب، تم تنفيذ بنية شبكة "سواب" (SWAP-network).
التحسين (Optimization): تم تطبيق سير عمل تحسين متعدد المراحل لتقليل عمق الدارة وعدد البوابات:
ضغط الشبكة الموترية (Tensor-Network Compression): استُخدم روتين AQC-Tensor لضغط دارة "تروتر" عبر تحسين نموذج (Ansatz) معلمي لزيادة الدقة (Fidelity) مع الحالة المستهدفة.
الترجمة (Transpilation): خضعت الدارات المضغوطة لمسار تحسين قياسي يتضمن Qiskit Level-3، وعمليات pytket، ثم تحسين نهائي من نوع Qiskj Level-3.
تخفيف الخطأ (Error Suppression): تم تطبيق إلغاء الاقتران الديناميكي (DD) وتدوير "باولي" (PT) للتخفيف من فقدان الترابط والأخطاء المتماسكة.
التنفيذ على الأجهزة: أُجريت التجارب على معالج IBM_Marrakesh (بمعمارية Heron r2) باستخدام 12 كيوبت. تم تهيئة النظام في حالتين من حالات الضرب (Product States): موجة كثافة الشحنة (CDW) وجدار النطاق (DW).
المساهمات الرئيسية
كفاءة الموارد: أدى دمج الضغط القائم على الشبكة الموترية مع الترجمة القياسية إلى تقليل كبير في متطلبات الموارد. سجل المؤلفون انخفاضاً بنسبة تقارب 88% في عدد بوابات الكيوبت الثنائية وانخفاضاً بنسبة 87% في عمق الدارة مقارنة بدارات "تروتر" غير المحسنة.
إثبات وجود Stark MBL: نجحت الدراسة في إظهار الانتقال من الديناميكيات الحرارية (عند ميول ضعيفة، Δ0=1,3) إلى سلوك التموضع القوي (عند ميول قوية، Δ0=15) على أجهزة كمومية فعلية.
المقارنة المرجعية (Benchmarking): تمت مقارنة نتائج الأجهزة بدقة مع التقطيع الدقيق (Exact Diagonalization) والمحاكاة الخالية من الضجيج، مما أظهر توافقاً نوعياً وكمياً جيداً مع الملاحظات المقاسة.
النتائج تم توصيف الديناميكيات باستخدام عدة ملاحظات:
عدم التوازن الشحني (Ic): في حالات الميل الضعيفة والمتوسطة، تضاءل عدم التوازن بسرعة، مما يشير إلى التوازن الحراري. أما عند الميل القوي (\Delta__0 = 15)، فقد ظل عدم التوازن مرتفعاً، مما حافظ على ذاكرة حالة الـ CDW الأولية.
معلومات فيشر الكمومية (QFI) وإنتروبيا المشاركة (S2P): في نظام التموضع، تم كبح نمو QFI وإنتروبيا المشاركة بقوة، مما يشير إلى محدودية انتشار الارتباطات وحصر الدالة الموجية في فضاء هيلبرت.
إنتروبيا التشابك: ظلت إنتروبيا التشابك الجزئي منخفضة وأظهرت اعتماداً ضعيفاً على حجم الجزء الفرعي عند الميول القوية، وهو ما يتوافق مع مقياس قانون المساحة (Area-law scaling) النموذجي لـ MBL.
تطور الكثافة: أظهرت قياسات الكثافة الموضعية أن الجسيمات ظلت محصورة في مواقعها الأولية تحت تأثير الميول القوية، بينما انتشرت عبر الشبكة في حالات المجالات الأضعف.
نقل الجسيمات: بالنسبة لحالة البداية (DW)، ظل عدد الجسيمات في النصف المشغول أولياً (Nhalf) ثابتاً عند الميول القوية، مما يؤكد كبح النقل.
الأهمية والادعاءات يدعي البحث أن الأجهزة الكمومية الحالية، عندما تقترن باستراتيجيات تجميع وتحسين فعالة للموارد، قادرة على دراسة الظواهر الديناميكية مثل MBL في الأنظمة الكمومية المتعددة الأجسام المترابطة. ويؤكد المؤلفون أن الطبيعة الخالية من الاضطراب لـ Stark MBL مفيدة بشكل خاص للمحاكاة الكمومية، حيث تلغي الحاجة إلى المتوسطات عبر حالات اضطراب متعددة، مما يقلل من الأعباء التجريبية. يعمل هذا العمل كمعيار لمدى فائدة أجهزة NISQ في استكشاف الديناميكيات الكمومية غير الإرغودية (Non-ergodic)، موضحاً أنه رغم ضجيج الأجهزة، يمكن التقاط السمات الأساسية لتموضع "ستارك" من خلال تنفيذ الدارات المحسنة وتقنيات تخفيف الخطأ. لا تدعي الدراسة حل المخطط الطوري الكامل للنموذج، بل تركز على إثبات كبح النقل والاحتفاظ بالذاكرة في نظام تفاعلي ثابت عند الميل القوي.