Approximate simulation of complex quantum circuits using sparse tensors

ملخص تقني: المحاكاة التقريبية للدوائر الكمومية المعقدة باستخدام الموترات المتفرقة

بيان المشكلة
تعد محاكاة الدوائر الكمومية على الحواسيب الكلاسيكية أمراً ضرورياً لتحديد حدود التفوق الكمومي القابل للتحقق، وحل مشكلات الأنظمة الكمومية متعددة الأجسام، وتوجيه تطوير الأجهزة. وبينما أصبحت طرق الشبكات الموترية (مثل حالات المصفوفة الناتجة عن منتج المصفوفات MPS، ومجموعة إعادة التطبيع لـ DMRG) هي الأحدث في محاكاة الدوائر المعقدة، إلا أنها غالباً ما تعتمد على تماثلات كامنة أو أبعاد منخفضة لتحقيق الكفاءة. علاوة على ذلك، فإن المحاكاة الكلاسيكية العامة مقيدة بنظرية التعقيد ($P \neq NP$)، مما يمنع وجود محاكيات عالمية ذات زمن تشغيل متعدد الحدود. حدد المؤلفون حاجة لأدوات محاكاة يمكنها تقريب الحالات الكمومية دون الحاجة إلى تفكيك القيم المفردة (SVD)، أو التماثلات الكامنة، أو الأبعاد المنخفضة، مستهدفةً بشكل خاص الدوائر التي قد تواجه فيها الشبكات الموترية التقليدية صعوبات.

المنهجية: محاكاة الموتر المتفرق المبتور (TruSTS)
يقدم البحث منهجية TruSTS، وهي طريقة للمحاكاة التقريبية للدوائر الكمومية باستخدام موترات متفرقة التوزيع. يعتمد النهج الجوهري على تمثيل الحالة الكمومية لـ $N$ من الكيوبتات $|\psi\rangle$ كحالة متفرقة $|\phi\rangle$، وهي عبارة عن تركيب خطي لـ $k$ من حالات الأساس ($k \le 2^N$).

• هيكل البيانات: يتم تخزين الحالة المتفرقة باستخدام مصفوفتين ثابتتي الطول بحجم $k$:
  1. $x[x_i]$: مصفوفة من الإحداثيات الصحيحة التي تمثل السلاسل الثنائية لحالات الأساس.
  2. $\alpha[\alpha_i]$: مصفوفة من سعات الاحتمال المركبة المقابلة لتلك الحالات.
  3. مؤشر صحيح nnz يتتبع عدد العناصر المملوءة حالياً.
• خوارزمية الانكماش بتنسيق البت (Bitwise Contraction): لتطبيق بوابة كمومية (ممثلة بموتر $\hat{U}$)، تستخدم الخوارزمية استراتيجية فرز بتات تهدف إلى محاكاة ضرب مصفوفة الصف المترابط (CSR)، وتتضمن الخطوات التالية:
  1. تحديد قناع بت (bitmask) لاختيار الكيوبتات المرتبطة بالبوابة.
  2. فرز مصفوفة الحالة بناءً على "البتات غير المرتبطة" (متمم القناع). يؤدي هذا إلى تجميع حالات الأساس في مصفوفات فرعية مستقلة حسابياً حيث تعمل البوابة بشكل متطابق.
  3. تطبيق البوابة عبر عمليات انكماض مصفوفة $4\times4$ مستقلة لكل مصفوفة فرعية.
  4. يمكن لهذه العملية أن تزيد عدد الحدود غير الصفرية (nnz) بمعامل يصل إلى أربعة لكل تطبيق للبوابة.
• البتور (Truncation): لمنع النمو الأسي لـ nnz، تقوم الخوارزمية ببتور الحالة للعودة إلى الحجم $k$ بعد كل بوابة. قارن البحث بين استراتيجيتين للبتور:
  • Top-k: تحتفظ بـ $k$ من الحدود ذات أكبر سعات احتمالية.
  • Random-k: تحتفظ بـ $k$ من الحدود المختارة عشوائياً.
  • يتم تطبيق عامل إعادة تطبيع $\gamma$ على $|\phi\rangle$ لمعالجة الكتلة الاحتمالية المفقودة أثناء البتور.

المساهمات الرئيسية والاستنتاجات النظرية
اشتق المؤلفون علاقات بين معامل إعادة التطبيع $\gamma$ ودقة المحاكاة $f = |\langle \psi | \phi \rangle|^2$.

