Lindblad-driven quarkonium production in heavy-ion collisions

ملخص تقني: إنتاج الكوارك-أونيوم المدفوع بمعادلة ليندبلاد في تصادمات الأيونات الثقيلة

بيان المشكلة
تُعد حالات الكوارك-أونيوم الثقيلة مجسات دقيقة لـ "بلازما الكوارك-غلوون" (QGP) المتكونة في تصادمات الأيونات الثقيلة فائقة النسبية. ويحكم إنتاجها التفاعل المتبادل بين حجب "ديباي" للجهد الثقيل، وفقدان الترابط (decoherence) داخل الوسط، مما يؤدي إلى تثبيط هرمي يعتمد على طاقة الترابط. وبينما يُعرف أن الجهد داخل الوسط يكون ذا قيمة مركبة (بجزء حقيقي يصف الحجب وجزء تخيلي يرمز لظاهرة "تخميد لانداو")، إلا أنه كان يفتقر البحث إلى إطار نظري موحد يصف في آن واحد كلاً من تثبيط الحالات ما قبل التكون، وإعادة التوليد (الاندماج/coalescence) للكواركات الثقيلة المتوازنة حرارياً من المبادئ الأولى. وغالباً ما كانت المقاربات السابقة تعامل التثبيط وإعادة التوليد كآليتين ظاهرتين (phenomenological) مستقلتين.

المنهجية
يستخدم المؤلفون إطار نظام كمي مفتوح يعتمد على معادلة "ليندبلاد" لوصف تطور مصفوفة الكثافة لزوج الكوارك-كوارك مضاد ($Q\bar{Q}$) المقترنة بحمام بلازما الكوارك-غلوون (QGP). وتتقدم المنهجية عبر الخطوات التالية:

