Amplitude Encoding of Slater-Type Orbitals via Matrix Product States: Efficient State Preparation and Integral Evaluation on Quantum Hardware

ملخص تقني لورقة بحثية بعنوان: "ترميز سعة مدارات من نوع سلاتر عبر حالات نواتج المصفوفات: إعداد الحالة الفعال وتقييم التكامل على الأجهزة الكمومية"

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

في الكيمياء الكمومية الحسابية، تُعد مدارات من نوع سلاتر (STOs) الوصف الفيزيائي الصحيح للدوال الموجية الذرية، حيث تحقق شرط "النتوء النووي" (nuclear cusp condition) وتُظهر اضمحلالاً أسياً صحيحاً. ومع ذلك، تم استبدالها إلى حد كبير بـ مدارات من نوع غاوس (GTOs) لأن التكاملات متعددة المراكز فوق مدارات GTOs لها حلول مغلقة، بينما تتطلب مدارات STO تكاملاً عددياً مكلفاً.

تعتمد النهج الحالية للكيمياء الكمومية على الحواسيب الكمومية غالباً على مدارات GTOs أو تعاني من أخطاء منهجية في مجموعة القواعد الأساسية عند تقريب مدارات STO. تعالج هذه الورقة تحدي ترميز مدارات STO مباشرة في الحالات الكمومية لتمكين مجموعات قواعد STO دقيقة. العقبة الرئيسية هي كفاءة إعداد الحالة: حيث تتطلب الحالات العامة المكونة من $n$ من الكيوبتات $O(2^n)$ من البوابات، مما يلغي الميزة الكمومية. الهدف هو تحديد ما إذا كان يمكن إعداد مدارات STO بكفاءة ($O(\text{poly}(n))$) باستخدام حالات نواتج المصفوفات (MPS) وتقييم تكاملات الإلكترون الواحد (التداخل، الطاقة الحركية، والجذب النووي) مباشرة على الأجهزة الكمومية.

2. المنهجية

يستخدم المؤلفون ترميز السعة (Amplitude Encoding)، حيث يتم تخزين دالة $f(x)$ على شبكة مكونة من $N=2^n$ من النقاط كـ "سعات" لحالة كمومية مكونة من $n$ من الكيوبتات:
$$|f\rangle = \frac{1}{\mathcal{N}} \sum_{j=0}^{N-1} f(x_j) |j\rangle$$

لضمان الإعداد الفعال، يتم تفكيك الحالة إلى حالة نواتج مصفوفات (MPS). ويحكم كفاءة ذلك بُعد الرابط (bond dimension $\chi$)، الذي يحدد مقدار التشابك عبر التقسيمات الثنائية للكيوبتات.

• البناء التحليلي: بالنسبة للدوال أحادية البعد من الشكل $P_d(x)e^{-\zeta x}$ (متعدد حدود $\times$ أسي)، اشتق المؤلفون موتّرات MPS تحليلية مباشرة من معاملات المدار، وهذا يتجنب أخذ عينات الشبكة ويوفر بُعد رابط ثابت $\chi = d + 1$.
• البناء العددي: بالنسبة لإحداثيات ديكارتية ثلاثية الأبعاد والدوال الموزونة بالجهد ($V(x)\psi(x)$)، حيث يكون التفكيك التحليلي مستحيلاً، استخدم المؤلفون تفكيك القيم المفردة العددي (SVD) لبناء الـ MPS.
• تقييم التكاملات:
  • التداخل (Overlap): يُحسب عبر طريقة الحساب/إلغاء الحساب ($\langle 0| U_A^\dagger U_B |0\rangle$).
  • الطاقة الحركية: تُختزل إلى تداخل حالات مشتقة باستخدام التكامل بالأجزاء ($\langle \partial_x \psi_A | \partial_x \psi_B \rangle$).
  • الجذب النووي: يتم ترميزه عن طريق إعداد حالة جديدة للدالة الناتجة عن الضرب $\phi(x) = V(x)\psi(x)$ وحساب تداخلها مع المدار الأصلي.

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

1. بناءات MPS التحليلية: توفر الورقة موترات MPS ذات صيغة مغلقة للمدارات أحادية البعد (1s, 2s، والمشتقات) مع بُعد رابط ثابت $\chi = d+1$. وهذا يتيح توسع الموارد الكلاسيكية والكمومية بمقدار $O(n)$.
2. مسار كامل للإلكترون الواحد: تم عرض مسار كامل لتكاملات التداخل، والطاقة الحركية، والجذب النووي في بُعد واحد وتم التحقق من صحته على الأجهزة.
3. التوسع ثلاثي الأبعاد وتحليل التشابك: توسع الدراسة إلى الإحداثيات الديكارتية ثلاثية الأبعاد. ومن الأهمية بمكان أن الورقة كشفت أن بُعد الرابط لمدارات STO ثلاثية الأبعاد يتشبع بدلاً من النمو بشكل أسي مع دقة الشبكة.
4. التحقق من الأجهزة: تم تنفيذ ناجح لتكاملات التداخل والطاقة الحركية على معالجات IBM Heron (5 كيوبت)، مع تحقيق خطأ أقل من 1% باستخدام استخراج الضجيج الصفري (ZNE).
5. تحسين الموارد: تحديد عتبة قطع SVD كمعامل عملي لتقليل بُعد الرابط (وبالتالي تقليل عمق الدائرة) مع تأثير ضئيل على الدقة.

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

