Efficient, precise DFT calculations of NMR shieldings: Revisiting the finite field approach

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

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

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

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

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

اختبر الفريق طريقتهم الجديدة على مجموعة متنوعة من الجزيئات، بدءًا من الحلقات البسيطة مثل البنزين وصولاً إلى الأدوية الكبيرة والمعقدة مثل أتورفاستاتين (atorvastatin)، المستخدم لعلاج الكوليسترول المرتفع. وقارنوا نتائجهم بالمعيار الذهبي للطرق التحليلية، المعروفة بدقتها ولكنها صعبة التشغيل للغاية على الأنظمة الكبيرة. أظهرت النتائج أن طريقتهم الجديدة أنتجت أرقامًا مطابقة تقريبًا للمعيار الذهبي، مع أخطاء صغيرة جدًا لدرجة أنها لا تُذكر مقارنة بعدم اليقين الطبيعي في القياسات التجريبية. وبالنسبة للجزيئات الأصغر، كانت الطريقة الجديدة أسرع من النهج التقليدي. أما بالنسبة لأكبر الجزيئات، التي تحتوي على آلاف الذرات، فقد أصبحت الطريقة الجديدة أسرع بشكل كبير، حيث تفوقت على النهج التقليدي بفارق واسع عند تشغيلها على مجموعة من الحواسيب (cluster). أثبت هذا أن طريقتهم تتوسع بكفاءة، مما يعني أنها تتحسن مع زيادة القدرة الحاسوبية بدلاً من الاصطدام بسقف الأداء.

ولإثبات القيمة العملية لعملهم، طبق الباحثون الطريقة على لغز تاريخي شهير في الكيمياء: هيكل ناتج طبيعي يسمى هيكساسيكلينول (hexacyclinol). تم عزل هذا الجزيء من فطر وُجد في خشب البتولا السيبيري، وكان شكله الدقيق موضوع نقاش حاد لسنوات. اقترحت مجموعات بحثية مختلفة هياكل مختلفة، ولم ينتهِ الجدل إلا عندما أكد التخليق الكلي للجزيء الشكل الصحيح في عام 2006. استخدم الباحثون خوارزميتهم الجديدة لحساب إشارات NMR لكل من الهيكل المقترح أصلاً والهيكل المعدل والصحيح. أظهرت النتائج بوضوح أن الهيكل المعدل يتوافق مع البيانات التجريبية بشكل أفضل بكثير من الاقتراح الأصلي. أكد هذا التطبيق الناجح أن طريقتهم ليست مجرد تحسين نظري، بل هي أداة عملية قادرة على حل المشكلات الهيكلية الواقعية التي استعصت على العلماء في الماضي.

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

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