The Role of Magnetic Reconnection in Energizing Protons and Heavier Ions at the Heliospheric Current Sheet

ملخص تقني: دور إعادة الاتصال المغناطيسي في طاقة البروتونات والأيونات الثقيلة عند مستوى التيار الهليوسفيري

بيان المشكلة
كشفت الملاحظات الحديثة في الموقع بواسطة مسبار باركر الشمسي (PSP)، خلال عبوراته بالقرب من الشمس لمستوى التيار الهليوسفيري (HCS)، عن وجود مجموعات من البروتونات فوق الحرارية والأيونات الثقيلة (He, O, Fe) التي تم تسريعها إلى طاقات تترا-وحوالي 10 إلى 100 كيلو إلكترون فولت لكل نوية (keV nucleon$^{-1}$). تشير هذه الملاحظات إلى أن إعادة الاتصال المغناطيسي هي آلية تسريع متكررة وفعالة، تنتج أطياف طاقة ذات قوانين قوة (power-law) بمعاملات طيفية $\delta \sim 4-6$ وانقطاعات طاقة عالية ($E_{max}$) تتناسب مع نسبة الشحنة إلى الكتلة للأيون ($Q/M$) وفق العلاقة $E_{max} \propto (Q/M)^\alpha$، حيث $\alpha \sim 0.6 - 1.7$.

ومع ذلك، فإن جهود النمذجة السابقة باستخدام محاكاة الجسيمات في الخلايا (PIC) الحركية محدودة بسبب القيود الحسابية في حل التباين الهائل في المقاييس بين أطوال القصور الذاتي للأيونات ($\sim 10$ كم) ومناطق إعادة الاتصال الكلية ($\sim 10^6$ كم) عند مستوى التيار الهليوسفيري (HCS). ومن ناحية أخرى، واجهت نماذج الديناميكا المائية المغناطيسية (MHD) واسعة النطاق المرتبطة بمعادلات النقل صعوبة في إعادة إنتاج التناسب الملحوظ لـ $E_{max}$ مع $Q/M$. وتحديداً، وجدت دراسة سابقة (Murtas et al. 2024) استخدمت نهج النقل قيمة معامل تناسب $\alpha \sim 0.44$، وهي قيمة تقل عن الحد الأدنى لملاحظات مسبار باركر (PSP). وقد فُرض أن هذا التباين ينبع من نموذج حقن جسيمات مبسط للغاية افترض طاقة حقن ثابتة لكل نوية لجميع الأنواع، متجاهلاً الفيزياء المعتمدة على الكتلة لحقن الأيونات والمستمدة من المحاكاة الحركية.

المنهجية
لمعالجة هذه القيود، يستخدم المؤلفون نهج نمذجة هجين يجمع بين محاكاة MHD ثنائية الأبعاد (2D) وحلال نقل حركي.

1. محاكاة MHD: تم إجراء محاكاة MHD مقاومة ثنائية الأبعاد ذات عدد لوندكويست مرتفع ($S \sim 1.3 \times 10^5$) لورقة تيار تعيد الاتصال باستخدام كود Athena++. تحاكي المحاكاة ورقة تيار هاريس بمعلمات مطابقة لملاحظات مسبار باركر (لقاءات E07/E08)، بما في ذلك مجال مغناطيسي علوي قدره 200 نانو تسلا وسرعة ألففين قدرها 112 كم/ثانية. يتضمن الإعداد مجالاً موجهاً ضعيفاً ولزوجة لتحفيز عدم استقرار التمزق وتشكيل البلازمويد.
2. نقل الجسيمات: يتم حل معادلة باركر للنقل باستخدام كود (GPAT) العالمي لتسريع الجسيمات. يقوم هذا الكود بدمج المعادلات التفاضلية العشوائية للجسيمات النشطة باستخدام المجالات المغناطيسية والتدفقات البلازمية المعتمدة على الزمن من محاكاة MHD.
3. فيزياء الحقن: تدرس الدراسة أربعة سيناريوهات حقن متميزة (المسوحات A، B، C، وD) لتحديد كيفية تناسب طاقة الحقن لكل نوية ($E$) مع كتلة الأيون ($M$):
   • المسح A وB: طاقة حقن كلية ثابتة ($E_0 = 5$ keV و $10$ keV على التوالي)، مما يعني أن $E \propto 1/M$ (هيمنة السرعة الحرارية).
   • المسح C: طاقة حقن ثابتة لكل نوية ($E = \text{const}$)، مما يعني أن $E \propto 1$ (هيمنة سرعة التدفق الخارج).
   • المسح D: نظام وسيط حيث $E \propto 1/\sqrt{M}$.
     يتم حساب معاملات الانتشار باستخدام نظرية شبه خطية (QLT) مع طول ارتباط اضطرابي قدره $5 \times 10^4$ كم.

