2025-11-12T10:37:09.722784

Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states

Calcluth, Hahn, Bermejo-Vega et al.
Classically simulating circuits with bosonic codes is challenging due to the prohibitive cost of simulating quantum systems with many, possibly infinite, energy levels. We propose an algorithm to simulate circuits with encoded Gottesman-Kitaev-Preskill (GKP) states, specifically for odd-dimensional encoded qudits. Our approach is tailored to be especially effective in the most challenging but practically relevant regime, where the codeword states exhibit high (but finite) squeezing. Our algorithm leverages the Zak-Gross Wigner function introduced by J. Davis et al. [arXiv:2407.18394], which represents infinitely squeezed encoded stabilizer states positively. The runtime of the algorithm scales with the negativity of the Wigner function, allowing for efficient simulation of certain large-scale circuits - namely, input stabilizer GKP states undergoing generalized GKP-encoded Clifford operations followed by modular measurements - with a high degree of squeezing. For stabilizer GKP states exhibiting 12 dB of squeezing, our algorithm can simulate circuits with up to 1,000 modes with less than double the number of samples required for a single input mode, in stark contrast to existing simulators. Therefore, this approach holds significant potential for benchmarking early implementations of quantum computing architectures utilizing bosonic codes.
academic

المحاكاة الكلاسيكية للدوائر ذات حالات غوتسمان-كيتايف-بريسكيل الفعلية ذات الأبعاد الفردية

المعلومات الأساسية

  • معرّف الورقة: 2412.13136
  • العنوان: Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states
  • المؤلفون: Cameron Calcluth, Oliver Hahn, Juani Bermejo-Vega, Alessandro Ferraro, Giulia Ferrini
  • التصنيف: quant-ph (الفيزياء الكمية)
  • تاريخ النشر: 17 ديسمبر 2024 (نسخة arXiv التمهيدية)
  • رابط الورقة: https://arxiv.org/abs/2412.13136

الملخص

تقترح هذه الورقة خوارزمية محاكاة كلاسيكية لمحاكاة الدوائس الكمية ذات الأبعاد الفردية القائمة على ترميز غوتسمان-كيتايف-بريسكيل (GKP). تم تحسين الخوارزمية بشكل خاص للسيناريو الأكثر تحديًا لكن ذي الصلة العملية: حالات الكلمات الرمزية ذات الضغط العالي لكن المحدود. تستفيد الخوارزمية من دالة Zak-Gross Wigner، التي تمكن التمثيل الإيجابي لحالات المثبت المرمزة ذات الضغط اللانهائي. يتناسب وقت تشغيل الخوارزمية مع سلبية دالة Wigner، مما يتيح المحاكاة الفعالة لبعض الدوائر الكبيرة. بالنسبة لحالات المثبت GKP ذات ضغط 12dB، تستطيع الخوارزمية محاكاة دوائر تصل إلى 1000 وضع، مما يتطلب عدد عينات أقل من ضعف العدد المطلوب لوضع إدخال واحد، وهذا يتناقض بشكل حاد مع المحاكيات الموجودة.

الخلفية البحثية والدافع

تعريف المشكلة

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

تحليل الأهمية

  1. القيمة العملية: تم تحقيق ترميز GKP في الأنظمة التجريبية مثل الدوائر الفائقة الموصلية المقترنة بالتجاويف الميكروية والأيونات المحبوسة والمنصات الفوتونية، وهو مرشح مهم للحوسبة الكمية المتسامحة مع الأخطاء
  2. الأهمية النظرية: توفير أدوات معايرة للتحقق من تحقيق معالجات البوزون الكمية المبكرة
  3. التحديات التقنية: الطرق التقليدية لا تستطيع التعامل مع محاكاة الدوائس الكبيرة لحالات GKP عالية الضغط

قيود الطرق الموجودة

  1. طريقة توسيع أساس Fock: نظرًا لأن كلمات رمز GKP تتضمن عددًا لانهائيًا من مستويات الطاقة، تصبح طرق المحاكاة الغاشمة غير قابلة للتطبيق حتى على عدد قليل من كيوبتات GKP
  2. طرق التوزيع شبه الاحتمالي: تتمتع دالة Wigner لحالات GKP المرمزة بسلبية كبيرة جدًا، مما يؤدي إلى انفجار وقت التشغيل
  3. طرق تحليل الحالة الغاوسية: غير قابلة للتطبيق على حالات GKP عالية الضغط ذات الصلة العملية، لأن هذه الحالات تتمتع بغير-غاوسية كبيرة

