Two-dimensional quantum lattice gas algorithm for anisotropic Burger-like equations

إليك ملخص تقني مفصل للورقة البحثية بعنوان "خوارزمية غاز الشبكة الكمومية ثنائية الأبعاد لمعادلات بيرجرز متباينة الخواص" لنيكولو فونيو، بيير ساغات، وجوزيبي دي مولفيتا.

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

تواجه عملية تطبيق الحوسبة الكمومية (QC) على ديناميكا الموائع الحسابية (CFD) تحديات في الموازنة بين الميزة الكمومية والدقة الفيزيائية والاستقرار العددي. تعتمد النهج الحالية غالبًا على ترميز نماذج الغاز الشبكي الكلاسيكية (مثل طريقة لابتز بولتزمان - LBM) في حالات كمومية، والتي قد تعاني من معاملات تصادم غير وحدوية (non-unitary)، مما يتطلب تنفيذات معقدة متعددة الخطوات الزمنية أو تقنيات خطية تعاني مع الاضطراب (turbulence).

على وجه التحديد، يعيد المؤلفون فحص نموذج ييبيز للغاز الشبكي الكمومي (QLG) (Phys. Rev. A, 2006). وبينما يعد هذا النموذج وحدويًا ومكتملًا من الناحية المفاهيمية، إلا أن الصيغ السابقة احتوت على تقريب في اشتقاق اللزوجة العيانية (macroscopic viscosity). علاوة على ذلك، فإن توسيع هذا النموذج إلى بُعدين (2D) مع الحفاظ على القابلية للتحليل الرياضي وتقليل استهلاك الموارد (الكيوبتات) يظل تحديًا مفتوحًا. الهدف هو اشتقاق لزوجة مصححة، وتعميم النموذج إلى بُعدين بأقل قدر من الموارد (2 كيوبت لكل موقع)، وتحديد المعادلات التفاضلية الجزئية (PDEs) العيانية الناتجة.

2. المنهجية

أ. الإطار النظري: نموذج Q-D1Q2

يحلل المؤلفون نموذج غاز شبكي كمومي أحادي البعد بـ 2 كيوبت لكل موقع (Q-D1Q2).

• تمثيل الحالة: يحمل كل موقع شبكي $x$ حالة $|\psi(x,t)\rangle = |q_0\rangle \otimes |q_1\rangle$، حيث تمثل الكيوبتات تجمعات الجسيمات التي تتحرك بسرعات منفصلة.
• الخوارزمية: يتكون التطور من ثلاث خطوات:
  1. التهيئة: يتم تهيئة الكيوبتات بناءً على دالات التوزيع الكلاسيكية $f_i$.
  2. التصادم: يتم تطبيق مؤثر وحدوي $\hat{C}$ (مُعامل بزوايا أويلر $\theta, \zeta, \xi$) للحفاظ على الكتلة.
  3. التدفق والقياس: يتم قياس التجمعات (مما يؤدي إلى انهيار الحالة) لتحديث التوزيعات، يليه تدفق كلاسيكي (streaming).
• الاشتقاق: يقوم المؤلفون بإجراء توسيع تشابمان-إنسكوغ (Chapman-Enskog expansion) على معادلة لابتز بولتزمان المشتقة من التطور الكمومي. يقومون بتوسيع معادلة الفرق المحدود حتى الدرجة الثانية في الفضاء ($\delta x$) والدرجة الأولى في الزمن ($\delta t$)، بافتراض نظام القياس الانتشاري ($\delta x \sim \epsilon, \delta t \sim \epsilon^2$).

ب. تصحيح اللزوجة

خطوة منهجية رئيسية هي الحساب الدقيق للمصطلح $J_1 - J_0$ (مشتقة حد التصادم بالنسبة للتجمعات عند التوازن).

• العمل السابق (ييبيز) قدّم تقريبًا لهذا المصطلح، مما أدى إلى تعبير عن اللزوجة $\nu \propto \cot^2(\theta)$.
• يشتق المؤلفون حلاً مغلقًا دقيقًا لتجمعات التوازن وحد التصادم، مما ينتج عنه صيغة لزوجة مصححة تتضمن حدودًا من رتب أعلى في حد منخفض لعدد ماخ (low-Mach limit).

ج. التعميم إلى بُعدين (Q-D2Q2)

تم توسيع النموذج إلى بُعدين باستخدام شبكة مربعة مع 2 كيوبت فقط لكل موقع (أقل استخدام للموارد).

• التدفق (Streaming): تتدفق التجمعات في اتجاهي $x$ و $y$.
• التصادم: يتم تطبيق نفس مؤثر التصادم أحادي البعد محليًا.
• التحليل: يشتق المؤلفون المعادلة التفاضلية الجزئية العيانية لمجموعات مختلفة من السرعات المنفصلة (مثل الشبكة المربعة، الشبكة المثلثية) لملاحظة كيفية ظهور متباينة الخواص (anisotropy) وحدود الحمل (advection).

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

