Molecular Dynamics with Nuclear Effects on Quantum Computers
تقدم هذه الورقة البحثية وتتحقق من صحة خوارزمية هجينة جديدة بين الكم والتقليدية تدمج التأثيرات الكمومية النووية في محاكاة الديناميكا الجزيئية من المبادئ الأولية باستخدام حل القيم الذاتية الكمومية التبايني، مما يظهر دقة محسنة في الأطياف الاهتزازية ونقل البروتونات عديم الحاجز لأنظمة مثل أيون زوندل باستخدام دوائر كفؤة عتاديًا مناسبة للأجهزة الكمومية الحالية.
المؤلفون الأصليون:Lukas Haßfurth, Juliane Heitkämper, Elias Walter, Birger Horstmann
تُعد الذرات هي لبنات بناء المادة، لكنها ليست كرات صلبة وساكنة كما تُصورها الرسوم التوضيحية المدرسية غالباً. ففي قلب كل ذرة يكمن النواة، وتحيط بها سحب من الإلكترونات. وفي معظم الحسابات الكيميائية، يعامل العلماء النوى الثقيلة كنقاط ثابتة بينما تتحرك الإلكترونات الأخف وزناً حولها بسرعة، وهو تبسيط يعمل بشكل جيد مع العديد من المواد. ومع ذلك، عندما يتعلق الأمر بالهيدروجين، أخف العناصر، ينهار هذا التبسيط؛ فبسبب خفة وزن نوى الهيدروجين، لا تتصرف كنماذج رخامية صغيرة ثابتة، بل تهتز وتنفذ عبر حواجز الطاقة بطريقة لا يمكن وصفها إلا من خلال الفيزياء الكمومية. هذا السلوك، المعروف باسم التأثيرات الكمومية النووية، أمر بالغ الأهمية لفهم كيفية تحرك البروتونات في العمليات البيولوجية مثل التمثيل الضوئي وكيفية تماسك جزيئات الماء. ومع ذلك، فإن محاكاة هذه الطبيعة الكمومية المهتزة للنوى الذرية على الحواسيب الفائقة التقليدية غالباً ما يكون مستحيلاً، لأن التعقيد الرياضي ينمو بسرعة كبيرة يصعب التعامل معها.
لقد اتخذ فريق من الباحثين في مركز الفضاء الألماني ومعهد هلمهولتز أولم خطوة كبيرة نحو حل هذه المشكلة من خلال نقل المحاكاة إلى حاسوب كمومي. فقد طوروا طريقة هجينة جديدة تجمع بين قوة المعالج الكمومي وموثوقية الحاسوب التقليدي. في نهجهم، يعمل الحاسوب الكمومي كآلة حاسبة متخصصة للغاية تحدد الطاقة والقوى المؤثرة على الذرات في أي لحظة، بينما يتولى الحاسوب التقليدي عملية حركة الذرات نفسها. والأهم من ذلك، أن هذه الطريقة لا تعامل بعض النوى الذرية كنقاط ثابتة، بل كأجسام كمومية يمكن أن توجد في أماكن متعددة في آن واحد، مما يسمح بصورة أكثر واقعية لكيفية اهتزاز الجزيئات وتفاعلها.
ولاختبار خوارزميتهم، قام الباحثون بمحاكاة ثلاثة أنظمة محددة: جزيء هيدروجين بسيط، وجزيء ماء، وأيون زوندل (Zundel ion)، وهو عبارة عن عنقود من جزيئي ماء يتشاركان بروتوناً واحداً. ويعد أيون زوندل مهماً بشكل خاص لأنه يعمل كنموذج أدنى لنقل البروتون، وهي عملية أساسية في الكيمياء والبيولوجيا. ووجد الفريق أنه عندما تجاهلوا الطبيعة الكمومية للنوى، أنتجت عمليات المحاكاة نتائج غير دقيقة في كثير من الأحيان. ومع ذلك، بمجرد تضمين هذه التأثيرات الكمومية النووية، تغير سلوك الجزيئات بشكل كبير. ففي حالة أيون زوندل، أدى إدراج هذه التأثيرات إلى إزالة حاجز الطاقة الذي يمنع البروتون عادةً من التحرك بحرية بين جزيئي الماء. وبدلاً من أن يكون عالقاً في "وادي"، أصبح تحرك البروتون بلا عائق فعلياً، مما سمح له بالتنقل ذهاباً وإياباً بسهولة. وقد تماشى هذا التحول بشكل أكبر مع البيانات التجريبية، مما قلل بشكل كبير من الأخطاء في الترددات الاهتزازية المتوقعة للجزيئات.
كما استكشف الباحثون كيفية جعل هذه المحاكاة عملية للحواسيب الكمومية المتاحة اليوم، والتي لا تزال في مراحلها الأولى والمشوبة بالضجيج. فقد اختبروا طرقاً مختلفة لترميز المشكلة، وقارنوا بين طريقة معقدة وعالية الدقة وبين مناهج أبسط وأكثر إيجازاً مصممة للأجهزة الحالية. واكتشفوا أنه بينما وفرت الطريقة المعقدة النتائج الأكثر دقة، كانت المناهج المبسطة والموجزة فعالة بشكل مدهش. فبالنسبة لجزيء الماء وجزيء الهيدروجين، أنتجت هذه المناهج الانسيابية نتائج تقارب تلك التي قدمتها الطرق المعقدة، مع تطلب موارد حوسبة أقل بكثير. وهذا يشير إلى إمكانية تشغيل الخوارزمية على الأجهزة الكمومية القريبة، حتى تلك ذات القدرة المحدودة.
ومع ذلك، لم يكن الطريق نحو الدقة المثالية خالياً من التحديات. فقد وجد الفريق أن نجاح المحاكاة يعتمد بشدة على كيفية اختيارهم لتمثيل الحالات الكمومية للنوى. فبالنسبة لجزيء الماء، كان التمثيل البسيط كافياً، ولكن بالنسبة لأيون زوندل الأكثر تعقيداً، كانوا بحاجة إلى تضمين أوصاف أكثر تفصيلاً للمدارات النووية لالتقاط التفاعلات الصحيحة. وعندما استخدموا وصفاً أكثر تفصيلاً، نجحت المحاكاة في إعادة إنتاج نمط "القمة المزدوجة" المحدد في طيف الاهتزاز الخاص بالأيون، والذي يتطابق مع ما يُرى في التجارب الواقعية. وبدون هذا التفصيل الإضافي، افتقرت المحاكاة لهذه الميزات الحاسمة. وتخلص الدراسة إلى أنه على الرغم من أن الحواسيب الكمومية الحالية ليست قوية بما يكفي لمحاكاة كل تفاصيل الأنظمة الكيميائية الكبيرة، إلا أن هذا النهج الهجين يوفر مساراً قابلاً للتطبيق للمضي قدماً. فمن خلال الاختيار الدقيق للأجزاء التي سيتم محاكاتها على الحاسوب الكمومي وتلك التي سيتم التعامل معها تقليدياً، يمكن للعلماء البدء في نمذجة السلوك الكمومي للنوى الذرية بمستوى من الدقة كان بعيد المنال في السابق.
ملخص تقني: الديناميكا الجزيئية مع التأثيرات النووية على الحواسيب الكمومية
بيان المشكلة تعد التأثيرات الكمومية النووية (NQEs)، مثل نفق البروتون وطاقة نقطة الصفر، بالغة الأهمية للوصف الدقيق لعمليات مثل انتقال البروتون والترابط الهيدروجيني. ومع ذلك، فإن دمج هذه التأثيرات في حسابات الكيمياء الكمومية على الأجهزة الكلاسيكية غالبًا ما يكون مكلفًا حوسبيًا، حيث يمكن أن يتسبب ذلك في نمو فضاء التفاعل الكامل (FCI) بشكل أسي. وبينما توفر الحواسيب الكمومية حلًا محتملاً نظرًا لتوسعها الخطي في متطلبات الكيوبت بالنسبة لحجم النظام، إلا أن إثبات ميزة كمومية عملية للمحاكاة الديناميكية التي تتضمن التأثيرات الكمومية النووية لا يزال يمثل تحديًا. تفتقر الخوارزميات الحالية غالبًا إلى تضمين التأثيرات الكمومية النووية، أو محاكاة الديناميكا، أو تعتمد على مسارات محددة مسبقًا.
المنهجية يقترح المؤلفون خوارزمية هجينة (كمومية-كلاسيكية) جديدة للديناميكا الجزيئية من المبادئ الأولى (ab-initio) تدمج التأثيرات الكمومية النووية عبر إطار عمل المدارات النووية-الإلكترونية (NEO).
إطار عمل NEO: يتم التعامل مع النوى الخفيفة المختارة (تحديدًا البروتونات في هذا العمل) كميكانيكا كمومية جنبًا إلى جنب مع الإلكترونات باستخدام مدارات نووية مماثلة للمدارات الذرية الإلكترونية. وهذا يتجنب الفصل القياسي بين بورن-أوبنهايمر لهذه النوى المحددة.
خوارزمية الديناميكا الهجينة: تتبع الخوارزمية نهج التكميم الثاني بناءً على ديناميكا بورن-أوبنهايمر الكلاسيكية.
المكون الكمومي: يتم استخدام الخوارزمية التباينية لليجين (VQE) على الحاسوب الكمومي لتحديد طاقة الحالة الأرضية اللحظية (E0) والقوى المؤثرة على النوى. يتم حساب القوى باستخدام مبرهنة هيلمان-فايجمان من خلال تقييم القيمة المتوقعة لمؤثر القوة ∇IH(R) على الحالة الأرضية ∣ψ0(R)⟩.
المكون الكلاسيكي: يتم تحديث الهندسة الجزيئية كلاسيكيًا باستخدام خوارزمية "فيرليت" للسرعة (velocity verlet) بناءً على القوى المحسوبة بواسطة المعالج الكمومي.
تقنيات الفضاء النشط: لمعالجة عدد الكيوبتات المحدود في الأجهزة ذات الكم المتوسطة المليئة بالضجيج (NISQ)، طبق المؤلفون تضمين "هارتري-فوك". يتم اختيار فضاء نشط للحاسوب الكمومي بينما تُعامل المدارات غير النشطة على مستوى هارتري-فوك. تشمل معايير الاختيار الإنتروبيا للمدار المنفرد، والتباين على طول مسار التفاعل، والحفاظ على التماثل.
اختيار الـ Ansatz (الدالة التجريبية): تستخدم الدراسة كلاً من الـ Ansatz المستوحى كيميائيًا (UCCSD/NEOUCCSD) والـ Ansatz الفعال من حيث الأجهزة (Real Amplitudes و Tiled Unitary Product Ansatz - tUPS) لتقييم الجدوى على الأجهزة القريبة من العصر الحالي.
مجموعات القاعدة (Basis Sets): تستخدم المحاكاة مجموعة القاعدة الإلكترونية 6-31G ومجموعات قاعدة نووية محددة مشتقة من DZSPDN (وهي DZSNB للنوع s فقط؛ و DZSPNB للنوعين s و p).
المساهمات الرئيسية
تنفيذ الخوارزمية: يقدم العمل أول تنفيذ لخوارزمية ديناميكا جزيئية من المبادئ الأولى (AIMD) هجينة (كمومية-كلاسيكية) تدمج صراحةً التأثيرات الكمومية النووية باستخدام طريقة NEO على حاسوب كمومي.
استراتيجية الفضاء النشط: طور المؤلفون بروتوكولًا محددًا لاختيار الفضاءات النشطة في ديناميكا NEO التي تأخذ في الاعتبار تباين المدارات والتماثل، مما يضمن الاتساق الفيزيائي أثناء العمليات الديناميكية.
الكفاءة في الأجهزة: توضح الدراسة أن الـ Ansatzs المدمجة والفعالة من حيث الأجهزة يمكن أن تحقق نتائج مقاربة لـ الـ Ansatzs الأكثر دقة (UCCSD)، مما يشير إلى جدوى التنفيذ على الأجهزة الكمومية القريبة من العصر الحالي.
حساسية مجموعة القاعدة: يسلط البحث الض {الضوء على الاعتماد الحرج لدقة المحاكاة على اختيار مجموعات القاعدة النووية (تحديدًا ضرورة وجود مدارات النوع p لنمذجة اقتران الأنماط المعقدة) والفضاءات النشطة.
النتائج تم التحقق من صحة الخوارزمية من خلال عمليات محاكاة لثلاثة أنظمة: جزيء الهيدروجين (H2)، والماء (H2O)، وأيون زوندل (H5O2+).
الهيدروجين (H2): أدى إدراج التأثيرات الكمومية النووية عبر NEO إلى تصحيح المبالغة في تقدير الترددات الاهتزازية التي لوحظت في الطرق الإلكترونية البحتة، رغم أنه أدى إلى نقص طفيف في التقدير. حققت الـ Ansatzs الفعالة من حيث الأجهزة (tUPS و Real Amplitudes) أخطاءً أقل من 1% في حالة NEO، وهو ما يُعزى إلى إلغاء الأخطاء المصادفي بين المبالغة في التقدير للـ Ansatz والنقص في التقدير لـ NEO.
الماء (H2O): أعطت محاكاةات NEO أخطاء نسبية أقل من الطرق الإلكترونية البحتة لجميع الأنماط الطبيعية (الانحناء، التمدد المتماثل، والتمدد غير المتماثل). ظل نمط الانحناء، الذي يتضمن تفاعلًا وثيقًا بين البروتون والبروتون، يمثل تحديًا ولكنه أظهر تحسنًا مع NEO.
أيون زوندل (H5O2+):
سطح طاقة الجهد: في الحالة الإلكترونية البحتة، يظهر جهد انتقال البروتون بنية ذات بئر مزدوج مع حاجز. عند إدراج التأثيرات الكمومية النووية، تصبح حركة انتقال البروتون فعليًا بدون حاجز، مما يحول الجهد إلى بئر واحد. يشير هذا إلى أن حركة نقطة الصفر للبروتون المشترك تمارس تأثيرًا ملحوظًا، مما يضع الحالة الاهتزازية الأرضية فوق الحاجز الإلكتروني البحت.
الأطياف: فشلت المحاكاة الإلكترونية البحتة في التقاط الثنائي الطيفي المميز حول 1000 سم−1 والاقتران بين انتقال البروتون وأنماط تمايل الماء.
نجاح NEO: نجحت محاكاةات NEO مع المدارات النووية من النوع p في إعادة إنتاج الثنائي الطيفي والاقتران بين نمط انتقال البروتون وتمايل/تمدد الماء. كان من الضروري استخدام نسخة مختزلة من NEOUCCSD (باختيار 17 استثارة من أصل 63) لحل مشكلة الأمثلة للفضاء النشط الأكبر.
الدقة: أظهرت نتائج NEO حدًا أقصى للخطأ يقل عن 40 سم−1 لكل قمة مقارنة بحل NEOCASCI الدقيق ضمن الفضاء النشط، مع مطابقة نوعية لأطياف ما قبل التفكك التجريبية.
الأهمية والادعاءات يدعي البحث أن خوارزميتهم الهجينة تنجح في إثبات جدوى محاكاة الديناميكا الجزيئية مع التأثيرات الكمومية النووية على الحواسيب الكمومية. ويؤكد المؤلفون أن:
التأثيرات الكمومية النووية ضرورية لنمذجة انتقال البروتون في أيون زوندل بشكل صحيح، وتحديدًا لإزالة حاجز الطاقة الاصطناعي وتمكين الاقتران الصحيح للأنماط.
بينما يوفر UCCSD أعلى دقة، فإن الـ Ansatzs الفعالة من حيث الأجهزة توفر مسارًا قابلًا للتطبيق للأجهزة القريبة من العصر الحالي، بشرًا بشرط اختيار الفضاء النشط ومجموعات القاعدة النووية بعناية.
اختيار الفضاء النشط ومجموعة القاعدة النووية (تحديدًا تضم مدارات النوع p لأيون زوندل) أكثر أهمية للدقة من الـ Ansatz المحدد المستخدم.
يضع هذا العمل أساسًا للدراسات المستقبلية على الأنظمة الأكبر، مع ملاحظة أن المزيد من التحسينات في استخدام الموارد أو الدقة يمكن تحقيقها من خلال الانتقال من القواعد الكانونية إلى قواعد المدارات الطبيعية.
يظل المؤلفون متواضعين فيما يتعلق بالحالة الراهنة للتكنولوجيا، مع الإقرار بأن نتائجهم تعتمد على محاكيات الحالة المتجهة (statevector) الخالية من الضجيج، وأن تحقيق ميزة كمومية كبيرة لم يتم إثباته بعد على الأجهزة الفيزيائية. يعمل هذا العمل كإثبات مفهوم لتوسيع طرق الكيمياء الكمومية لتشمل درجات الحرية النووية ديناميكيًا.