Clifft: Fast Exact Simulation of Near-Clifford Quantum Circuits

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

تواجه المحاكاة الكلاسيكية للدوائر الكمومية المتسامحة مع الأخطاء (fault-tolerant) مقايضة جوهرية بين القابلية للتوسع والدقة:

• دوائر كليفورد النقية (Pure Clifford Circuits): يمكن محاكاتها بكفاءة باستخدام أدوات مثل Stim باستخدام صيغة الموحد (stabilizer formalism) (مبرهنة غوتسمان-كنيل)، مما يسمح بمليارات الجولات (shots). ومع ذلك، فإن الحوسبة الكمومية الشاملة تتطلب بوابات غير كليفورد (مثل بوابات $T$)، مما يكسر هذه الكفاءة.
• طرق متجه الحالة الكثيفة (Dense State Vector Methods): دقيقة ولكنها تتوسع أسياً مع إجمالي عدد الكيوبتات ($N$)، مما يجعلها غير مجدية للدوائر الضخمة المتسامحة مع الأخطاء (مئات الكيوبتات).
• محاكيات "قرب كليفورد" الموجودة (Existing Near-Clifford Simulators): تحاول طرق مثل SOFT (المعتمدة على وحدات معالجة الرسومات لتعميم الموحدات) و Tsim (رتبة الموحد - stabilizer-rank) سد هذه الفجوة. ومع ذلك، فإنها غالباً ما تعاني من تكلفة عالية لكل جولة بسبب التحديثات الديناميكية لجدول (tableau) الموحد أو التوسع الأسي مع إجمالي عدد العمليات غير الكليفوردية ($T$-count)، مما يحد من قدرتها على محاكاة البروتوكولات واسعة النطاق مثل زراعة الحالة السحرية (Magic State Cultivation - MSC) بشكل كامل.

تحديداً، كانت محاكاة بروتوكول MSC الكامل (بما في ذلك "مرحلة الهروب" إلى كود سطح كبير) مستحيلة بالنسبة للمحاكيات الدقيقة بسبب الحجم الهائل (463 كيوبت فيزيائي) والحاجة إلى تريليونات الجولات لتقدير معدلات الخطأ المنطقية المنخفضة.

2. المنهجية: تمثيل الحالة المجزأ للإطار (Frame-Factored State Representation)

الابتكار الجوهري في Clifft هو تمثيل حالة هجين يفصل بين تطور الإحداثيات الحتمي وتطور السعة العشوائي. يقوم بتجزئة الحالة الكمومية $|\psi(t)\rangle$ إلى ثلاثة مكونات:

$$|\psi(t)\rangle = \gamma(t) U_C^{(t)} \tilde{P}^{(t)} \left( |\phi(t)\rangle_A \otimes |0\rangle_D \right)$$

حيث:

1. إطار كليفورد غير المتصل بالزمن (Offline Clifford Frame - $U_C^{(t)}$): وحدة (unitary) حتمية تمثل التأثير التراكمي لجميع عمليات كليفورد. يتم حسابها مسبقاً (وقت التجميع) وتظل ثابتة عبر جميع جولات المحاكاة.
2. إطار باولي الافتراضي (Virtual Pauli Frame - $\tilde{P}^{(t)}$): مؤثر باولي خفيف الوزن يعتمد على الجولة، يتتبع تقلبات الطور وتقلبات البت. يتم تحديثه عبر عمليات بتية (bitwise) أثناء وقت التشغيل.
3. متجه الحالة النشط (Active State Vector - $|\phi(t)\rangle_A$): متجه حالة كثيف بأبعاد $2^k$، حيث $k$ هو البعد الافتراضي النشط. يتتبع هذا المتجه فقط درجات الحرية غير الكليفوردية.
   • الكيوبتات الخاملة (Dormant Qubits - $D$): الكيوبتات في حالة $|0\rangle$ في القاعدة الافتراضية، ولا تتطلب تخزيناً كثيفاً.
   • الكيوبتات النشطة (Active Qubits - $A$): الكيوبتات المشاركة في التشابك غير الكليفوردي أو التراكب.

الآليات الرئيسية:

• رسم خريطة هايزنبرغ (Heisenberg Mapping): يتم استيعاب بوابات كليفورد الفيزيائية داخل $U_C$. ويتم رسم العمليات غير الكليفوردية في القاعدة الافتراضية.
• توطين باولي (Pauli Localization): تقوم خوارزمية جشعة بتحويل مولدات باولي متعددة الكيوبتات إلى مؤثرات أحادية الكيوبت. إذا كان الكيوبت المستهدف خاملاً، يتم ترقيته إلى المجموعة النشطة (مما يزيد $k$)؛ وإذا كان نشطاً، يتم تدويره مباشرة.
• البعد الديناميكي ($k$): يتوسع البعد $k$ عندما تخلق البوابات غير الكليفوردية تشابكاً، وينكمش عندما تؤدي القياسات إلى انهيار الحالة. بالنسبة لبروتوكولات قرب كليفورد، يكون $k_{max}$ (ذروة البعد النشط) غالباً أصغر بكثير من إجمالي عدد الكيوبتات $N$.

3. نموذج التنفيذ: تجميع لمرة واحدة، أخذ عينات لمرات عديدة

يتبنى Clifft خط أنابيب تنفيذ يتكون من مرحلتين يشبه Stim ولكن مع توسيع ليشمل البوابات غير الكليفوردية:

