Extracting conserved operators from a projected entangled pair state

إليك ملخص تقني مفصل للورقة البحثية بعنوان "استخراج المشغلات المحفوظة من حالة الأزواج المتشابكة المسقطة" (Extracting conserved operators from a projected entangled pair state) للباحثين Xu و Frías Pérez و Yang.

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

تتناول الورقة البحثية المشكلة العكسية لـ تعلم الهاملتوني (Hamiltonian learning): بالنظر إلى حالة كمومية متعددة الأجسام (تحديداً حالة شبكة تنسور مثل حالة الأزواج المتشابكة المسقطة - PEPS)، كيف يمكن إعادة بناء الهاملتوني الأساسي أو غيره من المشغلات المحفوظة التي تكون هذه الحالة حالة ذاتية لها؟

• السياق: بينما توجد خوارزميات فعالة لأنظمة (MPS) أحادية البعد، فإن التوسع إلى أنظمة (PEPS) ثنائية الأبعاد أصعب بكثير بسبب الطبيعة المعقدة (P-sharp complete) لعملية اختزال (contraction) الـ PEPS الدقيقة.
• التحدي: تعتمد الحسابات العملية في 2D على مخططات اختزال تقريبية (مثل Corner Transfer Matrix Renormalization Group - CTMRG)، والتي تُدخل أخطاءً في دالات الارتباط. كان من غير الواضح سابقاً ما إذا كان بإمكاننا استخراج مشغلات محفوظة بشكل موثوق من هذه الحالات التقريبية، خاصة بالنسبة للهاملتونيات غير الخالية من الإحباط (non-frustration-free) أو الأنظمة الحرجة (gapless).
• الهدف: تطوير طريقة لاستخراج مشغلات محفوظة ذات محلية هندسية $k$ (بما في ذلك الهاملتونيات) من حالة iPEPS لانهائية، حتى عندما تكون الحالة مجرد حالة ذاتية تقريبية أو ناتجة عن تحسين تبايني (variational optimization).

2. المنهجية

يقترح المؤلفون طريقة تعتمد على عوامل البنية الساكنة (static structure factors) للمشغلات متعددة المواقع. المنطق الجوهري هو أنه إذا كانت الحالة $|\Psi\rangle$ حالة ذاتية لمشغل محفوظ $\hat{H}$، فإن التقلبات الكمومية لـ $\hat{H}$ في تلك الحالة يجب أن تتلاشى (أو تكون في حدها الأدنى في حالة الحالات ذاتية التقريب).

أ. الإطار النظري

1. مصفوفة عامل البنية الساكنة ($S$):
   تُعرف الطريقة مصفوفة $S^{\alpha, \beta}(q)$ التي تمثل تحويل فورييه لدالات الارتباط المتصلة لمجموعة أساس من المشغلات الهيرميتية المحلية هندسياً $\{\hat{o}^\alpha_x\}$:
   $$S^{\alpha, \beta}(q) = \sum_x e^{-iq\cdot x} \left( \frac{\langle \Psi | \hat{o}^\alpha_x \hat{o}^\beta_0 | \Psi \rangle}{\langle \Psi | \Psi \rangle} - \frac{\langle \Psi | \hat{o}^\alpha_x | \Psi \rangle \langle \Psi | \hat{o}^\beta_0 | \Psi \rangle}{\langle \Psi | \Psi \rangle^2} \right)$$
2. استخراج النواة (Kernel Extraction):
   المشغلات المحفوظة $\hat{H} = \sum_x e^{iq\cdot x} \hat{h}_x$ تقابل النواة (null space) للجزء المتماثل من $S$. إذا كان $\hat{H}|\Psi\rangle = \lambda|\Psi\rangle$، فإن القيمة الذاتية لـ $S$ المرتبطة بالمعاملات المتجهة لـ $\hat{H}$ تكون صفراً (مما يمثل صفراً في التباين/التقلب).
