Particle number projected energies at finite temperature

إليك ملخص تقني مفصل للورقة البحثية "طاقات إسقاط عدد الجسيمات عند درجات الحرارة المحدودة" (Particle number projected energies at finite temperature) بقلم تشين، كيانغ، وبي.

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

تعد نظرية الكثافة الوظيفية النووية (DFT) أداة قوية لوصف النوى الثقيلة وفائقة الثقل، لا سيما فيما يتعلق بديناميكيات الانشطار والاندماج. ومع ذلك، تعاني حسابات نظرية الكثافة الوظيفية القياسية (مثل هارتري-فوك-بوغوليوبوف أو HF+BCS) من كسر التناظر عند مستوى المجال المتوسط:

• كسر تناظر القياس (Gauge Symmetry Breaking): تكسر نظرية BCS وتحويل بوغوليوبوف حفظ عدد الجسيمات، مما يؤدي إلى تقلبات في عدد الجسيمات.
• تعقيدات درجة الحرارة المحدودة: عند درجات الحرارة المحدودة، يتم وصف الاستثارات الحرارية عادةً باستخدام المجموعة الكبرى (grand canonical ensemble)، والتي تؤدي بطبيعتها إلى تقلبات في عدد الجسيمات. مع ارتفاع درجة الحرارة، تتلاشى ارتباطات الازدواج (pairing correlations) والآثار الكمومية، ويصبح غياب حفظ عدد الجسيمات قضية هامة بالنسبة للنوى المحدودة.
• محدودية الطرق الحالية: بينما يعد إسقاط عدد الجسيمات (PNP) أمراً راسخاً عند درجة الصفر، فإن تنفيذه عند درجات حرارة محدودة يعد معقداً. غالباً ما تعتمد مناهج PNP السابقة عند درجات الحرارة المحدودة على تقريبات (مثل تقريب نقطة السرج أو تقريب غاوس المتقطع)، وتفتقر إلى الحسابات الدقيقة للطاقات المسقطة ضمن أطر نظرية الكثافة الوظيفية ذاتية الاتساق. علاوة على ذلك، فإن حساب الإنتروبيا بدقة عبر $Tr(\hat{D}_p \ln \hat{D}_p)$ أمر مكلف حاسوبياً بالنسبة للنوى الثقيلة.

المشكلة الجوهرية التي تعالجها الورقة هي غياب إطار عمل صارم وذاتي الاتساق لحساب الطاقات وحرية الطاقة المسقطة بدقة لعدد الجسيمات للنوى المركبة عند درجات حراریة محددة ضمن إطار Skyme DFT.

2. المنهجية

طوّر المؤلفون صيغة لدمج إسقاط عدد الجسيمات بدقة في حسابات Skyrme-Hartree-Fock+BCS ذاتية الاتساق عند درجات حرارة محددة.

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

  • مصفوفة الكثافة المسقطة: يوصف النظام بمصفوفة كثافة مسقطة $\hat{D}_p = \frac{1}{Z_p} \hat{P} \hat{D} \hat{P}$، حيث $\hat{P}$ هو مؤثر الإسقاط المعرف عبر تكامل على زاوية القياس $\theta$.
  • حساب الطاقة الدقيق: يتم حساب الطاقة المسقطة $E$ باستخدام صيغة الأثر (trace formula):
    $$E = \frac{1}{Z_p} \frac{1}{2\pi} \int d\theta e^{-i\theta N} \text{Tr}(e^{i\theta \hat{N}} \hat{D} \hat{H})$$
  • مبرهنة ويك المعممة: اشتق المؤلفون تعبيرات لمصفوفات أحادية الجسيم ($\tilde{\rho}, \tilde{\lambda}, \tilde{\kappa}$) بدلالة زاوية القياس $\theta$، وطاقات شبه الجسيمات ($E_p$)، وأعداد إشغال BCS ($u_p, v_p$). يسمح هذا ببناء دالة الطاقة $E(\theta)$ دون الاعتماد على تقريب نقطة السرج للطاقة نفسها.
  • تقريب الإنتروبيا: بما أن حساب الإنتروبيا الدقيق مكلف للغاية، فقد اعتمدوا الطريقة المستخدمة من قبل Fanto وآخرون، حيث يتم حساب الإنتروبيا عبر دالة التجزئة: $S_p = \ln(Z_p) - \beta \frac{\partial \ln(Z_p)}{\partial \beta}$. ويُلاحظ أن هذه الطريقة دقيقة فوق درجة الحرارة الحرجة ($T_c$) حيث يتلاشى فجوة الازدواج.
  • كثافة المستويات: يتم حساب كثافة المستويات $\rho_N(E^*)$ باستخدام تقريب نقطة السرج بناءً على دالة التجزئة المسقطة.

