A Study of Entanglement and Ansatz Expressivity for the Transverse-Field Ising Model using Variational Quantum Eigensolver

إليك ملخص تقني مفصل للورقة البحثية بعنوان "دراسة حول التشابك وتعبيرية الـ Ansatz لنموذج إيزينج في المجال المستعرض باستخدام خوارزمية القيم الذاتية الكمومية التباينية (VQE)" من إعداد أشوتوش ب. تريباثي، ونيلماني ماثور، وفيكرام تريباثي.

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

تعد خوارزمية القيم الذاتية الكمومية التباينية (VQE) خوارزمية هجينة رائدة بين الكم والكلاسيكية لمحاكاة الأنظمة متعددة الأجسام في عصر التقنيات الكمومية متوسطة الحجم ذات الضجيج (NISQ). ومع ذلك، تعتمد فعاليتها بشكل كبير على عاملين:

1. تعبيرية الـ Ansatz (Ansatz Expressivity): قدرة الدائرة الكمومية ذات المعلمات (PQC) على تمثيل الحالة الأرضية الحقيقية، لا سيما في الأنظمة التي تتميز بتشابك قوي أو تدهور (degeneracy).
2. استقرار التحسين (Optimization Stability): قدرة المحسن الكلاسيكي على التنقل في مشهد الخسارة (loss landscape) لإيجاد المعلمات المثلى.

يعتبر نموذج إيزينج في المجال المستعرض (TFIM) منصة اختبار مثالية لهذه المشكلة لأن حالته الأرضية تمر بتحول طوري كمي. في نظام المجال المستعرض المنخفض ($h_x \ll h_c$)، تكون الحالة الأرضية عالية التشابك (تراكب للحالات $|00\dots0\rangle$ و $|11\dots1\rangle$)، بينما في نظام المجال العالي ($h_x \gg h_c$)، تكون حالة ناتج ضرب (product state) ذات تشابك منخفض. يهدف المؤلفون إلى اختبار بنيات (architectures) مختلفة للـ ansatz لتحديد مدى قدرتها على التقاط خصائص التشابك والسلوكيات الحرجة هذه عبر شبكات أحادية وثنائية وثلاثية الأبعاد.

2. المنهجية

النموذج والإعداد

• الهاملتوني (Hamiltonian): يتم تعريف هاملتوني TFIM كالتالي: $H = J_z \sum_{\langle ij \rangle} \sigma^z_i \sigma^z_j + h_x \sum_i \sigma^x_i$، حيث $J_z = -1.0$ مع شروط حدودية دورية.
• أحجام الأنظمة: أُجريت عمليات المحاكاة على أنظمة تصل إلى 27 كيوبت (تحديداً 10 مواقع في 1D، و $4\times4$ في 2D، وأحجام مختلفة في 3D).
• إطار المحاكاة: أُجريت محاكاة VQE دقيقة باستخدام محاكي NVIDIA CUDA-Q على وحدات معالجة الرسومات (GPUs)، مما سمح بمحاكاة متجهة الحالة (state-vector) بدقة لمقارنتها بالتشخيص الدقيق (exact diagonalization).

بنيات الـ Ansatz

تمت مقارنة ثلاث فئات متميزة من الـ ansatz:

1. الـ Ansatz الفعال من الناحية العتادية (HEA): وتحديداً EfficientSU2 من Qiskit. يستخدم بوابات عتادية أصلية ($R_z, R_x, CNOT$) وهو مصمم للتعبيرية العالية.
2. الـ Ansatz المتغير الهاملتوني (HVA): يتم بناؤه عن طريق عملية "Trotterization" لمصطلحات TFIM. يستخدم بوابات $R_{ZZ}$ للارتباط وبوابات $R_x$ للمجال المستعرض. هو مستوحى فيزيائياً ولكنه مقيد بفضاء فرعي محدد.
3. الـ HVA مع كسر التماثل (HVA-SB): امتداد لـ HVA يضيف طبقة $R_z$ لكسر تماثل التكافؤ $Z_2$ صراحةً، مما يسمح للدائرة بالوصول إلى حالات هي عبارة عن تراكبات لقطاعات التكافؤ.

استراتيجيات التحسين

• HEA: تم تحسينه باستخدام L-BFGS نظراً لمشهد التحسين السلس لديه.
• HVA و HVA-SB: تم تحسينهما باستخدام COBYLA (الخالي من المشتقات) لأن مشاهدهما أكثر وعورة وعرضة للحد الأدنى المحلي.
• التمدد ثلاثي الأبعاد: استُخدم ansatz "السعة الحقيقية" (Real Amplitude) -عن طريق إزالة دورات $R_z$ للتقيد بالسعات الحقيقية- لتقليل المعلمات والحفاظ على مشهد سلس.

مقاييس التقييم

