Real-time Krylov Diagonalisation for Open Quantum Systems

ملخص تقني: قطرية كريفلوف في الوقت الحقيقي للأنظمة الكمومية المفتوحة

بيان المشكلة
تتناول الورقة التحدي المتمثل في محاكاة الأنظمة الكمومية المفتوحة التي توصف بمعادلات ليندبلاد الرئيسية (Lindblad master equations)، مع التركيز بشكل خاص على تقدير فجوة ليوفيليان (Liouvillian gap). وبينما برزت طريقة قطرية الفضاء الجزئي الكمومي (QSD) كخوارزمية رائدة للأنظمة المغلقة (الهاملتونية) غير المتوافقة مع الخطأ، فإن تطبيقها على الأنظمة المبددة (dissipative systems) ليس بالأمر الهين. فالأنظمة المفتوحة تحكمها مؤثرات فائقة (superoperators) من نوع ليوفيليان غير هيرميتية وغير طبيعية (non-normal)، مما يؤدي إلى اختلافات هيكلية مثل المتجهات الذاتية ثنائية التناظر (bi-orthogonal eigenvectors)، والقيم الذاتية المركبة، والديناميكيات غير الوحدوية. هذه السمات تعقد طرق فضاء كريفلو (Krylov subspace) القياسية، لا سيما فيما يتعلق بالاستقرار العددي (سوء التكييف) وتفسير المرشحات الطيفية. يهدف المؤلفون إلى تكييف طريقة QSD في الوقت الحقيقي لتقدير فجوة ليوفيليان، وهي معلمة حرجة تحدد زمن الاسترخاء ووجود الأنماط البطيئة في التحولات الطورية المبددة.

المنهجية
يقترح المؤلفون تعديلاً لطريقة قطرية الفضاء الجزئي الكمومي (QSD) في الوقت الحقيقي للتعامل مع ديناميكيات ليوفيليان. وتتضمن المنهجية الجوهرية ما يلي:

1. بناء فضاء كريفلو الجزئي: بدلاً من تطوير حالة تحت هاملتوني وحدوي $H$، يتطور النظام تحت مؤثر ليوفيليان $\mathcal{L}$. يتم إنشاء فضاء كريفلو الجزئي عن طريق تطبيق قوى مؤثر التطور الزمني $e^{\mathcal{L}\tau}$ على حالة ابتدائية $|\rho_0)$:
   $$\text{span}\{|\rho_0), e^{\mathcal{L}\tau}|\rho_0), \dots, e^{\mathcal{L}(D-1)\tau}|\rho_0)\}$$
2. مسألة القيم الذاتية المعممة (GEVP): بما أن الدارات الكمومية لا يمكنها إجراء إعادة تعامد آني، فإن الطريقة تحل مسألة GEVP $\bar{H}v = \lambda Sv$. بالنسبة لليوفيليانات، يتم بناء مصفوفة التداخل $S$ ومصفوفة الهاملتوني $\bar{H}$ باستخدام المتجهات الذاتية اليسرى واليمنى ثنائية التناظر.
3. تفسير المرشح الطيفي: يتم صياغة الطريقة كعملية تحسين لمرشح طيفي. وخلافاً للأنظمة الهاملتونية حيث تكون الحالات الذاتية حقيقية والتطور الزمني تذبذبياً (مما يتطلب أخذ عينات نيويست صارمة)، فإن ديناميكيات ليوفيليان هي ديناميكيات مبددة. الأنماط سريعة الاضمحلال تلتف نحو الداخل في المستوى المركب، مما يسمح بأوقات تطور $\tau$ أطول دون حدوث آثار التعرج (aliasing). وهذا يتيح عزل الأنماط البطيئة (القيم الذاتية ذات الأجزاء الحقيقية الصغيرة) حتى مع بُعد كريفلو $D$ منخفض.
4. التعامل مع عدم الطبيعية (Non-Normality): يقر المؤلفون بأن المؤثرات غير الطبيعية تؤدي إلى حساسية عالية للاضطرابات (مبرهنة باوير-فيكي). تصبح مصفوفة التداخل $S$ سيئة التكييف بشكل أسي مع زيادة $D$. وللتخفيف من ذلك، تعتمد الطريقة على الطبيعة المبددة للنظام لتثبيط الأنماط السريعة، مما يسمح بقيم $\tau$ أكبر و $D$ أصغر، وبالتالي تقليل عدم الاستقرار العددي المرتبط بأبعاد كريفلو الكبيرة.
5. اختيار الحالة الابتدائية: من الضروري أن يكون للحالة الابتدائية إسقاط غير صفري على قطاع تتبع الصفر (التماسك أو عدم توازن الجمهرة) لإثارة الأنماط البطيئة. قطاع التتبع-1 (الحالة المستقرة) لا يساهم في الإشارة لتقدير الفجوة في هذا الإطار.

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