• التنفيذ الحسابي:

  • البرمجية: أُجريت الحسابات باستخدام كود SkyrmeAx في فضاءات إحداثيات متماثلة محورياً.
  • التفاعلات: استُخدمت قوة SkM* لقناة ثقب-جسيم (مناسبة للتشوهات الكبيرة)، واستُخدم تفاعل $\delta$ المعتمد على الكثافة لقناة الازدواج.
  • المنهج التبايني: تطبق الدراسة نهج "الإسقاط بعد التباين" (PAV). وبينما يعد "التباين بعد الإسقاط" (VAP) أكثر دقة من الناحية النظرية، إلا أنه مستحيل عملياً عند درجات الحرارة المحددة في DFT ذاتي الاتساق بسبب عدم القدرة على تطبيق مبرهنة ويك على عناصر مصفوفة الإنتروبيا.

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

1. اشتقاق صيغة PNP دقيقة عند درجات الحرارة المحددة: توفر الورقة أول اشتقاق وتنفيذ صارم لطاقات إسقاط عدد الجسيمات الدقيقة للنوى المركبة ضمن إطار Skyrme DFT ذاتي الاتساق عند درجات حرارة محددة.
2. تحليل نسب دالات التجزئة: تدرس الدراسة نسبة دالة التجزئة المسقطة إلى دالة التجزئة الكبرى ($\text{Z}_p/\Xi$)، مما يوضح الانتقال من التذبذب الفردي-الزوجي عند درجات الحرارة المنخفضة إلى توزيع غاوسي عند درجات الحرارة العالية.
3. مقارنة التقريبات: تقارن الدراسة بشكل منهجي بين طاقات PNP الدقيقة والطاقات الكانونية التقريبية (المشتقة من المشتقات الجزئية لدالة التجزئة) ونتائج FT-BCS غير المسقطة القياسية.
4. معاملات كثافة المستويات: استخرج المؤلفون معاملات كثافة المستويات ($a_0$ عند الحالة الأرضية و $a_f$ عند حاجز الانشطار) لـ $^{238}\text{U}$ و $^{292}\text{Fl}$، ومقارنة طريقة PNP مقابل تقريب غاوس المتقطع (DG).

4. النتائج

• دالة التجزئة والتذبذب:

  • عند درجات الحرارة المنخفضة ($T < 0.4$ MeV)، يهيمن مكون العدد الزوجي على النظام، مما يظهر تذبذباً فردياً-زوجياً واضحاً.
  • مع ارتفاع درجة الحرارة، تزدل مساهمة المكونات الفردية، ويتلاشى التذبذب.
  • فوق $T \approx 0.5$ MeV، يصبح توزيع مكونات عدد الجسيمات شبيهاً بالتوزيع الغاوسي، مدفوعاً بالآثار الحرارية.

• حواجز الانشطار:

  • اتجاهات الطاقة: طاقات PNP الدقيقة أقل بشكل منهجي (بمقدار $\sim 2.0$ MeV) من طاقات FT-BCS غير المسقطة بسبب تضمين ارتباطات إضافية.
  • ارتفاعات الحواجز: رغم الإزاحة في الطاقة المطلقة، فإن ارتفاعات حاجز الانشطار عند درجات الحرارة العالية ($T=1.0$ MeV) متشابهة جداً بين الحسابات المسقطة وغير المسقطة (5.47 MeV مقابل 5.20 MeV). يشير هذا إلى أن FT-BCS غير المسقط كافٍ لتقدير ارتفاعات الحواجز عند درجات الحرارة العالية.
  • درجة الصفر: عند $T=0$، يكون الفرق أكثر وضوحاً (7.38 MeV مع PNP مقابل 7.69 MeV بدونها)، مما يسلط الضال على أهمية PNP عندما تكون ارتباطات الازدواج قوية.
  • الطاقات التقريبية: تنحرف الطاقات الكانونية التقريبية (المعادلة 13) بشكل كبير عن طاقات PNP الدقيقة، خاصة عند حاجز الانشطار، مما يؤكد ضرورة الحسابات الدقيقة.

• كثافات المستويات:

  • تتفق كثافات المستويات المحسوبة لـ $^{238}\text{U}$ باستخدام طريقة PNP جيداً مع البيانات التجريبية عند تضمين عوامل التعزيز الجماعي.
  • يُظهر منهج غاوس المتقطع (DG) انحرافات عند طاقات الاستثارة المنخفضة ($E_x < 1$ MeV) لأنه يفشل في مراعاة التذبذب الفردي-الزوجي الناتج عن تفاعلات الازدواج عند درجات الحرارة المنخفضة.
  • المعاملات: وُجد أن نسبة معامل كثافة المستوى $a_f/a_0$ (الحاجز إلى الحالة الأرضية) هي 1.0718 لـ $^{238}\text{U}$ و 1.0625 لـ $^{292}\text{Fl}$ عند درجات الحرارة العالية. هذه القيم قريبة من القيمة التجريبية الشائعة وهي 1.1.
  • وُجد أن المعامل $a_0$ للحالة الأرضية هو تقريباً $A/9.97$ و $A/10.38$ للنوى المعنية، مما يوفر توجيهاً مجهرياً للنماذج الإحصائية.

5. الأهمية

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