Sensitivity of fission barrier heights to nuclear matter properties
تُظهر هذه الدراسة، باستخدام نظرية الدالة الوظيفية للكثافة المتغيرة تباينياً والمحاكيات القائمة على العملية الغاوسية، وجود ارتباط قوي بين معامل طاقة السطح وارتفاعات حواجز الانشطار في 240Pu و236U، مع تحديد ثابت اقتران مشتق المتجه-السكلي (isoscalar-vector) كعامل أساسي يحكم ارتفاع الحاجز الخارجي.
في قلب النواة الذرية يكمن توازن دقيق للقوى يحدد ما إذا كانت الذرة ستظل مستقرة أم ستنفصم. فداخل النواة، تتماسك البروتونات والنيوترونات معاً بفعل قوة جذب قوية، لكن البروتونات تدفع بعضها البعض أيضاً لأنها جميعاً تحمل شحنة كهربائية موجبة. هذا التجاذب والتنافر يخلق نوعاً من المشهد الطاقي. ولكي تنقسم النواة، أو تخضع لعملية الانشطار، يجب عليها تسلق تلة طاقة تُعرف بالحاجز. فإذا كانت التلة عالية، تظل النواة متماسكة؛ وإذا كانت التلة منخفضة، يمكن للنواة أن تتفكك بسهولة. إن فهم ارتفاع هذه التلال أمر بالغ الأهمية لكل شيء، بدءاً من توليد الطاقة النووية وصولاً إلى تفسير كيفية صهر العناصر الثقيلة في الانفجارات العنيفة للنجوم. ومع ذلك، فإن التنبؤ بدقة بمدى ارتفاع هذه الحواجز أمر صعب لأنه يتطلب حساب سلوك آلاف الجسيمات التي تتفاعل بطرق معقدة.
لقد اتخذت دراسة جديدة أجراها باحثون من جامعة كيوتو وجامعة سون يات سين نظرة جديدة على هذه المشكلة من خلال ربط سلوك النوى الفردية بخصائص المادة النووية بشكل عام. ركز الفريق على عنصرين ثقيلين، هما البلوتونيوم-240 واليورانيوم-236، وهما مثالان كلاسيكيان لذرات يمكن أن تنقسم. أراد الباحثون معرفة كيف ترتبط "طاقة السطح" للنواة — وهي أساساً تكلفة إنشاء مساحة سطح جديدة عندما تتمدد النواة وتتشوه — بارتفاع حاجز الانشطار. في "نموذج قطرة السائل"، وهو طريقة تفكير راسخة منذ فترة طويلة حول النوى، تُعد طاقة السطح هذه عاملاً رئيسياً. وقد استخدم الباحثون إطار عمل حاسوبي متطور يُسمى "نظرية الكثافة الوظيفية المتغيرة" (covariant density functional theory)، والتي تعامل النواة كأنها سحابة من الجسيمات التي تحكمها قواعد رياضية محددة، لمحاكاة كيفية سلوك هذه الذرات.
ولاستيعاب العدد الهائل من النتائج المحتملة، استخدم العلماء أداة إحصائية تُعرف باسم "المُحاكي ذو العملية الغاوسية" (Gaussian-process emulator). فبدلاً من إجراء عملية حسابية واحدة، قاموا بتوليد آلاف الاختلافات لنموذجهم النظري، حيث يتضمن كل منها إعدادات مختلفة قليلاً للقوى الأساسية داخل النواة. ثم استخدموا الخصائص المعروفة للبلوتونيوم واليورانيوم، مثل حجمهما وطاقة الترابط الخاصة بهما، لاستبعاد الاختلافات غير الواقعية والتركيز على الأكثر احتمالاً. سمح هذا لهم برسم النطاق المرجح لارتفاعات حاجز الانشطار ومعرفة مدى ارتباطها بمعامل طاقة السطح. وكشفت النتائج عن وجود رابط مباشر وقوي: فكلما زاد معامل طاقة السطح، ارتفعت حواجز الانشطار. وهذا منطقي فيزيائياً، حيث يتطلب السطح الأكثر صلابة طاقة أكبر للتمدد والتشوه، مما يرفع فعلياً التلة التي يجب على النواة تسلقها لتنقسم.
ذهبت الدراسة إلى أبعد من ذلك لتحديد أي الإعدادات المجهرية في نموذجهم كانت مسؤولة عن هذه النتائج. ووجدوا أن ارتفاع حاجز الانشطار الخارجي — وهي التلة الثانية والأكبر التي يجب على النواة عبورها — محكوم بالكامل تقريباً بإعداد محدد واحد في نموذجهم، يُعرف باسم "الاقتران الاشتقاقي المتجهي اللامتماثل" (isoscalar-vector derivative coupling). يتحكم هذا الإعداد في كيفية تغير القوة النووية مع تغير كثافة النواة. وفي المقابل، فإن الحاجز الداخلي، وهو التلة الأولى التي تواجهها النواة، يتأثر بمجموعة من الإعدادات المختلفة التي تعمل معاً. هذا التمييز مهم لأنه يخبر علماء الفيزياء أن أجزاء مختلفة من عملية الانشطار حساسة لجوانب مختلفة من القوة النووية. فبينما يعد الحاجز الخارجي انعكاساً مباشراً لتفاعل محدد، فإن الحاجز الداخلي هو نتيجة أكثر تعقيداً لعوامل متعددة.
كما فحص الباحثون كيفية ارتباط طاقة السطح بخصائص أخرى، مثل "معامل القابلية للانشطار"، الذي يقيس مدى احتمالية انقسام النواة. ووجدوا أن معامل طاقة السطح هو المحرك المهيمن لهذه الاحتمالية، وليس طاقة كولوم أو تأثيرات التماثل الأخرى. وهذا يؤكد أن قدرة النواة على مقاومة التشوه هي مسأ matter أساساً تتعلق بالتوتر السطحي. وتخلص الدراسة إلى أنه من خلال قياس خصائص الانشطار في المختبر، يمكن للعلماء الآن وضع قيود أكثر صرامة على معامل طاقة السطح، وهي قيمة يصعب قياسها مباشرة. ومن شأن هذا الفهم المحسن أن يصقل نماذجنا في الفيزياء الفلكية النووية، لا سيما في البيئات القاسية للنجوم النيوترونية، حيث تحدد الخصائص السطحية للمادة كيفية سلوك قشرة النجم وتطورها.
ملخص تقني: حساسية ارتفاعات حواجز الانشطار لخصائص المادة النووية
بيان المشكلة يُعد الانشطار النووي عملية حرجة في إنتاج الطاقة، وتخليق العناصر فائقة الثقل، والنشوء النووي الفيزيائي الفلكي. ومن الكميات المركزية في وصف الانشطار هو ارتفاع حاجز الانشطار، الذي يتطلب معالجة مجهرية متسقة للتشوه الجماعي واسع السعة وبنية القشرة أحادية الجسيم. وبينما نجحت نظرية الدالة الوظيفية للكثافة النووية (DFT) في وصف الانشطار عبر الخريطة النووية، فإن الربط بين خصائص المادة النووية الكلية وملاحظات الانشطار لا يزال أقل استكشافاً من الارتباطات الأخرى (مثل سمك القشرة النيوترونية وطاقة التماثل). وقد أثبتت الدراسات السابقة وجود ارتباط خطي بين معامل طاقة السطح (asurf) وارتفاعات الحواجز باستخدام بارامتريات "سكيرمي" (Skyrme)، ولكن كان ينقص البحث تحديد كمي دقيق لعدم اليقين البارامتري وتحديد المحركات المجهرية المحددة لهذا الارتباط.
المنهجية يستخدم المؤلفون نظرية الدالة الوظيفية للكثافة المتغايرة (CDFT) باستخدام دالة الاقتران النقطي PC-PK1 لاستقصاء الأصل المجهري للارتباطات بين طاقة السطح وخصائص الانشطار. تركز الدراسة على نواتين نموذجيتين من الأكتينيدات: 240Pu و 236U.
الإطار النظري:
حواجز الانشطار: يتم حساب أسطح طاقة الجهد (PES) ضمن نظرية المجال المتوسط النسبي (RMF)، بما في ذلك ارتباطات الازدواج (تقريب BCS)، وتصحيحات مركز الكتلة المجهرية، والتصحيحات الدورانية. يتم توليد سطح طاقة الجهد كدالة للتشوه الرباعي (β2) باستخدام قواعد التذبذب التوافقي المحوري.
استخراج طاقة السطح: يتم استخراج معاملات طاقة السطح المتماثل (asurf) وطاقة السطح المتباين (asym) من المادة النووية شبه اللانهائية (SINM) باستخدام تقريب توماس-فيرمي النسبي. وهذا يسمح بتوصيف خصائص السطح بشكل منفصل عن خصائص المادة النووية اللانهائية (INM).
التحليل الإحصائي: لاستكشاف فضاء البارامترات، قام المؤلفون بتغيير تسعة من ثوابت اقتران الـ EDF وهي (C=[αS,βS,γS,δS,αV,γV,δV,αTV,δTV]) باستخدام أخذ عينات "لاتين هايبركوب" (Latin hypercube sampling).
المحاكلات والاستدلال البايزي: تم تدريب محاكلات العمليات الغاوسية (GP)، وتحديداً PCA-GP لمنحنيات طاقة الجهد، على 200 مجموعة من بارامترات الـ EDF للتنبؤ بالملاحظات لمجموعة أكبر تضم 6000 مجموعة. ويُستخدم الاستدلال البايزي لبناء التوزيعات البعدية لهذه الكميات، المقيدة بخصائص الحالة الأرضية (طاقة الربط ونصف قطر الشحنة) لـ 240Pu و 236U.
تحليل الحساسية: تم إجراء تحليل حساسية عالمي (GSA) يعتمد على تأثير "شابلي" (Shapley effect) لتحديد مساهمة كل ثابت اقتران في تباين الملاحظات. وتأخذ هذه الطريقة في الاعتبار الارتباطات بين بارامترات الـ EDF البعدية.
المساهمات والنتائج الرئيسية
الارتباط بين asurf وحواجز الانشطار: يؤكد التحليل وجود ارتباط إيجابي قوي بين معامل طاقة السطح asurf وكل من ارتفاع الحاجز الداخلي (EA) والخارجي (EB). وهذا متسق فيزيائياً: فكلما زاد asurf، زادت صلابة السطح النووي، مما يزيد من تكلفة الطاقة اللازمة للتشوه.
الحساسية التفاضلية للحواجز:
ارتفاع الحاجز الداخلي (EA): يرتبط بارتفاع الحاجز الدالي ارتباطاً وثيقاً بـ asurf أكثر من الحاجز الخارجي. ويعزو المؤلفون ذلك إلى أن الحاجز الداخلي يحدث عند تشوهات أصغر حيث تهيمن المنافسة بين طاقة السطح وطاقة كولوم.
ارتفاع الحاجز الخارجي (EB): يُظهر ارتفاع الحاجز الخارجي ارتباطاً أضعف مع asurf مقارنة بالحاجز الداخلي. ويُعزى ذلك إلى التشوه الكبير وتكون العنق عند الحاجز الخارجي، حيث تصبح تأثيرات الانحناء التي لا تغطيها حد aبسيط ذات أهمية كبيرة.
الأصل المجهري (تحليل الحساسية):
الحاجز الخارجي (EB): يحكم ارتفاع الحاجز الخارجي بشكل أساسي ثابت الاقتران المشتق المتغاير-السكلي (δV). ويظهر تحليل حصة "شابلي" أن δV يفسر الغالبية العظمى من التباين في EB، بغض النظر عما إذا كانت الارتباطات بين ثوابت الاقتران الأخرى قد تم الاحتفاظ بها أم لا.
الحاجز الداخلي (EA) و asurf: يتم حكم هاتين الكميتين بشكل مشترك بواسطة δVومجموعةمنثوابتالاقترانالمترابطةتبادلياً(\beta_S, \gamma_S, \alpha_V$). ويُفسر التباين في هاتين الكميتين من خلال التأثير المشترك لهذه البارامترات الأربعة.
الملاحظات الأخرى: تُظهر طاقة كولوم (ECoul) حساسية ضئيلة تجاه δV، بينما يهيمن ثابت الاقتران المشتق المتباين-السكلي (δTV) على معامل تناظر السطح (asym).
تحديد عدم اليقين: تحدد الدراسة عدم اليقين البارامتري لملاحظات الانشطار. وُجد أن عدم اليقين النسبي للحاجز الخارجي (EB) يبلغ حوالي 25%، وهو أكبر بكثير من عدم اليقين في الحاجز الداخلي أو خصائص الحالة الأرضية، مما يشير إلى أن تقييد الـ EDF بخصائص الحالة الأرضية وحدها يوفر معلومات إضافية محدودة على الكميات المتعلقة بالانشطار مقارنة بالمعلومات السابقة، مما يعني أن بيانات الانشطار التجريبية (مثل أعمار النصف) ستكون ضرورية لتقييد هذه البارامترات بشكل أكبر.
الأهمية والادعاءات يزعم البحث أنه يقدم أول تحليل إحصائي شامل يربط ارتفاعات حواجز الانشطار بمعامل طاقة السطح وبالاقترانات المجهرية المحددة ضمن إطار اقتران نقطي متغاير.
التحديد المجهري: يحدد العمل δV كمحرك رئيسي للحاجز الخارجي للانشطار، مما يوفر هدفاً مجهرياً محدداً لتحسين دوال الـ EDF.
استراتيجية التقييد: يشير الارتباط القوي بين خصائص الانشطار و asurf إلى إمكانية استخدام بيانات الانشطار التجريبية كقيد لتحديد معامل طاقة السطح بدقة أكبر.
الأهمية الفيزيائية الفلكية: يشير المؤلفون إلى أن التحديد الأكثر دقة لـ asurf ذو صلة بالفيزياء الفلكية النووية، لا سيما في نمذجة تحولات الشكل في القشرة الداخلية للنجرات النيوترونية حيث قد تظهر أطوار "المعكرونة النووية" (nuclear pasta).
الاتجاهات المستقبلية: تقتصر الدراسة على التشوهات المحورية التناظرية. ويقترح المؤلفون أن توسيع الإطار ليشمل التشوهات ثلاثية المحاور ومجموعة أوسع من النوى المنشطرة يعد خطوة مستقبلية ضرورية.
يحافظ المؤلفون على نبرة متواضعة فيما يتعلق بنطاق العمل، مقرين بأنه بينما تقلل قيود الحالة الأرضية من عدم اليقين في الخصائص الكلية، إلا أنها غير كافية لتقييد حواجز الانشطار بشكل كامل دون دمج ملاحظات انشطار محددة.