3. الارتباط بالهندسة الكمومية:
   تؤسس الورقة أن مصفوفة عامل البنية الساكنة تتناسب مع تنسور المترية الكمومية (حساسية الوفاء - fidelity susceptibility) لمنطويد الحالات المشوهة بواسطة المشغلات. وتتوافق المشغلات المحفوظة مع اتجاهات تلاشي حساسية الوفاء.

ب. التنفيذ الحسابي: طريقة دالة التوليد (Generating Function Method)

لتقييم $S$ للمشغلات متعددة المواقع بكفاءة في 2D، يستخدم المؤلفون ويطورون طريقة دالة التوليد:

1. تعريف دالة التوليد:
   بدلاً من جمع السلاسل اللانهائية مباشرة، يعرّفون دالة توليد للحالة:
   $$\langle G_\alpha(\mu, q) | = \langle \Psi | \prod_x (1 + \mu e^{-iq\cdot x} \hat{o}^\alpha_x)$$
   حيث $\mu$ هو معامل قياسي.
2. الاشتقاق:
   يتم الحصول على عامل البنية الساكنة عن طريق اشتقاق القيمة المتوقعة بالنسبة لـ $\mu$ عند $\mu=0$:
   $$S^{\alpha, \beta}(q) = \left. \frac{d}{d\mu} \frac{\langle G_\alpha(\mu, q) | \hat{o}^\beta_0 | \Psi \rangle}{\langle G_\alpha(\mu, q) | \Psi \rangle} \right|_{\mu=0}$$
