Quantum-Classical Auxiliary Field Quantum Monte Carlo with Matchgate Shadows on Trapped Ion Quantum Computers
تستعرض هذه الورقة البحثية سير عمل متكاملًا، معززًا بوحدات معالجة الرسومات ومحطمًا للأرقام القياسية، لمحاكاة حواجز التفاعلات الكيميائية باستخدام طريقة مونت كارلو للمجال المساعد الكمي-الكلاسيكي مع ظلال البوابة المطابقة (matchgate shadows) على حاسوب "IonQ Forte" الكمي، محققةً تسريعًا كبيرًا ونتائج ضمن نطاق 10 كيلو كالوري/مول من القيم المرجعية.
المؤلفون الأصليون:Luning Zhao, Joshua J. Goings, Willie Aboumrad, Andrew Arrasmith, Lazaro Calderin, Spencer Churchill, Dor Gabay, Thea Harvey-Brown, Melanie Hiles, Magda Kaja, Matthew Keesan, Karolina Kulesz, Andrii MLuning Zhao, Joshua J. Goings, Willie Aboumrad, Andrew Arrasmith, Lazaro Calderin, Spencer Churchill, Dor Gabay, Thea Harvey-Brown, Melanie Hiles, Magda Kaja, Matthew Keesan, Karolina Kulesz, Andrii Maksymov, Mei Maruo, Mauricio Muñoz, Bas Nijholt, Rebekah Schiller, Yvette de Sereville, Amy Smidutz, Felix Tripier, Grace Yao, Trishal Zaveri, Coleman Collins, Martin Roetteler, Evgeny Epifanovsky, Arseny Kovyrshin, Lars Tornberg, Anders Broo, Jeff R. Hammond, Zohim Chandani, Pradnya Khalate, Elica Kyoseva, Yi-Ting Chen, Eric M. Kessler, Cedric Yen-Yu Lin, Gandhi Ramu, Ryan Shaffer, Michael Brett, Benchen Huang, Maxime R. Hugues, Tyler Y. Takeshita
المؤلفون الأصليون: Luning Zhao, Joshua J. Goings, Willie Aboumrad, Andrew Arrasmith, Lazaro Calderin, Spencer Churchill, Dor Gabay, Thea Harvey-Brown, Melanie Hiles, Magda Kaja, Matthew Keesan, Karolina Kulesz, Andrii Maksymov, Mei Maruo, Mauricio Muñoz, Bas Nijholt, Rebekah Schiller, Yvette de Sereville, Amy Smidutz, Felix Tripier, Grace Yao, Trishal Zaveri, Coleman Collins, Martin Roetteler, Evgeny Epifanovsky, Arseny Kovyrshin, Lars Tornberg, Anders Broo, Jeff R. Hammond, Zohim Chandani, Pradnya Khalate, Elica Kyoseva, Yi-Ting Chen, Eric M. Kessler, Cedric Yen-Yu Lin, Gandhi Ramu, Ryan Shaffer, Michael Brett, Benchen Huang, Maxime R. Hugues, Tyler Y. Takeshita
تخيل أنك تحاول التنبؤ بكيفية سلوك آلة معقدة بمجرد النظر إلى مخططاتها. في عالم الكيمياء، هذه الآلة هي الجزيء، والمخطط هو مجموعة من المعادلات التي تصف رقص الإلكترونات حول الذرات. لعقود من الزمن، استخدم العلماء أجهزة كمبيوتر كلاسيكية قوية لحل هذه المعادلات، ولكن عندما تصبح الجزيئات معقدة — مثل تلك التي تتضمن معادن انتقالية تُستخدم في الأدوية المنقذة للحياة أو المحفزات الصناعية — يصبح الحساب معقدًا للغاية لدرجة أن حتى أسرع الحواسيب الفائقة في العالم تتعثر. الأمر يشبه محاولة فك عقدة من مليار خيط وأنت معصوب العينين؛ فكلما أضفت المزيد من الخيوط، زادت الصعوبة، حتى تصبح المهمة مستح المستحيل. هنا يأتي دور الحوسبة الكمومية في القصة. فبدلاً من محاولة فك العقدة بيد واحدة، تستخدم الحواسيب الكمومية القواعد الغريبة للعالم الكمومي لاستكشاف العديد من الاحتمالات في وقت واحد، مما يوفر طريقًا مختصرًا محتملاً لفهم هذه التفاعلات الكيميائية المعقدة. السؤال الكبير هو: هل يمكننا حقًا استخدام هذه الآلات الجديدة الهشة والمليئة بالضجيج لحل مشكلات كيميائية واقعية اليوم، أم أننا لا نزال ننتظر مستقبلًا مثاليًا وخاليًا من الأخطاء؟
تحكي هذه الورقة قصة فريق قرر التوقف عن الانتظار والبدء في بناء جسر بين العالمين الكمومي والكلاسيكي لحل لغز كيميائي صعب ومحدد. لقد ركزوا على تفاعل يسمى "تفاعل سوزوكي-مياورا" (Suzuki–Miyaura cross-coupling)، وهو عملية حيوية تُستخدم لربط ذرات الكربون معًا لإنشاء الأدوية. وبشكل أكثر تحديدًا، أرادوا حساب حاجز الطاقة لخطوة حيث تقوم ذرة النيكل بالإمساك بجزء من جزيء لبدء التفاعل. هذه مشكلة صعبة للغاية لأن ذرة النيكل تمتلك إلكترونات "مرتبطة" بقوة، مما يعني أنها تعمل بطريقة منسقة للغاية وفوضوية، وهو ما يكسر النماذج الحاسوبية القياسية. لم يكتفِ الباحثون بمحاكاة هذا على جهاز كمبيوتر فحسب؛ بل أجروا تجربة حقيقية باستخدام حاسوب كمومي من شركة "IonQ"، والذي يستخدم الأيونات المحاصرة (الذرات المشحونة) كوحدات معالجة، مدمجًا مع جيش هائل من وحدات معالجة الرسومات (GPUs) من "NVIDIA" للقيام بالعمل الشاق.
طوّر الفريق سير عمل ذكيًا يسمى "مونت كارلو الكمومي المساعد للمجال الكلاسيكي" (QC-AFQMC). فكر في هذا كسباق تتابع. أولاً، يقوم الحاسوب الكمومي بإعداد "حالة تجريبية"، وهي تخمين تقريبي لشكل إلكترونات الجزيء. وبدلاً من مطالبة الحاسوب الكمومي بالقيام بكل الرياضيات الصعبة، استخدم الفريق تقنية تسمى "ظلال بوابة الماتش" (matchgate shadows). تخيل التقاط صورة لجسم ثلاثي الأبعاد من زوايا عشوائية مختلفة لإنشاء "ظل" يلتقط جوهره دون الحاجة لرؤية كل تفصيل فيه. يتم بعد ذلك إرسال هذه الظلال إلى الحواسيب الفائقة الكلاسيكية، والتي تستخدم حيلًا رياضية متقدمة (مثل التفاضل الخوارزمي) لإكمال الحساب والتنبؤ بحاجز الطاقة.
كانت النتائج مزيجًا من الانتصار والواقع المرير. فعندما استخدموا "المحاكي المثالي" للحاسوب الكمومي (نسخة مثالية خالية من الضجيج)، كانت نتائجهم قريبة جدًا من المرجع الذهبي، حيث وقعت ضمن نطاق ±4 كيلو كالوري/مول من القيمة المتوقعة. ومع ذلك، عندما أجروا التجربة على أجهزة "IonQ Forte" الفعلية المليئة بالضجيج، انحرفت النتائج قليًا، لتنتهي ضمن نطاق 10 كيلو كالوري/مول من المرجع. ورغم أن هذا ليس مثاليًا بعد، إلا أنه يمثل قفزة هائلة للأمام. وتستبعد الورقة صراحةً فكرة أن هذه الطريقة بطيئة جدًا أو مكلفة لدرجة تجعلها غير مفيدة؛ بل في الواقع، حققوا تسريعًا قدره 656 ضعفًا في الوقت اللازم لحل المشكلة مقارنة بالطرق الرائدة السابقة، بفضل خوارزمياتهم الجديدة وتسريع وحدات معالجة الرسومات. كما وجدوا أنه بينما أدخلت الأجهزة الكمومية بعض الضجيج الذي أدى إلى إزاحة قيم الطاقة، إلا أن الطريقة ظلت قوية بما يكفي لتقديم رؤى كيميائية مفيدة.
اكتشف الفريق أيضًا أن ضبط الحاسوب الكمومي خصيصًا لهذه المهمة أحدث فرقًا كبيرًا. فمن خلال تحسين كيفية إرسال وتنفيذ الدوائر، سرعوا عملية جمع البيانات بمقدار 9 مرات. استخدموا 24 كيوبت (qubit) في المجمل: 16 لتمثيل الحالة التجريبية للجزيء و8 كيوبت إضافية "مساعدة" (ancilla) للمساعدة في رصد وتصحيح الأخطاء. الخلاصة النهائية هي أنه على الرغم من أننا لم نصل بعد إلى مرحلة الكيمياء المثالية والخالية من الأخطاء على الأجهزة الكمومية، إلا أن سير العمل المتكامل هذا يثبت أنه يمكننا الآن التعامل مع تفاعلات حفزية حقيقية ومعقدة كانت بعيدة المنال سابقًا. إنها خطوة مهمة نحو جعل الكيمياء الكمومية أداة عملية لاكتشاف الأدوية، مما يظهر أنه مع المزيج الصحيح من الأجهزة الكمومية، وقوة الحوسبة الفائقة الكلاسيكية، والخوارزميات الذكية، يمكننا البدء في حل العقد المتشابكة للكيمياء الجزيئية.
ملخص تقني: طريقة مونت كارلو للمجال المساعد الكمي الكلاسيكي-الكمي باستخدام ظلال "ماتش-غيت" (Matchgate Shadows) على الحواسيب الكمية ذات الأيونات المحاصرة
بيان المشكلة لا تزال النمذجة الدقيقة لحواجز التفاعلات الكيميائية التي تشمل محفزات المعادن الانتقالية (TM)، مثل تفاعل اقتران سوزوكي-مياورا، تشكل تحدياً كبيًرا للكيمياء الحاسوبية الكلاسيكية. تظهر هذه الأنظمة ارتباطاً إلكترونياً قوياً وطابعاً متعدد المراجع بسبب الإلكترونات غير المزدوجة في مدارات d أو f، مما يؤدي إلى فشل أو انهيار طرق المرجع الواحد القياسية مثل نظرية الدالة الوظيفية للكثافة (DFT) وحتى الطريقة التي تُعد "المعيار الذهبي" وهي [CCSD(T)]. وبينما توفر خوارزمية تقدير الطور الكمي (QPE) حلولاً دقيقة، إلا أنها تتطلب أجهزة مقاومة للأخطاء تتجاوز القدرات الحالية. أما الخوارزميات التباينية (VQE)، فهي مقاومة للضوضاء ولكنها تعاني من تكاليف قياس باهظة وحساسية تجاه ضوضاء الأجهزة. وبناءً على ذلك، هناك حاجة إلى خوارزميات هجينة كمية-كلاسيكية يمكنها الاستفادة من الأجهزة الكمية القريبة من المدى (near-term) لإعداد حالات تجريبية مرتبطة مع التخفيف من الاختناقات الحسابية للمعالجة الكلاسيكية اللاحقة.
المنهجية يقدم المؤلفون سير عمل متكامل لخوارزمية مونت كارلو للمجال المساعد الكمي الكلاسيكي-الكمي (QC-AFQMC)، والتي يتم تنفيذها ضمن بيئة حوسبة فائقة كمية-كلاسيكية مسرعة. يدمج سير العمل حاسوب IonQ Forte الكمي ذو الأيونات المحاصرة مع وحدات معالجة الرسومات NVIDIA على خدمة أمازون ويب سيرفيسز (AWS).
الإطار الخوارزمي: تستخدم الطريقة نهج مونت كارلو للمسقط (projector Monte Carlo) الكلاسيكي لـ AFQMC، والذي يعمل على تطوير دالة موجية في زمن تخيلي لإيجاد الحالة الأرضية. ولمعالجة مشكلة الطور الفرميوني، يلزم وجود دالة موجية تجريبية (∣ΨT⟩). في هذا التنفيذ، يتم إعداد الحالة التجريبية على حاسوب كمي باستخدام الخوارزمية التباينية (VQE) مع نموذج (upCCD) المزدوج المزدوج الموحد.
ظلال ماتش-غيت (Matchgate Shadows): بدلاً من استخدام اختبار هادامارد القياسي أو ظلال كليفورد الكلاسيكية (التي تتوسع أسياً لحسابات التداخل في AFQMC)، يستخدم سير العمل ظلال ماتش-غيت. تسمح هذه التقنية بالتقييم الفعال للتداخلات ⟨ΨT∣ϕ⟩ بين الحالة التجريبية الكمية وحالات السائر (walker states) الكلاسيكية. ويتوسع حساب التداخل بشكل متعدد الحدود (O(N4.5)) مع حجم النظام N.
طاقة الارتباط الافتراضية (VCE): للتعامل مع الأنظمة الكيميائية الكبيرة حيث يتجاوز فضاء المدارات الكامل عدد الكيوبتات المتاحة، يتم تقسيم النظام إلى فضاء نشط (يتم التعامل معه بواسطة الحالة التجريبية الكمية) وفضاء افتراضي. تتيح تقنية VCE التقاط الارتباطات الإلكترونية خارج الفضاء النشط باستخدام محدد سلايتر واحد، مما يسمح للموارد الكمية بالتركيز على الفضاء النشط شديد الارتباط.
التحسينات الخوارزمية: قدم المؤلفون تحسينات خوارزمية كبيرة لخطوة المعالجة الكلاسيكية اللاحقة. فمن خلال استخدام التفاضل الخوارزمي (تحديداً مشتقات مصفوفات بافافيان) بدلاً من الفروق المحدودة، انخفض التعقيد الحسابي لحساب الطاقات المحلية ومصطلحات انحياز القوة من O(N8.5) و O(N7.5) إلى O(N5.5) و O(N4.5) على التوالي.
التنفيذ على الأجهزة:
الكمي: يتم تحسين الحالة التجريبية على محاكي مثالي، ويتم تنفيذ دوائر ظلال ماتش-غيت على معالج IonQ Forte. تم تنفيذ ضبط الأداء لبرمجيات التحكم في QPU والبرامج الثابتة لزيادة الإنتاجية.
الكلاسيكي: تتم المعالجة اللاحقة، بما في ذلك الانتشار في الزمن التخيلي وتقييم الطاقة، على AWS ParallelCluster باستخدام وحدات معالجة الرسومات NVIDIA H100 وH200، مع الاستفادة من مكتبات cuBLAS وcuSOLVER وcuTENSOR.
المساهمات الرئيسية
أول عرض تجريبي واسع النطى: يسجل هذا العمل أكبر تجربة QC-AFQMC باستخدام ظلال ماتش-غيت تم إجراؤها على أجهزة كمية حتى الآن، باستخدام 24 كيوبت (16 للحالة التجريبية، و8 للمساعدة في تخفيف الخطأ) على معالج IonQ Forte.
تسريع خوارزمي: يحقق هذا التنفيذ تحسناً بمقدار 656 ضعفاً في "الوقت اللازم للحل" (time-to-solution) لخطوة المعالجة الكلاسيكية اللاحقة مقارنة بأحدث تنفيذ من قِبل Huang وآخرون [41]، مدفوعاً بشكل أساسي بتقليل التعقيد التقاربي عبر التفاضل الخوارزمي وتسريع وحدات معالجة الرسومات.
تحسين إنتاجية الأجهزة: أدى ضبط الأداء الخاص ببرمجيات التحكم في IonQ Forte إلى تسريع بمقدار 9 أضعاف في جمع قياسات دوائر ماتش-غيت.
سير عمل متكامل: يضع البحث حجر الأساس لسير عمل متكامل يجمع بين إعداد الحالة الكمية، وتوموجرافيا الظل، وانتشار مونت كارلو الكلاسيكي المسرع بوحدات معالجة الرسومات للمشكلات الكيميائية المعقدة.
النتائج طُبقت الخوارزمية لمحاكاة خطوة الإضافة التأكسدية لتفاعل سوزوكي-مياورا المحفز بالنيكل (يتضمن النوع B، والحالة الانتقالية [B-C]‡، والمنتج C) باستخدام فضاء نشط (8,8).
الدقة: عندما تم أخذ عينات الماتش-غيت على محاكي مثالي، كانت حواجز التفاعل المحسوبة ضمن نطاق عدم اليقين البالغ ±4 كيلو كالوري/مول من نتيجة مرجع CCSD(T).
تأثير الضوضاء: عند قياسها على QPU الفعلي، انحرفت النتائج بنحو 10 كيلو كالوري/مول عن المرجع. تسببت الضوضاء في إزاحة نوعية في ترتيب الطاقة المتوقع للمتفاعل والمنتج، مما يسلط الضوء على القيود الحالية لتخفيف الخطأ لهذا التطبيق المحدد.
عدد الجسيمات: بينما أعطت المحاكة المثالية أعداد جسيمات قريبة من الحقيقة الأرضية (8 إلكترونات)، أظهرت قياسات QPU انحرافاً إلى حوالي 10.5 إلكترون، مما يشير إلى ضوضاء كبيرة في الأجهزة، رغم أن نتائج الطاقة ظلت أكثر صموداً أمام هذه الضوضاء من تقديرات عدد الجسيمات.
الأهمية يزعم هذا البحث أن هذا العمل يمثل "خطوة نحو عمليات محاكاة الكيمياء الكمية العملية على الأجهزة الكمية". إنه يتحدى الإجماع السابق بأن QC-AFQMC مكلف للغاية للأنظمة ذات الأهمية العلمية الحقيقية، من خلال إثبات أنه مع التحسينات الخوارخوارزمية وتسريع الأجهزة، يمكن تطبيق الطريقة على معقدات المحفزات المعدنية الانتقالية. تحدد الدراسة أنه بينما تعتبر الطريقة قابلة للتطبيق، إلا أن هناك حاجة لمزيد من التطوير في تقنيات تخفيف الخطأ واختيار حالات تجريبية قوية لاستغلال إمكاناتها بالكامل على أجهزة NISQ (الكمية ذات المدى المتوسط والمشوبة بالضوضاء). ويؤكد العمل على ضرورة التصميم المشترك لبرمجيات التحكم في الأجهزة الكمية مع أعباء العمل الخوارزمية المحددة لتحقيق تسريع ملموس.