R-matrix calculations for opacities: V. Temperature-density dependence of photoabsorption cross sections and opacity spectra of oxygen ions O VI and O VII

ملخص تقني: حسابات مصفوفة آر (R-Matrix) لدرجات العتامة: الجزء الخامس. اعتماد التعتامة على درجة الحرارة والكثافة للمقاطع العرضية للامتصاص الضوئي لأيونات الأكسجين O VI و O VII

بيان المشكلة
يتناول هذا العمل الحاجة إلى بيانات ذرية دقيقة لحل التناقضات في النماذج الشمسية الداخلية، وتحديداً "مشكلة الوفرة الشمسية" حيث تشير التحليلات الطيفية الحديثة ثلاثية الأبعاد (3D NLTE) إلى وفرة أقل للكربون والنيتروجين والأكسجين عما كان يُعتقد سابقاً. إن دمج هذه الوفرة المنخفضة في النماذج الشمسية القياسية (SSMs) يخلق صراعات كبيرة مع البيانات الهيليوسيزمية (علم الزلازل الشمسية)، بما في ذلك الانحرافات في ملفات سرعة الصوت وعمق قاعدة منطقة الحمل الحراري (BCZ). وبينما ركزت الأوراق السابقة من سلسلة مصفوفة آر للعتامة (RMOP) على أيونات الحديد، فإن الأكسجين أكثر وفرة من الحديد في الشمس بـ 17 مرة ويساهم بنحو 25% من إجمالي عتامة روسلاند المتوسطة عند قاعدة منطقة الحمل الحراري (BCZ). لذا، فإن النمذجة الدقيقة لأيونات الأكسجين، وخاصة O VI (شبيه الليثيوم) و O VII (شبيه الهيليوم)، أمر ضروري لحل هذه التناقضات. علاوة على ذلك، تفتقر بيانات مشروع العتامة (OP) الحالي إلى نطاق الطاقة اللازم وتفاصيل الرنين المطلوبة لظروف البلازما عالية كثافة الطاقة (HED).

المنهجية
يستخدم المؤلفون طريقة مصفوفة آر لـ "بريت-باولي" (Breit–Pauli R-matrix - BPRM) لحساب المقاطع العرضية للتأين الضوئي والعتامات أحادية اللون لـ O VI و O VII. يتضمن الإطار الحسابي ما يلي:

1. البنية الذرية: استخدام توسيعات دالة الموجة للارتباط الوثيق (Close-Coupling - CC) التي تتضمن مستويات البنية الدقيقة للأيون الأعلى التالي كحالات لُب (حتى $n \le 10, \ell \le 9$). يلتقط هذا النهج التداخل الكمي بين القنوات المفتوحة والمغلقة، مما يولد سلاسل كثيفة من رنينات التأين الذاتي (AI)، بما في ذلك ميزات التحفيز الضوئي للُّب (PEC) القوية.
2. توسيع البلازما: يتم تعديل مقاطع BPRM العرضية الأصلية غير الموسعة لتأخذ في الاعتبار تأثيرات بيئة البلازما (التصادمات الإلكترونية، مجالات الأيونات الدقيقة/تأثير ستارك، تأثير دوبلر الحراري، والانتقالات من نوع free-free). ويتم تحقيق ذلك عبر عملية تلافيف لورنتز (Lorentzian convolution) للمقاطع العرضية غير الموسعة، حيث يكون العرض الإجمالي $\Gamma$ هو مجموع مساهمات التصادم، وستارك، ودوبلر، وfree-free. عند ظروف قاعدة منطقة الحمل الحراري (BCZ) ($T \sim 10^6$ K، $N_e = 10^{22}–10^{23}$ سم$^{-3}$)، تم تحديد التوسيع الناتج عن اصطدام الإلكترونات كآلية مهيمنة.
3. الشبكة والظروف: تُجرى الحسابات على طول خطوط تساوي حرارة (isotherms) ممثلة لقاعدة منطقة الحمل الحراري الشمسية: $T = 1 \times 10^6$ K و $T = 2 \times 10^6$ K. وتتراوح كثافة الإلكترونات ($N_e$) من $10^{18}$ إلى $10^{23}$ سم$^{-3}$.
4. حساب العتامة: تُشتق العتامات أحادية اللون ($\kappa_\nu$) من خلال الجمع بين المقاطع العرضية للارتباط الحر (bound-free) الموسعة ومقاطع الخطوط (bound-bound) (التي تُعالج باستخدام نماذج التوسيع القياسية) وتوزيعات المستويات التي توفرها معادلة الحالة "ميالاس-هومر-دابن" (MHD).

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

