Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

ملخص تقني: بنية وديناميكيات النوى ذات الغلاف المفتوح من نظرية كابلد-كلستر الكروية

بيان المشكلة
بينما حققت نظرية النواة القائمة على المبادئ الأولية (ab initio) نجاحاً ملحوظاً مع النوى ذات الغلاف المغلق باستخدام طريقة كابلد-كلستر (CC)، فإن توسيع هذه القدرات لتشمل الأنظمة ذات الغلاف المفتوح يظل تحدياً قائماً. تعتمد النهج التقليدية غالباً على تقنيات كسر التناظر (مثل Bogoliubov CC أو Deformed CC) متبوعة باستعادة التناظر، مما يزيد من التعقيد الحسابي بشكل كبير. بدلاً من ذلك، توفر طريقة معادلة الحركة (EOM) لـ CC إطاراً يحافظ على التناظر، حيث تُعامل النوى ذات الغلاف المفتوح كحالات مثارة لنواة ذات غلاف مغلق. وبينما تم إرساء منهجية EOM-CC القائمة على إضافة جسيمين (2PA) للأنوية التي تمتلك جسيمين إضافيين، وتم ربطها مؤخراً بتحويل لورنتز المتكامل (LIT) لحساب دالات الاستجابة، فقد كان هناك نقص في صياغة مقابلة قائمة على إزالة جسيمين (2PR) للأنوية التي تمتلك جسيمين أقل. تعالج هذه الورقة الحاجة إلى توسيع إطار كابلد-كلستر الكروي لوصف بنية واستجابة ثنائي القطب للأنوية ذات الغلاف المفتوح الناتجة عن إزالة جسيمين من إغلاق الغلاف.

المنهجية
قام المؤلفون بتطوير وتنفيذ صياغة إزالة جسيمين (2PR) لـ EOM-CC ضمن إطار كابلد-كلستر الكروي، مقترنة بتقنية تحويل لورنتز المتكامل (LIT) لحساب دالات الاستجابة.

1. الصياغة:

   • الحالة الأرضية: تستخدم الطريقة نهج مرجع واحد (CCSD و CCSDT-1) لتحديد الحالة الأرضية لنواة مرجعية ذات غلاف مغلق (كتلتها $A$). يتم تطبيق إزاحة للكتلة بحيث تقابل طاقة المرجع نواة الهدف (كتلتها $A-2$).
   • 2PR-EOM: يتم وصف نواة الهدف ذات الغلاف المفتوح عبر مؤثر استثارة يقوم بإزالة فرميونين من المرجع. تتضمن الصيغة (ansatz) تكوينات $0p-2h$ (ثقبين) في الرتبة الأساسية، مكملة بمساهمات $1p-3h$. يتم بناء المؤثرات اليمنى ($\hat{R}$) واليسرى ($\hat{L}$) لحل مسألة القيم الذاتية غير الهيرميتية لطاقات الاستثارة بالنسبة للمرجع المزاح الكتلة.
   • دالات الاستجابة: لحساب الاستجابات الكهرومغناطيسية، يتم دمج صيغة 2PR-EOM مع تقنية LIT. يتضمن ذلك حل معادلات تشبه معادلة شرودنجر غير متجانسة لحالات مساعدة باستخدام طريقة لانزوس (Lanczos). ويتم استخراج قطبية ثنائي القطب الكهربائي ($\alpha_D$) من عزمات LIT.
   • التفاعلات: تستخدم الحسابات تفاعلات نظر المجال الفعال الكيرالي ($\Delta$NNLOGO(394)) بما في ذلك قوى الجسيمين والثلاثة جسيمات.

2. مخططات التقطيع (Truncation Schemes):

   • تُجرى حسابات الحالة الأرضية عند مستويات CCSD (يُرمز لها بـ D) و CCSDT-1 (يُرمز لها بـ T-1).
   • تختلف حسابات EOM باختلاف تقطيع مؤثرات الاستثارة (على سبيل المثال، $0p-2h$ مقابل $1p-3h$ لـ 2PR).
   • تُقدر أوجه عدم اليقين من خلال مقارنة مخططات التقطيع المختلفة (مثل D مقابل T-1 للحالات الأرضية؛ وتقطيعات EOM المختلفة للأنظمة ذات الغلاف المفتوح).

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