المساهمات الأساسية

  1. اقتراح دالة Zak-Gross Wigner متعددة الأوضاع: تعميم دالة Zak-Gross أحادية الوضع إلى حالة n-وضع، وإثبات أنها تحقق البديهيات المعدلة لـ Stratonovich-Weyl
  2. إنشاء العلاقة بين دالة ZGW ودالة Gross Wigner: إثبات المراسلات بين دالة ZGW لحالات CV GKP ودالة Gross Wigner للحالات المنطقية المرمزة
  3. تطوير خوارزمية محاكاة فعالة:
    • توفير خوارزمية محاكاة ضعيفة لحالات المثبت GKP المثالية
    • توفير خوارزمية تقدير احتمالي لحالات GKP ذات الضغط المحدود
  4. تحقيق تحسن أداء كبير: بالنسبة لحالات المثبت GKP ذات ضغط 12dB، تتطلب محاكاة دوائس 1000 وضع عددًا من العينات أقل من ضعف العدد المطلوب لوضع واحد

شرح الطريقة

تعريف المهمة

محاكاة توزيع احتمالية نتائج القياس لدائرة كمية:

  • الإدخال: كيوبتات مرمزة بـ GKP عشوائية، مع معامل فك الضغط Δ
  • التطور: عمليات أحادية موصوفة بمصفوفات辛 صحيحة وإزاحات عشوائية (تشمل جميع عمليات Clifford المرمزة)
  • القياس: قياس معياري في أساس الموضع MZ(s)=1dℓ∑ne−iℓs⋅nT(0,n)M_Z(s) = \frac{1}{d\ell}\sum_n e^{-i\ell s \cdot n}T_{(0,n)}

تعريف دالة ZGW

التعريف 1: تُعرّف دالة Zak-Gross Wigner متعددة الأوضاع ذات الأبعاد الفردية كالتالي: Wρ^(η)=Tr(ρ^A^η)W_{\hat{\rho}}(\eta) = \text{Tr}(\hat{\rho}\hat{A}_\eta)

حيث يكون معامل نقطة الطور: A^η=1(2π)n∑a∈Z2neiℓ[a,η]+iπaXTaZT^a\hat{A}_\eta = \frac{1}{(2\pi)^n}\sum_{a \in \mathbb{Z}^{2n}} e^{i\ell[a,\eta]+i\pi a_X^T a_Z}\hat{T}_a

حيث η∈[0,dℓ)2n\eta \in [0,d\ell)^{2n}، وdd عدد صحيح فردي موجب، وT^a=eiπaXTaZ/dT^aXT^aZ\hat{T}_a = e^{i\pi a_X^T a_Z/d}\hat{T}_{a_X}\hat{T}_{a_Z}.

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

الليما 1: العلاقة بين دالة ZGW ودالة Gross Wigner Wρ^(η)=Wˉρˉ(u)(t)W_{\hat{\rho}}(\eta) = \bar{W}_{\bar{\rho}(u)}(t) حيث η=ℓ(u+t)\eta = \ell(u+t)، وu=1ℓη mod 1u = \frac{1}{\ell}\eta \bmod 1، وt=1ℓη−ut = \frac{1}{\ell}\eta - u.

النظرية 1: تحويل دالة ZGW تحت التطور辛 بالنسبة للتطور الأحادي الغاوسي الموصوف بمصفوفة辛 صحيحة SS، تكون دالة ZGW للحالة المتطورة: WU^Sρ^U^S†(η)=Wρ^(Sη−t)W_{\hat{U}_S\hat{\rho}\hat{U}_S^\dagger}(\eta) = W_{\hat{\rho}}(S\eta - t)

النظرية 2: حساب احتمالية القياس توزيع احتمالية قياس معامل المنطق Z^L⊗n\hat{Z}_L^{\otimes n} هو: Tr(ρ^M^Z(s))=∫dηXWρ^((ηXs))\text{Tr}(\hat{\rho}\hat{M}_Z(s)) = \int d\eta_X W_{\hat{\rho}}\left(\begin{pmatrix}\eta_X\\s\end{pmatrix}\right)

خوارزمية المحاكاة

