Normalizing-flow-based density of states for (1+1)D U(1) lattice gauge theory with a $\theta$-term

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

تواجه محاكاة نظرية الحقل الشبكي، ولا سيما طريقة مونت كارلو الهجينة (HMC)، تحديين رئيسيين في أنظمة معينة:

• التباطؤ الحرج (Critical Slowing Down): بالقرب من حد الاستمرارية، تزدัง أوقات الارتباط الذاتي بشكل كبير، مما يجعل عملية أخذ العينات غير فعالة.
• مشكلة الإشارة (The Sign Problem): تعاني النظريات التي تتضمن أفعالاً مركبة (مثل تلك التي تتضمن حد $\theta$ طوبولوجي) من مشكلة إشارة عددية، حيث يصبح مقياس تكامل المسار مركباً، مما يمنع أخذ العينات بالاعتماد على الأهمية (importance sampling) التقليدية.

يقدم منهج كثافة الحالات (Density of States - DoS) حلاً محتملاً عن طريق تحليل مسار التكامل. فبدلاً من أخذ عينات من وزن بولتزمان الكامل مباشرة، يقوم منهج DoS بحساب كثافة الحالات $\rho(c)$ لشرائح ثابتة من متغير الفعل $c$، مما يسمح بإعادة إدخال اعتماد الاقتران عبر تكامل أحادي البعد. ومع ذلك، تعتمد طرق DoS التقليدية غالباً على إعادة بناء $\rho(c)$ عن طريق تكامل مشتق لوغاريتمه، وهو ما قد يؤدي إلى تراكم الأخطاء.

2. المنهجية

يقترح المؤلفون تطبيقاً يعتمد على التدفق الطبيعي (Normalizing Flow - NF) لمنهج DoS، مع توسيع العمل السابق على نظريات الحقل القياسي ليشمل نظريات الحقل الشبكي.

• التدفقات الطبيعية المتوافقة مع القياس (Gauge-Equivariant Normalizing Flows): يستخدم المؤلفون تدفقات طبيعية (خرائط قابلة للعكس مُعلمة بشبكات عصبية) تحترم تناظر القياس الخاص بالنظرية. وتحديداً، يستخدمون شبكات تلافيفية متوافقة مع القياس وقائمة على الاقتران (بناءً على المراجع [3, 4]) لضمان أن التكوينات المولدة تحقق ثبات القياس لمجموعة U(1).
• القيد المنظم (Regularized Constraint): بما أن تعريف DoS يتضمن دالة ديراك دلتا ($\delta[c - S(\phi)]$)، والتي يستحيل أخذ عينات منها مباشرة، استبدل المؤلفون دالة ديراك دلتا بـ توزيع غاوسي ذي عرض محدود $P$:
  $$\delta(c-x) \approx \frac{1}{\mathcal{N}} e^{-\frac{P}{2}(c-x)^2}$$
  وهذا يخلق فعلاً فعالاً $S_{c,P}(\phi)$ يعتمد على قيد $c$ والمنظم $P$.
• المقدر المباشر (Direct Estimator): بخلاف الطرق التقليدية التي تقيس المشتقات، يسمح إطار عمل NF بـ التقدير المباشر لكثافة الحالات $\rho_P(c)$ عبر القيمة المتوقعة:
  $$\rho_P(c) = \langle e^{-S_{c,P}(\phi) - \log q_w(\phi)} \rangle_{\phi \sim q_w}$$
  هنا، يتم تدريب تدفق طبيعي منفصل لكل قيمة قيد $c$ لتقريب التوزيع المقيد.
• كثافة الحالات المعممة (gDoS): بالنسبة للنظرية ذات حد $\theta$، يصبح الفعل مركباً ($S = \beta S_R + i \theta S_I$). يطبق المؤلفون صيغة gDoS، حيث يتم أخذ العينات بالنسبة للجزء الحقيقي من الفعل المقيد على الجزء التخيلي (الشحنة الطوبولوجية)، ويتم استعادة الطور المركب عبر تحويل فوريه.

3. حالة الاختبار: نظرية قياس U(1) في (1+1) بُعد

