Digital Quantum Simulation of Spin Transport

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

يعد فهم ظواهر النقل في أنظمة الغزل الكمومية (quantum spin systems) أمراً حاسماً لتطوير أجهزة السبنترونيكس (spintronic devices) والكيوبتات القائمة على الغزل. وبينما استُخدمت المحاكاة الكمومية لـ دوال الارتباط الذاتي للغزل-الغزل (spin-spin ACF) بنجاح لاستقصاء نقل الغزل، إلا أن المحاكاة القائمة على دالة الارتباط الذاتي لتيار الغزل (spin-current ACF) ظلت غير مستكشفة إلى حد كبير على الأجهزة الكمومية.

• التحدي: مشغل تيار الغزل أكثر تعقيداً في التنفيذ من مشغل الغزل القياسي، مما يؤدي إلى تكلفة بوابات أعلى بكثير.
• محدودية الطرق الحالية: تتطلب الطرق التقليدية مثل اختبار هادامارد (Hadamard test) كيوتبات مساعدة (ancilla qubits) وبوابات وحدوية متحكم بها (مثل بوابات controlled-$J$). بالنسبة لنظام مكون من $N$ جسيم، يتوسع اختبار هادامارد بمعدل $O(N^2)$ من حيث عدد الدارات وعمق الدارة، مما يجعله مكلفاً من حيث الموارد للأجهزة المتاحة حالياً.
• المحدودية الكلاسيكية: تواجه الطرق الكلاسيكية صعوبة في التعامل مع الدوال الموجية متعددة الأجسام بسبب قيود الذاكرة، مما يستدعي حلولاً كمومية للأنظمة الأكبر.

2. المنهجية

يقترح المؤلفون مخطط قياس مباشر يستخدم القياسات في منتصف الدارة (MCMs) والعمليات غير الوحدوية لتقييم دالة الارتباط الذاتي لتيار الغزل دون الحاجة إلى كيوتبات مساعدة.

أ. الإطار النظري

• النموذج: تركز الدراسة على نموذج XXZ هايزنبرغ أحادي البعد ذو الهاملتوني:
  $$H = J \sum_{i=1}^{L} (S^X_i S^X_{i+1} + S^Y_i S^Y_{i+1} + \Delta S^Z_i S^Z_{i+1})$$
• المراقب المستهدف: دالة الارتباط الذاتي لتيار الغزل، وتُعرف بـ $C_J(t) = \sum_{i,j} \langle J_i(t) J_j \rangle / N$، حيث $J_r$ هو تيار الغزل المحلي المشتق من معادلة الاستمرارية: $J_r = i[S^Z_r, H_{r,r+1}]$.
• بروتوكول القياس المباشر: بناءً على المرجع [28]، يقيس المؤلفون $\langle U^\dagger W U G \rangle$ مباشرة.
  • الجزء الحقيقي: يتم تقييمه باستخدام المؤثرات الإسقاطية $(I \pm G)/2$. وبما أن هذه ليست وحدوية، يتم تنفيذها عبر القياسات في منتصف الدارة (MCMs).
  • الجزء التخيلي: يتم تقييمه باستخدام بوابات وحدوية $e^{\pm i\pi G/4}$.
  • التحسين: للتعامل مع حقيقة أن مشغل تيار الغزل $J$ ليس في قاعدة Z، يقدم البروتوكول تحويلات وحدوية ($U_B$) لتدوير القاعدة، مما يسمح باستخدام قياسات قاعدة Z القياسية.

ب. كفاءة الموارد

• التوسع (Scaling): يقلل الأسلوب المقترح من توسع الدارة من $O(N^2)$ (اختبار هادامارد) إلى $O(N)$.
• عدد الدارات:
  • الجزء الحقيقي: يتطلب دارتين فقط (باستخدام نوعين مختلفين من MCMs) لقياس جميع الملاحظات المحلية في وقت واحد.
  • الجزء التخيلي: يتطلب 4 دارات (باستخدام أربع بوابات وحدوية مختلفة).
• تخفيف الضجيج: تستخدم التجارب تقنية إلغاء التماثل الديناميكي (تسلسل XY4) و100 دائرة من نوع Pauli-twirled (PT). يتم إعادة تطبيع النتائج باستخدام نموذج قناة إزالة الاستقطاب لتصحيح الضجيج، من خلال تعلم عامل إزالة الاستقطاب $(1-p)$ من الدارات المرجعية.