محاكاة ضعيفة لحالات المثبت GKP المثالية

  1. أخذ عينة من المتجه η\eta من التوزيع الاحتمالي الفعال Wρ^0(η)W_{\hat{\rho}_0}(\eta)
  2. تحويل تحت تأثير مصفوفة辛 SS ومتجه إزاحة tt: Sη−tS\eta - t
  3. تطبيق معامل الإزاحة الخطي: Sη−t+cS\eta - t + c
  4. يتم إعطاء نتيجة القياس بواسطة النصف الخلفي من المتجه

تقدير الاحتمالية لحالات GKP ذات الضغط المحدود

استخدام طريقة أخذ العينات من التوزيع شبه الاحتمالي: Pr(η)=1MU^S,cρ^0U^S,c†∣WU^S,cρ^0U^S,c†(η)∣\text{Pr}(\eta) = \frac{1}{M_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}}|W_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}(\eta)|

يتم إعطاء عدد العينات المطلوبة بواسطة: N=2ϵ2Mρ^02log⁡(2/δ)N = \frac{2}{\epsilon^2}M_{\hat{\rho}_0}^2\log(2/\delta)

حيث Mρ^M_{\hat{\rho}} هي سلبية دالة ZGW (حجم السلبية).

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

إطار التحليل النظري

تركز هذه الورقة بشكل أساسي على التحليل النظري والحسابات العددية، مع التركيز على:

  1. حساب سلبية دالة ZGW: لحالات GKP بمستويات ضغط مختلفة
  2. تحليل تعقيد العينة: بناءً على سلوك تحجيم السلبية
  3. المقارنة مع الطرق الموجودة: مقارنة التعقيد النظري

مؤشرات التقييم

  • السلبية اللوغاريتمية لـ ZGW: log⁡Mρ^\log M_{\hat{\rho}}، تقيس تعقيد المحاكاة
  • تحجيم عدد العينات: نمو عدد العينات بالنسبة لحالة وضع واحد
  • سعة عدد الأوضاع: الحد الأقصى لعدد الأوضاع التي يمكن محاكاتها بفعالية

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

دالة ZGW لحالات GKP ذات الضغط المحدود

من خلال الحساب الصريح، تم الحصول على دالة ZGW لحالة المنطق 0 الفعلية: WGKP((uv))∝ϑ(Γ;z)W_{\text{GKP}}\left(\begin{pmatrix}u\\v\end{pmatrix}\right) \propto \vartheta(\Gamma; z)

حيث z=(v/(dℓ),−u/(dℓ),0,0)Tz = (v/(d\ell), -u/(d\ell), 0, 0)^T، وΓ\Gamma مصفوفة معقدة يحددها معامل الضغط Δ\Delta.

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

بالنسبة لـ Δ=0.25\Delta = 0.25 (يقابل ضغط 12dB):

  • السلبية أحادية الوضع: Mρ^0Δ≈e3×10−4M_{\hat{\rho}_0^\Delta} \approx e^{3 \times 10^{-4}}
  • محاكاة بألف وضع: نظرًا للطبيعة الضربية لسلبية، تتطلب محاكاة 1000 وضع إدخال عددًا من العينات أقل من ضعف العدد المطلوب لوضع واحد
  • تحسن الأداء: تحسن بعدة رتب من حيث الحجم مقارنة بالطرق الموجودة

تحليل الحالات السحرية

بالنسبة لحالة π ∣ψπ⟩=13(∣0L⟩+∣1L⟩−∣2L⟩)|\psi_\pi\rangle = \frac{1}{\sqrt{3}}(|0_L\rangle + |1_L\rangle - |2_L\rangle):

  • تظهر سلبية ملحوظة عند جميع قيم Δ\Delta
  • تكون التكاليف كبيرة عند محاكاة الدوائس الكبيرة، وهذا متوقع

الأعمال ذات الصلة

تطور خوارزميات المحاكاة الكلاسيكية

  1. نظرية Gottesman-Knill: محاكاة دوائس المثبت في وقت متعدد الحدود
  2. طرق شبه الاحتمالية: خوارزميات أخذ العينات بناءً على دالة Gross Wigner
  3. محاكاة الأنظمة البوزونية: استخدام إيجابية التوزيعات شبه الاحتمالية وشبكات الموتر والتمثيلات النجمية

الحالة الحالية لمحاكاة ترميز GKP

تقتصر الطرق الموجودة بشكل أساسي على:

  • حالات GKP المثالية (ذات الضغط اللانهائي)
  • العمليات الغاوسية وقياس الضربة الصفرية
  • عدم القدرة على التعامل مع الحالات ذات الضغط المحدود ذات الصلة العملية