أ. نتائج البعد الواحد

• أبعاد الرابط:
  • مدار 1s ($e^{-\zeta|x-R|}$): $\chi = 2$.
  • مدار 2s ($x e^{-\zeta x}$): $\chi = 2$.
  • هيدروجين 2s ($(2-r/a)e^{-r/2a}$): $\chi = 3$.
  • ناتج الجذب النووي ($V \cdot \psi$): $\chi = 11$ (يتشبع).
• الدقة:
  • تقترب تكاملات التداخل والطاقة الحركية من دقة الآلة مع زيادة عدد الكيوبتات (خطأ التجزئة $<0.01\%$ عند 10 كيوبت).
  • أداء الأجهزة: على جهاز IBM Kingston (5 كيوبت)، حقق تكامل التداخل خطأً ناتجاً عن الجهاز بنسبة 0.67% باستخدام ZNE. أظهرت الطاقة الحركية خطأً نسبياً أعلى (9.6%) بسبب صغر حجم الإشارة ($P(0) \approx 0.004$) التي تتنافس مع مستوى ضجيج الأجهزة.

ب. نتائج البعد الثلاثي

• الترميز الديكارتي: تم حساب تكاملات التداخل متعددة المراكز بين مدارات 1s و 2s.
  • عند 6 كيوبت لكل إحداثي (18 كيوبت إجمالاً)، كان خطأ التجزئة لتداخل 1s-1s هو 0.02%.
  • تم التحقق من التعامد بين مدارات 1s و 2s ببواقي متسقة مع دقة الشبكة.
• تشبع التشابك:
  • تشبع أقصى بُعد رابط ($\chi_{max}$) لمدار الهيدروجين 1s في الإحداثيات الديكارتية عند $\sim 138$ (العتبة $10^{-12}$) مع زيادة دقة الشبكة إلى 12 كيوبت لكل إحداثي.
  • تشابك الإحداثيات المتقاطعة ($\chi_{x|y}$) يتشبع عند $\sim 41$.
  • الأهمية: يثبت هذا أن الترميز السعة في الإحداثيات الديكارتية ثلاثية الأبعاد له تعقيد محدود (مكلف ولكنه ليس مستعصياً)، خلافاً للمخاوف الأولية من التوسع الأسي.
• تأثير القطع: تقليل عتبة SVD من $10^{-12}$ إلى $10^{-6}$ يقلل $\chi_{max}$ من 138 إلى 39 (عامل قدره 3.4) مع تأثير ضئيل على دقة التكامل.

ج. مقاييس الدائرة

• دوائر البعد الواحد: فعالة، تتطلب حوالي 100-250 من بوابات الكيوبت الثنائي لـ 5-8 كيوبت.
• دوائر البعد الثلاثي: عميقة جداً حالياً بالنسبة لأجهزة NISQ. إعداد 1s ثلاثي الأبعاد يتطلب 3 كيوبت/إحداثي حوالي 9,500 بوابة ECR؛ وعند 4 كيوبت/إحداثي، يتجاوز 220,000 بوابة. تم التحقق من هذه الدوائر في المحاكاة فقط.

5. الأهمية والآثار المترتبة

• جدوى مدارات STO الدقيقة: تثبت هذه الدراسة أن ترميز السعة القائم على MPS هو مسار قابل للتطبيق لاستخدام مدارات من نوع سلاتر دقيقة على الحواسيب الكمومية، متجاوزةً تقريب غاوس.
• التعقيد المحدود: إن اكتشاف تشبع أبعاد الرابط في الفضاء ثلاثي الأبعاد هو اكتشاف نظري حاسم. وهذا يعني أنه بينما يعد الترميز ثلاثي الأبعاد كثيف الموارد، إلا أنه ليس صعباً بشكل أسي، مما يجعله قابلاً للتنفيذ في الحواسيب الكمومية المستقبلية القادرة على تصحيح الأخطاء.
• الإدارة العملية للموارد: يسمح تحديد عتبة SVD كمعامل قابل للضبط للباحثين بالمقايضة بين دقة طفيفة وتقليل كبير في عمق الدائرة (عدد الكيوبتات وبوابات الدائرة).
• جاهزية الأجهزة: يؤكد نجاح العرض التجريبي المكون من 5 كيوبت صحة طريقة الحساب/إلغاء الحساب لتقييم التكاملات، رغم أن نسب الإشارة إلى الضجيج لا تزال تشكل تحدياً للعناصر المصفوفية الصغيرة (مثل الطاقة الحركية) على أجهزة NISQ الحالية.
• المسار المستقبلي: يضع هذا الإطار المسرح لحساب تكاملات الجذب النووي ثلاثية الأبعاد، والتوسعات متعددة المراكز (عبر قنوات الزخم الزاوي)، وفي النهاية تكاملات الإلكترونين (عبر توسيع متعدد الأقطاب)، بشرط بقاء بُعد الرابط تحت السيطرة.

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