1. المُجمع (Offline - وقت التجميع):

   • يقبل دوائر متوافقة مع Stim (موسعة ببوابات غير كليفوردية).
   • ينفذ رسم خريطة هايزنبرغ لاستيعاب بوابات كليفورد في الإطار.
   • ينفذ توطين باولي لتحديد جدول المجموعة النشطة.
   • يولد كود بايت (bytecode) مُحسّن لآلة شرويدنجر الافتراضية (Schrödinger Virtual Machine - SVM).
   • النتيجة: هندسة كليفورد وجدول المجموعة النشطة ثابتان قبل بدء أخذ العينات.

2. وقت التشغيل (Online):

   • ينفذ كود البايت المُجمع مسبقاً لكل جولة.
   • تُختزل العمليات إلى:
     • تحديثات بتية لإطار باولي.
     • أخذ عينات متفرقة للضوضاء.
     • عمليات جبر خطي كثيفة فقط على متجه الحالة النشط بحجم $2^{k_{max}}$.
   • يستخدم تقنية SIMD (تعليمات واحدة لبيانات متعددة) لعمليات المصفوفة النشطة و OpenMP للتوازي عندما يكون $k$ كبيراً.

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

• بنية مبتكرة: قدم تمثيل الحالة المجزأ للإطار، مما نقل التكلفة الأسية من إجمالي الكيوبتات ($N$) إلى ذروة البعد الافتراضي النشط ($k_{max}$).
• محاكي مفتوح المصدر (Clifft): حزمة Python/C++ بواجهة برمجة تطبيقات تشبه Stim تدعم الضوضاء، والقياسات وسط الدائرة، والتحكم الكلاسيكي.
• أول محاكاة شاملة لـ MSC: نجح في محاكاة بروتوكول زراعة الحالة السحرية (MSC) الكامل (بما في ذلك مرحلة الهروب) باستخدام 463 كيوبت فيزيائي و $k_{max}=10$، مع تشغيل أكثر من مئات المليارات من الجولات على معالجات عادية.
• تحسين الأداء: أثبت أنه من خلال التجميع المسبق لتحويلات كليفورد، فإن تكلفة الجولة الواحدة تهيمن عليها عمليات $O(2^{k_{max}})$ بدلاً من تحديثات الجدول $O(N^2)$ المطلوبة في محاكيات الموحد المتفرقة.

5. النتائج والاختبارات المعيارية

• نطاق كليفورد النقي: Clifft أبطأ بنحو 10 أضعاف من Stim (بسبب التكاليف الإضافية) ولكنه يظل منافساً.
• نطاق قرب كليفورد (السحر المنخفض): يتفوق Clifft بشكل كبير على الأدوات الموجودة.
  • في دائرة زراعة الحالة السحرية $d=3$، حقق Clifft إنتاجية أعلى بمقدار 370 ضعفاً من Tsim.
  • في دائرة $d=5$، فشل Tsim في التجميع ضمن ميزانية دقيقتين، بينما حافظ Clifft على سرعة 314,000 جولة/ثانية.
• النطاق الكثيف: في الحالة الأسوأ ($k_{max} = N$)، يعمل Clofft ضمن عامل ثابت من محاكيات متجه الحالة الكثيفة الرائدة (Qiskit-Aer, Qulacs, qsim) في اختبارات حجم الكم (Quantum Volume).
• نتائج زراعة الحالة السحرية (MSC):
  • كفاءة التكلفة: حقق Clifft تقديرات لمعدل الخطأ المنخفض مقاربة لعنقود مكون من 16 وحدة معالجة رسومات (SOFT) باستخدام مثيل CPU واحد، مما قلل ساعات العمل الآلي بنسبة ~32 ضعفاً.
  • تحليل فجوة T/S: كشفت عمليات المحاكاة أن التباين بين دائرة بوابة $T$ الحقيقية ووكيل بوابة $S$ (المستخدم في العمل السابق) يتم حجبه عند عتبات فجوة فك التشفير المنخفضة بواسطة إخفاقات فك تشفير مرحلة الهروب. ومع ذلك، عند العتبات العالية (التي تستبعد إخفاقات فك التشفير)، يقترب سلوك البروتوكول الكامل من التباين الكبير الملاحظ في مراحل الزراعة وحدها (يصل إلى نسبة خطأ 30 ضعفاً).

6. الأهمية

• سد الفجوة: يشغل Clifft "نقطة مثالية" بين محاكيات الموحد السريعة ولكن التقريبية، ومحاكيات متجه الحالة الكثيفة الدقيقة ولكن البطيئة. إنه يتيح المحاكاة الدقيقة للدوائر الضخمة المتسامحة مع الأخطاء التي كانت غير قابلة للمعالجة سابقاً.
• التحقق من البروتوكولات: من خلال تمكين المحاكاة الدقيقة لـ MSC من البداية إلى النهاية، يوفر Clifft تحققاً حاسماً للبروتوكولات المتسامحة مع الأخطاء، كاشفاً أن الدوائر الوكيلة (مثل تقريبات بوابة $S$) قد تقلل بشكل كبير من تقدير معدلات الخطأ في سيناريوهات معينة.
• القابلية للتوسع: القدرة على تشغيل المحاكاة على معالجات عادية بدلاً من الحاجة إلى عناقيد GPU مكلفة تساهم في دمقرطة الوصول إلى محاكاة الدوائر الكمومية عالية الدقة، مما يسهل تصميم وتصحيح أكواد تصحيح الخطأ الكمومي المستقبلية.
• إطار المجمع: يوفر تمثيل هايزنبرج الوسيط (HIR) الذي تم تطويره لـ Clifft أساساً جديداً لتحسين وتجميع البرامج الكمومية المبكرة المتسامحة مع الأخطاء بما يتجاوز مجرد المحاكاة.