1. تصحيح تحليلي للزوجة:
   يشتق المؤلفون تعبيرًا مصححًا للزوجة الكينماتيكية $\nu$ في النموذج أحادي البعد. الصيغة الجديدة هي:
   $$\nu = -\frac{\delta x^2}{\delta t} \frac{1}{2} \left( 1 - \frac{\sqrt{\alpha^2 + 1}}{\sin^2(\theta)} \right)$$
   حيث $\alpha = \cot(\theta)\cos(\zeta - \xi)$. وهذا يصحح تقريب $\cot^2(\theta)$ السابق، مما يمد صلاحية معادلة تشبه معادلة بيرجرز إلى الرتبة $O(\epsilon^4)$.

2. تعميم أدنى للبُعدين (Q-D2Q2):
   يقترحون نموذج Q-D2Q2 (بُعدان، 2 كيوبت لكل موقع). على عكس غازات الشبكة التقليدية ثنائية الأبعاد (مثل FHP) التي تتطلب 6 أو 7 سرعات، يستخدم هذا النموذج 2 فقط.

• النتيجة: يعيد النموذج إنتاج ديناميكيات تشبه معادلة بيرجرز متباينة الخواص.
• المعادلة: المعادلة العيانية المشتقة هي:
  $$\partial_t \rho + \mathbf{d} \cdot \nabla \rho - \nabla \cdot (\mathbf{D} \nabla \rho) = 0$$
  حيث $\mathbf{d}$ هو متجه حمل، و $\mathbf{D}$ هو موتر انتشار متباين الخواص يعتمد على مجموعة السرعات المختارة.

3. لزوجة صغيرة للغاية مع الاستقرار:
   تثبت الورقة أن نموذج QLG يمكنه تحقيق لزوجات صغيرة للغاية دون الاستقرار العددي الذي يوجد عادة في نماذج LBM ذات زمن الاسترخاء الأحادي (والتي تصبح غير مستقرة عندما تقترب $\tau \to 0.5$). ويرجع ذلك إلى البنية الوحدوية المتأصلة في النموذج وتحقيقه لنظرية H (H-theorem).

4. النتائج

• التحقق العددي من اللزوجة:
  أكدت المحاكاة في بُعد واحد صيغة اللزوجة المصححة. طابقت اللزوجة التجريبية التنبؤ النظري (المعادلة 22) بشكل أفضل بكثير من التقريب السابق (المعادلة 19)، خاصة في نظام اللزوجة المنخفضة.
• تكون الصدمة والاستقرار:
  • نجحت المحاكاة في تمثيل تكون الصدمة في معادلة بيرجرز.
  • المقايضة: يتطلب تحقيق لزوجة منخفضة جدًا (صدمات حادة) أوقات محاكاة أطول للترمالية (الاسترخاء) و/أو دقة شبكة أعلى، لكن المخطط يظل مستقرًا بشكل مطلق.
• السلوك متباين الخواص في بُعدين:
  • أنتجت مجموعات السرعات المختلفة سلوكيات عيانية مختلفة.
  • المجموعات المتماثلة (مثل $c_0 = -c_1$) نتج عنها ديناميكيات تشبه البعد الواحد مسقطة على بُعدين.
  • المجموعات غير المتماثلة (مثل $c_0$ أفقي، $c_1$ رأسي) أدخلت حدود المشتقات المختلطة ($\partial_{xy}\rho$) وحدود حمل إضافية، مما يؤكد الطبيعة متباينة الخواص للنموذج الأدنى ثنائي الأبعاد.
• المقارنة مع الطرق الكلاسيكية:
  أظهرت المقارنات العددية بين محاكاة QLG وطريقة الفروق المحدودة (FDM) للمعادلة التفاضلية الجزئية المستمدة توافقًا نوعيًا وكميًا ممتازًا (تحليل خطأ معيار $L_2$)، مما يثبت صحة المعادلة المشتقة.

5. الأهمية والآفاق المستقبلية

• المحاكاة الأصلية الكمومية (Quantum-Native): يرسخ هذا العمل مسارًا لمحاكاة الموائع "الأصلية الكمومية" التي لا تكتفي بمجرد ترميز الخوارزميات الكلاسيكية، بل تستخدم الديناميكيات الوحدوية الكمومية لدمج الظواهر التبددية (اللزوجة) بشكل طبيعي ضمن قواعد مجهرية عكسية.
• كفاءة الموارد: من خلال تحقيق ديناميكيات متباينة الخواص في بُعدين باستخدام 2 كيوبت فقط لكل موقع، يقدم النموذج بديلًا فعالًا من حيث الموارد لغازات الشبكة الكلاسيكية (مثل FHP) ونماذج LBM القياسية، والتي تتطلب حالات أكثر لالتقاط فيزياء مماثلة.
• الجسر إلى معادلات نافيير-ستوكس: بينما يلتقط النموذج الحالي ديناميكيات تشبه معادلة بيرجرز (حفظ قياسي)، يحدد المؤلفون أنه بإضافة سرعة ثالثة (3 كيوبت)، يمكن تمكين حفظ الزخم، مما قد يؤدي إلى حل كامل لمعادلات نافيير-ستوكس على شبكة كمومية.
• الصلة بالأجهزة (Hardware): إن الحد الأدنى من متطلبات الكيوبت والاستقرار المطلق يجعل هذا النهج واعدًا للغاية للتنفيذ على الأجهزة الكمومية في عصر NISQ (الأجهزة الكمومية ذات الضجيج المتوسط).

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