3. بناء شبكة التنسور:
   • ينتج عن حاصل الضرب حاصلاً لـ (iPEPO) لانهائي.
   • يوضح المؤلفون أنه بالنسبة للمشغلات متعددة المواقع، يمكن بناء هذا الـ iPEPO ببعد رابط ($D'$) يمكن التحكم فيه.
   • استراتيجية الاشتقاق: وجدوا أن الاشتقاق التلقائي (Automatic Differentiation - AD) غير مستقر للمشغلات متعددة المواقع بسبب القيم المفردة الصفرية في تنسورات iPEPO عند $\mu=0$. بدلاً من ذلك، استخدموا صيغة الفرق المحدود المركزية الخماسية النقاط لتحقيق دقة عالية.
4. التعامل مع التقريبات:
   تستخدم الطريقة CTMRG لاختزال شبكات التنسور الناتجة. أثبت المؤلفون أنه رغم أخطاء CTMRG (التي يتم التحكم فيها عبر بعد الرابط للبيئة $\chi$)، فإن أدنى القيم الذاتية لـ $S$ تظل قوية بما يكفي لتحديد المشغلات المحفوظة.

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

• التعميم إلى 2D: نجاح في توسيع تعلم الهاملتوني من 1D MPS إلى 2D iPEPS، وتجاوز تعقيد تقييم المشغلات متعددة المواقع.
• المتانة تجاه التقريب: إثبات إمكانية استخراج الهاملتونيات المحفوظة حتى عندما تكون الحالة المدخلة حالة تقريبية تباينية أو عندما يُدخل مخطط الاختزال (CTMRG) أخطاءً.
• ما وراء الخلو من الإحباط: على عكس بناءات الهاملتوني "الأب" القياسية التي غالباً ما تنتج هاملتونيات خالية من الإحباط (حيث الحالة هي الحالة الأرضية)، يمكن لهذه الطة تحديد مشغلات غير خالية من الإحباط تكون فيها الحالة حالة ذاتية مثارة.
• تحسين المحلية: الهاملتونيات المستخرجة غالباً ما تمتلك محلية أفضل (عدد أقل من المواقع المعنية) مقارنة بالبناءات القياسية.

4. النتائج الرئيسية والاختبارات المعيارية

قام المؤلفون بالتحقق من صحة الطريقة على عدة أنظمة فيزيائية متميزة:

أ. حالة 2D AKLT (اختبار معياري دقيق)

• النظام: حالة Spin-2 Affleck-Kennedy-Lieb-Tasaki على شبكة مربعة.
• النتيجة: نجحت الطريقة في استعادة مولدات تناظر SO(3) الدقيقة (موقع واحد) وحددت 81 هاملتوني "أب" غير خالٍ من الإحباط مكون من موقعين (projectors onto the fusion channel 4).
• الأهمية: تؤكد أن الطريقة تعمل بشكل مثالي مع PEPS دقيقة ذات طول ارتباط محدود.

ب. نموذج 2D XX الكمومي (اختبار تبايني)

• النظام: نموذج Spin-1/2 XX في مجال Z متدرج، تم تحسينه عبر iPEPS تبايني.
• النتيجة:
  • استعادة مولد تناظر U(1) (وهو $\sum Z_i$) حتى في الطور الفيرومغناطيسي حيث ينكسر التناظر تلقائياً.
  • استخراج الهاملتوني "الأب" غير الخالي من الإحبط $H(J)$ بدقة عالية (مطابقة المعاملات حتى المرتبة الرابعة بعد الفاصلة العشرية) رغم أن الحالة المدخلة كانت تقريباً تباينياً.
• الأهمية: توضح القابلية للتطبيق على الحالات التباينية والأنظمة ذات كسر التناظر التلقائي.

ج. حالة RVB قصيرة المدى (نظام حرج)

• النظام: حالة الرنين لروابط فالنس (Resonating Valence Bond - RVB) على شبكة مربعة (نظام فجوي/gapless، ذو طول ارتباط متباعد).
• النتيجة: استخراج هاملتوني تقريبي مكون من 4 مواقع (Plaquette local Hamiltonian) يحتوي على حدود هايزنبرغ وتفاعلات رباعية الأجسام.
  • أظهرت حالة RVB أنها حالة ذاتية تقريبية جيدة (تباين منخفض).
  • طاقة التوقع لحالة RVB في هذا الهاملتوني تكاد تكون مطابقة لطاقة الحالة الأرضية للهاملتوني نفسه.
• الأهمية: تظهر أن الطريقة تعمل مع الأنظمة الحرجة ذات أطوال الارتباط المتباعدة، حيث تجد هاملتوني "أب" محلي أكثر كفاءة من بناءات الـ 8 مواقع القياسية.

د. كود توريك/إيزينج المشوه (ندوب كمومية متعددة الأجسام - QMBS)

• النظام: عائلة من دالات موجة إيزينج المشوهة $|\Psi_{Ising}(\beta)\rangle$ وحالات كود توريك المزدوجة المشوهة.
• النتيجة: وجدت هاملتونياً تكون فيه هذه الحالات حالات ذاتية مثارة بطاقة صفر (تقع في منتصف الطيف)، وليست حالات أرضية.
• الأهمية: يحدد هذا احتمال وجود ندوب كمومية متعددة الأجسام (Quantum Many-Body Scars)، حيث توجد حالات ذات تشابك منخفض كحالات مثارة في نماذج غير متكاملة. الهاملتوني يتوافق مع معامل التشويه، مما يجعل عائلة الحالات بأكملها حالات ذاتية.

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

• الفيزياء القائمة على البيانات: توفر طريقة منهجية ومؤتمتة لاكتشاف الهاملتونيات والتناظرات مباشرة من دالات موجة (ansatzes) أو من البيانات التجريبية (عبر عوامل البنية)، مما يقلل الاعتماد على الحدس الفيزيائي.
• القابلية للتوسع: تتناسب التكلفة الحسابية حدودياً مع بعد الرابط، مما يجعلها قابلة للتطبيق على الأنظمة الكبيرة.
• الاتجاهات المستقبلية:
  • تطبيق الطريقة لاستخراج التناظرات المستمرة الناشئة عند نقاط الحرجة الكمومية منزوعة التمركز (deconfined quantum critical points).
  • تعميم الطريقة لاستخراج الهاملتونيات من حالات PEPO التي تمثل حالات غيبس (Gibbs states) عند درجات حرارة نهائية.
  • تشديد الحدود النظرية لتعقيد التعلم لحالات شبكة التنسور.

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