Scattering Faddeev calculations in the double continuum

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

تتناول الورقة البحثية تحدي إجراء حسابات تشتت ميكانيكا الكم لثلاثة جسيمات في نظام المتصل المزدوج (double continuum)، حيث تكون الجسيمات الثلاثة حرة ($1+1+1$ state).

• السياق: بينما تُعد صياغة فاديف (Faddeev formalism) راسخة في الحالات المقيدة وتشتت المتصل المنفرد (جسيم واحد حر، وجسيمان مرتبطان، حالة $1+2$)، فإن المتصل المزدوج يطرح صعوبات فريدة.
• الصعوبة الجوهرية: في المتصل المزدوج، تحتوي الدالة الموجية المحسوبة على تراكب لسلوكين فيزيائيين متميزين:
  1. قنوات ثنائية الجسيمات: حيث يشكل جسيمان عنقوداً مرتبطاً والجسيم الثالث يكون حراً.
  2. قنوات التكسر ثلاثية الجسيمات: حيث تكون الجسيمات الثلاثة حرة تماماً.
• هذه القنوات ليست متعامدة بالمعنى القياسي، مما يجعل من الصعب فصل سعات التشتت (scattering amplitudes) للتشتت المرن، والتكسر، وإعادة الاتحاد من حل عددي واحد. علاوة على ذلك، يختلف السلوك التقاربي للدالة الموجية بشكل كبير بين هذين النظامين (الموجات المستوية المرتبطة مقابل الموجات الأسطوانية)، مما يعقد استخراج القيم الفيزيائية الملحوظة.

2. المنهجية

يستخدم المؤلف صياغة فاديف في فضاء التشكيل (configuration-space Faddeev formalism) باستخدام إحداثيات ج Jacobi ذات المقياس الكتلي. تنقسم المنهجية إلى عدة مكونات رئيسية:

أ. الإطار الرياضي

• الإحداثيات: يستخدم النظام متجهات ج Jacobi سداسية الأبعاد $(\mathbf{x}_i, \mathbf{y}_i)$. يتم تفكيك الدالة الموجية الكلية $\Psi$ إلى ثلاثة مكونات فاديف $\phi_i$.
• توسع الموجة الجزئية: يتم توسيع الجزء الزاوي باستخدام التوافقات الكروية ثنائية القطب، مما يحول المشكلة إلى معادلات تكاملية تفاضلية مترابطة للدوال الموجية الشعاعية $f_{i,a}(x_i, y_i)$.
• تحويل ج Jacobi: أداة رياضية مركزية هي تحويل التكامل $J_{ij}^{ab}$، الذي يربط مكونات فاديف المختلفة عن طريق التحويل بين ترتيبات إحداثيات ج Jacobi المختلفة.

ب. التحليل التقاربي

تحدد الورقة بدقة السلوك التقاربي للنظامين:

1. المتصل المنفرد ($1+2$): عندما تؤول $y \to \infty$ (مع بقاء $x$ محدوداً)، يسلك الحل سلوك "موجة مستوية مرتبطة" (حاصل ضرب دالة موجية لحالة مرتبطة ودالة ريكاتي-هانكل خارجة).
2. المتصل المزدوج ($1+1+1$): عندما تؤول كل من $x, y \to \infty$، يُوصف الحل بشكل أفضل في الإحداثيات القطبية $(\rho, \alpha)$. يتكون الشكل التقاربي من "موجات أسطوانية" (دوال هانكل من الرتبة $\nu$) معدلة بواسطة دوال ديلفز (Delves functions) (وهي الدوال الذاتية للمؤثر الزاوي).

ج. طريقة الاستخراج المبتكرة (إعادة أخذ العينات - Resampling)

الابتكار الجوهري لهذا العمل هو طريقة لاستخراج مساهمات ثنائية الجسيمات وثلاثية الجسيمات من الحل العددي:

• استراتيجية الشبكة: تُجرى الحسابات على شبكة قطبية $(\rho, \alpha)$، وهي الأمثل لوصف سلوك الموجة الأسطوانية لقناة التكسر.
• إعادة أخذ العينات (Resampling): لاستخراج مساهمة الحالة المرتبطة (ثنائية الجسيمات)، يتم إعادة أخذ عينات الدالة الموجية القطبية المحسوبة على شبكة ديكارتية $(x, y)$.
• إجراء الملاءمة (Fitting):
  • يتم إسقاط الدالة الموجية على الحالة المرتبطة للديوترون لعزل الاعتماد الشعاعي.
  • يتم إجراء ملاءمة غير خطية لفصل معاملات الموجات المستوية المنتظمة/الخارجة (قناة التشتت المرن) والموجات الأسطوانية المنتظمة/الخارجة (قناة التكسر).
  • يسمح هذا بالاستخراج المتزامن لـ عنصر مصفوفة التشتت المرن ($S_{11}$) و سعة التكسر ($S_{12}, S_{13}$) من حساب واحد.

