VQE Validation on the Mononuclear T1 (Blue-Copper) Site: A Stepping Stone to the Multi-Copper Laccase Cluster

لقد أمضت الطبيعة مليارات السنين في صقل الإنزيمات للقيام بمهام كيميائية يعجز المهندسون البشر عن محاكاتها. ومن بين هذه الآلات البيولوجية، يبرز إنزيم "اللاكاز" (laccase) لقدرته على تفكيك الألياف النباتية الصلبة وحتى البلاستيك الاصطناعي، وهي عملية مدفوعة بعنقود من ذرات النحاس في قلبه. ولفهم كيفية عمل هذا الإنزيم، يجب على العلماء رسم خريطة لسلوك الإلكترونات التي تحوم حول ذرات النحاس هذه. وهذه مهمة بالغة الصعوبة لأن الإلكترونات لا تعمل بشكل مستقل؛ بل هي متشابكة بعمق، وتتفاعل مع بعضها البعض بطرق تخلق شبكة معقدة من التفاعلات. ولعقود من الزمن، اصطدمت أقوى الحواسيب الفائقة بحائط مسدود عند محاولة محاكاة هذه العناقيد النحاسية المحددة. فالحسابات الرياضية المطلوبة لتتبع كل تفاعل محتمل تصبح ضخمة للغاية لدرجة أن أسرع الأجهزة تنفد ذاكرتها ووقتها، مما يترك فجوة في فهمنا لكيفية عمل هذه الإنزيمات.

ولسد هذه الفجوة، يتجه الباحثون إلى نوع جديد من الحواسيب: الحاسوب الكمي. وخلافاً للآلات التقليدية التي تعالج المعلومات في شكل بتات ثنائية، تستخدم الحواسيب الكمية بتات كمية، أو "كيوبتات" (qubits)، والتي يمكن أن توجد في حالات متعددة في آن واحد، مما يسمح لها نظرياً بالتنقل عبر شبكات الإلكترونات المعقدة هذه بسهولة. ومع ذلك، وقبل أن تتمكن هذه الآلات من معالجة إنزيم "اللاكاز" بالكامل، يجب على العلماء إثبات نجاح أساليبهم على قطع أصغر وأكثر قابلية للإدارة من اللغز. وهذا هو بالضبط ما قامت به الباحثة لوسيا ماليكوفا (Lucia Malícková)؛ حيث ركزت على موقع نحاسي واحد داخل الإنزيم، يُعرف باسم موقع T1، واختبرت طريقة محددة تسمى "المحلل الذاتي المتغير الكمي" (Variational Quantum Eigensolver أو VQE). لم يكن هدفها حل مشكلة الإنزيم بأكمها على الفور، بل التحقق من صحة مسار العمل (pipeline)—التأكد من أن البرمجيات والخوارزميات يمكنها محاكاة موقع النحاس بدقة على معالج كمي قبل التوسع نحو الهدف الأكبر والأكثر تعقيداً والذي يستعصي على الحواسيب التقليدية.

بدأ الباحثون ببناء نموذج رقمي دقيق لموقع T1، مع التركيز على ذرة النحاس وجيرانها المباشرين. واستخدموا تقنية لعزل الإلكترونات الأكثر أهمية، مما قلص المشكلة إلى نظام يمكن رسم خرائطته على 18 كيوبت. وللتحقق من إعدادهم، قاموا أولاً بتشغيل محاكاة مثالية خالية من الضجيج على "محاكي" (emulator) قوي، وهو برنامج متطور يحاكي سلوك حاسوب كمي حقيقي. سمح لهم ذلك بحساب الطاقة الدقيقة للنظام باستخدام طريقة تسمى "القطرية الدقيقة" (exact diagonalization)، والتي تعمل كنقطة مرجعية نهائية. ووجدوا أن طاقة موقع النحاس هذا كانت تقريباً -2518.995 Eh. ومع وضع هذا الخط المرجعي، اختبروا نهجين مختلفين لإيجاد هذه الطاقة باستخدام طريقة VQE.

