Extracting average properties of disordered spin chains with translationally invariant tensor networks

ملخص تقني: استخراج الخصائص المتوسطة لسلاسل الغزل المضطربة باستخدام الشبكات التنسورية المترابطة تباعيًا

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

المنهجية
يقترح المؤلفون خوارزمية تعتمد على الشبكة التنسورية تستغل الثبات الإحصائي الانتقالي—وهو الشرط الذي يكون فيه توزيع متغيرات الاضطراب متطابقاً عبر جميع مواقع الشبكة. وتتمثل الاستراتيجية الجوهرية في تمثيل مصفوفة الكثافة لمتوسط الاضطراب، $\rho = \sum_{\{R_n\}} [\prod_n P(R_n)] \rho_G[\{R_n\}]$، كـ "مؤثر منتج مصفوفي" (MPO) واحد وثابت انتقالياً.

تتبع الخوارزمية الخطوات التالية:

1. إدخال الـ "أنسيلا" (Ancilla): للتعامل مع الاضطراب، يتم إدخال "كوديتات" أنسيلا $|R_n\rangle$ ذات بُعد $N_D$ (عدد قيم الاضطراب الممكنة) على كل موقع. يتم بناء هاميلتوني ثابت انتقالياً $H = \sum_n h_n$، حيث تعمل الحدود المحلية $h_n$ كعمليات محكومة: $h_n = \sum_{R_n} \tilde{h}_n[R_n] \otimes |R_n\rangle\langle R_n|$. تعمل هذه الـ "أنسيلا" كأدوات تحكم، تختار حد الهاميلتوني المحلي المحدد بناءً على قيمة الاضطراب.
2. التطور الزمني التخيلي المعدل: تبني الخوارزمية حالة "غيبس" عبر تطور زمني تخيلي باستخدام نهج "تطوير بلوك لزمن التطور" (TEBD) معدل. بدلاً من مجرد تطوير $e^{-\tau H}$، تقوم الطريقة بتطوير مؤثر غير معير $\tilde{\rho}(\tau) = N(\tau)e^{-\tau H}$. حيث $N(\tau)$ هو مؤثر قطري يعمل على الـ "أنسيلا" ويحتوي على مقلوبات دالات التجزئة $1/Z[\{R_n\}]$ لكل تكوين اضطراب. يضمن هذا أن تتبع (trace) دوران (spins) الفيزياء يؤدي إلى مصفوفة الوحدة على كوديتات الاضطراب، مما يحافظ على المعايرة الصحيحة لمزيج حالات "غيبس".
3. المعايرة والقلب (Inversion): تتضمن الخطوة الحرجة تحديث مؤثر المعايرة $N(\tau)$ عند كل خطوة زمنية. يتطلب هذا قلب MPO يُرمز له بـ $\Lambda$ مشتق من أثر (trace) الـ "spins" الفيزيائية. وبما أن $\Lambda$ تكون قريبة من مصفوفة الوحدة عند الخطوات الزمنية الصغيرة، يستخدم المؤلفون تحسيناً تباينياً (تعميماً لخوارهازمية VOMPS) لإيجاد تقريب MPO ذي بُعد رابط منخفض $\Lambda_M^{-1}$ يعظم التطابق (fidelity) مع مؤثر الوحدة.
4. القطع (Truncation): بعد تحديث المعايرة، يتم قطع الـ MPO مرة أخرى إلى بُعد رابط مستهدف $D$ باستخدام قطع قيم "شميدت" القياسي. تتكرر دورة التطور، وتعديل المعايرة، والقطع حتى الوصول إلى درجة الحرارة العكسية المستهدفة $\beta$. ويتم الحصول على متوسط كثافة الاضطراب النهائي عن طريق تتبع (tracing out) كوديتات الـ "أنسيلا".

النتائج الرئيسية
تم اختبار الطريقة على نموذج إيزينج للمجال المستعرض العشوائي (RTFIM) بالقرب من نقطة الحرجة ذات العشوائية اللانهائية ($\delta=0$).

• الكفاءة: نجحت الخوارما في إعادة إنتاج الخصائص المتوسطة لنموذج RTFIM باستخدام أبعاد رابطة صغيرة نسبياً (رتبة $\sim 100$) وصولاً إلى درجات حرارة منخفضة.
• التحجيم الحرج: يُظهر طول الارتباط $\xi$ المستخرج "التحجيم الديناميكي المنشط" المتوقع نظرياً $\xi \sim (\ln \beta)^2$ والذي يميز نقطة الثبات ذات العشوائية اللانهائية. أظهرت النتائج تقارباً فيما يتعلق بالبُعد الرابط وحجم الخطوة الزمنية.
• توزيع الاضطراب: من خلال تغيير عدد قيم الاضطراب المنفصلة $N_D$، أكد المؤلفون أن تحجيم $\xi$ يعتمد على تباين توزيع الاضطراب، بما يتوافق مع التوقعات النظرية.
• المناطق النادرة: من خلال تحليل مصفوفات النقل لتكوينات اضطراب فردية (دون تتبع الـ "أنسيلا")، استخرج المؤلفون توزيع أطوال الارتباط. وقد لاحظوا وجود "ذيل ثقيل" في التوزيع، مما يؤكد وجود المناطق النادرة. وُجد أن طول الارتباط النموذجي ($\xi_{typ} \approx 35.4$) أصغر بكثير من متوسط طول الارتباط ($\xi_{av} \approx 47.5$)، وهو ما يتفق مع التوقعات النظرية لفيزياء العشوائية اللانهائية.
• سلوك الطور: أكدت نتائج إضافية النمو الجبري لطول الارتباط $\xi \sim 1/T^{\alpha(\delta)}$ على الجانب المرتب من الانتقال، وهو ما يختلف عن النمو الأسي في الأنظمة النقية.

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

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