• نطاق طاقة ممتد وبنية رنينية: على عكس قاعدة بيانات OP، التي تنتهي عند طاقات منخفضة وتفتقر إلى رنينات التأين الذاتي، فإن حسابات RMOP تمتد باستمرار إلى طاقات عالية ($\sim 250$ Ry لـ O VI و $\sim 700$ Ry لـ O VII). وهي تحل بوضوح هياكل رنينات التأين الذاتي (AI) الغنية التي تتقارب مع عتبات اللب المختلفة.
• نظميات توسيع البلازما: ترسم الدراسة خريطة اعتماد الرنين على درجة الحرارة والكثافة فيما يتعلق بالتحلل. وقد حددت عتبات كثافة محددة تبدأ عندها رنينات AI في التوسع وتبدأ بعدها في التحلل تماماً في المستمر (continuum)، مما يؤدي إلى تسطيح المقطع العرضي الخلفي.
• المقارنة مع OP: يوفر العمل مقارنة مباشرة بين بيانات BPRM الجديدة ونتائج OP الموروثة، مع تسليط الضوء على الاختلافات الكمية في كل من المقاطع العرضية والعتامات الناتجة.

النتائج

• المقاطع العرضية للتأين الضوئي:
  • O VI (شبيه الليثي): تظهر مقاطع RMOP العرضية بنية غنية من رنينات AI تمتد لعدة رتب مقدارية فوق المستمر، وهي غائبة في بيانات OP. ومع زيادة كثافة الإلكترونات، تتلاشى هذه الرنينات تدريجياً. عند كثافات BCZ ($N_e = 10^{23}$ سم$^{-3}$)، يكون هيكل الرنين قد تحلل إلى حد كبير، مما ينتج عنه مقطع عرضي شبه مسطح في منطقة 40–60 Ry.
  • O VII (شبيه الهيليوم): لوحظ سلوك مماثل، مع وجود عناقيد رنين كثيفة في نافذة 40–63 Ry. تشير الدراسة إلى استجابة تفاضلية بين عناقيد الرنين (على سبيل المثال، 47–50 Ry مقابل 57–60 Ry)، حيث تتحلل العناقيد ذات الرنينات الأكثر تقارباً عند كثافات أقل.
• العتامات أحادية اللون:
  • O VI: تعطي RMOP عتامات أعلى بشكل منهجي من OP في منطقة الطاقة المنخفضة (0.02–0.2 keV)، وهو أمر بالغ الأهمية لعتامة روسلاند المتوسطة عند $T \sim 10^6$ K. تبلغ نسبة RMOP/OP المتوسطة حوالي 7 (الوسيط $\sim 4$)، مع تعزيزات محلية تتجاوز عشرة أضعاف عند ميزات الرنين.
  • O VII: عتامات RMOP أيضاً أعلى من OP، وإن كان التفاوت أقل (المتوسط $\sim 1.4$، الوسيط $\sim 1.2$). ويُعزى هذا التعزيز إلى أن تحلل هياكل الرنين يرفع العتامة الخلفية وإلى تضمين المزيد من تأثيرات الارتباط في معالجة BPRM.
• حدود التوسيع: يلخص الجدول 2 حدود الكثافة للتوسيع الرنيني. بالنسبة لـ O VI وعنقود O VII (47–50 Ry)، يبدأ التوسع عند $N_e \sim 10^{20}$ سم$^{-3}$، ويحدث التحلل بالقرب من $10^{22}–10^{23}$ سم$^{-3}$. أما عنقود O VII (57–60 Ry) فيستمر حتى كثافات أعلى ($\sim 10^{23}$ سم$^{-3}$).

الأهمية
تؤكد الورقة أن هذه النتائج المحددة للأيونات هي لبنات بناء أساسية لحسابات عتامة الأكسجين الكلية والخلطات الشمسية المستقبلية. وتشير النتائج إلى أن معالجات R-Matrix المحسنة، والتي تتضمن البنية الدقيقة النسبية، والارتباط الواسع، وتوسيع البلازما، تؤدي إلى عتامات أعلى بشكل منهجي مما كان متوقعاً سابقاً. يُقدم هذا العمل كعمل قابل للتطبيق بشكل عام لنمذجة بلازما HED في الفيزياء الفلكية وتحليل أطياف الإرسال في التجارب المختبرية، مثل تلك التي تُجرى على أجهزة الاندماج بالحصر العطالي (ICF). ويشير المؤلفون إلى أنه بينما تم قياس إجمالي عتامة الأكسجين في مرافق مثل Sandia Z، إلا أنه لم يكن من الممكن استنتاج مساهمات الأيونات الفردية من تلك التجارب، مما يؤكد ضرورة البيانات النظرية المقدمة هنا. ويعد هذا العمل جزءاً من جهد أوسع لمعالجة التناقضات في النماذج النجمية القياسية باستخدام بيانات ذرية محدثة.