On the role of positivity preservation for high order approximations of the Dean--Kawasaki equation

ملخص تقني: الحفاظ على الإيجابية في التقريبات عالية الرتبة لمعادلة دين-كاوساكي (Dean–Kawaki)

بيان المشكلة
تصف معادلة دين-كاوساكي التطور للقياس التجريبي لـ $N$ من الجسيمات البراونية المستقلة. وبينما تعد هذه المعادلة دقيقة تماماً على المستوى المجهري كحل مارينجال (martingale solution) لقياس عشوائي ذري، إلا أنها شديدة التفرد بحيث لا يمكن تفسيرها كمعادلة تفاضلية جزئية عشوائية (SPDE) ذات كثافة، وذلك بسبب معامل ضجيج الجذر التربيعي غير لبيسي (non-Lipschitz) والنمو الفائق للتحجيم (scaling supercriticality). ولتقريب نظام الجسيمات هذا إلى نموذج مستمر، يلزم إجراء عملية تنظيم (regularization).

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

1. مناهج الفروق المحدودة (finite-difference approaches) (مثل [CF23, CFIR26]): تحقق تقارباً ضعيفاً عالي الرتبة في أنظمة الكثافة العالية، لكنها تفشل في الحفاظ على الإيجابية بسبب الطبيعة غير الموضعية لمؤثر التباعد (divergence operator)، مما يسمح للحلول بأخذ قيم سالبة.
2. المناهج الحافظة للإيجابية (positivity-preserving approaches) (مثل [DKP24, DJP26]): تقوم بقطع الضجيج أو تنظيم التفرد لضمان عدم السالبية، لكنها تعاني من معدلات تقارب أضعف أو تتطلب قطعاً قوياً للضجيج مما يحد من جودة التقريب.

يبحث هذا البحث في التأثير المحدد لـ انتهاك الحفاظ على الإيجابية على خطأ التقريب الضعيف، وتحديداً في أنظمة الكثافة المنخفضة حيث قد تتلاشى الكثافة الأولية محلياً. ويتساءل البحث عما إذا كانت التقريبات التي لا تحافظ على الإيجابية يمكن أن تحقق معدلات تقارب فوق-متعددة الحدود (superpolynomial convergence rates)، وكيف يتصرف الخطأ عندما لا تكون الكثافة الأولية بعيدة عن الصفر.

المنهجية
يحلل المؤلفون تنظيماً طيفياً (spectral regularization) لمعادلة دين-كاوساكي:
$$du_t = \frac{1}{2}\Delta u_t dt + \frac{1}{\sqrt{N}} \nabla \cdot \left( K_\varepsilon * (f(u_t^+) dW_t) \right)$$
حيث $K_\varepsilon$ هو ممهد (mollifier) ذو تحويل فورييه مدعوم بشكل مضغوط، و $u_t^+ = \max(u_t, 0)$، و $f$ يحاكي ضجيج الجذر التربيعي (على سبيل المثال $f(x)=\sqrt{x}$). وخلافاً للمخططات السابقة الحافظة للإيجبية، لا يتم قطع عملية ويينر (Wiener process) $W$، ويعمل التنظيم الطيفي على معامل الضجيج بأكمله. هذا الخيار يسمح للحل $u_t$ بأخذ قيم سالبة، مما يكسر الحفاظ على الإيجابية، ولكنه يوفر ميزة نقل التنظيم إلى دالة الاختبار عبر الالتفاف (convolution)، مستغلاً انتظام دالة الاختبار.

يعتمد التحليل على:

1. الثنائية مع معادلة هاميلتون-جاكوبي-بيلمان (HJB): توسيع نتيجة الثنائية لـ [KLvR19] لتشمل دوال اختبار معقدة (complex-valued) لربط دالة المولد للعزوم لحل الـ SPDE بمعادلة كولموغوروف الخلفية ذات الأبعاد اللانهائية.
2. تحليل تحويل لابلاس: قياس الخطأ الضعيف عبر الفرق في تحويلات لابلاس، $|E[e^{\langle u_t, \phi \rangle}] - E[e^{\langle \mu_t, \phi \rangle}]|$، والذي يتحكم في العزوم من جميع الرتب عبر صيغة كوشي التكاملية.
3. تفكيك الخطأ: يتم تفكيك الخطأ الكلي إلى خطأ التقريب الأولي، وخطأ القطع الفورييه، وحد مأخوذ من الخطأ الناجم تحديداً عن الجزء السالب من الحل ($u_t^-$).
4. بناء الحد الأدنى: تم إنشاء مثال محدد أحادي البعد حيث تتلاشى الكثافة الأولية في فترات فرعية لاستخلاص حد أدنى لخطأ التقريب، مما يوضح قيود المخططات التي لا تحافظ على الإيجابية في أنظمة الكثافة المنخفضة.

المساهمات والنتائج الرئيسية

