Probability distribution of observables from a Bogoliubov vacuum projected onto good particle number: application to scission configurations of an actinide

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "توزيع الاحتمالات للملحوظات من فراغ بوغوليوبوف المسقط على عدد جسيمات جيد: تطبيق على تكوينات الانشطار في نويدات الأكتينيدات."

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

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

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

2. المنهجية

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

أ. الإطار النظري

• تمثيل الحالة: يُوصف النظام كحاصل ضرب تينسوري لغازات النيوترونات والبروتونات، $|\psi\rangle = |\psi_n\rangle \otimes |\psi_p\rangle$، حيث يتكون كل مكون من فراغ بوغوليوبوف مسقط على عدد ثابت من الجسيمات ($N$ و $Z$).
• فضاء التكوين: تعتمد الطريقة على أخذ عينات من التكوينات في تمثيل الموضع/اللف المغزلي الجوهري، حيث يتم تعريف التكوين بمجموعة من الإحداثيات $\{ \mathbf{r}_1, \dots, \mathbf{r}_N \}$، حيث $\mathbf{r}_i = (x, y, z, \sigma)$.
• كثافة الاحتمال: يتناسب احتمال تكوين معين مع مربع معيار تداخل الدالة الموجية: $p(\mathbf{r}_1 \dots \mathbf{r}_N) \propto |\langle \mathbf{r}_1 \dots \mathbf{r}_N | \psi \rangle|^2$.
• صياغة المحدد: بالنسبة لحالات بوغوليوبوف المسقطة، يمكن التعبير عن مربع هذا التداخل كمحدد لمصفوفة مائلة التناظر $Z$، تعتمد عناصرها على الدوال أحادية الجسيم الكنسية وسعات الازدواج ($u_k, v_k$).

ب. خوارزمية أخذ العينات (MCMC)

للتعامل مع الأبعاد العالية لفضاء التكوين (أكثر من 600 بُعد للأكتينيدات)، يستخدم المؤلفون نهج مونت كارلو لسلاسل ماركوف (MCMC) باستخدام خوارزمية ميتروبوليس:

1. التهيئة: يتم وضع النوكليونات في مناطق الكثافة العالية للدالة الكثافية أحادية الجسم مع لفات مغزلية عشوائية.
2. خطوة الاقتراح: يتم إنشاء تكوين مرشح عن طريق تغيير موضع نوكليون واحد (توزيع غاوسي بعرض $\sigma_{space}$) و/أو قلب لفه المغزلي (احتمالية $P_{flip}$).
3. القبول: يتم قبول الحركة بناءً على نسبة الاحتمال الجديد إلى القديم.
4. فترة الاحتراق والترقيق: لضمان التقارب وتقليل الارتباط الذاتي، يتم استبعاد "فترة احتراق" (burn-in) ويتم تسجيل التكوينات فقط كل $N_{jump}$ خطوة.

ج. تقدير الملحوظات

تقتصر الطريقة على الملحوظات التي تكون قطرية في تمثيل الموضع/اللف (مثل عدد الجسيمات، وعزوم متعدد الأقطاب، والتفاعلات ثنائية الجسم المحلية).

• بالنسبة لملحوظة قطرية $\hat{O}$، يكون القيمة المتوقعة هي متوسط النواة $O(\mathbf{r}_1 \dots \mathbf{r}_N)$ عبر التكوينات التي تم أخذ عينات منها.
• والأهم من ذلك، يتم الحصول على دالة توزيع الاحتمالات الكاملة للملحوظة مباشرة عن طريق رسم المدرج التكراري للقيم التي تحسبها النواة لكل تكوين تم أخذ عينة منه.

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

1. تقنية أخذ عينات مبتكرة: تطوير طريقة لأخذ عينات من مواضع النوكليونات ولفاتها المغزلية من فراغ بوغوليوبوف مسقط، مما يتيح حساب دوال توزيع الاحتمالات لملحوظات متعددة الأجسام التي كان يصعب الوصول إليها سابقاً في أطر دالة كثافة الطاقة النووية.
2. تنفيذ مفتوح المصدر: أصدر المؤلفون الكود NucleoScope (بلغة ++C) لتسهيل هذه الحسابات.
3. التحقق من الصحة: إجراء اختبارات مرجعية صارمة مقابل الحسابات الحتمية (في أساس المتذبذب التوافقي) للأنوية الخفيفة ($^{20}\text{Ne}$) والثقيلة ($^{252}\text{Cf}$)، مما أظهر التقارب والدقة ضمن حدود عدم اليقين الإحصائي.
4. التطبيق على الانشطار: تطبيق الطريقة على تكويني انشطار محددين لـ $^{252}\text{Cf}$ (نمطي Standard II و Super-Long) لتحليل التقلبات في الهندسة، وطاقات التفاعل، وعزوم الدوران.