استخدم النهج الأول دارة مصممة مسبقاً تُعرف باسم "الأنسا" (ansatz) الكفؤة للأجهزة. تشبه هذه الطريقة استخدام أداة عامة تناسب العديد من الوظائف ولكنها قد لا تكون مثالية لأي منها بمفردها. وعندما قام الباحثون بتشغيل هذه الطريقة على المحاكي الخاص بهم، كانت النتائج جيدة ولكنها لم تكن دقيقة بما يكفي للوصول إلى الدقة الكيميائية. فقد انحرفت الطاقة التي حسبوها بنحو 173 إلى 278 مللي-هارتي، وهي فجوة كبيرة في عالم الكيمياء الكمية. أظهر هذا أنه رغم نجاح الطريقة الأساسية، إلا أنها افتقرت إلى المرونة اللازمة لاستيعاب التعقيد الكامل لسلوك إلكترونات النحاس. ثم جرب الفريق نهجاً ثانياً أكثر تطوراً يسمى ADAPT-VQE. فبدلاً من استخدام دارة ثابتة، تقوم هذه الطريقة ببناء الحل خطوة بخطوة، عبر إضافة مؤثرات رياضية محددة فقط عندما تكون هناك حاية إليها لتحسين النتيجة. إنها عملية ديناميكية تُفصّل الدارة لتناسب المشكلة المطروحة تحديداً.

كانت نتائج هذا النهج التكيفي مذهلة. فمن خلال الاختيار المتكرر للمؤثرات الأكثر فائدة من مجموعة كبيرة من الاحتمالات، تمكنت الخوارزمية من إزالة أكثر من 99 بالمائة من الخطأ الموجود في الطريقة الأبسط. وكانت حسابات الطاقة النهائية منحرفة بنسبة 4.03 مللي-هارتي فقط، وهو تحسن هائل جعل المحاكاة ضمن النطاق الذي يعتبر دقيقاً كيميائياً. مثل هذا النجاح تحسناً بمقدار أربعين ضعفاً مقارنة بالنهج القياسي. وأكد الباحثون أن هذا النجاح لم يكن مجرد صدفة برمجية؛ فقد تحققوا بدقة من أن المحاكاة تحترم القوانين الفيزيائية للنظام، مثل حفظ عدد الإلكترونات، مما يضمن أن النتائج حقيقية فيزيائياً.

ومع ذلك، فإن المحاكاة على محاكي مثالي لا تماثل التشغيل على آلة حقيقية معرضة للأخطاء والضجيج. وللاستعداد للمستقبل، قام الفريق بمحاكاة ما سيحدث إذا قاموا بتشغيل أفضل دارة لديهم على حاسوب كمي حقيقي بمستويات الضجيج الحالية. ووجدوا أن قراءة الطاقة ستنزاح بشكل كبير، لتصبح أقل دقة بسبب تداخل أخطاء البوابات ثنائية الكيوبت. وقد أكد ذلك أن مجرد تشغيل الكود على الأجهزة لن يكون كافياً؛ إذ سنحتاج إلى تطبيق استراتيجيات تخفيف الخطأ. ولحسن الحظ، وبما أن طريقتهم تحافظ بصرامة على عدد الإلكترونات، يمكنهم استخدام تقنية تسمى "الاختيار اللاحق للتماثل" (symmetry post-selection) لاستبعاد أي نتائج تنتهك هذه القاعدة الفيزيائية. وبالاقتران مع طرق تصحيح الخطأ الأخرى، يوفر هذا مساراً واضحاً للحصول على نتائج دقيقة حتى على الأجهزة الصاخبة وغير المثالية.

الهدف النهائي من هذا العمل ليس مجرد محاكاة ذرة نحاس واحدة، بل فك أسرار عنقود إنزيم "اللاكاز" بأكمله، والذي يتطلب نظاماً أكبر بكثير يتراوح بين 50 إلى 70 كيوبت. عند هذا المقياس، تصبح المشكلة معقدة للغاية لدرجة أن أي حاسوب فائق تقليدي لن يتمكن أبداً من حلها. إن التحقق من صحة مسار العمل هذا على موقع T1 المكون من 18 كيوبت يثبت أن الطريقة تعمل وأن الأجهزة قادرة على التعامل مع المهمة. إنه يؤسس قاعدة متينة، تظهر أنه مع الخوارزميات المناسبة وإدارة الأخطاء، يمكن للحواسيب الكمية في النهاية معالجة المشكلة الكاملة المستعصية كلاسيكياً والمتمثلة في كيفية تفكيك هذه الإنزيمات للمواد. وتعد هذه الدراسة بمثابة خطوة حاسمة، تنقل المجال من مجرد إمكانية نظرية إلى مسار عملي ومتحقق منه لحل أحد أكثر الألغاز الكيميائية تعقيداً في الطبيعة.