Tailored coupled cluster method with sample-based quantum diagonalization: Application to titanium-based metallocene catalytic reactions for 1-hexene production

ملخص تقني: طريقة كلوستر المزدوجة المصممة مع التشخيص الكمي القائم على العينات

بيان المشكلة
يتطلب التنبؤ الدقيق بالانتقائية في تفاعلات الحفز من نوع الميتالوسين القائم على التيتانيوم لإنتاج 1-هكسين، حساب فروق الطاقة الحرة بين حالات الانتقال المتنافسة (إدراج الإيثيلين مقابل انتقال بيتا-هيدريد) بدقة تبلغ 1 سعرة حرارية/مول (kcal/mol). تواجه طرق البنية الإلكترونية الحالية ثنائية: حيث تتسم نظرية الدالة الكثافية (DFT) بالكفاءة الحسابية ولكنها غالباً ما تفتقر إلى الدقة اللازمة بسبب اعتمادها على الدوال الوظيفية، بينما تعد طرق الدالة الموجية عالية المستوى (مثل CCSD(T)) مكلفة حسابياً لأنظمة المعادن الانتقالية الواقعية. علاوة على ذلك، فإن مقاربات الحوسبة الكمية مثل التشخيص الكمي القائم على العينات (SQD) محدودة حالياً بقيود الأجهزة التي تحصر الفضاء النشط في أحجام صغيرة جداً لا يمكنها استيعاب الارتباط الديناميكي الكامل المطلوب لنتائج دقيقة كيميائياً في الأنظمة الحفزية. وبناءً على ذلك، يظل التعامل مع كل من الارتباط الساكن (داخل الفضاء النشط) والارتباط الديناميكي (من المدارات الخارجية) في آن واحد تحدياً كبيراً.

المنهجية
تقترح هذه الدراسة إطار عمل هجين (كمي-كلاسيكي) يدمج التشخيص الكمي القائم على العينات (SQD) مع نظرية كلوستر المزدوجة المصممة (TCC) لمعالجة قيود حجم الفضاء النشط ومعالجة الارتباط.

1. النظام المستهدف: تم استخدام نموذج جزيئي أدنى لمحفز تيتانيوم من نوع (ansa-cyclopentadienyl-arene) كاتيوني لاستقصاء الطاقات النسبية لحالتين انتقاليتين، T3 (إدراج الإيثيلين) وT4 (انتقال بيتا-هيدريد). ولضمان إطار مداري متسق، تم تعزيز حالة T3 بجزيء إيثيلين صريح (T3(+C2H4))، ليتطابق مع عدد الذرات والإلكترونات في T4.
2. SQD للارتباط الساكن: تم استخدام SQD لمعالجة البنية الإلكترونية للفضاء النشط. وباستخدام أنساب (LUCJ - Local-Unitary Cluster Jastrow)، تم تطبيق SQD على معالج "هيرون" (Heron) الكمي من IBM، حيث يقوم SQD بأخذ عينات عشوائية من السلاسل البتية (bitstrings) لبناء هاميلتوني مسقط ضمن فضاء جزئي. يتيح ذلك معالجة الفضاءات النشطة (التي تصل إلى 30 مداراً مكانياً) والتي تعد ضخمة تركيبياً (combinatorially) بالنسبة لطريقة (CASCI) الكلاسيكية. وتستخدم الطريقة استعادة التكوين ذاتي الاتساق للتخفيف من ضوضاء الأجهزة وضمان صحة أعداد الإلكترونات واللف المغزلي.
3. TCC للارتباط الديناميكي: عملت الدوال الموجية المستمدة من SQD كحالات مرجعية لحسابات كلوستر المزدوجة المصممة (TCC). في هذا الإطار، يتم "تصميم" (تثبيت) سعات الكلوستر للفضاء النشط ($\hat{T}_{AS}$) بناءً على نتائج SQD، مما يلتقط الارتباط الساكن والإثارات من الرتب العليا (الثلاثيات، والرباعيات) داخل الفضاء النشط. أما السعات المتبقية للفضاء الخارجي ($\hat{T}_{Ext}$)، فيتم تحسينها عبر معادلات (CCSD) القياسية لالتقاط الارتباط الديناميكي.
4. التوسعات: طبقت الدراسة أيضاً توسعة الثلاثيات الاضطرابية، TCC(T)، لتحسين الدقة. كما تم استخدام مخطط تصحيح طاقة لضمان أن تعكس طاقة الفضاء النشط مستوى SQD بدلاً من تقريب CCSD الأدنى.
5. الإعداد الحسابي: أجريت جميع حسابات النقطة الواحدة (single-point) اللاحقة لـ DFT في الطور الغازي لعزل تأثيرات الارتباط الإلكتروني، رغم أن عمليات تحسين الهندسة شملت تأثيرات المذيب (التولوين) عبر نموذج الاستمرارية المستقطب (PCM).

