Quantum algorithms for viscosity solutions to nonlinear Hamilton-Jacobi equations based on an entropy penalisation method

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

تتناول الورقة البحثية تحدي محاكاة المعادلات التفاضلية الجزئية (PDEs) غير الخطية على الحواسيب الكمومية، مع التركيز بشكل خاص على معادلات هاملتون-جاكوبي (HJ).

• الصعوبة الجوهرية: تقتصر معظم الخوارزميات الكمومية الموجودة على الأنظمة الخطية. أما المعادلات التفاضلية الجزئية غير الخطية، وخاصة معادلات هاملتون-جاكوبي، فتطرح عقبتين رئيسيتين:
  1. اللاخطية والتفردات: غالبًا ما تطور معادلات هاملتون-جاكوبي تفردات في زمن محدد (caustics) حيث تنهار الحلول الكلاسيكية. وتكون الحلول ذات الصلة الفيزيائية في هذه الأنظمة هي "حلول اللزوجة" (viscosity solutions)، وهي حلول ضعيفة تُعرف لما بعد مرحلة التفرد.
  2. لعنة الأبعاد: تعاني الطرق العددية الكلاسيكية لحلول اللزوجة (التي غالبًا ما تتضمن مخططات عالية الرتبة مع لزوجة عددية أو محددات ميل) من لعنة الأبعاد، مما يجعلها غير قابلة للتنفيذ في المسائل عالية الأبعاد الشائعة في التحكم الأمثل، والألعاب ذات المجال المتوسط، والتعلم الآلي.
• الفجوة: بينما توجد بعض الخوارزميات الكمومية للتمثيلات الخطية للمعادلات التفاضلية الجزئية غير الخطية (مثل طرق مجموعة المستويات)، إلا أنها غالبًا ما تفشل في التقاط الميزات غير الخطية الجوهرية مثل الصدمات أو الـ caustics، أو تكون محدودة باللاخطية الضعيفة. هناك نقص في الخوارزميات الكمومية الفعالة التي يمكنها التعامل مع اللاخطية القوية والتبدد الضعيف (الذي يميز حلول اللزوجة) مع البقاء صالحة لأزمنة طويلة للغاية.

2. المنهجية

يقترح المؤلفون إطار عمل موحد يحول معادلة هاملتون-جاكوبي اللزجة غير الخطية إلى معادلة تفاضلية جزئية معلمية (parabolic PDE) خطية، والتي يمكن معالجتها بكفاءة عبر المحاكاة الكمومية. تتكون المنهجية من ثلاث مراحل:

أ. الصيغة الخطية عبر جزاء الإنتروبيا (Entropy Penalisation)

بدلاً من استخدام اللزوجة العددية القياسية (وهي عملية غير خطية ويصعب محاكاتها كموميًا)، يستخدم المؤلفون نهج اللزوجة المعنوية الاصطناعية مقترنًا بـ طريقة جزاء الإنتروبيا (Gomes and Valdinoci).

1. معادلة HJ اللزجة: يبدأون بالمعادلة المنتظمة:
   $$\frac{\partial S_\nu}{\partial t} + H(t, x, \nabla S_\nu) = 2a\nu \Delta S_\nu$$
   حيث $\nu$ هو معامل لزوجة صغير.
2. تعميم كول-هوف (Cole-Hopf Generalization):
   • بالنسبة لـ هاملتونيان التربيعي، يقوم تحويل كول-هوف الكلاسيكي ($S_\nu = -2\nu \ln u$) بنقل المعادلة غير الخطية بدقة إلى معادلة حرارية خطية.
   • بالنسبة لـ هاملتونيان محدب عام، يعمم المؤلفون ذلك باستخدام طريقة جزاء الإنتروبيا. يقدمون مؤثر تطور في زمن منفصل يقرب الديناميكيات غير الخطية. ومن خلال عملية زمنية منفصلة تتضمن لاغرانجيان $L(x,v)$، يبنون مؤثرًا خطيًا $\mathcal{L}$ بحيث يتطور الحل $u(t,x)$ خطيًا:
     $$u(t, x) \approx \text{solution to a linear parabolic PDE}$$
   • يتم الحفاظ على العلاقة عبر التحويل: $S_\nu(t, x) = -2\nu \ln u(t, x)$.
3. الحد الزمني المستمر: يُظهر المخطط الخطي ذو الزمن المنفصل تقاربه نحو معادلة تفاضلية جزئية معلمية خطية في الزمن المستمر على شكل:
   $$\frac{\partial u}{\partial t} = \sum \mu_i \frac{\partial u}{\partial x_i} + \sum \nu_{jk} \frac{\partial^2 u}{\partial x_j \partial x_k}$$
   يمكن محاكاة هذه المعادلة الخطية بكفاءة على الأجهزة الكمومية.

ب. المحاكاة الكمومية (Schrödingerisation)

لمحاكاة المعادلة التفاضلية الجزئية المعلمية الخطية الناتجة، يستخدم المؤلفون بروتوكول الشرويدنجرية (Schrödingerisation).

• الآلية: تقوم هذه الطريقة بدمج حل المعادلة المعلمية الخطية غير الوحدوية (non-unitary) داخل تطور كمومي وحدوي (unitary) عن طريق إدخال بُعد مساعد (ancilla).
• التنفيذ:
  • النموذج التناظري (Analog): يستخدم أنظمة كمومية ذات متغيرات مستمرة (qumodes) حيث المتغير المكاني $x$ مستمر. يتم رسم المعادلة التفريقية الجزئية كمعادلة تشبه معادلة شرويدنجر $du/dt = -iA(t)u$.
  • النموذج الرقمي (Digital): يتم تقسيم الفضاء والزمن، مما يحول المسألة إلى نظام قائم على الكيوبت (qubit). تستخدم الخوارزمية تقنيات ترميز الكتلة (block-encoding) أو محاكاة هاملتونيان الوصول المتفرق (sparse-access).
