Classical Simulations of Low Magic Quantum Dynamics

ملخص تقني: المحاكاة الكلاسيكية لديناميكيات الكم ذات السحر المنخفض

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

المنهجية
يطور المؤلفون خوارزمية محاكاة كلاسيكية حتمية تعتمد على تفكيك المثبت منخفض الرتبة (LRSD). وبخلاف الأساليب القائمة على أخذ العينات، تحافظ هذه الطريقة على تمثيل دقيق لمصفوفة الكثافة $\rho(t) = |\psi(t)\rangle\langle\psi(t)|$ طوال تطور الدائرة.

1. تمثيل LRSD: يتم تمثيل الحالة كمجموع فوق مؤثرات منطقية تعمل على حالة مثبت واحدة:
   $$|\psi\rangle\langle\psi| = \sum_{l \in \mathcal{L}_{LRSD}} \lambda_l \sigma_l \rho_S$$
   حيث $\rho_S$ هي حالة مثبت معرفة بواسطة مجموعة المثبت $S$، و $\sigma_l$ هي مؤثرات باولي المنطقية التي تتبادل مع $S$، و $\lambda_l$ هي معاملات حقيقية. يعزل هذا الهيكل المحتوى غير المثبت في مجموعة مدمجة من المؤثرات المنطقية.

2. التحديثات الديناميكية: تقوم الخوارزمية بتحديث هذا التفكيك تحت ثلاثة أنواع من العمليات:

   • وحدات كليفورد (Clifford Unitaries): تقوم برسم خرائط لمؤثرات باولي إلى مؤثرات باولي، مما يحدث مجموعة المثبت والمؤثرات المنطقية دون زيادة عدد الحدود.
   • قياسات باولي (Pauli Measurements): اعتمادًا على علاقات التبادل بين مؤثر القياس $P$ ومجموعة المثبت/المؤثرات المنطقية، تقوم الخوارزمية إما بتحديث الحالة مباشرة، أو استبعاد المؤثرات المنطقية المتضادة، أو إجراء تفكيك للمثبت للتعامل مع الحالات المتضادة. يمكن للقياسات تقليل عدد المؤثرات المنطقية، مما يعادل النمو الناجم عن بوابات غير كليفوردية.
   • بوابات T (غير كليفوردية): يتم تفكيك بوابة T إلى أجزاء تتبادل أو تتضاد مع مؤثرات باولي. يمكن لتطبيق بوابة T أن يؤدي إلى "تفرع" المؤثرات المنطقية، مما قد يضاعف عدد الحدود ثلاث مرات في أسوأ الحالات (الحالتان II و IV) أو يضاعفها مرتين (الحالة III).

3. القطع (Truncation) وقياس السحر: لإدارة التكلفة الحسابية، يمكن إدخال معامل قطع $\epsilon$ لاستبعاد الحدود التي يكون فيها $|\lambda_l| < \epsilon$. يعرّف المؤلفون العدمية المثبتة ($M$) كمقياس للسحر: $M(|\psi\rangle) = L - \log_2 |S(|\psi\rangle)|$. ولحساب $M$ للأنظمة الكبيرة، قاموا بتنفيذ خوارزمية أخذ عينات بيل (Bell sampling) قابلة للتوسع (مقتبسة من المرجع [48]) تقوم بأخذ عينات من سلاسل باولي من الحالة $|\psi\rangle \otimes |\psi\rangle$ لتقدير أبعاد الفضاء الذي تشغله فروق بيل.

4. نماذج الدوائر: تم اختبار الخوارزمية على نماذج الدوائر العشوائية شاملة الاتجاه (single-pair all-to-all)، حيث تعمل بوابات ثنائية الكيوبت (CZ) والقياسات أحادية الكيوبت على كيوبتات مختارة عشوائيًا. تمت دراسة متغيرين:

   • قياسات أساس Z: قطرية في قاعدة الحوسبة، مما يسمح بنقل بوابات T إلى نهاية الدائرة.
   • قياسات أساس X: غير قطرية، مما يمنع نقل بوابات T ويعمل على تعديل التشابك بشكل نشط.

النتائج الرئيسية

