Determination of nuclear deformations with an emulator for sub-barrier fusion reactions

ملخص تقني لورقة بحثية بعنوان: "تحديد التشوهات النووية باستخدام محاكي لتفاعلات الاندماج تحت حاجز كولوم"

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

• السياق: يؤثر التشوه النووي بشكل كبير على تفاعلات اندماج الأيونات الثقيلة منخفضة الطاقة، لا سيما بالقرب من حاجز كولوم. تُظهر النوى المشوهة توزيعاً لارتفاعات الحواجز بناءً على اتجاهها، مما يؤدي إلى تعزيز مقاطع الاندماج العرضية.
• التحدي: يتطلب استخراج معاملات التشوه النووي (مثل التربيعي $\beta_2$، والكساسيليبول $\beta_4$، والأوكتوبول $\beta_3$) من مقاطع الاندماج التجريبية العرضية حل معادلات القنوات المزدوجة (CC) بشكل متكرر.
• العقبة: حسابات القنوات المزدوجة القياسية مكلفة حوسبياً. ولإيجاد معاملات التشوه المثلى، يجب إجراء بحث متعدد الأبعاد في المعاملات (تقليل $\chi^2$)، وهو ما يتضمن آلاف الحسابات المتكررة للقنوات المزدوجة. وهذا يخلق تكلفة حوسبية باهظة، مما يحد من الاستكشاف المنهجي للأشكال النووية.

2. المنهجية

يقترح المؤلفون نموذجاً بديلاً (محاكياً/Emulator) يعتمد على استمرار المتجهات الذاتية (Eigenvector Continuation - EC)، وهو تطبيق محدد لـ طريقة القاعدة المختزلة (Reduced Basis Method - RBM)، لتسريع حسابات القنوات المزدوجة مع الحفاظ على دقة عالية.

• الإطار النظري:

  • صيغة القنوات المزدوجة (Coupled-Channels Formalism): يستخدم إطار عمل كود CCFULL. يتضمن الهاملتوني الحركة النسبية، والاستثارة الجوهرية، وجهود الاقتران (النووية وكولوم) التي تعتمد على معاملات التشوه ($\beta_\lambda$).
  • طريقة القاعدة المنفصلة (Discrete Basis Method - DBM): يتم تمثيل معادلة شرودنجر المستمرة باستخدام شبكة محدودة. يستخدم المؤلفون مبدأ كوهن التبايني (Kohn variational principle) بالتزامن مع طريقة القاعدة المنفصلة لتحويل مشكلة التشتت إلى مسألة قيم ذاتية مصفوفية.
  • طريقة نوميروف (Numerov Method): لتعزيز الاستقرار والتقارب دون الحاجة إلى صقل مفرط للشبكة، طبق المؤلفون طريقة نوميروف المعدلة لمؤثر الطاقة الحركية بدلاً من صيغة الفروق المحدودة ثلاثية النقاط القياسية.

• بناء المحاكي (استمرار المتجهات الذاتية - EC):

  1. توليد اللقطات (مرحلة الـ Offline): يتم حساب الحلول عالية الدقة (الدوال الموجية للتشتت $\Psi_i$) لمجموعة صغيرة من معاملات التدريب المختارة (مثل قيم محددة لـ $\beta_2, \beta_4$).
  2. بناء الفضاء الجزئي: تشكل هذه الحلول فضاءً جزئياً منخفض الأبعاد.
  3. الإسقاط: يتم تقريب الدالة الموجية لأي مجموعة جديدة من المعاملات كتراكب خطي للقطات التدريب:
     $$\Psi_{EC} = \sum c_i \phi_i + \text{boundary terms}$$
  4. الحساب الفعال: يتم تحديد المعاملات $c_i$ عن طريق حل نظام خطي صغير مشتق من إسقاط الهاملتوني على القاعدة المختزلة. وهذا يتجنب حل المعادلة المصفوفية كاملة الأبعاد لكل مجموعة جديدة من المعاملات.

• سير العمل:

  1. تحديد نقاط التدريب في فضاء المعاملات.
  2. توليد لقطات عالية الدقة.
  3. بناء محاكي (EC).
  4. إجراء تقليل سريع لـ $\chi^2$ مقابل البيانات التجريبية لاستخراج معاملات التشوه المثلى.

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

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

4. النتائج

طُبق المحاكي على ثلاثة أنظمة تفاعلية محددة:

