نظرية الكثافة الوظيفية

نظرية الدالة الوظيفية للكثافة  (DFT, Density Functional Theory) هي أحد أهم الطرق المستعملة في الفيزياء و الكيمياء النظريين و بواسطتها نستطيع أن نحدد خصائص نظام متعدد الجسيمات (الطاقة الكلية للنظام, الكثافة الإلكترونية للمدارات, المعاملات الفيزيائية والضوئية للمادة....)، وهي واحدة من أكثر الطرق استخداما في العمليات الحسابية الكمومية بسبب إمكانية تطبيقها على أنظمة متنوعة و بتكلفة و سرعة العالية.
تعتمد الطرق التقليدية في حلها معادلة شرودنغر لنظام مكون من عدة ذرات ولا سيما نظرية هارتري-فوك و الطرق المستنبطة منها على الدالة الموجية ذات 3N متغير (حيث N هو العدد الإجمالي لجزيئات النظام)، لذلك تكون لمعادلات المراد حلها جد معقدة و تتطلب جهدا كبيرا· الهدف الرئيسي نظرية الرئيسي نظرية الدالة الوظيفية للكثافة هو استبدال الدالة الموجية بالدالة الوظيفية للكثافة ذات ثلاثة متغيرات فقط و جعلها كقاعدة للحساب، لذلكفالتعامل معها كمفهوم رياضي أو فيزيائي هو أسهل بكثير ، فمبدأ DFT هو إعادة صياغة للمسألة الكمومية وتحويلها من مسألة لنظام متعددالجسيمات إلى مسألة أحادية الجسيمة .
تعود جذور DFT للنموذج الذي وضعه لويلين توماس و انريكو فيرمي في أواخر 1920. ومع ذلك فإنه لم يكن ب الإمكان إستعمالها حتى منتصف 1960. ومنذ ذلك الوقت و هي في تطور متصاعد. مع تطور وقوة البرامج المعلوماتية وسرعة تطبيق العمليات الحسابية. جاءت نتائج (DFT) النظرية قريبة على نحو مرض تماما مع البيانات التجريبية وبتكاليف منخفضة نسبيا مع الطرق التقليدية التي تستهلك المال والوقت معا, ولكن على الرغم من التحسن الكبير الذي طرأ مؤخرا والتحسين المستمر للبرامج، لا تزال هناك بعض الصعوبات في استخدام نظرية الكثافة الوظيفية (DFT) لوصف التفاعلات بين الجزيئات وخاصة القوى الضعيفة (فان دير وال), الحالات الانتقالية للإلكترونات, حساب الفجوة (band Gap) في أشباه الموصلات, هذه الثغرات يمكن أن تؤثر على نتائج (DFT) "على الأقل عندما تستعمل وحدها دون تصحيح" حيث يتم تطوير طرق جديدة ل (DFT) للتغلب على هذه المشكلة, من خلال إحداث تعديلات على الدالة الوظيفية أو إدراج شروط في برامج الحساب وهو موضوع البحث الحالي في هذا المجال.

أساسيات

"إن القوانين الفيزيائية الأساسية اللازمة للتعبير الرياضي لجزء مهم من الفيزياء و الكيمياء كلها معروفة سلفا، وتكمن الصعوبة فقط في التطبيق الدقيق لهذه القوانين الذي يؤدي إلى المعادلات جد معقدة يتعين حلها " عن - بول ديراك 1929

معادلة شرودنجر