النتائج الرئيسية
تم تشغيل المحاكاة حتى الوصول إلى حالة شبه مستقرة ($t \approx 16.5 \tau_A$)، مما سمح بتحليل توزيعات تدفق الطاقة ($J$) والخصائص الطيفية.

• المعاملات الطيفية ($\delta$): أنتجت جميع المسوحات توزيعات قانون القوة بمعاملات طيفية تتراوح بين $\delta \sim 2.8$ و $5.4$.

  • في المسح C (حقن مستقل عن الكتلة)، زاد $\delta$ مع كتلة الأيون (الأيونات الأثقل كانت تمتلك أطيافاً أكثر ليونة)، وهو ما يتعارض مع ملاحظات مسبار باركر.
  • في المسوحات A وB وD (حقن معتمد على الكتلة)، انخفض $\delta$ مع الكتلة (البروتونات كانت تمتلك الأطياف الأكثر ليونة)، وهو ما يتماشى مع اتجاه "الأيون الأخف، الطيف الأكثر ليونة" الملحوظ بواسطة مسبار باركر.
  • كان للطاقة الكلية للحقن ($E_0$) تأثير هامشي على ميل الطيف ولكنه أثر على التقييس ومدى الوصول إلى الطاقة القصوى.

• انقطاعات الطاقة العالية ($E_{max}$): حققت البروتونات باستمرار أعلى طاقات انقطاع ($E_{max} \sim 93-97$ keV)، بينما أظهرت الأيونات الأثقل انقطاعات أقل. وقعت قيم $E_{max}$ الناتجة لجميع الأنواع ضمن نطاق $\sim 20-90$ keV لكل نوية، وهو ما يتوافق مع قياسات مسبر باركر.

• تناسب الشحنة إلى الكتلة ($\alpha$): حسبت الدراسة معامل التناسب $\alpha$ لـ $E_{max} \propto (Q/M)^\alpha$:

  • المسح C: $\alpha \approx 0.30$، وهو ما يتوافق مع نتيجة Murtas et al. (2024) السابقة، ولكنه غير متوافق مع بيانات مسبار باركر.
  • المسح A ($E \propto 1/M$): $\alpha \approx 1.13 \pm 0.11$.
  • المسح D ($E \propto 1/\sqrt{M}$): $\alpha \approx 0.81 \pm 0.13$.
    أعطى كل من المسح A والمسح D قيم $\alpha$ ضمن النطاق الواسع الملحوظ في مسبار باركر ($\alpha \sim 0.6 - 1.7$).

الأهمية والادعاءات
يزعم البحث أن التباين بين نماذج النطاق الكبير السابقة وملاحظات مسبر باركر فيما يتعلق بتناسب $Q/M$ لتسريع الأيونات مدفوع بشكل أساسي بالافتراضات الواردة في نموذج حقن الجسيمات. من خلال دمج فيزياء الحقن التي تأخذ في الاعتبار اعتماد طاقة الحقن على كتلة الأيون (تحديداً الأنظمة التي يساهم فيها كل من السرعة الحرارية وسرعة التدفق الخارج)، نجح نموذج النقل في إعادة إنتاج:

1. معاملات القوة الطيفية الملحوظة واعتمادها على كتلة الأيون.
2. انقطاعات الطاقة العالية المتوافقة مع البيانات في الموقع.
3. معامل التناسب الصحيح لنسبة الشحنة إلى الكتلة.

يخلص المؤلفون إلى أن إعادة الاتصال المغناطيسي عند مستوى التيار الهليوسفيري (HCS) هي آلية معقولة لإنتاج الأيونات الثقيلة فوق الحرارية، بشرط ألا تكون عملية الحقن مهيمنة تماماً على التدفق الخارج. ويقترحون أن المسح D (الاعتماد الوسيط على الكتلة) قد يكون الأكثر تمثيلاً فيزيائياً، لأنه ينتج معاملات طيفية تنخفض مع الكتلة (مطابقة للملاحظات) مع الحفاظ على معامل تناسب متوافق مع البيانات. ومع ذلك، يشير المؤلفون بتواضع إلى أن نماذجهم لا تزال تميل إلى إنتاج معاملات طيفية عند الحد الأدنى من القيم المرصودة، مما يشير إلى أن عوامل مثل نسبة معاملات الانتشار العمودي إلى الموازي ($\kappa_\perp/\kappa_\parallel$) والعمليات المهملة مثل تشتت زاوية الميل أو اضطراب الموجات السيكلوترونية قد تتطلب مزيداً من الصقل لاستيعاب تعقيد عملية تزويد الطاقة في مستوى التيار الهليوسفيري بشكل كامل.