• التطبيع (Normalization): الميزة الرئيسية لبروتوكول الشرويدنجرية المستخدم هنا هي قدرته على تقدير ثابت التطبيع $\|u(t)\|$ كمنتج ثانوي للمحاكاة، وهو أمر بالغ الأهمية لاستعادة الحل الفيزيائي $S_\nu$.

ج. استخراج الكميات الفيزيائية

بما أن التصوير المقطعي الكامل للحالة (full state tomography) غير فعال، طورت الورقة بروتوكولات كمومية محددة لاستخراج الملاحظات الفيزيائية ذات الصلة مباشرة من الحالة الكمومية $|u(t)\rangle$ دون إعادة بناء الدالة الموجية بالكامل:

1. القيمة عند نقطة $S(t, x_a)$: تُقدر عبر سعة الاحتمال $|\langle x_a | u(t) \rangle|^2$ وثابت التطبيع.
2. التدرج $\nabla S(t, x_a)$: يُقدر باستخدام تقنيات القياس الضعيف (weak measurement). يتفاعل حالة مساعدة مع النظام عبر بوابة وحدوية متحكم بها، ويكشف انزياح تربيعات الحالة المساعدة (أو حالة الكيوبت) عن التدرج.
3. الحد الأدنى العالمي $S_{min}$: يُقدر باستخدام ثابت التطبيع $\|u(t)\|^2$. في حد $\nu \to 0$، يقترب $-\nu \ln \|u(t)\|^2$ من القيمة الدنيا للحل (مرتبط بدالة LogSumExp).
4. قيمة الدالة عند الحد الأدنى $f(x^*)$: تُقدر عبر قياس القيمة المتوقعة $\langle u(t) | f(\hat{x}) | u(t) \rangle$. نظرًا لتركيز كثافة الاحتمال حول الحد الأدنى مع $\nu \to 0$، فإن هذه القيمة المتوقعة تتقارب نحو $f(x^*)$.

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

• الصيغة الخطية المعممة: توسع الورقة تحويل كول-هوف ليشمل هاملتونيان محدب عام (وليس فقط التربيعي) باستخدام طريقة جزاء الإنتروبيا. يسمح هذا بتمثيل خطي لمعادلات هاملتون-جاكوبي اللزجة التي تلتقط الميزات غير الخطية الجوهرية (مثل الصدمات) عالميًا عبر الزمن.
• خوارزميات كمومية لحلول اللزوجة: تقدم الورقة أول مجموعة شاملة من الخوارزميات الكمومية (التناظرية والرقمية) المصممة خصيصًا لاستخراج حلول اللزوجة ومشتقاتها، متجاوزة حاجز اللاخطية الذي يعيق معظم خوارכות المعادلات التفاضلية الجزئية الكمومية الأخرى.
• استخراج الملاحظات بكفاءة: صمم المؤلفون بروتوكولات محددة وفعالة لاستخراج التدرجات، والحدود الدنيا، وقيم الدوال عند الحدود الدنيا دون الحاجة إلى التصوير المقطعي الكامل للحالة. تعتمد هذه البروتوكولات على قياس ثوابت التطبيع والقياسات الضعيفة، مما يوفر تسريعًا حدوديًا (polynomial speedups) مقارنة بالطرق الكلاسيكية في الأبعاد العالية.
• تحليل الخطأ: تقدم الورقة حدود خطأ صارمة، حيث تُظهر أن خطأ التقريب يعتمد على معامل اللزوجة $\nu$ وخطوة التجزئة $h$. وتثبت أنه للوصقة المستهدفة $\epsilon$، فإن الموارد المطلوبة تتوسع حدوديًا مع البعد $d$ (لتجنب لعنة الأبعاد) تحت شروط محددة على $\nu$ و $h$.

4. النتائج

• الصلاحية النظرية: تم إثبات صحة الإطار لأي لاخطية مقابلة لهاملتونيان محدب ولأزمنة طويلة للغاية.
• التعقيد:
  • تحقق المحاكاة الكمومية للمعادلة المعلمية الخطية تعقيد استعلام قدره $\tilde{O}(\alpha_A t \ln(1/\epsilon) / \|u(t)\|)$، حيث تعتمد $\alpha_A$ على البعد وخصائص الهاملتونيان.
  • لتحقيق دقة $\epsilon$ للحل $S$، يجب اختيار معامل اللزوجة بحيث يكون $\nu = O((\epsilon/d)^2)$.
  • لتقدير التدرج، يجب أن يتوسع حجم الشبكة كـ $\Delta x = O(\epsilon/d^{3/2})$.
• التطبيق على معادلة بورغرز (Burgers' Equation): أظهرت الدراسة أن الطريقة تنطبق مباشرة على معادلة بورغرز القسرية (تدرج معادلة هاملتون-جاكوبي)، وهي نموذج أساسي في ديناميكا السوائل.
• الاستقرار العددي: يتجنب استخدام اللزوجة المعنوية الخطية التذبذبات العددية المرتبطة بمحددات الميل غير الخطية، مما يضمن التقارب نحو حل اللزوجة الفريد.

5. الأهمية

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