النتائج الرئيسية

• مشكلات التقارب مع SQD بمفرده: عند تطبيق SQD بمعزل عن غيره، لم تتقارب الطاقات النسبية لحالات الانتقال ضمن الفضاءات النشطة المتاحة حالياً. كانت النتائج حساسة لحجم الفضاء النشط، مما يشير إلى أن الارتباط الساكن وحده غير كفي لتمثيل هذا النظام.
• نجاح SQD-TCC: حقق النهج الهجين SQD-TCC نتائج متقاربة بشكل معقول للطاقات النسبية. ومن خلال دمج الارتباط الديناميكي عبر إطار عمل TCC، نجحت الطريقة في موازنة معالجة الارتباط الإلكتروني.
• تأثير الارتباط الساكن: لوحظ وجود تباين كبير (يتجاوز 1 سعرة حرارية/مول) بين نتائج CCSD(T) القياسية ونتائج SQD-TCC(T). وهذا يسلط الضوء على أن تأثيرات الارتباط الساكن حاسمة في نظام التيتانيوم هذا، وأن تجاهلها (كما في حالة CCSD(T) أحادي المرجع) يؤدي إلى أخطاء كبيرة في التنبؤ بالانتقائية.
• قابلية توسع الفضاء النشط: أظهرت الدراسة أن SQD يتيح معالجة فضاءات نشطة (مثل 30 مداراً) يصعب الوصول إليها عملياً باستخدام طرق CASCI الكلاسيكية في ظل الموارد الحسابية المتاحة (حتى 192 نواة CPU و4 تيرابايت ذاكرة).

الأهمية والادعاءات
يزعم البحث أن دمج الدوال الموجية متعددة المراجع المستمدة من الحوسبة الكمية مع نظرية كلوستر الكلاسيكية يوفر مساراً قابلاً للتطبيق نحو محاكاة دقيقة كيميائياً للأنظمة الحفزية المعقدة. وتحديداً، يوضح العمل أن:

1. الضرورة الهجينة: كلا الارتباطين الإلكترونيين، الساكن والديناميكي، ضروريان للوصف الموثوق للكيمياء الأساسية في حفز الميتالوسين القائم على التيتانيوم.
2. التقدم المنهجي: نجح إطار عمل SQD-TCC في سد الفجوة بين قدرات الحوسبة الكمية الحالية (المحدودة بالفضاءات النشطة الصغيرة) ومتطلبات المحاكاة الحفزية الواقعية (الفضاءات النشطة الكبيرة + الارتباط الديناميكي).
3. التنبؤ بالانتقائية: تؤثر معالجة الارتباط الإلكتروني بشكل كبير على الانتقائية المتوقعة للتفاعل، حيث توفر الطريقة الهجينة نهجاً أكثر قوة من الطرق الكلاسيكية عالية المستوى للأنظمة ذات الارتباط الساكن القوي.

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