ج. الإعداد التجريبي

• الأجهزة: معالج IBM Heron (ibm_kingston) المزود بـ 156 كيوتب ترانسمون ثابت التردد واتصال heavy-hex.
• حجم النظام: تم إجراء عمليات المحاكاة على نظام مكون من 40 كيوتب.
• تحضير الحالة:
  • التحقق: حالة Néel ($|0101...\rangle$).
  • النقل: حالة جدار النطاق (Domain wall state) ($|0\rangle^{\otimes 20} \otimes |1\rangle^{\otimes 20}$) لدراسة ديناميكيات نقل الغزل.
• تحسين الدارة: يستفيد المؤلفون من خاصية المخروط الضوئي (light-cone property) لانتشار المعلومات. حيث يتم حذف البوابات والقياسات خارج المخروط الضوئي الفعال لتقليل عمق الدارة.

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

1. أول عرض لـ ACF تيار الغزل: يقدم هذا العمل أول محاكاة كمومية رقمية لـ ACF تيار الغزل على كيوتبات فائقة التوصيل، متجاوزاً دراسات ارتباط الغزل-الغزل التي كانت مهيمنة سابقاً.
2. بروتوكول قياس فعال: تقديم مخطط قياس مباشر باستخدام MCMs يلغي الحاجة إلى كيوتبات مساعدة ويقلل من تعقيد الدارة من $O(N^2)$ إلى $O(N)$.
3. القدرة على قياس الارتباط في زمنين: تم تعميم البروتوكول لقياس دوال الارتباط في زمنين $\langle J(t_1) J(t_2) \rangle$، وهو أمر ضروري لدراسة الديناميكيات غير المتوازنة.
4. القابلية للتوسع: تم إثبات قابلية الطريقة للتوسع لأنظمة يصل حجمها إلى (40 كيوتب)، وهي أنظمة يصعب محاكاتها كلاسيكياً بدقة عالية باستخدام التقطيع الدقيق (exact diagonalization)، وإن كانت قابلة للمعالجة عبر طريقة MPS للتحقق.

** 4. النتائج**

تحقق المؤلفون من صحة بروتوكولهم مقابل محاكاة Matrix Product State (MPS) واستقصوا ثلاثة أنظمة نقل متميزة عن طريق تغيير معامل التباين $\Delta$:

• التحقق:

  • تم العثور على توافق ممتاز بين النتائج التجريبية على ibm_kingston ومحاكاة MPS لكل من الجزء الحقيقي والتخيلي لـ ACF تيار الغزل المحلي.
  • نجح البروتوكول في رصد خصائص انتشار المخروط الضوئي.

• أنظمة النقل:

  • القريب من البالستي ($\Delta = 0.2$):
    • لوحظ وزن درودي غير متلاشٍ ($\tilde{D}_w \approx 23.7 \times 10^{-4}$).
    • أس معامل انتشار $\alpha \approx 0.827$ (التوقع النظري: 1).
    • انتشر تيار الغزل بفعالية عبر مسافات طويلة.
  • الانتشار الفائق ($\Delta = 1$):
    • كان الوزن الدرودي صغيراً ولكنه موجب ضمن نافذة المحاكاة ($\approx 5.5 \times 10^{-4}$).
    • أس التوسع $\alpha \approx 0.325$، وهو ما يتوافق مع فئة عالمية Kardar–Parisi–Zhang (KPZ) حيث ($\alpha = 1/3$).
  • الانتشار المنتظم ($\Delta = 1.6$):
    • تلاشى الوزن الدرودي ($\approx 0.9 \times 10^{-4}$).
    • أس التوسع $\alpha \approx 0$، مما يشير إلى انتشار طبيعي.
    • تلاشى تيار الغزل بسرعة بعيداً عن جدار النطاق.

• عمق الدارة: بلغ أقصى عمق للدارات حوالي 99 مع حوالي 1,880 بوابة ثنائية الكيوتب، مما يثبت جدوى تشغيل محاكاة نقل معقدة على الأجهزة الحالية قبل مرحلة مقاومة الخطأ (pre-fault-tolerant).

5. الأهمية

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