Point-group symmetry analysis of many-electron wavefunctions on a quantum computer

ملخص تقني: تحليل تناظر زمرة النقطة للدوال الموجية متعددة الإلكترونات على حاسوب كمي

بيان المشكلة
تعد زمر النقطة (Point groups)، التي تمثل عمليات التناظر المكاني في الأنظمة الجزيئية، أدوات أساسية في الكيمياء لتحليل المدارات الجزيئية والتحليل الطيفي. وبينما اقتُرحت خوارزميات كمية لاستغلال هذه التناظرات — مثل تقليل عدد الكيوبتات أو تبسيط دوائر الـ (ansatz) — إلا أن التطبيقات العملية لعمليات زمرة النقطة والتحليلات التفصيلية للتناظر للدوال الموجية الواقعية متعددة الإلكترونات لا تزال غير مستكشفة إلى حد كبير. وتحديداً، هناك نقص في الطرق التي يمكنها تنفيذ عمليات زمرة النقطة بشكل ملموس وتحليل خصائص التناظر (مثل أوزان التمثيلات غير القابلة للاختزال) لحالات واقعية متعددة الإلكترونات ومرتبطة (correlated) على أجهزة كمية كمية قريبة من المدى (near-term).

المنهجية
يقترح المؤلفون طريقة هجينة (كمية-كلاسيكية) لا تتطلب كيوبتات مساعدة (ancilla-free) لتحليل تناظرات زمرة النقطة للحالات متعددة الإلكترونات، وهي قابلة للتطبيق على الزمر الآبلية (abelian) وغير الآبلية (non-abelian). يعتمد جوهر الطريقة على حساب "الوزن" ($w_\Gamma$) للتمثيلات غير القابلة للاختزال ($\Gamma$) داخل دالة موجية معينة $|\Psi\rangle$. ويُعرف هذا الوزن بأنه مربع معامل الإسقاط لـ $|\Psi\rangle$ على الفضاء الجزئي لـ $\Gamma$.

يعتمد الإطار النظري على صيغة مؤثر الإسقاط:
$$w_\Gamma = \frac{d_\Gamma}{|G|} \sum_C \chi^*_\Gamma(C) \sum_{g \in C} \langle\Psi|\hat{g}|\Psi\rangle$$
حيث $d_\Gamma$ هو بُعد التمثيل، و $|G|$ هو رتبة الزمرة، و $\chi_\Gamma(C)$ هو رمز الشخصية (character) للفئة المرافقة $C$، و $\hat{g}$ هي عملية التناظر.

لتقييم $\langle\Psi|\hat{g}|\Psi\rangle$ على حاسوب كمي دون استخدام كيوبتات مساعدة، يستخدم المؤلفون عملية قطرية لمصفوفات التمثيل. بالنسبة للمدارات المتكيفة مع التناظر (مثل هارتري-فوك)، يتم تفكيك عملية التناظر $U(g)$ إلى دورات مدارية. ومن خلال قطرية مصفوفة التمثيل $\bar{D}_\gamma(g)$ لمجموعات المدارات المتعددة، يتم إعادة كتابة الوحدة $U(g)$ كحاصل ضرب دوران قاعدة $\tilde{U}(g)$ ومؤثر طور قطري. يتم الحصول على القيمة المتوقعة عن طريق تحضير الحالة $|\tilde{\Psi}(g)\rangle = \tilde{U}(g)|\Psi\rangle$، وقياسها في القاعدة الحسابية، ودمج النتائج مع عوامل الطور المناسبة. يتجنب هذا النهج الدوائر العميقة ومتطلبات الكيوبتات المساعدة المستخدمة في اختبارات هادامارد (Hadamard tests) أو تقنية الجمع الخطي للوحدات (LCU).

تستخدم الدراسة استراتيجيتين رئيسيتين للتحقق من الصحة:

1. المحاكاة العددية الكلاسيكية: باستخدام مكتبات ffsim و pyscf ، تم تطبيق الطريقة على البنزين ($D_{6h}$) والفروسين المتدرج ($D_{5d}$). يحلل المؤلفون محددات سلايستر (Slater determinants) الفردية (الحالة الأرضية لهارتري-فوك والحالات المثارة المفردة) والدوال الموجية المرتبطة (دوال Unitary Cluster Jastrow و Local Unitary Cluster Jastrow).
2. العرض التجريبي على الأجهزة: تم تنفيذ الطريقة على جهاز ibm_kawasaki فائق التوصيل التابع لشركة IBM (يصل إلى 32 كيوبت). النظام المستهدف هو البنزين ذو تناظر $D_{2h}$ (وهو زمرة فرعية من $D_{6h}$). تم تحضير الدوال الموجية التجريبية (الحالة الأرضية والحالة المثارة الأولى) باستخدام حسابات كثافة مصفوفة التناقص (DMRG)، وضغطها في قالب "جدار الطوب" (brick-wall ansatz) عبر تقنيات الشبكة التنسورية، ثم رسم خرائطها إلى دوائر كمية. ولمعالجة الضجيج، تمت مقارنة تقنيتين لتخفيف الخطأ: الاستقراء صفر الضجيج (ZNE) عبر طي البوابة (gate-folding)، ونهج التخفيف شبه الاحتمالي (QESEM).

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

