Kinetic approach of light-nuclei production in intermediate-energy heavy-ion collisions

خاتمة تقنية: النهج الحركي لإنتاج النوى الخفيفة في تصادمات الأيونات الثقيلة ذات الطاقة المتوسطة

بيان المشكلة
تعد تصادمات الأيونات الثقيلة ذات الطاقة المتوسطة (من طاقة فيرمي حتى نطاق جيجا إلكترون فولت) بالغة الأهمية لاستقصاء التفاعلات الفعالة بين النوكليونات ومعادلة الحالة النووية. وبينما تم تحليل ملاحظات النوكليونات والبيونات (مثل التدفق الجمعي، ونسب الأطياف) بشكل مكثف، إلا أن النوى الخفيفة ($A \le 4$)، التي تُنتج بكثرة في هذه التصادمات، لم تحظَ تاريخياً بالاهتمام النظري الكافي. إن معظم مقاربات النظرية الحركية القياسية، مثل تلك القائمة على معادلة بولتزمان-أوهلينج-أولينج (BUU)، تقتصر تفاعلاتها عند التصادمات ثنائية الجسم وتعتبر فقط درجات الحرية النووية. وبناءً على ذلك، فإنها تفشل في الوصف الديناميكي لإنتاج وتفكك النوى الخفيفة (الديوترون، التريتون، الهيليوم-3، وجسيمات ألفا) وكذلك تأثير موت (Mott effect) — وهي الظاهرة التي تذوب فيها النوى الخفيفة في وسط نووي كثيف عندما يصبح كثافة فضاء الطور للنوكليونات المحيطة مرتفعة جداً.

تظهر البيانات التجريبية الجديدة من تعاوني INDRA وFOPI، لا سيما في تصادمات Au+Au المركزية، زيادة كبيرة في عوائد جسيمات ألفا عند طاقات الحزمة المنخفضة. تشير هذه الملاحظة إلى تأثير "موت" أضعف لجسيمات ألفا مقارنة بالعناقيد الأخف، ويرجع ذلك على الأرجح إلى طاقة ربطها الأكبر. ومع ذلك، تفتقر النماذج النظرية الحالية إلى إطار ديناميكي موحد يتضمن جميع النوى الخفيفة حتى $A=4$ لتفسير هذه العوائد وتأثيرات الوسط الكامنة.

المنهجية
طوّر المؤلفون نهجاً حركياً يعتمد على صيغة دالة غرين للمجموعات الكبيرة في الزمن الحقيقي. يعامل هذا الإطار النوى الخفيفة كدرجات حرية ديناميكية جنباً إلى جنب مع النوكليونات، البيونات، ورنينات دلتا ($\Delta$).

1. المعادلات الحركية: يتم تحديد التطور الزمني لدوال ويغنر (توزيعات فضاء الطور) $f_\tau(\vec{r}, \vec{p}, t)$ لأنواع الجسيمات $\tau$ (بما في ذلك $n, p, d, t, h, \alpha, \pi, \Delta$) عبر معادلات حركية مقترنة:
   $$(\partial_t + \vec{\nabla}_p \epsilon_\tau \cdot \vec{\nabla}_r - \vec{\nabla}_r \epsilon_\tau \cdot \vec{\nabla}_p)f_\tau = I^{coll}_\tau$$
   حيث $\epsilon_\tau$ هي طاقة الجسيم الواحد المستمدة من جهد سكيرمي (Skyrme pseudopotential)، و $I^{coll}_\tau$ يمثل تكامل التصادم.

2. تكاملات التصادم والتفاعلات: يتضمن تكامل التصادم حدود الربح والفقد المستمدة من توسعات دالة غرين للمجموعات الكبيرة. يدمج النموذج صراحةً التفاعلات التحفيزية المستحثة بالنوكليونات لتكوين وتحلل النوى الخفيفة:

   • $NNN \leftrightarrow Nd$
   • $NNNN \leftrightarrow Nt(h)$
   • $NNNNN \leftrightarrow N\alpha$
   • $NNt(h) \leftrightarrow N\alpha$
   • $N\alpha \leftrightarrow dt(h)$
     يتم حساب سعات الانتقال باستخدام تقريب الدفع (impulse approximation)، حيث يتم تفكيك مربع القيمة المطلقة لعنصر المصفوفة إلى حاصل ضرب دالة الموجة في فضاء الزخم الداخلي للنواة الخفيفة (المُنمذجة كدالة غاوسية) وسعة التشتت المرن بين النوكليونات. يتم إدخال عامل يعتمد على الطاقة $F(\sqrt{s})$ لتصحيح إهمال التشتت المرن وأوجه القصور في التقريب؛ ويتم ضبطه لإعادة إنتاج مقاطع عرض اضمحلال النوكليون-النواة التجريبية.