اختار المؤلفون نظرية قياس U(1) في (1+1) بُعد كنموذج مرجعي لأنها:

• قابلة للحل تماماً في الاستمرارية ولها دالة تقسيم شبكية معروفة تحليلياً، مما يوفر "الحقيقة الأرضية" للتحقق من الصحة.
• تُظهر قطاعات طوبولوجية وشحنة طوبولوجية ($Q_{top}$) محددة جيداً.
• إدخال حد $\theta$ يسبب مشكلة إشارة، مما يجعلها بيئة اختبار مناسبة للخوارزميات المصممة للتعامل مع الأفعال المركبة.

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

أ. بدون حد $\theta$ (فعل ويلسون النقي)

• التحقق من الصحة: نجح منهج DoS القائم على NF في إعادة إنتاج النتائج التحليلية الدقيقة لدالة التقسيم وكثافة الحالات.
• دقة القيد: حلل المؤلفون العلاقة بين قوة القيد $c$ والقيمة المتوقعة لفعل ويلسون.
  • عند قوة القيد المنخفضة ($P=100$)، كانت الانحرافات كبيرة، خاصة بالقرب من حدود مجال الفعل.
  • عند قوة القيد العالية ($P=2000$)، تحسن الاتفاق مع النتائج الدقيقة بشكل كبير عبر المجال.
• حجم العينة الفعال (ESS): كان الـ ESS مرتفعاً (يصل إلى 72%) في المنطقة المركزية للفعل، لكنه انخفض بشكل كبير (إلى 3.5%) بالقرب من الحدود (الأذيال). يشير هذا إلى أن التدفق يواجه صعوبة في التقاط التكوينات النادرة، وهي ضرورية للتكامل الدقيق.

ب. مع حد $\theta$ (الفعل المركب)

• أخذ عينات الشحنة الطوبولوجية: نجحت الطريقة في توليد تكوينات حقل القياس عند قيم صحيحة ثابتة من الشحنة الطوبولوجية $Q_{top}$، حتى للقيم النادرة نسبياً.
• تحقيق القيد: تم تحقيق التعريف الهندسي لـ $Q_{top}$ بدقة لقوى القيد المتوسطة ($P=5, 50$).
• القيود: إن إعادة بناء دالة التقسيم كدالة لـ $\theta$ محدودة حالياً بـ قدرة التعبير (expressivity) لبنية NF. مع زيادة الاقتران $\beta$ أو عند التقييد بقيم $Q_{top}$ صحيحة كبيرة، ينخفض الـ ESS، مما يؤدي إلى تقليل دقة إعادة البناء النهائية.

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

1. التوسع إلى نظريات القياس: التوسع بنجاء في منهج DoS القائم على NF من نظريات الحقل القياسي إلى نظرية حقل شبكي (U(1))، باستخدام بنيات متوافقة مع القياس.
2. التقدير المباشر لـ DoS: إثبات جدوى التقدير المباشر لكثافة الحالات باستخدام NFs دون الحاجة لخطوة وسيطة لدمج المشتقات اللوغاريتمية.
3. أخذ عينات الشحنة الطوبلية الثابتة: إظهار أن إطار عمل gDoS المدمج مع NFs يمكنه توليد تكوينات عند شحنة طوبولوجية ثابتة، مما يتجاك بفعالية مشكلة الإشارة لقطاعات محددة.
4. التحليل المنهجي: تقديم تحليل مفصل للمقايضة بين قوة القيد ($P$) وكفاءة أخذ العينات (ESS)، مع تسليط الضوء على صعوبة أخذ عينات التكوينات النادرة (أذيال التوزيع).

6. الأهمية والآفاق المستقبلية

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

باختملال، يوضح البحث أن التدفقات الطبيعية يمكنها بفعالية إعادة بناء كثافة الحالات لنظرية قياس U(1) في (1+1) بُعد، مما يؤكد صحة المنهج مقابل النتائج الدقيقة ويظهر وعداً في التعامل مع مشكلات الإشارة، بشرما كانت بنيات الشبكة العصبية معبرة بما يكفي لالتقاط التكوينات النادرة.