Symplectic H2 Model Reduction for High-Dimensional Linear Quantum Systems

ملخص تقني: اختزال النموذج المتماثل (Symplectic) من نوع $H_2$ للأنظمة الكمومية الخطية عالية الأبعاد

بيان المشكلة
يتناول هذا البحث مسألة اختزال النموذج من نوع $H_2$ للأنظمة الكمومية الخطية عالية الأبعاد. ويتمثل التحدي الرئيسي في هذا المجال في قيد التحقق الفيزيائي (Physical Realizability - PR). فخلافاً للأنظمة الكلاسيكية، يجب أن تحقق الأنظمة الكمومية الخطية هويات جبرية محددة تربط مصفوفات النظام $(A, B, C, D)$ للحفاظ على علاقات التبديل القانونية (canonical commutation relations) وبنية المدخلات والمخرجات الكمومية. وغالباً ما تفشل طرق الاختزال القائمة على الإسقاط التقليدية في الحفاظ على قيود التحقق الفيزيائي غير الخطية هذه، مما يؤدي إلى نماذج مختزلة ذات رتبة أقل غير صالحة فيزيائياً. كما أن فرض هذه القيود مباشرة كشروط تحسين يؤدي إلى معضلات جبرية غير خطية صعبة وغير قابلة للتوسع للأنظمة الكبيرة.

المنهجية
يقترح المؤلفون إطار عمل "بتروف-غاليركين متماثل" (Symplectic Petrov-Galerkin framework) يضمن تحقيق شرط التحقق الفيزيائي (PR) عن طريق البناء وليس عبر الفرض اللاحق.

1. الإسقاط المتماثل (Symplectic Projection): تعتمد الطريقة على بناء نموذج مختزل الرتبة عبر الإسقاط على فضاء جزئي متماثل. إذا كان الأساس التجريبي $V \in \mathbb{R}^{2n \times 2r}$ يحقق الشرط المتماثل $V^\top J_n V = J_r$ (حيث $J_k$ هي المصفوفة المتماثلة القانونية)، وتم اختيار أساس الاختبار ليكون $U = J_r^{-1} V^\top J_n$، فإن المصفوفات المختزلة الناتجة $(A_r, B_r, C_r, D_r)$ تحقق تلقائياً هويات التحقق الفيزيائي (PR).
2. خوارزمية Q-IRKA: طور المؤلفون نسخة متماثلة من خوارزمية "إيراكا" (IRKA) - وهي خوارزمية التكرار العقلاني (Iterative Rational Krylov Algorithm). تعمل الخوارزمية كعملية تكرار لنقطة ثابتة:
   • تحديثات الإزاحة (Shift Updates): يتم تحديث نقاط الاستيفاء (الإزاحات) بشكل تكراري بناءً على أقطاب النموذج الحالي المختزل. وتحديداً، يتم ضبط الإزاحات لتكون سالبة الأقطاب المختارة (مرآة القطب - pole-mirroring)، مع احترام التماثل الطيفي للنظام الحقيقي المختزل.
   • الاستيفاء المماسي (Tangential Interpolation-): في كل تكرار، يتم توليد مجموعة مرشحة عبر حل أنظمة خطية مزاحفة $(A - \sigma_i I)^{-1} B t_{i,\ell}$ لإزاحات محددة واتجاهات مماسية.
   • استخراج الأساس المتماثل: من المجموعة المرشحة، يتم استخراج أساس متماثل باستخدام إجراء يشبه عملية "غرام-شميت" لربط المتجهات مع مرافقاتها المتماثلة ($J_n^\top v$). وتضمن خطوة تطبيع لاحقة أن الأساس يحقق بدقة الشرط $V^\top J_n V = J_r$.
   • الحفاظ على البنية (Structure Preservation): يتم الحصول على المصفوفات المختزلة حصرياً من خلال الإسقاط $A_r = UAV$، و $B_r = UB$، و $C_r = CV$، و $D_r = D$. ولا يتم إجراء أي تعديل مستقل على المصفوفات المختزلة، مما يضمن الحفاظ على التحقق الفيزيائي (PR) بدقة الآلة طوال عملية التكرار.

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

• خوارزمية Q-IRKA: تقديم خوارزمية تكرارية قابلة للتوسع وتحافظ على البنية لاختزال النموذج الأمثل من نوع $H_2$ للأنظمة الكمومية الخطية. وهي تجمع بين كفاءة خوارزمية IRKA (باستخدام حلول خطية مزاحفة) والقيود المتماثلة الصارمة.
• التحقق الفيزيائي القائم على البناء: يلغي هذا الإطار الحاجة إلى حل مسائل التحسين غير الخطية المقيدة. فالتحقق الفيزيائي هو خاصية متأصلة في هندسة الإسقاط، مما يضمن بقاء النموذج المختزل نظاماً كمومياً صالحاً.
• التحليل العددي المنهجي: يقدم البحث دراسة شاملة لكيفية اعتماد جودة الاختزال على السمات الهيكلية الفيزيائية، وتحديداً:
  • هندسة التبديد (Dissipation Geometry): التخميد المتجانس مقابل غير المتجانس.
  • توزيع القنوات (Channel Placement): تكوينات القنوات المنخفضة حيث تتشارك عدة مذبذبات في مسارات تبديد واحدة.
  • عدم تجانس النظام (System Heterogeneity): التباينات في ترددات الموقع وقوى التخميد.

النتائج
أُجريت تجارب عددية على مجموعتين مرجعيتين:

1. أنظمة بورت-هاميلتونيان ذات القنوات المنخفضة: سلاسل من المذبذبات ذات قنوات خارجية مجمعة.
2. أنظمة مستوحاة من سلسلة كيتايف البوزونية (Bosonic Kitaev-Chain): أنظمة مهيمنة على النقل ومنظمة ذات اقتران الجوار القريب.

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

• الدقة والاستقرار: تنتج خوارزمية Q-IRKA نماذج مختزلة الرتبة ذات تقريب دقيق من نوع $H_2$. وتظل التماثلات وبواقي التحقق الفيزيائي ($\Delta_{symp}$ و $\Delta_{PR}$) عند مستوى دقة الآلة ($\approx 10^{-14}$ إلى $10^{-16}$) طوال التكرارات.
• الاعتماد على عدم التجانس: تتأثر جودة الاختزال بشدة بهندسة التبديد للنظام. فالأنظمة المتجانسة ذات الاضمحلال السريع لقيم "هانكل" المنفردة (Hankel singular value) تعطي تقريبات من رتبة منخفضة بدقة عالية (خطأ $H_2 \approx 10^{-5}$). أما الأنظمة غير المتجانسة فتظهر اضمحلالاً طيفياً أبطأ، مما يؤدي إلى رتب ديناميكية فعالة أعلى وأخطاء أكبر (خطأ $H_2 \approx 10^{-3}$)، حتى مع نفس الرتبة المختزلة.
• التكلفة الحسابية: تتناسب التكلفة الحسابية بشكل أساسي مع أبعاد النظام $n$ والرتبة المختزلة $r$. أما الاعتماد على عدد القنوات $m$ فهو طفيف. وتتطلب الحالات غير المتجانسة عموماً عدداً أكبر من التكرارات للوصول إلى التقارب وتتكبد تكاليف حسابية أعلى بسبب الاضمحلال الطيفي الأبطأ.
• المقارنة: في مقارنة صغيرة النطاق، حققت Q-IRKA خطأ $H_2$ أقل (1.2130) مقارنة بطريقة التوازن شبه المتوازن (1.86) وطريقة التحسين غير الخطية (1.25).

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