4. النتائج الرئيسية

طُبقت الدراسة على تكويني انشطار لـ $^{252}\text{Cf}$: نمط Standard II (SII) (غير متماثل) ونمط Super-Long (SL) (متماثل).

أ. التقلبات الهندسية

• النظام المركب: تظهر عزوم الاستطالة ($\beta_{20}$) والأوكتوبول ($\beta_{30}$) للنواة المركبة توزيعات قريبة من التوزيع الطبيعي (الغاوسي) مع تقلبات نسبية صغيرة ($\sim 2\%$ لـ $\beta_{20}$).
• منطقة العنق: يظهر عدد الجسيمات في العنق ($Q_{neck}$) تقلبات كبيرة (تقلب نسبي $\approx 100\%$) وتوزيعاً غير غاوسي، حيث يصل ذروته غالباً عند صفر جسيمات.
• الشظايا: تتقلب أشكال ما قبل الشظايا بمقدار $0.05$إلى$0.1$وحدة. تشير تقلبات الشكل هذه إلى أن تقلب طاقة الإثارة يتراوح بين$1\text{--}2$ ميجا إلكترون فولت فقط، وهو ما لا يكفي لتفسير كامل التقلب التجريبي في تعدد النيوترونات الفورية.

ب. طاقات التفاعل بين الشظايا

• تقلب الطاقة الحركية الكلية (TKE): وجدت الدراسة أن جزءاً كبيراً من التقلب التجريبي في الطاقة الحركية الكلية (TKE) لشظايا الانشطار ينشأ من تقلب جهد التفاعل بين الشظايا عند نقطة الانشطار.
• النووي مقابل الكولومبي: من المثير للدهشة أن التفاعل النووي المتبقي (قصير المدى) هو المصدر المهيمن لهذا التقلب، وليس التنافر الكولومبي.
  • التفاعل النووي يرتبط بقوة ($\rho \approx 0.75$) مع عدد الجسيمات في العنق.
  • التفاعل الكولومبي أقل حساسية لجسيمات العنق ($\rho \approx 0.1$).
• الآلية: ينشأ التقلب لأن مواضع عدد قليل من النوكليونات في منطقة العنق تتغير بشكل كبير، مما يغير القوة النووية قصيرة المدى بين الشظايا المتباعدة.

ج. عزوم الدوران والزخم الزاوي

• عزوم الدوران المتبقية: حسبت الدراسة دوال توزيع الاحتمالات لعزوم الدوران التي تؤثر على الشظايا بسبب التفاعلات بين الشظايا.
• الأصل: على غرار الطاقة، فإن التفاعل النووي هو المحرك الأساسي لتقلبات عزم الدوران، بينما تتبع عزوم الدوران الكولومبية التوزيع الطبيعي.
• الآثار المترتبة: يكون عزم الدوران المتبقي موجهاً في الغالب بشكل عمودي على محور الانشطار. وبينما قدر المؤلفون أن هذا قد يساهم بنحو $8\hbar$ في الزخم الزاوي، إلا أنهم أشاروا إلى أن هذا قد يكون مبالغاً فيه بسبب الافتراضات الكلاسيكية، مما يشير إلى أن صورة المجال المتوسط تلتقط بالفعل جزءاً كبيراً من توليد الدوران.

5. الأهمية والاستنتاج

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

باختصار، يجسد هذا العمل الفجوة بين الأوصاف الساكنة للمجال المتوسط والطبيعة الإحصائية لنتائج الانشطار، موضحاً أن التقلبات الكمومية في مواضع النوكليونات عند نقطة الانشطار هي المحرك الأساسي للتباين الملحوظ في خصائص شظايا الانشطار.