Grow, Compress, and Perturb: A Hybrid Quantum-Classical Pathway to Scalable Quantum Chemistry on Quantum Computers
تقدم الورقة البحثية VULCAN، وهو إطار عمل هجين كمي-كلاسيكي يجمع بين توليد العمليات الخالي من تدرج الهاميلتوني عبر تدفق القوس المزدوج التبايني، واستراتيجية تحسين الدوائر القائمة على النمو ثم الضغط، ونظرية إبستين-نيسبيت الاضطرابية العشوائية، لتحقيق الدقة الكيميائية في محاكاة الكيمياء الكمية بتكاليف قياس ومتطلبات كيوبت أقل بكثير مقارنة بالطرق الحالية مثل batched ADAPT-VQE.
إن البحث عن فهم كيفية ترابط الذرات لتكوين الجزيئات هو أحد التحديات الأكثر جوهرية في العلم. بالنسبة للجزيئات البسيطة، يمكن للحواسيب الكلاسيكية حساب طاقة إلكتروناتها بدقة عالية. ومع ذلك، كلما كبرت الجزيئات أو تفاعلت إلكتروناتها بطرق معقدة، يتفاقم التعقيد الرياضي بشكل هائل. إذ ينمو عدد ترتيبات الإلكترونات الممكنة بسرعة كبيرة لدرجة أن أقوى الحواسيب الفائقة في العالم لا تستطيع حل المعادلات بدقة تامة للعديد من الأنظمة الكيميائية الهامة. وقد أعاق هذا القصور التقدم في تصميم أدوية ومواد ووقود جديد لفترة طويلة.
تقدم الحواسيب الكمومية طريقة محتملة لتجاوز هذا الجدار. فبدلاً من محاكاة الإلكترونات باستخدام بتات المعلومات، تستخدم هذه الآلات الكيوبتات (البتات الكمومية)، والتي يمكنها محاكاة سلوك الإلكترونات بشكل طبيعي. ومن الناحية النظرية، يمكن للحاسوب الكمومي أن ينمذج طاقة الجزيء بشكل مثالي. ومع ذلك، فإن الحواسيب الكمومية المتاحة اليوم صغيرة وعرضة للأخطاء، وهي مرحلة تُعرف باسم عصر "الكم المتوسط الضجيج". فهي لا تستطيع إجراء الحسابات العميقة والمعقدة المطلوبة للدقة المثالية. وللتعامل مع هذه الآلات غير المكتملة، طور العلماء طرقاً هجينة تقسم العمل بين حاسوب كلاسيكي وآخر كمومي. وأكثر هذه الطرق شعبية، وتسمى "المحلل الذاتي الكمومي المتغير" (Variational Quantum Eigensolver)، تحاول إيجاد أدنى حالة طاقة من خلال بناء نموذج رياضي خطوة بخطوة. ومع ذلك، أصبحت هذه العملية عائقاً؛ إذ تتطلب أفضل الطرق الحالية من الحاسوب إجراء آلاف القياسات لتحديد الخطوة التالية، وهي عملية بطيئة ومليئة بالضجيج وغالباً ما تتعثر في أخطاء محلية قبل الوصول إلى الإجابة الحقيقية.
قدم فريق من الباحثين نهجاً جديداً يسمى "VULCAN" لتجاوز هذه العقبات. تغير طريقتهم استراتيجية بناء النموذج الكمومي؛ فبدلاً من إضافة قطعة صغيرة واحدة في كل مرة والتوقف لقياس تدرج الطاقة لتحديد القطعة التالية، يقوم "VULCAN" بتنمية النموذج في طبقات كاملة وكبيرة. وهو يستخدم معادلة تدفق رياضية لتوليد مجموعة محددة من لبنات البناء التي تضمن صلتها بالجزيء قيد الدراسة. ومن خلال إلحاق طبقة كاملة من هذه اللبنات دفعة واحدة، تلغي الخوارزمية الحاجة إلى قياس تدرجات الطاقة لاختيار المؤثرات، مما يقلل بشكل كبير من عدد القياسات المطلوبة. هذا التحول يسمح للنظام بالوصو إلى إجابة دقيقة كيميائياً باستخدام جزء ضئيل فقط من موارد القياس التي تتطلبها الطرق السابقة.
اختبر الباحثون هذه الاستراتيجية على مجموعة متنوعة من الجزيئات، بما في ذلك سلاسل من ذرات الهيدروجين ومركبات أكثر تعقيداً مثل الماء والبروبين. وفي عمليات المحاكاة التي أجروها، وصلت الطريقة الجديدة إلى مستوى الدقة المطلوب باستخدام 5 إلى 20 بالمائة فقط من موارد القياس التي تتطلبها التقنية الرائدة الحالية. وبالنسبة لبعض الأشكال الجزيئية الأكثر صعوبة والممتدة، كان التحسن أكثر دراماتيكية. كما طور الفريق طريقة لضغط الدوائر الكمومية، عبر إزالة الخطوات غير الضرورية لإبقاء الحسابات قصيرة بما يكفي للآلات الحالية المعرضة للخطأ. علاوة على ذلك، أظهروا أنه بالنسبة للأنظمة الأكبر التي تتطلب كيوبتات أكثر مما توفره الأجهزة الحالية، يمكنهم إسقاط المشكلة على نسخة أصغر وأبسط من الجزيء دون فقدان الفيزياء الأساسية. وقد سمح لهم ذلك بحل مشكلة مكونة من 26 كيوبت على سجل محاكى بـ 20 كيوبت، بل وتناول نظام مكون من 32 كيوبت حيث فشلت الطرق القديمة في الوصول إلى الدقة اللازمة.
وللتعامل مع الأخطاء الصغيرة المتبقية دون زيادة عمق الدائرة الكمومية، أضاف الفريق خطوة نهائية يتم تنفيذها بالكامل على حاسوب كلاسيكي. تطبق هذه الخطوة تصحيحاً بناءً على "نظرية الاضطراب"، وهي تقنية قياسية في الفيزياء لتقدير تأثير التفاعلات الصغيرة المهملة. ومن خلال الجمع بين الحساب الكمومي وهذا التصحيح الكلاسيكي، حققوا دقة كيميائية لنظام مكون من 32 كيوبت في طبقة واحدة من الدائرة الكمومية. تشير هذه النتيجة إلى أن الجزء الأكثر صعوبة من الحساب يمكن إسناده إلى المعالجات الكلاسيكية، مما يسمح للحاسوب الكمومي بالقيام بعمل أقل مع تقديم نتائج دقيقة.
تُظهر النتائج، المستمدة من عمليات محاكاة عالية الأداء على ما يصل إلى 32 كيوبت منطقي، أن بنية مختلفة للكيمياء الكمومية هي أمر ممكن. يتبنى إطار عمل "VULCAN" استراتيجية تستبدل عبء القياس الثقيل للطرق الحالية باستراتيجية تنمو الحل في طبقات، وتضغط النتيجة، وتنتهي بتصحيح كلاسيكي. وبينما تستند هذه النتائج حالياً إلى عمليات محاكاة خالية من الضجيج ولم يتم نشرها بالكامل بعد على الأجهزة المادية، إلا أنها توفر مساراً واضحاً لكيفية قيام الحواسيب الكمومية المستقبلية بحل المشكلات الكيميائية التي لا يمكن الوصول إليها حالياً. ويشير العمل إلى أنه من خلال إعادة التفكير في كيفية بناء الخوارزميات الكمومية، يمكن للعلماء تحقيق استفادة أفضل من الموارد المحدودة المتاحة في أجهزتنا اليوم، مما يقرب وعد الكيمياء الكمومية من الواقع.
ملخص تقني: النمو، الضغط، والاضطراب: مسار هجين بين الكم والكلاسيكية نحو كيمياء كمية قابلة للتوسع على الحواسيب الكمية
بيان المشكلة تُعد الخوارزميات التباينية للمحلل الذاتي الكمي (VQE) من المرشحين الرائدين لحسابات البنية الإلكترونية على أجهزة الكم الكمي ذات المقياس المتوسط المليئة بالضجيج (NISQ). ومع ذلك، فإن قابليتها للتوسع تعوقها اختناقات جوهرية في تصميم الـ "ansatz" (المتخذات) والتحسين. تعاني المتخذات القياسية الفعالة للأجهزة من مشكلة "الهضاب القاحلة" (barren plateaus) والحساسية للحالات الأولية. أما النهج التكيفي مثل ADAPT-VQE، فبالرغم من معالجته لقابلية التعبير عبر اختيار المؤثرات ديناميكيًا من مجمع (pool)، إلا أنه يفرض عبئًا هائلًا في عمليات القياس؛ حيث يتطلب اختيار كل مؤثر تقييم تدرجات الطاقة لكامل المجمع، وهي عملية تتوسع بشكل سيئ وتكون عرضة لضجيج العينات (shot noise). علاوة على ذلك، فإن استراتيجيات الفرز المسبق الثابتة لتقليل أحجام المجمعات غالبًا ما تتطلب عمليات تشغيل أولية لـ VQE، مما يخلق تبعية دائرية غير مرغوب فيها للنشر القابل للتوسع.
المنهجية: إطار عمل VULCAN قدم المؤلفون إطار عمل VULCAN (النمو التبايني، الضغط الطبقي، الاضطراب التكيفي متعدد الأجسام)، وهو إطار هجين يجمع بين الكم والكلاسيكية مصمم لتخفيف هذه الاختناقات من خلال بنية "النمو، الضغط، والاضطراب".
النمو الطبقي الخالي من التدرج (النمو التبايني):
توليد مجمع المؤثرات: بدلاً من استخدام مجمعات ضخمة قائمة على الحدس (مثل مجمعات QEB أو مجمعات الإثارة القائمة على الكيوبت)، يقوم VULCAN بتوليد مجمع ديناميكي. فهو يستخدم معادلة التدفق المزدوج المتغير (VDBF) لتحديد الإثارات ذات الصلة فيزيائيًا مباشرة من الهاملتوني (Hamiltonian).
المجمعات الكاملة الدنيا (MCP): يتم تصفية مؤثرات VDBF ودمجها مع مجمع كامل أدنى (MCP) بحجم 2N−4 (حيث N هو عدد الكيوبتات). يضمن هذا أن المجمع يولد جبر لي (Lie algebra) المستهدف مع بقائه أصغر بمقدار 4 إلى 12 مرة من مجمعات QEB القياسية.
الإلحاق الطبقي: بدلاً من اختيار المؤثرات واحدًا تلو الآخر عبر فحص التدرج، يقوم VULCAN بإلحاق طبقة كاملة من المجمع بالمتخذ (ansatz) في وقت واحد. يتم تحسين معاملات الطبقة الجديدة بينما تظل الطبقات السابقة مجمدة، مما يلغي الحاجة تمامًا إلى قياسات تدرج الهاملتوني، مع الاعتماد فقط على تقييمات الطاقة.
إدارة عمق الدائرة (الضغط الطبقي):
يؤدي إلحاق الطبقات الكاملة إلى زيادة عمق الدائرة. ولمواجهة ذلك، يستخدم VULCAN مخطط ضغط موجه بالتداخل (overlap-guided compression).
بعد نمو وتحسين طبقة ما، تقوم حلقة ضغط فرعية بإعادة بناء دائرة أقل عمقًا. يتم إعادة اختيار المولدات بناءً على تدرج تداخل الحالة مع الهدف (الحالة الناتجة من مرحلة النمو) بدلاً من تدرجات الطاقة.
نظرًا لأن المجمع يتكون من سلاسل باولي مفردة، فإن مقدار تدرج التداخل يمكن الحصول عليه من تقييم تداخل واحد، مما يقلل تكاليف القياس بشكل كبير مقارنة بتدرجات الطاقة. يتم تقليم الدائرة حتى تصبح مكاسب عدم الدقة (infidelity) ضئيلة أو عند الوصول إلى كسر ضغط مستهدف.
تقليل عدد الكيوبتات (إسقاط الفضاء الجزئي):
بالنسبة للأنظمة التي تتجاوز سجلات الكيوبت المتاحة، يتم إسقاط الهاملتوني على فضاء جزئي تشغله المحددات المهيمنة لمتخذ أحادي الطبقة.
يتم اختيار أنماط الإثارة عبر مسح جشع (greedy sweep) وبحث شعاعي (beam search). ويتم إعادة كتابة الهاملتوني ومجمع المؤثرات على سجل الكيوبت المختزل.
يقود تقدير "إبستين-نيسبيت" (ENPT2) من الدرجة الثانية النمو التكيفي لحجم الفضاء الجزئي (k)، مما يوفر معيارًا للتوقف بناءً على مساهمة الطاقة المقدرة للمحددات المستبعدة.
المعالجة الكلاسيكية اللاحقة (الاضطراب):
لتحقيق الدقة الكيميائية دون تعميق الدوائر الكمية، يتم تطبيق تنفيذ عشوائي لـ ENPT2 من الدرجة الثانية على أجهزة GPU.
يستخدم هذا التصحيح دالة موجة VQE المحسنة والطاقة كبذرة. بالنسبة للأنظمة الأكبر حيث يكون التقييم الدقيق مستحيلاً، يستخدم أخذ العينات لأهمية القيم (importance sampling) جزءًا صغيرًا فقط من الحدود الاضطرابية، مما يسمح بالإنهاء المبكر لحلقة التحسين الكمي.
المساهمات الرئيسية
إلغاء فحص التدرج: من خلال إلحاق طبقات كاملة واستخدام مجمعات مشتقة من VDBF، يزيل VULCAN الحاجة إلى قياسات تدرج الهاملتوني، وهي المحرك الرئيسي للتكلفة في ADAPT-VQE.
بناء المجمع الديناميكي: يقوم الإطار بتوليد المؤثرات البادئة والمجمعات الكاملة في ثوانٍ على كمبيوتر كلاسيكي دون الحاجة لعمليات تشغيل أولية لـ VQE، مما يحل تحدي "اختيار البادئ" للمجمعات (MCPs).
قابلية ضبط مقايضة الموارد: يوفر الإطار مقايضة قابلة للضبط بين تكلفة القياس، وعمق الدائرة، وعدد الكيوبتات، والمعالجة الكلاسيكية اللاحقة، مما يسمح للمستخدمين بالتحسين وفقًا لقيود أجهزتهم الخاصة.
تكامل ENPT2 العشوائي: إن إدخال تصحيح ENPT2 العشوائي المعزز بـ GPU يتيح الوصول إلى الدقة الكيميائية باستخدام موارد كمية أقل بكثير، خاصة للأنظمة الأكبر.
النتائج قام المؤلفون باختبار VULCAN مقابل ADAPT-VQE بنظام الدفعات (batched) مع مجمعات QEB باستخدام محاكاة الحالة المتجهة الخالية من الضجيج (Hyperion) لأنظمة تشمل سلاسل الهيدروجين (H6 إلى H10)، وBeH2، وH2O (أساس 6-31G، 26 كيوبت)، وC3H6 (32 كيوبت).
تكلفة القياس: حقق VULCAN الدقة الكيميائية عند 5-20% فقط من عدد تقييمات حدود باولي (تكلفة القياس) مقارنة بـ ADAPT-VQE بنظام الدفعات في جميع الأنظمة المختبرة. وبالنسبة لـ H2O (6-31G)، انخفضت التكلفة إلى 0.05% من الخط المرجعي عند دمجه مع ENPT2 العشوائي.
عمق الدائرة: بينما أنتجت تكرارات النمو الأولية دوائر أعمق، أدى الضغط الموجه بالتداخل إلى تقليل عدد بوابات CNOT عن خط الأساس لـ ADAPT-VQE في خمسة من أصل سبعة أنظمة مختبرة. وبالنسبة للهندسات الممتدة، كانت التحسينات أكثر وضوحًا (على سبيل المثال، 0.14× خط الأساس لـ H10).
عدد الكيوبتات: نجح إسقاط الفضاء الجزئي في تقليل حجم السجل لـ H2O (6-31G) من 26 إلى 20 كيوبت مع الحفاظ على الدقة الكيميائية في التقطيع الدقيق (exact diagonalization). كما قلل هذا من تعقيد الهاملتوني من 25,035 إلى 8,911 حد باولي.
الأنظمة الكبيرة: بالنسبة لنظام C3H6 المكون من 32 كيوبت، حيث فشل ADAPT-VQE بنظام الدفعات في الوصول إلى الدقة الكيميائية، حققت طبقة واحدة من VALCON مع تصحيح ENPT2 العشوائي الدقة المستهدفة.
التحقق من الأجهزة: كإثبات أولي للمفهوم، تم اختبار الطريقة على جهاز Trapped-ion من نوع Quantinuum System Model H2، لـ H2، حيث أعطت نتائج قريبة من الدقة الكيميائية (1.7×10−3 فوق FCI).
الأهمية والادعاءات يزعم البحث أن VULCAN يوفر مسارًا عمليًا للكيمياء الكمية القابلة للتوسع من خلال معالجة عبء القياس وقيود عمق الدائرة التي تواجهها الخوارزميات التكيفية الحالية. ويؤكد المؤلفون أن الطريقة ليست خوارزمية ثابتة واحدة، بل هي إطار عمل يسمح للمستخدمين بالتنقل بين المقايضات بناءً على المورد الأكثر ندرة (ميزانية القياس، عدد الكيوبتات، أو عمق الدائرة).
تظل ادعاءات الورقة متواضعة فيما يتعلق بالنشر الفوري على الأجهزة. يقر المؤلفون بأن جميع النتائج الأساسية استندت إلى محاكاة خالية من الضجيج، وأن سلوك VULCAN تحت تأثير ضجيج العينات وأخطاء الأجهزة، خاصة للطبقات التي تحتوي على حوالي 103 من المعاملات، لم يتم إثباته بعد. كما يشيرون إلى أن التنفيذ الحالي لا يحافظ بدقة على عدد الجسيمات أو اللف المغزلي (spin)، مما قد يؤدي إلى تلوث اللف (spin contamination) في الأنظمة شديدة الارتباط، وهي مشكلة معروفة في الطرق غير المتكيفة مع التماثل. وقد تم تحديد العمل المستقبلي في التركيز على النشر الفيزيائي في عصر NISQ، ودمج الإسقاط مع الحلقة الكاملة لمجموعات أساس أكبر، واستعادة تماثلات الحالة عبر حدود الجزاء (penalty terms).