The diffusion equation for non-Markovian Gaussian stochastic processes

ملخص تقني: معادلة الانتشار للعمليات العشوائية الغاوسية غير الماركوفية

بيان المشكلة
تتناول الورقة اشتقاق معادلة تطور دقيقة لدالة كثافة الاحتمال (PDF) لإزاحات الجسيمات، $P(x, t)$، الناتجة عن عمليات سرعة غاوسية تعسفية. وبينما يوفر نموذج فوكر–بلانك (FPE) الكلاسيكي ووصف لانجفان إطاراً قوياً للعمليات الماركوفية والمستقرة (وتحديداً عملية أورن-أولينبك)، إلا أنهما يفشلان في رصد الديناميكيات الدقيقة للأنظمة غير الماركوفية حيث تكون ارتباطات السرعة محدودة وغير أسية.

تعتمد المقاربات الحالية للعمليات الغاوسية غير الماركوفية، مثل معادلة زوانزيغ–باليسكو (ZBE) ومعادلة باتشيلور–هانجي (BHE)، على تقريبات. ويشير المؤلفون إلى أنه على الرغم من أن هذه المعادلات تنطبق على العمليات الغاوسية المستقرة، إلا أنها ليست دقيقة؛ وتحديداً، فإنها تفشل في إعادة إنتاج العزوات الموضعية (position moments) الغاوسية لما وراء الرتبة الثانية. علاوة على ذلك، تعتمد معادلة BHE صراحة على الزمن الأولي $t_0$ وقد تنتهك شرط الإيجابية. وتتمثل المشكلة المركزية في اشتقاق معادلة انتشار مغلقة ودقيقة تحافظ على الطبيعة الغاوسية للعملية دون افتراض الاستقرار أو الماركوفية.

المنهجية
يتبنى المؤلفون نهجاً منهجياً يعتمد على الدالة المميزة لكثافة الموضع، $\hat{P}(q, t)$، بدلاً من البدء من معادلة لانجف مثل محددة. وتتقدم المنهجية كما يلي:

1. توسيع الدالة المميزة: يتم توسيع الدالة المميزة في متسلسلة تايلور حول $q=0$، مما يربطها بالعزوات الزوجية للإزاحة، $\langle x^{2n}(t) \rangle$.
2. تطبيق مبرهنة ويك (Wick's Theorem): بما أن عملية السرعة $v(t)$ مفترضة بأنها غاوسية، فإن جميع الخصائص الإحصائية تتحدد من خلال دالة الارتباط بين زمنين، $\langle v(t_1)v(t_2) \rangle$. يستخدم المؤلفون مبرهنة ويك لتفكيك ارتباطات السرعة ذات الرتب العليا إلى مجموع نواتج ارتباطات ثنائية النقاط.
3. التسلسل الهرمي المبتور: يتم تعريف دالة مميزة مبتورة العزوات، $\hat{\rho}_N(q, t)$. ويتم التعبير عن التطور الزمني لهذه الدالة المبتورة كمجموع لجميع احتمالات اقترانات ويك (Wick pairings) لمتغيرات السرعة.
4. إعادة الترتيب المخططاتي (Diagrammatic Rearrangement): يقدم المؤلفون تمثيلاً مخططاتياً لاقترانات ويك. ومن خلال إعادة تجميع هذه المخططات، يميزون بين مخططات ويك "المتصلة هندسياً" والمخططات غير المتصلة. تحاكي هذه الخطوة مبرهنة العناقيد المرتبطة (linked-cluster theorem) في نظرية الحقل الكمي، حيث تساهم المخططات المتصلة فقط في لوغاريتم الدالة المولدة.
5. البنية التكرارية: تؤدي إعادة الترتيب إلى علاقة تكرارية هرمية تتضمن معاملات $R_k$، والتي تمثل مجموع جميع مخططات ويك المتصلة هندسياً من الرتبة $2k$.

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

1. معادلة التطور الدقيقة: النتيجة الأساسية هي اشتقاق معادلة انتشار غير ماركوفية دقيقة ومغلقة لـ $P(x, t)$:
   $$\frac{\partial}{\partial t}P(x, t) = \sum_{k=1}^{\infty} \int_0^t dt_1 \int_0^{t_1} dt_2 \cdots \int_0^{t_{2k-2}} dt_{2k-1} R_k \frac{\partial^{2k}}{\partial x^{2k}} P(x, t_{2k-1})$$
   حيث $R_k$ هي نوى (kernels) مشتقة من دالات ارتباط السرعة المتصلة.

2. تعميم النماذج القائمة:

   • المصطلح الأول من المتسلسلة ($k=1$) يستعيد معادلة زوانزيغ–باليسكو (ZBE).
   • يوضح المؤلفون أن أي بتر نهائي لهذه المتسلسلة يدمر الطابع الغاوسي للحل. فقط المتسلسلة اللانهائية الكاملة تحافظ بدقة على الغاوسية.
   • المعادلة صالحة لكل من الحركيات المستقرة وغير المستقرة، بشر قدم أن كثافة احتمال السرعة عند زمن واحد هي غاوسية.

3. التحليل التوليفي (Combinatorial Analysis): تثبت الورقة وجود علاقة تكرارية لعدد مخططات ويك المتصلة ($\bar{R}_k$)، حيث أن إجمالي اقترانات ويك $(2n-1)!!$ هو مجموع المساهمات من المخططات المتصلة من الرتب الأدنى. تشبه هذه البنية علاقات التكرار للمخططات الفاينمانية المتصلة في الكهروديناميكا الكمية (QED).

4. التطبيقات على عمليات محددة:

   • عملية أورن–أولينبك (OU): في حد التخميد المفرط ($\gamma \to \infty$)، تقل المعادلة بشكل صحيح إلى معادلة انتشار سموكولوفسكي القياسية، مع تلاشي الحدود ذات الرتب الأعلى.
   • الحركة البرونية الكسرية (fBm): بالنسبة لـ fBm مع أس هيرست $1/2 < H < 1$، يعطي المصطلح الرئيسي للمعادلة المشتقة معادلة تفاضلية كسرية تتضمن مشتقة ريمان–ليفيلي (Riemann–Liouville) الكسرية. يشير المؤلفون إلى أن هذا الشكل المحدد لم يُدرس من قبل في الأدبيات.

الأهمية والادعاءات
تدعي الورقة أنها تقدم أول معادلة انتشار دقيقة لكثافة الاحتمال للإزاحات الناتجة عن عمليات سرعة غاوسية غير ماركوفية تعسفية. تكمن الأهمية في:

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

يصرح المؤلفون صراحةً أنه بالنسبة لنظام التنافر (antipersistent regime) للحركة البرونية الكسرية ($0 < H < 1/2$)، فإن اشتقاق معادلة الرتبة الأولى يظل مشكلة مفتوحة بسبب عدم تكامل نواة ارتباط السرعة عند الأوقات المتطابقة، وهو ما يعتزمون معالجته في أعمال مستقبلية.