المعادلة الأساسية التي يتعين حلها لوصف نظام يتكون مثلا من العديد من الجزيئات أو الذرات هي معادلة إرفين شرودنغر (1887-1961) المعروفة بمعادلة شرودنجر و تكتب على الشكل التالي :
بحيث  هو المؤثر الهاملتوني (Hamiltonian operator)الكمي للمجموعة قيد الدرس والمكونة من  إلكترون ( ذات الأساس  ) و  لنوى الذرات ( ذات الأساس  ). الحدين الأولين من المؤثر الهاملتوني هما على التوالي هو مؤثر الطاقة الحركية للإلكترونات و نوى الذرات على التوالي، أما الحدود الثلاثة المتبقية فهي تمثل مؤثرات الطاقة الوضعية أو الجهد لمختلف التفاعلات البنية ( الإلكترون-النواة ) و ( الإلكترون -الإلكترون ) و ( النواة- النواة ).
حل هذه المعادلة على حالتها هذه معقدة جدا، ولكن وباستخدام تقريب بورن-أوبنهايمر الذي وضعه كل من ماكس بورن و روبرت أوبنهايمر أصبح ممكنا إلى حد ما. و هذا التقريب ينص على إمكانية إهمال حركة النوى مقارنة مع حركة الإلكترونات. و ذلك راجع للكتلة النواة العالية و سرعتها البطيئة مقارنة بكتلة و سرعة الإلكترونات. هذا يمكننا من عزل حركة الإلكترونات عن حركة النوى و تكون النتيجة الابقاء على الحدود التي تحتوي تمثل طاقة الالكترون الحركية والوضعية بالاضافة إلى طاقة الوضع الناتج عن تفاعل الالكترون مع النواة. وبهذا نكون قد اقتصرنا على حل معادلة شرودنغر للإلكترونات . ويمكننا إعادة كتابة المعادلة اعلاه على الشكل التالي :
سنقتصر في دراستنا المتبقية على دراسة  و الذي يحدد معادلة شرودنغر المتعددة الإلكترونات:
وقد تم تطوير العديد من طرق لحل هذه المعادلة لكنها كانت تعتمد على الدوال الموجية كطريقة هرتري فوك . توفر DFT طريقة بديلة تعتبر كثافة الإلكترونات هي الأساس لمعرفة خصائص المجموعة المدروسة.

لدالة الوظيفة

كائن رياضي يمثل علاقة رياضية نربط بها دالة بدالة أخرى عكس مفهوم الد آلة العادية التي نربط بها دالة بعدد ذو قيمة معينة. فمثلا الدالة الموجية  هي دالة عادية، أما الطاقة  فيمكن التعبير عنها بدلالة  فتسمى دالة وظيفية

الكثافة الإلكترونية لزوج من الإلكترونات

كثافة زوج من الإلكترونات هو احتمال العثور في وقت واحد على اثنين من بين  من الإلكترونات و المتموضعين في وحدتي الحجم و  هي:

نموذج توماس فيرمي

نظرية دالة الكثافة الوظيفية نشأت من نموذج توماس فيرمي، التي وضعتها لويلين توماس (1903-1992) و إنريكو فيرمي (1901-1954) في عام 1927. و اعتمدت طريقة توماس فيرمي على نموذج إحصائي لتقريب التوزيع الإلكتروني حول الذرات. هذا النموذج أعطى صورة عن إمكانية الإعتماد على الكثافة الإلكترونية لحساب الطاقة الحركية رغم أن هذا النموذج أعطى نتائج ضعيفة لعدم دقته.ويمكن اعتباره

نظرية هوهنبرج و كوهين

التقريب الذي وضعه بيير هوهنبرج و والتر كوهين مكن من إعادة صياغة نظرية الكثافة الدالّية (Density Functional Theory) المقترحة من قبل توماس فيرمي ووضع نظرية دقيقة لنظام متعدد الجسيمات. هذا التقريب ينطبق على أي نظام من الجسيمات المتفاعلة تتحرك في جهد خارجي و يقوم على اثنين من النظريات الأساسية التي تم إثباتها من طرف كوهين في بحثه لسنة 1964.

نظرية هوهنبرج و كوهين

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

التقريبات المستعملة

كما هو موضح لاحقا فان نظرية الكثافة الوظيفية نظرية دقيقة تماما تدور في معادلات Khon Sham، ; (وبصرف النظر عن تقريب Born-Oppenheimer) عندما نقيس الكثافة الإلكترونية فإنها تمثل بدقة كثافة نظام مكون من N إلكترون في حالة تفاعل . ومع ذلك ، تبقى DFT غير قابل للتطبيق لان كمون التبادل-الارتباط (الذي يحتوي أيضا على تصحيح للطاقة الحركية) ما زال مجهولا.ولذلك من الضروري إيجاد أحسن تقريب لكمون التبادل-الارتباط لكي نستطيع العمل لأن كمون التبادل-الارتباط يمثل قيمة ضعيفة حوالي 20mh/é ولكنه مسئول في الحساب الكمي.
يوجد بشكل عام نوعين من التقريب :( تقريب الكثافة المحلية Local Density Approximation LDA)و تقريب الانحدار المعمم Generalized gradient approximations