نقاط الابتكار في هذه الورقة

مقارنة بالأعمال الموجودة، تحقق هذه الورقة للمرة الأولى محاكاة فعالة لحالات GKP ذات الضغط المحدود، خاصة مع تحسن أداء كبير في نظام الضغط العالي.

الخلاصة والمناقشة

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

  1. المساهمات النظرية: إنشاء إطار نظري كامل لدالة ZGW، بما في ذلك التعميم متعدد الأوضاع والقوانين التطورية
  2. الابتكار الخوارزمي: اقتراح خوارزمية محاكاة فعالة قابلة للتطبيق على حالات GKP الفعلية
  3. القيمة العملية: توفير أداة معايرة قوية للتحقق من معالجات البوزون الكمية

القيود

  1. قيد الأبعاد الفردية: الطرق الحالية تنطبق فقط على كيوبتات GKP ذات الأبعاد الفردية، ولا يزال التوسع إلى الأبعاد الزوجية (مثل الكيوبتات) مسألة مفتوحة
  2. قيود العمليات: محدودة بالعمليات الموصوفة بمصفوفات辛 صحيحة، لا تستطيع التعامل مع الضغط أحادي الوضع العشوائي
  3. نموذج الضوضاء: لم يتم النظر في تأثير الضوضاء الفعلية مثل الخسارة على سلبية دالة ZGW

الاتجاهات المستقبلية

  1. التوسع إلى كيوبتات GKP: تطوير طرق قابلة للتطبيق على الحالة ثنائية الأبعاد
  2. ترميزات بوزونية أوسع: تعميم الإطار على ترميزات بوزونية أخرى خارج نطاق GKP
  3. إطار شبه الاحتمالية: تطوير نظرية شبه احتمالية شاملة متخصصة للترميزات البوزونية

التقييم المتعمق

المميزات

  1. الصرامة النظرية: الاشتقاقات الرياضية كاملة والإثباتات مفصلة والأساس النظري متين
  2. قوة التطبيق العملي: حل مشكلة محاكاة حالات GKP ذات الضغط المحدود ذات الصلة العملية
  3. الأداء المتميز: تحقيق تحسن في الأداء بعدة رتب من حيث الحجم في نظام الضغط العالي
  4. الابتكار الطريقي: التعميم متعدد الأوضاع لدالة ZGW وتطبيقاتها مساهمة نظرية مهمة

أوجه القصور

  1. نطاق التطبيق محدود: ينطبق فقط على الأنظمة ذات الأبعاد الفردية، مما يحد من الشمولية
  2. نقص التحقق التجريبي: يركز بشكل أساسي على العمل النظري، ويفتقر إلى مقارنات مع الأنظمة الكمية الفعلية
  3. تحليل التعقيد: تحليل التعقيد الحسابي الفعلي للدوائس الكبيرة ليس مفصلاً بما يكفي

التأثير

  1. القيمة الأكاديمية: توفير أدوات ووجهات نظر جديدة لنظرية المحاكاة الكمية
  2. الآفاق العملية: ذات أهمية كبيرة للتجارب الكمية المبكرة المتسامحة مع الأخطاء
  3. الدفع التقني: قد تلهم تطوير المزيد من طرق محاكاة الترميز البوزوني

السيناريوهات القابلة للتطبيق

  • معايرة والتحقق من معالجات البوزون الكمية
  • التحليل النظري لخطط الحوسبة الكمية المتسامحة مع الأخطاء بترميز GKP
  • أبحاث المحاكاة الكلاسيكية للحوسبة الكمية ذات المتغيرات المستمرة

المراجع

تستشهد هذه الورقة بـ 68 مرجعًا ذا صلة، تشمل بشكل أساسي:

  • الأعمال الأصلية لترميز Gottesman-Kitaev-Preskill
  • الأدبيات الكلاسيكية لخوارزميات المحاكاة الكمية
  • أحدث التطورات في الحوسبة الكمية البوزونية
  • الأساس النظري لدوال Wigner والتوزيعات شبه الاحتمالية

التقييم الشامل: هذه ورقة فيزياء نظرية عالية الجودة حققت اختراقًا مهمًا في المحاكاة الكلاسيكية لدوائس ترميز GKP الكمية. على الرغم من وجود بعض القيود، فإن مساهماتها النظرية وقيمتها العملية كبيرة جدًا، وتتمتع بأهمية كبيرة لمجال الحوسبة الكمية.