1. الجهد داخل الوسط: تستخدم الدراسة نموذج قانون غاوس (Lafferty and Rothkopf) لتمثيل الجهد المركب داخل الوسط $V(r, T) = \text{Re } V + i \text{Im } V$. يجمع هذا النموذج بين جهد "كورنيل" في الفراغ ونظرية الاضطراب "HTL" (Hard-Thermal-Loop)، ويحكمه بارامتر واحد يعتمد على درجة الحرارة وهو كتلة "ديباي" $m_D(T)$. وقد تم تنظيم الجزء التخيلي بدقة لإزالة التباينات غير الفيزيائية في حد الوتر (string term).
2. الدوال الطيفية وخصائص التفكك: بدلاً من حل معادلة شرودنجر المعتمدة على الزمن مع جهد تخيلي يعتمد على الموقع (وهو أمر صعب عددياً)، يقوم المؤلفون بحل معادلة شرودنجر في فضاء التردد. ويستخرجون الدوال الطيفية داخل الوسط عبر حساب الجزء التخيلي لدالة "جرين" الارتدادية (retarded Green's function) المتتبعة.
   • طاقات الترابط ($E_n$): تُحدد من مواقع القمم في الدوال الطلطية بعد طرح كتل الكواركات المكونة.
   • عرض التحلل الحراري ($\Gamma_n$): يُستخرج من العرض الكامل عند نصف القمة (FWHM) للقمم، مع استخدام ملف "برايت-ويغنر" مائل (skewed Breit-Wigner) لمراعاة عدم التماثل بالقرب من عتبة الاستمرار (continuum threshold).
   • درجات حرارة التفكك ($T_d$): تُحدد باعتبارها درجة الحرارة التي تندمج عندها القمة الطيفية في حالة الاستمرار (حيث تتلاشى طاقة الترابط).
3. احتمالية البقاء: بالنسبة لنظام يخضع لتمدد "بيوركن" (Bjorken expansion)، يتم حساب احتمالية بقاء حالة ما قبل التكون عن طريق تكامل عرض التحلل المعتمد على درجة الحرارة $\Gamma(T(t))$ عبر الزمن.
4. نموذج إعادة التوليد: تم توسيع الإطار ليشمل إعادة التوليد عبر إسقاط مؤثرات القفز العشوائي لمعادلة "ليندبلاد" على الفضاء الجزئي للحالة المرتبطة. وتحت تقريب "أدياباتيكي" (حيث يتطور الهاملتوني الفعال ببطء مقارنة بنطاقات زمن القفزة الكمية)، يتم اشتقاق معدل اندماج. يتناسب هذا المعدل مع مربع عدد الكواركات الثقيلة الحرة ($N^2$) وحجم الوسط ($V$)، ويتم تعديله بعامل استرخاء لمراعاة الوقت المحدود المطلوب لتوازن الكواركات الثقيلة حرارياً.

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

• إطار موحد: يشتق البحث كلاً من آليات التثبيط وإعادة التوليد من معادلة حركة واحدة (معادلة ليندبلاد) لمصفوفة الكثافة $Q\bar{Q}$، مما يضمن الاتساق الداخلي بين الجزء التخيلي للجهد (الذي يقود التحلل) ومؤثرات القفز (التي تقود إعادة التوليد).
• اشتقاق من المبادئ الأولى: يتجنب هذا النهج الجمع العشوائي بين النماذج المستقلة. فالجهد المركب يحدد عروض التحلل الحراري، والتي تحدد بدورها كلاً من احتمالية البقاء ومعايرة معدل إعادة التوليد عبر نظرية التقلب والتبدد (fluctuation-dissipation theorem).
• الاندماج الأدياباتيكي: يشتق المؤلفون نموذج اندماج محدد من معادلة "ليندبلاد"، يتضمن عامل استرخاء لتصحيح عدم التوازن الحراري اللحظي للكواركات الثقيلة في الوسط.

النتائج
تطبق الدراسة هذا الإطار على حالات "تشارمونيوم" ($J/\psi, \psi(2S)$) و"بوتومونيوم" ($\Upsilon(nS)$) في تصادمات رصاص-رصاص (Pb–Pb) عند طاقة $\sqrt{s_{NN}} = 5.02$ تريليون إلكترون فولت.

• درجات حرارة التفكك:
  • التشارمونيوم: $T_{J/\psi} \approx 372$ ميجا إلكترون فولت؛ $T_{\psi(2S)} \lesssim 155$ ميجا إلكترون فولت (تتفكك عند أو تحت درجة الحرارة الحرجة).
  • البوتومونيوم: $T_{\Upsilon(1S)} \approx 564$ ميجا إلكترون فولت؛ $T_{\Upsilon(2S)} \approx 223$ ميجا إلكترون فولت؛ $T_{\Upsilon(3S)} \approx 164$ ميجا إلكترون فولت؛ $T_{\Upsilon(4S)} < 155$ ميجا إلكترون فولت.
• عامل التعديل النووي ($R_{AA}$):
  • $J/\psi$: يتنبأ النموذج بمساهمة كبيرة لإعادة التوليد في التصادمات المركزية ($N_{part} \gtrsim 200$)، مدفوعة بالوفرة العالية لكواركات التشارم عند طاقات "LHC". هذه إعادة التوليد تعوض جزئياً التثبيط، مما يؤدي إلى $R_{AA}$ إجمالي قريب من الواحد الصحيح، وهو ما يتوافق مع بيانات تجارب ALICE وCMS.
  • $\Upsilon(1S)$: نظراً لأن وفرة كواركات البوتوم أقل بنحو رتبتين من كواركات التشارم، فإن حد إعادة التوليد يكون ضئيلاً. لذا، فإن $R_{AA}$ لـ $\Upsilon(1S)$ يهيمن عليه التثبيط بشكل شبه كامل، وتتفق نتائج النموذج بشكل معقول مع بيانات CMS دون الحاجة لمعاملات ضبط إضافية.

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