• التعبيرية (Expressivity): تُقاس عبر جهد الإطار (Frame Potential) ($F_t$)، حيث تشير القيم الأقل إلى تعبيرية أعلى (أقرب إلى t-design وحدوي).
• دقة الطاقة: تُقاس عبر تباين الطاقة (Energy Variance) ($Var(E) = \frac{\langle H^2 \rangle - \langle H \rangle^2}{\langle H \rangle^2}$).
• الملاحظات الفيزيائية:
  • المغنطة (Magnetization): $M = \frac{1}{N}\sum \langle \sigma^z_i \rangle$.
  • ارتباط السبين (Spin Correlation): $M_{Corr} = \frac{1}{N/2}\sum \langle \sigma^z_i \sigma^z_{i+N/2} \rangle$ (يُستخدم لاستقصاء الحرجة في الأنظمة محدودة الحجم).
  • إنتروبيا التشابك (Entanglement Entropy): إنتروبيا فون نيومان ($S_{EE}$) المحسوبة عبر مصفوفات الكثافة المختزلة.

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

1. التقييم المنهجي للتعبيرية مقابل التحسين: تكمكن الدراسة من قياس المقايضة بين تعبيرية الـ ansatz (مثل HEA) واستقرار مشهد التحسين الخاص بها (مثل HVA).
2. دور كسر التماثل: توضح الورقة أن إضافة طبقات كسر التماثل (HVA-SB) أمر بالغ الأهمية للوصول إلى الحالة الأرضية الصحيحة في الأنظمة محدودة الحجم، حيث تكون الحالة الأرضية الحقيقية عبارة عن تراكب لحالات مكسورة التماثل ومتدهورة.
3. تدرج الأبعاد: يوسع العمل نماذج VQE من 1D إلى 2D و3D، مسلطاً الضوء على كيفية تسبب زيادة الاتصال والتشابك في الأبعاد الأعلى في تفاقم صعوبات التحسين.
4. التقليل من تقدير التشابك: حدد المؤلفون قصوراً منهجياً حيث تميل خوارزمية VQE، خاصة مع الدوائر الضحلة، إلى التقليل من تقدير إنتروبيا التشابك.

4. النتائج

النظام أحادي البعد (10 مواقع)

• تباين الطاقة: أظهر HEA تحسناً تدريجياً مع زيادة عمق الدائرة. في المقابل، أظهر HVA و HVA-SB انتقالاً حاداً من التقدير الضعيف إلى الدقيق للطاقة بمجرد الوصول إلى العمق الكافي. هذا يشير إلى أن HVA يستفيد من كونه مقيداً بالفضاء الفرعي للهاملتوني.
• التشابك والارتباطات:
  • أنتج HVA إنتروبيا تشابك موقع واحد عالية ($S \approx 1$) في نظام المجال المنخفض لأنه فشل في كسر التماكل تلقائياً في الأنظمة محدودة الحجم.
  • نجح كل من HEA و HVA-SB في التقاط تغيرات الانحناء النوعية بالقرب من النقطة الحرجة، رغم أن HEA مال إلى اختيار الحالات المكسورة التماثل (تشابك منخفض) وهو أمر صحيح فيزيائياً فقط في الحد الديناميكي الحراري.

النظام ثنائي الأبعاد (شبكة $4\times4$)

• عدم استقرار التحسين: أصبح التحسين غير مستقر بشكل كبير بسبب زيادة الاتصال.
• الأداء: كان أداء HVA ضعيفاً في الأنظمة ذات التشابك المنخفض. قلل HEA من تقدير الإنتروبيا في المناطق عالية التشابك في المجال المنخفض، مائلاً مرة أخرى نحو الحالات المكسورة التماثل.
• الحرجة (Criticality): نجح كل من ارتباط السبين وإنتروبيا التشابك في الإشارة إلى التحول الطوري بالقرب من المجال الحرج $h_x \approx 3.0$.

النظام ثلاثي الأبعاد

• التحديات: اشتدت تحديات التحسين بشكل أكبر.
• الحل: استُخدم HEA معدل (السعة الحقيقية) مع تهيئة المعلمات من نقاط المجال المستعرض القريبة، مما سمح بدقة معقولة لطاقة الحالة الأرضية.
• النتائج: ظلت طاقة الحالة الأرضية لكل موقع مستقلة عن حجم الشبكة، بينما استمرت إنتروبيا التشابك في إظهار علامات واضحة للسلوك الحرج.

5. الأهمية والخاتمة

تخلص الورقة إلى أن التعبيرية واستقرار التحسين هما متطلبان متنافسان في VQE:

• يوفر HEA تعبيرية عالية ولكنه يعاني من مشاهد تحسين غير منتظمة ويمكن أن يعلق في الحد الأدنى المحلي أو يختار قطاعات تماثل غير صحيحة في الأنظمة محدودة الحجم.
• يوفر HVA مساراً أكثر سلاسة نحو الحالة الأرضية ضمن فضاءه الفرعي المقيد، ولكنه يفتقر إلى المرونة لالتقاط الحالات خارج ذلك الفضاء دون تعديلات كسر التماثل.
• يعمل HVA-SB على سد هذه الفجوة من خلال السماح بالوصول إلى حالات انتهاك التكافؤ، مما يحسن تقدير الطاقة.

الآثار المترتبة:

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

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