Low-Depth Initial-State Preparation for Ground-State Energy Estimation of Two-Dimensional Strongly Correlated Systems
بالاسترشاد بالمحاكاة الكلاسيكية، تقترح هذه الورقة استراتيجية دارة كمومية منخفضة العمق لإعداد الحالات الأرضية لنماذج هوبارد ثنائية الأبعاد عن طريق تضمين حالات هايزنبرغ تقريبية وتطبيق بوابات تقلب الشحنة المشتقة من تحويل شريفير-ولف، مما يتيح نقل المعلمات القابل للتوسع إلى الشبكات الكبيرة بتكلفة إعداد تُقدر بنحو 105 من بوابات T.
في السعي لفهم كيفية سلوك المادة عند مستواها الأكثر جوهرية، غالباً ما يتجه الفيزيائيون إلى نماذج تصف كيفية تفاعل الإلكترونات داخل المادة الصلبة. ويعد التنبؤ بأدنى حالة طاقة للنظام أحد التحديات المستمرة في هذا المجال، وهي قيمة تحدد استقرار وخصائص مواد مثل الموصلات الفائقة أو المغناطيسات. لعقود من الزمن، عانت الحواسيب الكلاسيكية في إجراء هذه الحسابات لأن عدد الترتيبات المحتملة للإلكترونات ينمو بسرعة هائلة بحيث يفوق قدرة حتى أقوى الحواسيب الفائقة. وتوفر الحواسيب الكمومية حلاً محتملاً، حيث يمكنها محاكاة هذه التفاعلات المعقدة بشكل طبيعي. ومع ذلك، لكي يجد الحاسوب الكمومي الإجابة الصحيحة، يجب أولاً تزويده بنقطة بداية — حالة ابتدائية — تكون قريبة نوعاً ما من الحل الحقيقي. فإذا كانت نقطة البداية هذه بعيدة جداً، فإن الحاسوب يهدر الوقت أو يفشل في إيجاد الإجابة على الإطلاق. وتكمن الصعوبة في إنشاء نقاط البداية هذه بكفاءة، خاصة بالنسبة للشبكات ثنائية الأبعاد الكبيرة من الذرات حيث تكون الإلكترونات مرتبطة ببعضها البعض بقوة وتتصرف بطرق شديدة الترابط.
لقد صمم فريق من الباحثين الآن طريقة لإنشاء نقاط البداية الحاسمة هذه باستخدام دوائر كمومية ضحلة، أو قصيرة، للغاية. ويركز عملهم على نموذج محدد لتفاعل الإلكترونات يُعرف باسم "نموذج هوبارد"، والذي يصف انتقال الإلكترونات بين المواقع أثناء تنافرها مع بعضها البعض. أدرك الباحثون أنه عندما يكون التنافر بين الإلكترونات قوياً جداً، يتصرف النظام تقريباً مثل نموذج أبسط يسمى "نموذج هايزنبرغ"، الذي يتعامل فقط مع لفات (spins) الإلكترونات بدلاً من حركتها. وقد قاموا ببناء دائرة كمومية تُجهز أولاً حالة تقريبية لنموذج اللف البسيط هذا. وتتكون هذه الحالة الابتدائية من ربط أزواج من الإلكترونات في "سينجليتات" (singlets)، وهي نوع محدد من الروابط الكمومية حيث تلغي لفات اثنين من الجسيمات بعضها البعض. ومن خلال ترتيب هذه الروابط في نمط محدد عبر الشبكة ثم تطبيق سلسلة من العمليات التي تسم تسمح للإلكترونات بالتقلب قلياً في مواقعها، أنشأ الفريق دائرة تحاكي عن كثب السلوك المعقد للنظام الكامل.
اختبر الفريق هذا النهج عن طريق محاكاة الدوائر على حواسيب كلاسيكية لمعرفة مدى جودة أدائها. بدأوا بتحسين إعدادات دائرتهم على شبكة صغيرة تبلغ أربعة في أربعة مواقع. ومن المثير للاهتمام، عندما طبقوا نفس الإعدادات بالضبط على شبكات أكبر بكثير، وصلت إلى عشرة في عشرة مواقع، استمرت الدوائر في الأداء بشكل جيد دون الحاجة إلى أي تعديل إضافي. ويشير هذا "النقل البارامتري" إلى أن المبادئ الفيزيائية التي توجه النظام الصغير تظل صالحة للأنظمة الأكبر، وهو اكتشاف حيوي لأن محاكاة الأنظمة الأكبر مباشرة هي أمر مستح المستحيل حالياً بالنسبة للحواسيب الكلاسيكية. وللتحقق من جودة نتائجهم، قارن الباحثون حالات دوائرهم بحالات مرجعية عالية الدقة تم إنشاؤها بواسطة تقنيات محاكاة متقدمة أخرى. ووجدوا أنه بالنسبة لشبكة عشرة في عشرة، أنتجت طريقتهم حالة ابتدائية قريبة بما يكفي من الحالة الأرضية الحقيقية لتكون مفيدة للحسابات الكمومية المستقبلية.
ولضمان عمل الطريقة لنموذج هوبارد الكامل، أضاف الباحثون مجموعة محددة من العمليات المصممة لالتقاط اللحظات النادرة التي تشغل فيها الإلكترونات نفس الموقع لفترة وجيزة أو تترك موقعاً فارغاً. وتعمل هذه العمليات، المستمدة من تحويل رياضي معروف باسم "تحويل شريفر-ولف" (Schrieffer–Wolff transformation)، كإجراء تصحيحي للحالة الأبسط القائمة على اللف. وعندما اختبروا هذه الدائرة المجمعة على شبكة صغيرة (أربعة في أربعة)، حيث يمكنهم حساب الإجابة الدقيقة، وجدوا أن إضافة خطوات التصحيح هذه حسنت دقة النتيجة بشكل كبير. كما أن الدائرة النهائية، التي تجمع بين تحضير اللف وتصحيحات تقلب الشحنة، مصممة لتكون "صديقة للأجهزة"، مما يعني أنها تستخدم حداً أدنى من الخطوات وتعتمد على اتصالات من المرجح أن تكون متاحة في الآلات الكمومية المستقبلية.
حسب الباحثون الموارد المطلوبة لتشغيل هذه الدائرة على حاسوب كمومي مقاوم للأخطاء (fault-tolerant)، وهو آلة نظرية قادرة على تصحيح أخطائها ذاتياً. وقدروا أن تجهيز الحالة لشبكة عشرة في عشرة سيتطلب حوالي 100,000 عملية منطقية محددة تُعرف باسم بوابات T. ويُعتبر هذا الرقم منخفضاً في سياق الحوسبة الكمومية، حيث يقع ضمن النطاق المتوقع تحقيقه في المستقبل القريب، وهو عصر يُطلق عليه أحياناً عصر "الميجاكوب" (megaquop). ومن خلال إثبات أن دائرة بسيطة وضحلة يمكنها توليد حالة ابتدائية عالية الجودة لنظام شديد الترابط ومعقد، قدم الفريق مساراً عملياً لاستخدام الحواسيب الكمومية لحل مشكلات طويلة الأمد في فيزياء المادة المكثفة. ويشير عملهم إلى أنه من خلال الاستفادة من العلاقة بين النماذج الأبسط والأكثر تعقيداً، يمكن للعلماء تجاوز الحاجة إلى الدوائر العميقة والمعرضة للأخطاء والاقتراب من محاكاة المواد الواقعية.
ملخص تقني: تحضير الحالة الأولية منخفضة العمق لتقدير طاقة الحالة الأرضية للأنظمة ثنائية الأبعاد شديدة الارتباط
بيان المشكلة تتطلب خوارزميات الحوسبة الكمومية لتقدير طاقة الحالة الأرضية، مثل خوارزمية تقدير الطور الكمومي (QPE)، حالة أولية ذات دقة (fidelity) غير مهملة مع الحالة الأرضية الحقيقية. وتتدرج احتمالية استخراج طاقة الحالة الأرضية بنجัง هذه الدقة (F0)؛ فإذا تضاءلت F0 أسياً مع حجم النظام، تصبح الخوارزمية غير قابلة للتطبيق بسبب عدد التكرارات المطلوب (1/F0). وبينما توفر خوارزميات الملحمة المتغيرة الكمومية (VQE) نهجاً هجيناً، إلا أنها تواجه تحديات في تكاليف تقييم التدرج (gradient) وتكلفة المحاكاة الكلاسيكية للدوائر الضحلة. كما أن تجميع التقريبات الكلاسيكية الدقيقة (مثل حالات المصفوفات ذات المنتجات المترابطة - MPS) في دوائر كمومية يؤدي غالباً إلى أعماق دوائر مفرطة، وهو أمر محظور في أجهزة الحوسبة الكمومية المتسامحة مع الخطأ في مرحلة الخطأ التصحيحي الكامل (FTQC) حيث تكون الموارد مثل الكيوبتات المنطقية وحالات السحر (magic states) محدودة. ويتمثل التحدي المركزي في تصميم دوائر كمومية ضحلة وفعالة في الموارد لتحضير حالات أولية عالية الدقة للأنظمة شديدة الارتباط، وتحديداً نموذج هوبارد (Hubbard model) ثنائي الأبعاد، دون الاعتماد على حلقات تحسين عميقة أو تعقيد كلاسيكي.
المنهجية يقترح المؤلفون استراتيجية لبناء دوائر كمومية ضحلة من خلال الاستففادة من الاتصال القوي في الاقتران بين نموذج هايزنبرج (Heisenberg model) ثنائي الأبعاد ونموذج هوبارد عند نصف الامتلاء عبر تحويل شريفر-ولف (Schrieffer–Wolff - SW). وتتقدم المنهجية عبر ثلاث مراحل رئيسية:
بناء دائرة هايزنبرج:
الأنزات (Ansatz): تستهدف الدائرة فضاء الحالة الفردية (S=0) لنموذج هايزنبرج مضاد للفيرومغناطيسية. الحالة الأولية هي نتاج ثنائيات (dimers) مفردة (∣Φ⟩) محددة على مطابقة تامة للشبكة.
الديناميكيات: يتم تطوير الحالة باستخدام طبقات من بوابات التبديل الأسية (eSWAP). تحترم هذه البوابات تناظر SU(2)، مما يضمن بقاء الحالة في فضاء الحالة المفردة (singlet subspace).
البنية: يتم تقسيم روابط الشبكة إلى أربعة مجموعات فرعية منفصلة (Bx0,Bx1,By0,By1). تطبق الدائرة بوابات eSWAP في تسلسل محدد (By1→By0→Bx1→Bx0) عبر عدة كتل (Nb). تحاكي هذه البنية فيزياء الروابط المتذبذبة (RVB) للحالة الأرضية لنموذج هايزنبرج.
الرسم الفيرميوني (Fermionic Mapping): يتم رسم حالات السبين (spin) إلى فضاء الامتلاء الفردي الفيرميوني باستخدام تضمين وحدوي محلي (V) يُدخل كيوبتات مساعدة.
توسيع هوبارد عبر بوابات SW:
لالتقاط تقلبات الشحنة (الدوبلونات والهولونات) الموجودة في نموذج هوبارد عند قيم U/t محدودة، يتم "تزيين" دائرة هايزنبرج ببوابات تقلب الشحنة.
تُشتق هذه البوابات مباشرة من مولد SW من الدرجة الأولى (Ω)، والتي تم تقريبها عبر تفكيك سوزوكي-تروتر (Suzuki–Trotter) إلى وحدات ثنائية الموقع (W⟨i,j⟩).
يجمع الأنزات النهائي بين تطور هايزنبرج (UB) وكتلة تقلب الشحنة (WB)، والمعلمات بها هي زوايا الدوران θ و ϕ.
تقدير الدقة والمعايرة:
نقل المعلمات: تُنقل المعلمات التي تم تحسينها على شبكة 4×4 إلى شبكات أكبر (6×6 إلى 10×10) دون الحاجة لمزيد من التحسين.
المحاكاة: يتم محاكاة حالات الدائرة باستخدام حالات الأزواج المتشابكة المسقطة (PEPS) مع تحديثات بسيطة موجهة بواسطة انتشار الاعتقاد (belief propagation - BP).
حدود الدقة: بما أن الحالة الأرضية الدقيقة غير معروفة للأنظمة الكبيرة، يقدر المؤلفون حداً أدنى لدقة الحالة الأرضية (F0). ويتم ذلك من خلال:
حساب الدقة (Fref) بين حالة الدائرة وحالة مرجعية من مونت كارلو المتغيرة (VMC).
وضع حد للدقة بين حالة VMC والحالة الأرضية الدقيقة (Fref,0) باستخدام الطاقة المرجعية، وطاقة الحالة الأرضية الدقيقة (من سلسلة ستوكاستيك السلسلة - SSE)، وفجوة الإثارة المفردة (Δ0).
تطبيق متباينة المثلث لزوايا فوبيني-فوبسودي (Fubini–Study) لاشتقاق حد أدنى صارم لـ F0.
المساهمات والنتائج الرئيسية
فعالية نقل المعلمات: بالنسبة لنموذج هايزنبرج، توفر المعلمات التي تم تحسينها على شبكة 4×4 حدوداً دنيا قوية للدقة لشبكات تصل إلى 10×10 دون إعادة تحسين. وتحديداً، بالنسبة لـ Nb=3، يظل الحد الأدنى المقدر لدقة الحالة الأرضية موجباً (≈0.028) حتى لشبكة 10×10، مما يشير إلى أن الحالة الأولية مناسبة لتضخيم السعة (amplitude amplification).
تحسين دقة هوبارد: تُظهر الحسابات الدقيقة على شبكة 4×4 أن إدراج بوابات تقلب الشحنة المشتقة من SW يحسن بشكل كبير الدقة تجاه الحالة الأرضية لهوبارد. بالنسبة لـ Nb=4، تزداد الدقة من ≈0.44 (هايزنبرج فقط) إلى ≈0.89 (مع زوايا الشحنة المحسنة).
كفاءة الموارد: يقدر المؤلفون تكلفة تحضير الحالة لشبكة 10×10. بافتراض تركيب Clifford+T مع خطأ لكل دوران قدره 10−5، فإن التكلفة الإجمالية تبلغ حوالي 105 من بوابات T. وهذا يضع تكلفة التحضير ضمن نطاق "الميجاكويب" (megaquop) لأجهزة FTQC المبكرة، وهو أقل بكثير من التكاليف المرتبطة بخوارزميات تقدير طاقة الحالة الأرضية اللاحقة.
القابلية للتوسع: تتدرج عمق الدائرة خطياً مع عدد روابط الشبكة (∣B∣)، وتسمح استراتيجية مشاركة المعلمات ببناء دوائر للأنظمة الكبيرة بناءً على تحسين الأنظمة الصغيرة.
الأهمية والادعاءات يزعم البحث تقديم مسار قابل للتطبيق لتحضير حالة أولية منخفضة التكلفة لتقدير طاقة الحالة الأرضية في الأنظمة شديدة الارتباط. ومن خلال الاستفادة من الاتصال الفيزيائي بين نموذجي هايزنبرج وهوبارد، يوضح المؤلفون أنه يمكن للدوائر الضحلة تحقيق دقة غير مهملة للحالة الأرضية على شبكات تتجاوز قدرة التقطيع الكلاسيكي الدقيق. وتشير النتائج إلى أن هذه الدوائر تعد مرشحة جذابة لعصر FTQC المبكر، حيث يكون تقليل عمق الدائرة وعدد بوابات T أمراً حاسماً. يدعم هذا العمل الفرضية القائلة بأن محاكاة الشبكات الصغيرة يمكن أن توجه بفعالية تصميم الدوائر الكمومية للأنظمة الأكبر والأكثر تعقيداً كلاسيكياً، بشرما كانت المعلمات قابلة للنقل وظل عمق الدائرة ضحلاً. ويشير المؤلفون صراحة إلى أنه بينما تعتبر حدودهم صارمة، فإن الدقة الفعلية قد تكون أعلى، وأن العمل المستقبلي قد يتضمن أنزات هاملتونية أكثر عمقاً أو تقنيات محاكاة أكثر دقة (مثل PEPS بالتحديث الكامل) لزيادة دقة هذه التقديرات.