• الحالة 1: $^{16}\text{O} + ^{144}\text{Sm}$ (اهتزاز أوكتوبول)

  • الهدف: استخراج معامل تشوه الأوكتوبول $\beta_3$.
  • النتيجة: حدد المحاكي الحد الأدنى لـ $\chi^2$ عند $\beta_3 = 0.21$، وهو ما يطابق تماماً القيمة المستمدة من بيانات $B(E3)$ التجريبية وحسابات CC الدقيقة.
  • الدقة: كانت مقاطع الاندماج العرضية وتوزيعات الحواجز المحسوبة بواسطة المحاكي غير قابلة للتمييز عن النتائج الدقيقة.

• الح الحالة 2: $^{16}\text{O} + ^{154}\text{Sm}$ (تشوه ساكن)

  • الهدف: استخراج معاملي التربيعي ($\beta_2$) والكساسيليبول ($\beta_4$) في فضاء بحث ثنائي الأبعاد.
  • النتيجة: وجد المحاكي المعاملات المثلى $\beta_2 = 0.31$ و $\beta_4 = 0.06$. هذه القيم متسقة مع الحسابات الدقيقة والمسبارات التجريبية المستقلة (مثل تشتت ألفا غير المرن).
  • الارتباط: أعاد المحاكي إنتاج الارتباط السلبي القوي ($r \approx -0.82$) بين $\beta_2$ و $\beta_4$ الموجود في الحسابات الدقيقة.

• الحالة 3: $^{16}\text{O} + ^{186}\text{W}$ (تشوه ساكن)

  • الهدف: استخراج $\beta_2$ و $\beta_4$ لنواة ذات تشوه كساسيليبول سالب.
  • النتيجة: أعطى المحاكي $\beta_2 = 0.29$ و $\beta_4 = -0.02$، وهي نتائج متسقة مع الحسابات الدقيقة وبيانات تشتت النيوترونات.
  • الملاحظة: بينما كانت المتوسطات المثلى دقيقة، بالغ المحاكي قليلاً في تقدير الارتباط العكسي بين المعاملات بالقرب من حدود شبكة التدريب، مما يشير إلى الحاجة لتحسين اختيار اللقطات في النسخ المستقبلية.

• الأداء الحوسبي:

  • سرعة التنفيذ: حقق المحاكي عامل تسريع يتراوح بين 200 إلى 400 مقارنة بحسابات DBM الدقيقة للأنظمة التي تحتوي على 2-4 قنوات.
  • التقارب: انخفضت الأخطاء العددية بنحو رتبة مقدار واحدة مع كل متجه قاعدة إضافي ($N_{EC}$).
  • نقطة التعادل: يصبح المحاكي ذا جدوى حوسبية بعد حوالي 5 تقييمات (مع مراعاة تكلفة الـ offline لتوليد اللقطات).

5. الأهمية والآفاق المستقبلية

• الكفاءة: يحول هذا الأسلوب مشكلة عكسية مستعصية حوسبياً (استخراج المعاملات) إلى مشكلة يمكن التعامل معها، مما يتيح الاستكشاف المنهجي للأشكال النووية.
• الموثوقية: يحافظ المحاكي على دقة عالية تجاه الفيزياء الأساسية، حيث يعيد إنتاج مقاطع الاندماج العرضية وتوزيعات الحواجز بدقة عبر كامل نطاق الطاقة.
• القابلية للتوسع: الإطار جاهز للتوسع ليشمل تشوهات أكثر تعقيداً (مثل التشوه ثلاثي المحاور $\gamma$، أو الرتب المتعددة الأعلى) والأنظمة الأثقل (مثل $^{238}\text{U}$).
• التأثير الأوسع: يربط هذا النهج بين فيزياء التفاعلات النووية منخفضة الطاقة ودراسات تصادم الأيونات الثقيلة عالية الطاقة (حيث تُستقصى الأشكال النووية أيضاً عبر ملاحظات التدفق)، مما يوفر أداة موحدة لفهم الخصائص النووية الجوهرية.
• العمل المستقبلي: يخطط المؤلفون لدمج التحسين الخوارزمي لاختيار اللقطات واستخدام تحليل المكونات الرئيسية (PCA) لتقليل حجم القاعدة بشكل أكبر والتعامل مع فضاءات المعاملات ذات الأبعاد الأعلى.