Lattice determination of the higher-order hadronic vacuum polarization contribution to the muon $g-2$

إليك ملخص تقني مفصل للورقة البحثية "تحديد الشبكة لمساهمة الاستقطاب الفراغي الهادروني من الرتبة الأعلى في العزم المغناطيسي للميون $g-2$" لـ أرناو بلتران وآخرون.

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

يُعد العزم المغناطيسي غير الطبيعي للميون ($a_\mu$) مراقبة رئيسية لاختبار النموذج القياسي (SM) والبحث عن فيزياء جديدة. وبينما وصلت الدقة التجريبية إلى 124 جزءاً في المليار (Fermilab E989)، فإن التنبؤ النظري محدود بمساهمة الاستقطاب الفراغي الهادروني (HVP).

• التناقض: يوجد توتر كبير بين تنبؤ النموذج القياسي والبيانات التجريبية. وتأتي معظم عدم اليقين في تنبؤ النموذج القياسي من مساهمة الـ HVP من الرتبة الرئيسية (LO) ($a_\mu^{\text{hvp, lo}}$).
• فجوة الرتبة التالية (NLO): يتم تقدير مساهمة الـ HVP من الرتبة التالية ($a_\mu^{\text{hvp, nlo}}$) حالياً باستخدام الأساليب التشتتية القائمة على البيانات (المعتمدة على مقاطع عرضية لعمليات $e^+e^- \to \text{hadrons}$). تعاني هذه التقديرات القائمة على البيانات من نفس التوترات التجريبية (مثل التناقضات بين CMD-3 والقياسات السابقة) التي تعاني منها تحديدات الرتبة الرئيسية (LO).
• الحاجة: لتوفير فحص متسق للنموذج القياسي من المبادئ الأولى، يجب حساب كل من مساهمات HVP للرتبة الرئيسية (LO) والرتبة التالية (NLO) باستخدام نفس الإطار (Lattice QCD). وقبل هذا العمل، لم يكن موجوداً حساب شبكي (Lattice) بدقة تقل عن واحد في المئة لمساهمة الرتبة التالية (NLO).

2. المنهجية

قام المؤلفون بإجراء أول حساب لـ $a_\mu^{\text{hvp, nlo}}$ باستخدام نظرية الكروموديناميكا الكمية الشبكية (Lattice QCD) بدقة تقل عن واحد في المئة باستخدام تمثيل الزمان-الزخم (TMR).

أ. الصيغة والمخططات

يغطي الحساب ثلاثة أصناف من مخططات الرتبة التالية (NLO) (كما هو موضح في الشكل 1):

1. NLOa: خطوط الفوتونات وحلقات الميون (مساهمة سالبة مهيمنة).
2. NLOb: حلقات الإلكترون والتاو (مساهمة موجبة، تلغي جزئياً NLOa).
3. NLOc: إدراجات QCD مزدوجة (ثانوية).

يتم التعبير عن المساهمة عبر موتر المتجه المسقط عند الزخم الصفري $G(t)$:
$$a_\mu^{\text{hvp, (i)}} = \left(\frac{\alpha}{\pi}\right)^3 \int_0^\infty dt_1 \dots dt_m \, \tilde{f}^{(i)}(\hat{t}_1, \dots) G(t_1) \dots G(t_m)$$
حيث $\tilde{f}^{(i)}$ هي دالات نواة الزمن مكانية للرتبة التالية (NLO).

ب. إعداد الشبكة (Lattice Setup)

• المجموعات (Ensembles): 35 مجموعة من نماذج CLS (محاكاة الشبكة المنسقة) مع 3 أنواع من النكهات ($N_f = 2+1$) من فيرميونات ويلسون المحسنة غير اضطرابياً من الدرجة $\mathcal{O}(a)$.
• المعلمات: ستة مسافات شبكية ($a \approx 0.039 - 0.097$ فرسخ فلكي) وكتل بيون تتراوح من $\sim 430$ ميجا إلكترون فولت نزولاً إلى النقطة الفيزيائية.
• تحديد المقياس: تم تحديده عبر الملحوظة المعتمدة على تدفق التدرج (gradient-flow) $\sqrt{t_0}$ حيث $\sqrt{t_0^{\text{ph}}} = 0.1440(7)$ فرسخ فلكي.
• التيارات: تم استخدام تنفيذين لمتجه الكهرومغناطيسية (المحلي والمنفصل نقطياً) للتحكم في أخطاء التجزئة (discretization errors).

ค. الاستراتيجيات التقنية الرئيسية