1. حدود الخطأ الضعيف الموحدة عبر الزمن:
   يضع البحث حداً موحداً للخطأ الضعيف عبر الزمن مقاساً عبر تحويل لابلاس. يتكون حد الخطأ من ثلاثة مكونات:

• خطأ التقريب الأولي.
• خطأ القطع الفورييه (يعتمد على الدقة الطيفية $\varepsilon$).
• خطأ الجزء السالب: حد يتناسب مع احتمال وحجم أخذ الحل لقيم سالبة.

2. نظام الكثافة العالية ($\rho_{\min} > 0$):
   عندما تكون الكثافة الأولية بعيدة عن الصفر ($\rho_0 \ge \rho_{\min} > 0$)، يتم كبت الجزء السالب من الحل أسياً. وبناءً على ذلك، يحقق التقريب الطيفي معدل تقارب فوق-متعدد الحدود للدوال الناعمة ($\phi \in C^\infty$). وتحديداً، مع اختيار أمثل لـ $\varepsilon \sim N^{-1/d}(\log N)^{3/d}$، يتلاشى الخطأ بمعدل $O(N^{-1-\alpha/d}(\log N)^{3\alpha/d})$ لأي $\alpha > 0$. وتتحقق هذه النتيجة رغم غياب الحفاظ الصارم على الإيجابية، بشرط أن تكون الكثافة عالية بما يكفي بالنسبة لقوة الضجيج.

3. نظام الكثافة المنخفضة ($\rho_{\min} = 0$):
   عندما يُسمح للكثافة الأولية بالتلاشي، يختفي الكبت الأسي للجزء السالب. ويتدهور الحد الأعلى إلى معدل متعدد الحدود: $O(N^{-(3\alpha+d)/(2\alpha+d)})$.
   ومن الأهمية بمكان أن المؤلفين قدموا حداً أدنى لمثال أحادي البعد حيث $\rho_{\min}=0$. لقد أثبتوا أن الخطأ الضعي لا يمكن أن يتلاشى فوق-متعدد الحدود في هذا النظام. يصبح الجزء السالب من الحل كبيراً نسبياً، مما يهيمن على التقلبات ويحد من جودة التقريب. وهذا يشير إلى أن التقريبات التي لا تحافظ على الإيجبية في مناطق الكثافة المنخفضة محدودة بطبيعتها من حيث الدقة.

4. المقارنة مع المناهج الموجودة:

• في الأبعاد $d=1, 2$، يكون الحد الأعلى للمخطط غير الحافظ للإيجبية في نظام الكثافة المنخفضة أسوأ من نظام تقريب الـ SPDE الحافظ للإيجبية في [DKP24].
• تشير النتائج إلى أن المخططات الهجينة (التي تجمع بين التقريبات الطيفية عالية الرتبة في مناطق الكثافة العالية مع طرق تعتمد على الجسيمات أو الحفاظ على الإيجبية في مناطق الكثافة المنخفضة) قد تكون الاستراتيجية الأكثر كفاءة.

5. التحقق العددي:
   تؤكد التجارب العددية النتائج النظرية. تظهر المحاكاة أن:

• في أنظمة الكثافة العالية، يظهر تقريب غاليركن الطيفي (spectral Galerkin) تقارباً فوق-متعدد الحدود للملاحظات الناعمة.
• في أنظمة الكثافة المنخفضة، تزداد احتمالية اتخاذ الحل لقيم سالبة بشكل كبير، مما يربطها بالحدود النظرية الدنيا للخطأ.
• فقدان الإيجابية يكون نادراً عندما تكون $\rho_{\min}$ كبيرة، ولكنه متكرر عندما تكون $\rho_{\min} \approx 0$.

الأهمية والادعاءات
يزعم البحث أنه يقدم أول تحليل للتقريبات غير الحافظة للإيجبية لمعادلة دين-كاوساكي التي تسمح بكثافات أولية منخفضة. وتكمن أهميته الأساسية في تحديد تكلفة انتهاك الحفاظ على الإيجبية:

• يثبت أن الحفاظ على الإيجبية ضروري لتحقيق معدلات تقارب فوق-متعدد الحدود عندما تتلاشى الكثافة الأولية.
• يوضح أنه بينما يمكن للتنظيم الطيفي تحقيق دقة عالية الرتبة في أنظمة الكثافة العالية دون حفظ صارم للإيجبية، فإن "الانحرافات السالبة" للحل تصبح مصدراً مهيمناً للخطأ في أنظمة الكثافة المنخفضة، مما يمنع التقارب فوق-متعدد الحدود.
• يميز العمل نفسه عن المناهج ذات الأبعاد المحدودة (مثل [CF23]) من خلال استخدام البنية الماركوفية ذات الأبعاد اللانهائية والثنائية مع معادلة HJB، مما يسمح بحدود خطأ موحدة عبر الزمن ومعالجة دوال اختبار ذات قيم مركبة.

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