نظرية الكثافة الوظيفية
نظرية الدالة الوظيفية للكثافة (DFT, Density Functional Theory) هي أحد أهم الطرق المستعملة في الفيزياء و الكيمياء النظريين و بواسطتها نستطيع أن نحدد خصائص نظام متعدد الجسيمات (الطاقة الكلية للنظام, الكثافة الإلكترونية للمدارات, المعاملات الفيزيائية والضوئية للمادة...)، وهي واحدة من أكثر الطرق استخداما في العمليات الحسابية الكمومية بسبب إمكانية تطبيقها على أنظمة متنوعة و بتكلفة و سرعة العالية.
تعتمد الطرق التقليدية في حلها معادلة شرودنغر لنظام مكون من عدة ذرات ولا سيما نظرية هارتري-فوك و الطرق المستنبطة منها على الدالة الموجية ذات 3N متغير (حيث N هو العدد الإجمالي لجزيئات النظام)، لذلك تكون لمعادلات المراد حلها جد معقدة و تتطلب جهدا كبيرا· الهدف الرئيسي نظرية الرئيسي نظرية الدالة الوظيفية للكثافة هو استبدال الدالة الموجية بالدالة الوظيفية للكثافة ذات ثلاثة متغيرات فقط و جعلها كقاعدة للحساب، لذلكفالتعامل معها كمفهوم رياضي أو فيزيائي هو أسهل بكثير ، فمبدأ DFT هو إعادة صياغة للمسألة الكمومية وتحويلها من مسألة لنظام متعددالجسيمات إلى مسألة أحادية الجسيمة .
تعود جذور DFT للنموذج الذي وضعه لويلين توماس و انريكو فيرمي في أواخر 1920. ومع ذلك فإنه لم يكن ب الإمكان إستعمالها حتى منتصف 1960. ومنذ ذلك الوقت و هي في تطور متصاعد. مع تطور وقوة البرامج المعلوماتية وسرعة تطبيق العمليات الحسابية. جاءت نتائج (DFT) النظرية قريبة على نحو مرض تماما مع البيانات التجريبية وبتكاليف منخفضة نسبيا مع الطرق التقليدية التي تستهلك المال والوقت معا, ولكن على الرغم من التحسن الكبير الذي طرأ مؤخرا والتحسين المستمر للبرامج، لا تزال هناك بعض الصعوبات في استخدام نظرية الكثافة الوظيفية (DFT) لوصف التفاعلات بين الجزيئات وخاصة القوى الضعيفة (فان دير وال), الحالات الانتقالية للإلكترونات, حساب الفجوة (band Gap) في أشباه الموصلات, هذه الثغرات يمكن أن تؤثر على نتائج (DFT) "على الأقل عندما تستعمل وحدها دون تصحيح" حيث يتم تطوير طرق جديدة ل (DFT) للتغلب على هذه المشكلة, من خلال إحداث تعديلات على الدالة الوظيفية أو إدراج شروط في برامج الحساب وهو موضوع البحث الحالي في هذا المجال.
"إن القوانين الفيزيائية الأساسية اللازمة للتعبير الرياضي لجزء مهم من الفيزياء و الكيمياء كلها معروفة سلفا، وتكمن الصعوبة فقط في التطبيق الدقيق لهذه القوانين الذي يؤدي إلى المعادلات جد معقدة يتعين حلها " عن - بول ديراك 1929
المعادلة الأساسية التي يتعين حلها لوصف نظام يتكون مثلا من العديد من الجزيئات أو الذرات هي معادلة إرفين شرودنغر (1887-1961) المعروفة بمعادلة شرودنجر و تكتب على الشكل التالي :
بحيث formula_1 هو المؤثر الهاملتوني (Hamiltonian operator)الكمي للمجموعة قيد الدرس والمكونة من formula_2 إلكترون ( ذات الأساس formula_3 ) و formula_4 لنوى الذرات ( ذات الأساس formula_5 ). الحدين الأولين من المؤثر الهاملتوني هما على التوالي هو مؤثر الطاقة الحركية للإلكترونات و نوى الذرات على التوالي، أما الحدود الثلاثة المتبقية فهي تمثل مؤثرات الطاقة الوضعية أو الجهد لمختلف التفاعلات البنية ( الإلكترون-النواة ) و ( الإلكترون -الإلكترون ) و ( النواة- النواة ).
حل هذه المعادلة على حالتها هذه معقدة جدا، ولكن وباستخدام تقريب بورن-أوبنهايمر الذي وضعه كل من ماكس بورن و روبرت أوبنهايمر أصبح ممكنا إلى حد ما. و هذا التقريب ينص على إمكانية إهمال حركة النوى مقارنة مع حركة الإلكترونات. و ذلك راجع للكتلة النواة العالية و سرعتها البطيئة مقارنة بكتلة و سرعة الإلكترونات. هذا يمكننا من عزل حركة الإلكترونات عن حركة النوى و تكون النتيجة الابقاء على الحدود التي تحتوي تمثل طاقة الالكترون الحركية والوضعية بالاضافة إلى طاقة الوضع الناتج عن تفاعل الالكترون مع النواة. وبهذا نكون قد اقتصرنا على حل معادلة شرودنغر للإلكترونات . ويمكننا إعادة كتابة المعادلة اعلاه على الشكل التالي :
سنقتصر في دراستنا المتبقية على دراسة formula_6 و الذي يحدد معادلة شرودنغر المتعددة الإلكترونات:
وقد تم تطوير العديد من طرق لحل هذه المعادلة لكنها كانت تعتمد على الدوال الموجية كطريقة هرتري فوك . توفر DFT طريقة بديلة تعتبر كثافة الإلكترونات هي الأساس لمعرفة خصائص المجموعة المدروسة.
كائن رياضي يمثل علاقة رياضية نربط بها دالة بدالة أخرى عكس مفهوم الد آلة العادية التي نربط بها دالة بعدد ذو قيمة معينة. فمثلا الدالة الموجية formula_7 هي دالة عادية، أما الطاقة formula_8 فيمكن التعبير عنها بدلالة formula_9 فتسمى دالة وظيفية
احتمال العثور على إلكترون متمركز في الموضع formula_10 من بين العدد الإجمالي للإلكترونات formula_2 في حيز حجمه formula_12 هو :
حيث formula_13 هو احتمال الكثافة الإلكترونية والذي يعرف بأنه :
كما هو موضح لاحقا فان نظرية الكثافة الوظيفية نظرية دقيقة تماما تدور في معادلات Khon Sham، ; (وبصرف النظر عن تقريب Born-Oppenheimer) عندما نقيس الكثافة الإلكترونية فإنها تمثل بدقة كثافة نظام مكون من N إلكترون في حالة تفاعل . ومع ذلك ، تبقى DFT غير قابل للتطبيق لان كمون التبادل-الارتباط (الذي يحتوي أيضا على تصحيح للطاقة الحركية) ما زال مجهولا.ولذلك من الضروري إيجاد أحسن تقريب لكمون التبادل-الارتباط لكي نستطيع العمل لأن كمون التبادل-الارتباط يمثل قيمة ضعيفة حوالي 20mh/é ولكنه مسئول في الحساب الكمي.
يوجد بشكل عام نوعين من التقريب :( تقريب الكثافة المحلية Local Density Approximation LDA)و تقريب الانحدار المعمم Generalized gradient approximations