3. تنفيذ تأثير موت (Mott Effect): يتم تنفيذ تأثير موت عبر قطع في فضاء الطور (phase-space cutoff) في حد الإنتاج في تكامل التصادم. يمكن للنواة الخفيفة ذات العدد الكتلي $A$ أن تتكون فقط إذا كانت متوسط كثافة فضاء الطور للنوكليونات $\langle f_N \rangle_A$ في الوسط تحت عتبة حرجة $f^{cut}_A$:
   $$\langle f_N \rangle_A \equiv \int f_N\left(\frac{\vec{P}}{A} + \vec{p}\right)\rho_A(\vec{p})d\vec{p} \le f^{cut}_A$$
   هنا، تعني قيمة $f^{cut}_A$ الأصغر تأثيراً أقوى لـ "موت" (أي سهولة الذوبان).

4. الحل العددي: يتم حل المعادلات باستخدام نهج الجسيمات الاختبارية (test-particle approach) مقترناً بـ طريقة الهاميلتوني الشبكي (lattice Hamiltonian method) لدرجات الانجراف وطريقة عشوائية لتكامل التصادم.

النتائج الرئيسية
طُبِّق النموذج على تصادمات Au+Au المركزية عند طاقات حزمة تتراوح من $0.25A$ إلى $1.0A$ جيجا إلكترون فولت.

• إعادة إنتاج العوائد: باستخدام معاملات العتبة المُحسّنة ($f^{cut}_{A=2} = 0.11$, $f^{cut}_{A=3} = 0.16$, و $f^{cut}_{A=4} = 0.25$)، نجح النهج الحركي في إعادة إنتاج عوائد النوى الخفيفة المقاسة من تعاون FOPI عبر نطاق الطاقة الكامل.
• التطور الزمني: تُظهر المحاكاة أن النوى الخفية تُنتج بغزارة خلال مرحلة الانضغاط المبكرة للتصادم بسبب الكثافة النووية العالية. وهذا يؤكد ضرورة معاملتها ديناميكياً طوال عملية التصادم بأكملها، بدلاً من الاكتفاء بمعاملتها عند مرحلة التجمد الحركي (كما في نماذج الاندماج/الالتحام - coalescence models).
• زيادة عوائد جسيمات ألفا: ينجح النموذج في تفسير الزيادة الملحوظة في عوائد جسيمات ألفا عند الطاقات المنخفضة، حيث تتجاوز عوائد ألفا عوائد الهيليوم-3. ويُعزى ذلك إلى تأثير "موت" الأضعف على جسيم ألفا. فنظراً لطاقة ربطها الأكبر بكثير، تتطلب جسيمات ألفا كثافة فضاء طور أعلى لتذوب (أي زخم "موت" أصغر) مقارنة بالديوترونات أو التريتونات أو الهيليوم-3. وبناءً على ذلك، يلزم استخدام معامل عتبة أكبر ($f^{cut}_{A=4} = 0.25$) لمطابقة البيانات، بينما تؤدي القيمة $f^{cut}_{A=4} = 0.15$ إلى نقص كبير في تقدير عائد ألفا.

الأهمية والادعاءات
يزعم هذا العمل تقديم أول وصف ديناميكي لإنتاج النوى الخفيفة في تصادمات الأيونات الثقيلة ذات الطاقة المتوسطة يتضمن جميع الأنواع حتى $A=4$ ويأخذ تأثير "موت" في الاعتبار صراحةً. وتثبت الدراسة أن:

1. يجب معاملة النوى الخفيفة كدرجات حرية ديناميكية طوال تطور التصادم لضمان وصف دقيق للعوائد.
2. الزيادة الملحوظة في عوائد جسيمات ألفا عند الطاقات المنخفضة هي نتيجة مباشرة لتأثير "موت"، وتحديداً مقاومة جسيم ألفا للذوبان في الوسط النووي بسبب طاقة ربطه العالية.
3. يتيح النهج الحركي تحديد معاملات العتبة $f^{cut}_A$، والتي تعمل كوكلاء لقوة تأثيرات "موت" في المادة النووية.

يشير المؤلفون إلى أنه بينما يعيد النهج الحالي إنتاج البيانات لـ $A \le 4$، فإن وصف الشظايا الأثقل ($A \ge 5$) المتكونة عبر التحلل التلقائي (spinodal decomposition) يتطلب توسيع الصيغة لتشمل التقلبات في توزيعات فضاء طور النوكليونات. ويقترحون تطبيقات مستقبلية لهذا الإطار في مجالات التدرج النظائري (isoscaling)، وتكون عناقيد ألفا على النوى الثقيلة، والسيناريوهات الفيزيائية الفلكية المتعلقة بانفجار النجوم العظمى (core-collapse supernovae) والنجوم النيوترونية المتراصة.