د. التنفيذ العددي

• الأساس: تُستخدم شرائح هيرميت التكعيبية (Cubic Hermite splines) للتوسع الشعاعي.
• المحلل (Solver): يتم حل النظام الخطي الكبير الناتج (بأبعاد تقارب $256,000$) باستخدام طريقة GMRES التكرارية مع مسبق شرط (preconditioner) يعتمد على مؤثر الطاقة الحركية $(T-E)^{-1}$.
• الشروط الحدية:
  • تحت مستوى التكسر: شروط الموجة المستوية المرتبطة الخارجة.
  • فوق مستوى التكسر: شروط الموجة الأسطوانية الخارجة.

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

1. مصفوفة تشتت موحدة: يبني المؤلف مصفوفة تشتت $S$ واحدة تشمل جميع العمليات: التشتت المرن $1+2$، والتكسر $1+2 \to 1+1+1$، وإعادة الاتحاد $1+1+1 \to 1+2$، والتشتت المرن $1+1+1$.
2. تقنية تحويل الإحداثيات: تطوير تقنية إعادة أخذ عينات قوية (قطبي $\leftarrow$ ديكارتي) لفصل القنوات التقاربية غير المتعامدة.
3. الوحدوية والتبادلية المعممة: تحدد الورقة عملية ثنائية مخصصة $\langle \cdot, \cdot \rangle$ وحاصل ضرب مصفوفة $\ast$ يدمج تحويل ج Jacobi. يسمح هذا بالتحقق من الوحدوية ($S \ast S^\dagger = I$) و التبادلية ($S = S^T$) لمصفوفة تحتوي على عناصر قياسية ودوال مستمرة للزاوية القطبية $\alpha$.
4. التحقق المرجعي: تم تطبيق الطريقة على نظام تشتت النيوترون-ديوترون (n-d)، وهو معيار مرجعي في الفيزياء النووية.

4. النتائج

تم اختبار الطريقة على نظام النيوترون-ديوترون ($n-d$) في الحالة الثنائية ($J^\Pi = 1/2^+$) باستخدام جهود يوکاوا (Yukawa potentials).

• التشتت المرن ($S_{11}$): تظهر عناصر مصفوفة التشتت المرن المحسوبة توافقاً ممتازاً مع البيانات المرجعية (سواء تحت أو فوق مستوى التكسر) عبر نطاق واسع من الطاقة ($E_{lab} = 4$ إلى $42$ MeV).
• سعات التكسر ($S_{12}, S_{13}$): تتوافق سعات التكسر المستخرجة للقنوات الفردية والثلاثية مع البيانات المرجعية من Friar et al. (1995) بشكل مرضٍ جداً.
  • ملاحظة: لوحظ انحراف طفيف بالقرب من $\alpha \approx \pi/2$ بسبب الحجم المحدود للشبكة الحسابية، وهو قصور معروف في حسابات المتصل المزدوج.
• **إعادة الاتحاد ثلاثي الجسيمات و $npp$المرن:** تقدم الورقة أول النتالج المحسوبة لإعادة الاتحاد ثلاثي الجسيمات ($n+n+p \to n+d$) والتشتت المرن$n+n+p$، والتي لا توجد لها بيانات تجريبية.
• الدقة: العيوب الناتجة عن الوحدوية والتبادلية هي في حدود $10^{-3}$ إلى $10^{-2}$، وهي مقاربة للحسابات عالية الدقة السابقة، مما يؤكد الاستقرار العددي للمنهج.

5. الأهمية

• التقدم النظري: يحل هذا العمل الصعوبة طويلة الأمد المتمثلة في استخراج سعات تشتت متعددة من دالة موجية واحدة في المتصل المزدوج دون الحاجة إلى تعامد القنوات (كما في الطرق السابقة).
• الكفاءة الحسابية: من خلال استخدام الشبكات القطبية لمنطقة التكسر وإعادة أخذ العينات الديكارتية للمنطقة المرتبطة، تعمل الطريقة على تحسين وصف كلا السلوكين التقاربيين بكفاءة.
• التحقق المرجعي: يؤكد النجاح في إعادة إنتاج البيانات المرجعية لـ $n-d$ موثوقية نهج فاديف في فضاء التشكيل لمشكلات تكسر الجسيمات الثلاثة المعقدة.
• القدرة التنبؤية: القدرة على حساب إعادة الاتحاد والتشتت المرن النقي لثلاثة جسيمات توفر بيانات نظرية جديدة للتفاعلات النووية حيث تكون البيانات التجريبية شحيحة أو غير موجودة.

ختاماً، يقدم Guérout إطاراً موحداً وقوياً لتشتت الأجسام الثلاثة يتعامل بسلاسة مع الانتقال بين الحالات المرتبطة (العناقيد) والحالات التقاربية الحرة تماماً، مما يوفر أداة قوية للفيزياء النووية والذرية قليلة الأجسام.