• المحاكاة الكلاسيكية:
  • بالنسبة للحالات الأرضية لهارتري-فوك للبنزين والفروسين، حددت الطريقة بشكل صحيح وزن التمثيل غير القابل للاختزال المتماثل كلياً ($A_{1g}$) ليكون قريباً من 1.0، وهو ما يتوافق مع التكوينات مغلقة الغلاف.
  • بالنسبة للحالات المثارة المفردة (SE)، نجحت الطريقة في تفكيك الدوال الموجية إلى تمثيلات متعددة غير قابلة للاختزال. على سبيل المثال، تم تفكيك حالة SE معينة في البنزين إلى $\approx 0.005 B_{1u} + 0.495 B_{2u} + 0.5 E_{1u}$. وقد أعطى جمع الأوزان عبر جميع تكوينات SE المستقلة نتائج ثابتة متوافقة مع صيغ اختزال نظرية الزمر.
  • في تحليل دوال LUCJ المرتبطة، كشفت الطريقة أن إجراءات التحسين التي تقلل خطأ الطاقة (بدون تنظيم/regularization) يمكن أن تؤدي بشكل كبير إلى تدهور وزن التناظر (الانحراف عن وزن $A_{1g}$ الصحيح)، بينما ساعد التنظيم في الحفاظ على خصائص التناظر، وإن كان مع مقايضات مختلفة في الطاقة. يشير هذا إلى أن المقياس مفيد لتقييم جودة الدالة الموجية التجريبية.
• العرض على الأجهزة:
  • على جهاز ibm_kawasaki ، أظهرت القياسات الخام تدهوراً كبيراً في أوزان التناظر مع زيادة عدد الكيوبتات (على سبيل المثال، انخفض وزن $A_g$ للحالة الأرضية إلى ~0.54 لـ 32 كيوبت بدون تخفيف الخطأ).
  • كان تطبيق تخفيف الخطأ أمراً حاسماً. وبينما حسن ZNE النتائج، حقق نهج QESEM شبه الاحتمالي أعلى دقة، حيث أعاد إنتاج أوزان التمثيلات غير القابلة للاختزال ضمن خطأ بنسبة مئوية قليلة (على سبيل المثال، $0.968 \pm 0.010$ للحالة الأرضية على 32 كيوبت، مقارنة بالقيمة الدقيقة 0.999).
  • أكدت التجربة أن الطريقة المقترحة الخالية من الكيوبتات المساعدة، جنباً إلى جنب مع تحضير حالة الشبكة التنسورية وتخفيف الخطأ المتقدم، يمكنها إعادة إنتاج خصائص التناظر بدقة على الأجهزة الحالية.

الأهمية والادعاءات
يزعم البحث أنه يسد الفجوة بين الإسقاط النظري للتناظر والتنفيذ العملي على الأجهزة الكمية. وتعتبر الطالة المقترحة أداة عملية لـ:

1. تحليل التناظر: تحليل خصائص التناظر للدوال الموجية متعددة الإلكترونات في محاكاة المواد الواقعية على الحواسيب الكمية في المدى القريب والمدى البعيد (fault-tolerant).
2. تقييم الدالة الموجية: تقييم جودة الدوال الموجية التجريبية (مثل الخوارزميات التباينية) من خلال قياس مدى التزامها بتناظرات محددة.
3. الاختبار المعياري (Benchmarking): توفير معيار للأجهزة الكمية الحقيقية وتقنيات تخفيف الخطأ، حيث لا تتطلب الطريقة كيوبتات مساعدة إضافية وتتطلب حداً أدنى من العبء الحسابي للمدارات المتكيفة مع التناظر.

يؤكد المؤلفون أنه بينما استخدم العرض التجريبي على الأجهزة زمرة فرعية آبلية ($D_{2h}$) وتقييماً قائماً على "باولي"، فإن المنهجية الأساسية عامة وقابلة للتطبيق على الزمر غير الآبلية ودوال الأساس التعسفية. ويشير العمل إلى أن تحليل التناظر هو مكون قابل للتطبيق وضروري لمحاكاة المواد المعقدة كمياً في المستقبل.