• حدود الدقة: أثبتوا أن الدقة المتوقعة $\bar{f}$ للدوائر الكمومية العشوائية تساوي تقريباً $\gamma^2$ (الكتلة الاحتمالية المحتفظ بها)، مع انحراف إيجابي طفيف بسبب حدود التداخل.
• الحدود العليا والدنيا:
  • تم اشتقاق حد أدنى $\bar{f}_{min} \approx 1/2^N$ للحالة التي يتم فيها بتور الحالة إلى حالة أساس عشوائية واحدة.
  • تم اشتقاق حد أقصى $f_{max}$ بافتراض أن توزيع المخرجات يتبع توزيع "بورتر-توماس" (Porter-Thomas). يمثل هذا الحد أفضل دقة ممكنة إذا كان بإمكان المرء تحديد والاحتفاظ بـ $k$ من أكثر السلاسل الثنائية احتمالاً بعد تنفيذ الدائرة، بدلاً من البتور في منتصف الدائرة. لاحظ البحث أن $f_{max} \approx d(1 - \ln d)$، حيث $d = k/2^N$ هو كسر البتور.

النتائج
قام المؤلفون باختبار أداء TruSTS على دوائر كمومية عشوائية تتكون من وحدات ثنائية الكيوبت (Haar-random) موزعة في طبقات.

• الدقة مقابل البتور: أظهرت عمليات المحاكاة على $N=24$ كيوبت أن دقة الحالة تتدرج تقريباً مع مربع عامل إعادة التطبيع $\gamma^2$. وقد تفوق بتر Top-k باستمرار على بتر Random-k من خلال تركيز الاحتمالية في الحدود المحتفظ بها.
• تدرج وقت التشغيل:
  • مقابل $k$: يتدرج وقت التشغيل لكل بوابة خطياً مع عدد الحدود المحتفظ بها $k$.
  • مقابل $N$: بالنسبة لـ $k \ll 2^N$، يكون وقت التشغيل لكل بوابة ثابتاً تقريباً بالنسبة لعدد الكيوبتات $N$ (تم اختباره حتى $N=64$). يشير هذا إلى أن الخوارزمية لا تفرض عبئاً إضافياً يعتمد على حجم فضاء هيلبرت الكامل، طالما ظلت الحالة متفرقة.
  • الذاكرة: تتجنب الطريقة تخصيص ذاكرة كاملة بحجم $2^N$، حيث تخزن $k$ من الحدود فقط.
• المقارنة مع MPS: أظهرت مقارنة أولية مع حالات المصفوفة الناتجة عن منتج المصفوفات (MPS) باستخدام مكتبة quimb على $N=24$ كيوبت (5 طبقات) أنه بالنسبة للمعايير المختبرة، توفر MPS مقايضة أكثر ملاءمة بين وقت التشغيل والدقة عبر معظم نطاق البيانات، حيث تحافظ على دقة أعلى لأوقات تشغيل مماثلة. ومع ذلك، تُظهر TruSTS اضمحلالاً ثابتاً ويمكن التنبؤ به في الدقة مع زيادة وقت التشغيل، لتصل إلى حد أدنى عند $f \sim 10^{-8}$.

الأهمية والادعاءات
يضع البحث TruSTS كأداة لتوسيع قدرات المحاكيات القائمة على الشبكات الموترية، لا سيما في السيناريوهات التي:

1. يغيب فيها التماثل: لا تعتمد الطريقة على تماثلات كامنة.
2. يتم تجنب SVD: تحقق البتور دون التكلفة الحسابية لتفكيك القيم المفردة.
3. توجد فيها تفرقة متأصلة: تعد جذابة للحالات الكمومية التي تكون متفرقة بطبيعتها.

يزعم المؤلفون أن TruST توفر مقايضة مفيدة بين الدقة ووقت التشغيل، مما قد يجد تطبيقاً فورياً في الدوائر التي لا تتطلب دقة عالية جداً. ويؤكدون أن هذا العمل يعد أساساً لأبحاث مستقبلية في تحسين الشبكات الموترية المتفرقة ودمج TruSTS مع استراتيجيات تحسين ترتيب انكماش الشبكة الموترية لدفع حدود المحاكاة الكلاسيكية للأنظمة الكمومية المعقدة. كما يظل المؤلفون متواضعين في طرحهم، مشيرين إلى أن المقارنة الصارمة مع الطرق البديلة والتحقيقات في الكميات القابلة للقياس (مثل دقة التقاطع في الإنتروبي) هي مجالات هامة للعمل المستقبلي.