• تطوير 2PR-EOM-CC: تقدم الورقة الصياغة التفصيلية لمنهج 2PR-EOM-CC، مما يسمح بدراسة النوى التي تحتوي على ثقبين بالنسبة لإغلاق الغلاف.
• التكامل مع LIT: وسع المؤلفون تقنية 2PA-LIT-CC إلى قطاع 2PR، مما يسمح بالحساب المبدئي (ab initio) لدالات استجابة ثنائي القطب الكهربائي والقطبية في النوى ذات الغلاف المفتوح.
• التحقق الشامل: تم التحقق من المنهج عبر سلاسل نظائر الأكسجين والكالسيوم، لتغطية طاقات الحالة الأرضية، وحالات مختارة مثارة، وقطبية ثنائي القطب الكهربائي.

النتائج

• طاقات الحالة الأرضية:
  • بالنسبة لـ $^{22}$O و $^{38}$Ca، يستعيد نهج 2PR-EOM مع استثارات $1p-3h$ (المشار إليها بـ D/1p-3h) حوالي 85% من طاقة الارتباط التجريبية.
  • يؤدي تضمين الثلاثيات (triples) في الحالة الأرضية المرجعية (T-1/1p-3h) إلى تحسين الدقة بشكل كبير، مما يجعل النتائج ضمن نطاق $\sim$1% من القيم التجريبية ويطابق دقة حسابات CCSDT-1 المستقرة للغلاف المغلق.
  • تُظهر فحوصات الاتساق بين مخططات 2PA و 2PR لـ $^{38}$Ca التقارب نحو طاقات مماثلة عند تضمين استثارات من رتب أعلى.
• الحالات المثارة:
  • في $^{14}$C، تصف طريقة 2PR-EOM الحالات ذات التكافؤ الموجب (مثل $2^+$) بدقة، لكنها تواجه صعوبة مع الحالات ذات التكافؤ السالب ($1^-$)، والتي تتطلب تكوينات عابرة للغلاف (cross-shell) وارتباطات من رتب أعلى لا تستوعبها التقطيعات الحالية بالكامل.
  • في $^{22}$O، تنجح طريقة 2PR-EOM في إعادة إنتاج ترتيب الحالات الخمس الدنيا، متفوقة على الصيغة التقليدية للغلاف المغلق في هذه الحالة المحددة.
• قطبية ثنائي القطب الكهربائي ($\alpha_D$):
  • تُظهر حسابات $^{22}$O و $^{38}$Ca أن نتائج 2PR-LIT-CC تقلل عموماً من قيم $\alpha_D$ التجريبية مقارنة بتنبؤات CC للغلاف المغلق.
  • يُعزى هذا النقص إلى انزياح قوة ثنائي القطب إلى طاقات أعلى في حسابات 2PR، مما يؤدي إلى فقدان مساهمات نمط ثنائي القطب "اللين" (soft dipole mode) منخفض الطاقة المشاهدة في حسابات الغلاف المغلق و SCGF.
  • بالنسبة لسلسلة الكالسيوم، تميل كل من مخططات 2PA و 2PR إلى التقليل من قيمة $\alpha_D$، مما يشير إلى ضرورة وجود مساهمات جسيمات متعددة من رتب أعلى للحصول على تنبؤات دقيقة لدالة الاستجابة.

الأهمية والادعاءات
تدعي الورقة توفير إطار موحد يحافظ على التناظر لدراسة النوى ذات الغلاف المفتوح بالقرب من الإغلاقات الفرعية للغلاف باستخدام طريقة كابلد-كلستر الكروية. ويؤكد المؤلفون أن نهج 2PR-EOM-CC يحقق دقة تقارب دقة CCSDT-1 للغلاف المغلق فيما يتعلق بطاقات الارتباط والحالات المثارة منخفضة الطاقة عند استخدام التقطيعات المناسبة. ومع ذلك، يشيرون بتواضع إلى أنه بينما تعد الطريقة قوية للخصائص البنيوية، فإنها تقلل حالياً من تقدير قطبية ثنائي القطب الكهربائي، خاصة في نظائر الكالسيوم. يسلط هذا التباين الضوء على الحاجة إلى ارتباطات جسيمات متعددة من رتب أعلى في حسابات دالة الاستجابة. يرسخ هذا العمل تقنية 2PR كأداة قابلة للتطبيق لتوسيع نطاق النظرية النووية المبدئية (ab initio) لتشمل الأنظمة ذات الغلاف المفتوح، مع خطط مستقبلية لتضمين استثارات من رتب أعلى (حتى 2p-4h) وتطبيق مخططات إضافة/إزالة جسيم واحد على النوى الفردية.