1. إلغاء NLOa و NLOb: ميزة حاسمة هي الإلغاء القوي بين نواتي (kernels) NLOa و NLOb عند الأزمان الإقليدية الكبيرة ($t \gtrsim 1.5$ فرسخ فلكي). من خلال حساب المجموع NLOa&b مباشرة، يقوم المؤلفون بكبح مساهمة المسافة الطويلة (LD)، والتي عادة ما تكون مصدر الضجيج الإحصائي الأكبر وعدم اليقين في الحجم المحدود (FV) في حسابات الرتبة الرئيسية (LO).
2. معالجة المسافة القصيرة (SD): تحتوي نواة NLOa على آثار اصطناعية معززة لوغاريتمياً من النوع $\mathcal{O}(a^2 \ln^2 a)$. تمت إزالة هذه الآثار عن طريق طرح سلوك النواة الرئيسي واستعادة الجزء اضطرابياً باستخدام الاستقطاب الفراغي لـ QCD حتى الرتبة $\mathcal{O}(\alpha_s^4)$.
3. التحكم في ضجيج المسافة الطويلة (LD): تمت معالجته عبر:
   • متوسط النمط المنخفض (LMA) للـ كواركات الخفيفة.
   • إعادة بناء الطيف لذيل الموتر باستخدام تحليل GEVP.
   • طرق التحديد القياسية (bounding methods).
4. تصحيحات الحجم المحدود (FV): تم تطبيقها باستخدام إجراء من خطوتين: صيغة Hansen–Patella للمسافات الصغيرة $t$ وطريقة Meyer–Lellouch–Lüscher للمسافات الكبيرة $t$، متبوعة باستكمال مستمر (continuum extrapolation) باستخدام نظرية الاضطراب الكيرال (chiral perturbation theory) من الرتبة الثانية (NNLO).
5. كسر التماثل الإيزوسبيني (IB): بما أن المجموعات متماثلة إيزوسبينياً، فقد تم تقدير تصحيحات IB (الكهرومغناطيسية والقوية) ظواهرتياً باستخدام التمثيل المكاني. ومن المهم أن إلغاء NLOa/NLOb يقلل أيضاً من إزاحة IB، مما يؤدي إلى تصحيح إجمالي صغير ($\Delta_{\text{IB}} \approx 0.06 \times 10^{-11}$).

3. النتائج الرئيسية

حصل المؤلفون على النتائج التالية في الحد المستمر (بوحدات $10^{-11}$):

• NLOa: $-216.83(81)(87)(99)(81)(34)(31)[1.81]$
• NLOb: $+111.51(55)(64)(64)(58)(16)(25)[1.24]$
• NLOc: $+3.75(6)(6)(3)(3)[10]$
• إجمالي NLO HVP:
  $$a_\mu^{\text{hvp, nlo}} = -101.57(26)_{\text{stat}}(29)_{\text{syst}} \times 10^{-11}$$
  • إجمالي الخطأ النسبي: 0.6% (بدمج عدم اليقين الإحصائي والمنهجي).

تفصيل عدم اليقين:
يهيمن على إجمالي عدم اليقين تقريباً بالتساوي كل من:

1. الضجيج الإحصائي للمسافة الطويلة.
2. انتشار متوسط النموذج (الاستكمال الكيرالي-المستمر).
3. عدم يقين تحديد المقياس ($t_0$).
4. تصحيحات كسر التماثل الإيزوسبيني (IB).

4. الأهمية والمقارنة

• الدقة: هذه النتيجة أكثر دقة بمرتين من تقدير تحديث الورقة البيضاء (WP25) لعام 2025، والذي اعتمد كلياً على الطرق التشتتية القائمة على البيانات.
• فحص الاتساق: توفر هذه النتيجة تحديداً متسقاً تماماً من الشبكة لكل من LO و NLO HVP، مما يزيل عدم الاتساق المنهجي حيث كانت الرتبة LO فقط هي المعتمدة على الشبكة في التنبؤات السابقة للنموذج القياسي.
• التوتر مع التقديرات القائمة على البيانات:
  • النتيجة أقل بـ 1.4$\sigma$ من متوسط WP25 ($-99.6 \pm 1.3$).
  • تظهر توتراً قدره 4.6$\sigma$ مع تقييم KNT19 القائم على البيانات (الذي يستبعد بيانات CMD-3).
  • يعكس هذا النمط التوتر الملاحظ في LO HVP بين حسابات الشبكة ونتائج التشتت السابقة لـ CMD-3، مما يشير إلى أن بيانات CMD-3 قد تكون قيمة شاذة أو أن نتائج الشبكة قوية ضد التشتت التجريبي.

5. الخاتمة

يمثل هذا العمل علامة فارقة في فيزياء العزم المغناطيسي للميون من خلال تقديم أول تحديد للرتبة التالية (NLO) لـ HVP بدقة تقل عن واحد في المئة باستخدام الشبكة. يعتمد نجاح هذا العمل بشكل كبير على الإلغاء بين نواتي NLOa و NLOb، مما يجعل الحساب بطبيعته أقل حساسية لمنهجيات المسافة الطويلة مقارنة بحسابات الرتبة الرئيسية (LO). تعيد هذه النتيجة التحدي للتقييمات القائمة على البيانات التي لا تدمج قياسات CMD-3، وتعزز فرضية وجود تناقض محتمل بين النموذج القياسي والبيانات التجريبية.