1. تعقيد المساحة والكفاءة:

   • تظهر الخوارزامية تدرجًا متعدد الحدود ($N \sim L^\alpha$ مع $\alpha \approx 2$) في عدد مدخلات الذاكرة المطلوبة لتحديد الحالة، بشرط أن يتدرج معدل بوابات T بشكل دون استثنائي (على سبيل المثال، $p_T \sim 1/L$ أو $1/L^2$).
   • يتناقض هذا مع التدرج الأسي لطرق متجه الحالة الكاملة ($2^L$).
   • تظل الكفاءة قائمة حتى مع استبدال بوابات CZ ببوابات كليفورد عشوائية، مما يثبت متانة هيكل الدائرة.
   • التكلفة المتوسطة للمجموعة لا تهيمن عليها "المسارات الصعبة" النادرة ذات التكلفة الأسية؛ حيث يمتلك توزيع المؤثرات المنطقية ذيلًا ذا قانون قوة ولكنه يظل قابلاً للإدارة للأنظمة المدروسة ($L \le 256$).

2. انتقالات الطور المستحثة بالقياس (MIPTs):

   • انتقال التنقية: في نموذج أساس X، حدد المؤلفون انتقال تنقية حيث ينتقل النظام من طور مختلط (زمن تنقية أسي) إلى طور نقي (زمن تنقية ثابت). وُجد أن معدل القياس الحرج هو $p_{cp}^m \approx 0.26(1)$، مع أس ديناميكي $z_p \approx 0.22(2)$ وأس طول الارتباط $\nu_p \approx 0.45(5)$. وُجد أن هذا الانتقال غير حساس لمعدل بوابات T.
   • انتقال السحر: في نموذج أساس Z، لوحظ انتقال طور متميز في "السحر" (العدمية المثبتة). ينتقل النظام بين أطوار سحر قانون المساحة (قابلة للمحاكة بكفاءة) وأطوار سحر قانون القوة/دون الاستثنائي اعتمادًا على أس التدرج لـ كثافة بوابات T وهو $\beta$ (حيث $p_T = \eta/L^\beta$) ومعدل القياس.
   • التشابك مقابل السحر: تكشف الدراسة عن انفصال بين انتقالات التشابك وانتقالات السحر. في نموذج أساس Z، يكون التشابك متدهورًا بالنسبة لمعدلات بوابات T (بسبب التبادل)، بينما يتدرج السحر بشكل غير خطي. تم تحديد أربعة أطوار متميزة بناءً على الجمع بين تدرج التشابك (قانون المساحة/الحجم) وتدرج السحر (قانون المساحة/القوة).

3. الأداء الخوارزمي:

   • تصل طريقة أخذ عينات بيل إلى العدمية المثبتة الدقيقة في وقت متعدد الحدود لأنظمة تصل إلى $L=256$.
   • يقلل القطع باستخدام $\epsilon$ متوسط ($\sim 10^{-2}$) عدد المؤثرات المنطقية بشكل كبير (بعدة مراتب قدرية) مع الحفاظ على دقة الملحوظات غير التافهة مثل العدمية المثبتة.

الأهمية والادعاءات
يزعم البحث تقديم منهجية جديدة لمحاكة الدوائر الكمومية القريبة من حد المثبت، وهي مكملة لأساليب MPS. تكمن أهميتها الأساسية في:

• تطور مصفوفة الكثافة الدقيق: بخلاف أساليب أخذ العينات، تحافظ الطريقة على مصفوفة الك densité كاملة، مما يتيح حساب الملحوظات غير الخطية (التشابك، التنقية) التي لا يمكن الوصول إليها عبر المحاكيات التي تعطي احتمالات المخرجات فقط.
• التعامل مع تشابك قانون الحجم: تنجح في محاكة الأنظمة ذات تشابك قانون الحجم، وهو النظام الذي تفشل فيه طرق MPS، بشرط بقاء "السحر" (موارد غير كليفوردية) منخفضًا.
• توصيف أطوار جديدة: تحدد وتصف "انتقال السحر" المتميز عن انتقال التشابك، مما يقدم منظورًا قائمًا على نظرية الموارد لصعوبة الحوسبة في الدوائر المراقبة.
• القابلية للتوسع: تثبت الخوارثمة أن المحاكاة التفصيلية للمسارات لمصفوفات الكثافة الكاملة أمر ممكن لأنظمة كبيرة ($L \sim 100-250$) في طور قانون المساحة للسحر، مما يسد الفجوة بين المحاكاة الدقيقة صغيرة النطاق والطرق التقريبية واسعة النطاق.

يشير المؤلفون إلى القيود: تقتصر المحاكة الفعالة على طور قانون المساحة للسحر أو بالقرب من معدل القياس الحرج. هم لا يدعون حل الحالة العامة لديناميكيات السحر العالي والتشابك العالي، بل يقدمون أداة لهذا النطاق المحدد من ديناميكيات السحر المنخفض والتشابك العالي ذات الصلة بتصحيح الخطأ الكمومي (QEC) وانتقالات الطور المستحثة بالقياس (MIPTs).