• توسيع نطاق QSD ليشمل الليوفيليانات: توفر الورقة إطاراً نظرياً لتطبيق QSD في الوقت الحقيقي على مؤثرات ليوفيليان غير هيرميتية وغير طبيعية في صيغة ليندبلاد.
• استراتيجية تقدير الفجوة: توضح كيفية تقدير فجوة ليوفيليان ($\Delta = -\max \Re(\lambda_k)$) من خلال استغلال الطبيعة المبددة للنظام لتصفية الأنماط سريعة الاضمحلال، مما يتجاوز بفعالية بعض قيود نيويست التي تنطبق على التطور الوحدوي.
• تحليل الاستقرار العددي: يسلط العمل الضوء على المقايضة بين زمن التطور $\tau$ وبُعد كريفلو $D$. ويجادل المؤلفون بأن زيادة $\tau$ (بدلاً من $D$) هي استراتيجية قابلة للتطبيق للأنظمة المفتوحة لتحسين دقة الفجوة مع الحفاظ على الاستقرار العددي، حيث تضمحل الأنماط السريعة طبيعياً.
• التطبيق على كيوبتات كات كير (Kerr Cat Qubits): تُطبق الطريقة على نموذج مبسط لكيوبت كات كير مع فقدان الفوتون الواحد (مذبذب كير فائق التوصيل مدفوع بثنائي الفوتون).

النتائج
طبق المؤلفون طريقتهم على نموذج مبسط لكيوبت كات كير مع فقدان الفوتون الواحد.

• بارامترات المحاكاة: أُجريت عمليات المحاكاة باستخدام حد أقصى قدره 30 فوتون وبُعد كريفلو $D=20$، مما نتج عنه مصفوفة كريفلو بحجم $900 \times 20$.
• إعادة بناء الفجوة: نجحت الط طريقة في إعادة بناء فجوة ليوفيليان عبر مستويات مختلفة من قوة الدفع. أظهرت النتائج أن زيادة زمن التطور $\tau$ حسنت دقة تقدير الفجوة، خاصة في المناطق التي يكون فيها مُبدل ليوفيليان (مقياس عدم الطبيعية) مرتفعاً.
• تحديد النمط البطيء: أكدت عمليات المحاكاة ظهور نمط بطيء يقابل اضمحلال التماسك وعدم توازن الجمهرة. في نظام "الضوضاء المنحازة"، تتقلص الفجوة أسياً مع سعة الحالة المتماسكة ($\propto e^{-\alpha^2}$)، مما يشير إلى تشكل فضاء حسابي محمي.
• تمثيل ويغنر: استُخدمت الحالات الذاتية المعاد بناؤها من مصفوفة كريفلو لتصور تمثيل ويغنر للنمط البطيء، مما أظهر الانتقال من نظام مخمد إلى نظام بوجود فصين متميزين (كسر التماثل) والعودة إلى نظام مضغوط عند الدفع العالي.

الأهمية والادعاءات
تدعي الورقة إمكانية استخدام QSD في الوقت الحقيقي لدراسة الأنظمة الكمومية المفتوحة، وتحديداً لتقدير فجوة ليوفيليان في سياق التحولات الطورية المبددة. يضع المؤلفون هذا العمل كخطوة ضرورية لتكييف خوارزميات عصر NISQ لتصميم ومحاكاة الأجهزة الكمومية، والتي تتضمن بطبيعتها التبديد.

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