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

نظرية الدالة الوظيفية للكثافة (DFT, Density Functional Theory) هي أحد أهم الطرق المستعملة في الفيزياء والكيمياء النظريين وبواسطتها نستطيع أن نحدد خصائص نظام متعدد الجسيمات (الطاقة الكلية للنظام، الكثافة الإلكترونية للمدارات، المعاملات الفيزيائية والضوئية للمادة....)، وهي واحدة من أكثر الطرق استخداما في العمليات الحسابية الكمومية بسبب إمكانية تطبيقها على أنظمة متنوعة وبتكلفة وسرعة العالية.

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

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

أساسيات

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

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

المعادلة الأساسية التي يتعين حلها لوصف نظام يتكون مثلا من العديد من الجزيئات أو الذرات هي معادلة إرفين شرودنغر (1887-1961) المعروفة بمعادلة شرودنجر وتكتب على الشكل التالي:

$H\mathrm{\Psi }=[-\sum _i^N\frac{{\hslash }^2}{2m}\underset{i}{\overset{2}{\mathrm{\nabla }}}-\sum _I^A\frac{{\hslash }^2}{2M}\underset{I}{\overset{2}{\mathrm{\nabla }}}-{\sum }_{i,I}\frac{Z_I{e}^2}{|{\overrightarrow{r}}_i-{\overrightarrow{R}}_I|}+{\sum }_{i<j}\frac{e^2}{|{\overrightarrow{r}}_i-{\overrightarrow{r}}_j|}+{\sum }_{I<J}\frac{Z_I{Z}_J{e}^2}{|{\overrightarrow{R}}_I-{\overrightarrow{R}}_J|}]\mathrm{\Psi }=E\mathrm{\Psi }$

بحيث $H$ هو المؤثر الهاملتوني (Hamiltonian operator)الكمي للمجموعة قيد الدرس والمكونة من $N$ إلكترون (ذات الأساس $i$) و $A$ لنوى الذرات (ذات الأساس $I$). الحدين الأولين من المؤثر الهاملتوني هما على التوالي هو مؤثر الطاقة الحركية للإلكترونات ونوى الذرات على التوالي، أما الحدود الثلاثة المتبقية فهي تمثل مؤثرات الطاقة الوضعية أو الجهد لمختلف التفاعلات البنية (الإلكترون-النواة) و (الإلكترون -الإلكترون) و (النواة- النواة).

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

$H\mathrm{\Psi }=[H^e+E_A]\mathrm{\Psi }=E\mathrm{\Psi }$

سنقتصر في دراستنا المتبقية على دراسة $H^e$ والذي يحدد معادلة شرودنغر المتعددة الإلكترونات:

$H^e\mathrm{\Psi }=[T+V_{ext}+U]\mathrm{\Psi }=E^e\mathrm{\Psi }$

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

تعاريف

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

كائن رياضي يمثل علاقة رياضية نربط بها دالة بدالة أخرى عكس مفهوم الد آلة العادية التي نربط بها دالة بعدد ذو قيمة معينة. فمثلا الدالة الموجية $\mathrm{\Psi }(\overrightarrow{r})$ هي دالة عادية، أما الطاقة $E$ فيمكن التعبير عنها بدلالة $E[\mathrm{\Psi }({\overrightarrow{r}}_i)]$ فتسمى دالة وظيفية

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

احتمال العثور على إلكترون متمركز في الموضع $\overrightarrow{r}$ من بين العدد الإجمالي للإلكترونات $N$ في حيز حجمه $d\overrightarrow{r}$ هو:

$n(r)d\overrightarrow{r}$

حيث $n(r)$ هو احتمال الكثافة الإلكترونية والذي يعرف بأنه:

$n(\overrightarrow{r})=N\int |\mathrm{\Psi }(\overrightarrow{r}s,{\overrightarrow{r}}_2{s}_2,\dots ,{\overrightarrow{r}}_N{s}_N){|}^{2}dsd{s}_{2}d\overrightarrow{r_2}\dots d{s}_{N}d\overrightarrow{r_N}$

احتمال الكثافة الإلكترونية له خاصيتين أساسيتين هما:

1. $n(\overrightarrow{r}\to \mathrm{\infty })=0$
2. $\int n(\overrightarrow{r})d\overrightarrow{r}=N$

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

كثافة زوج من الإلكترونات هو احتمال العثور في وقت واحد على اثنين من بين $N$ من الإلكترونات والمتموضعين في وحدتي الحجم$d\overrightarrow{r}$ و $d{{\overrightarrow{r}}^{\prime }}^{\prime }$ هي: $n(\overrightarrow{r},\overrightarrow{r^{\prime }})drd{r}^{\prime }$

$n(\overrightarrow{r},\overrightarrow{r^{\prime }})drd{r}^{\prime }=N(N-1)\int |\mathrm{\Psi }(\overrightarrow{r}s,\overrightarrow{r^{\prime }}{s}^{\prime },{\overrightarrow{r}}_3{s}_3,\dots ,{\overrightarrow{r}}_N{s}_N){|}^{2}dsd{s}^{\prime }d{s}_{3}d\overrightarrow{r_3}\dots d{s}_{N}d\overrightarrow{r_N}$

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

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

الأسس الرياضية

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

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

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

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

الشق الأول للنظرية

إذا ما اعتبرنا غازا مكونا من العديد من الإلكترونات. فإن الجهد خارجي $V_{ext}(r)$ الذي تخضع له هو الذي يحدد الحالة الدنيا (الأكثر استقرارا) والكثافة الإلكترونية $n_0(r)$ لهذا الغاز. البرهان الرياضي للبرهنة على هذه النظرية نستعمل طريقة الإستدلال بالترجع. حيث نفترض أنه يوجد لدينا جهدان خارجيان مختلفان $V_{ext}^1(r)$ و $V_{ext}^2(r)$ ( الإختلاف هاهنا نقصد به أن فرقهما يساوي عددا معينا ) يرتبطان بالكثافة الدنيا $n(r)$ ( أصغر كثافة ممكنة تكون عندها المجموعة المدروسة مستقرة ). الجهدان يقودان إلى تحديدعاملان هاملتونيان مختلفين $H^{(1)}$ و $H^{(2)}$، بحيث تصف كل من الدالتين الموجتين $|{\mathrm{\Psi }}^{(1)}⟩$ و $|{\mathrm{\Psi }}^{(1)}⟩$ حالتين دنيويتن مختلفتين. وبما أن $|{\mathrm{\Psi }}^{(2)}⟩$ لا تصف الحالة الدنيا ل $H^{(1)}$ فانه يمكننا كتابة: $E^{(1)}=⟨{\mathrm{\Psi }}^{(1)}\left|H^{(1)}\right|{\mathrm{\Psi }}^{(1)}⟩<⟨{\mathrm{\Psi }}^{(2)}\left|H^{(1)}\right|{\mathrm{\Psi }}^{(2)}⟩$ هذه المتراجحة لا تكون صحيحة إلا إذا كانت مستويات الطاقة للإلكترونات غير منحلة وهذا الشرط محقق في نظرية هوهنبرج و كوهين. الحد الثاني من المتراجحة يساوي: $⟨{\mathrm{\Psi }}^{(2)}\left|H^{(1)}\right|{\mathrm{\Psi }}^{(2)}⟩=⟨{\mathrm{\Psi }}^{(2)}\left|H^{(2)}\right|{\mathrm{\Psi }}^{(2)}⟩+⟨{\mathrm{\Psi }}^{(2)}|H^{(1)}-H^{(2)}|{\mathrm{\Psi }}^{(2)}⟩$ $=E^{(2)}+\int [V_{ext}^1(r)-V_{ext}^2(r)]n_0(r)d^{3}r$ هذا يعطينا: $E^{(1)}<E^{(2)}+\int [V_{ext}^1(r)-V_{ext}^2(r)]n_0(r)d^{3}r$ بنفس المنهج و بتعويضنا ل$E^2$ ب$E^1$ يمكننا الحصول على : $E^{(2)}<E^{(1)}+\int [V_{ext}^2(r)-V_{ext}^1(r)]n_0(r)d^{3}r$ نجمع حدود المتراجحتن على حدا فنجد : $E^{(1)}+E^{(2)}<E^{(1)}+E^{(2)}$ و هذا تعارض يوضح أن إفترضنا الأول خاطئ ، ومنه نستنتج أن جهدان مختلفان لايقودان بأي حال من الأحوال إلى كثافة واحدة. الرسم البياني أدناه يوضح لنا جليا إمكانية تحديد جميع خصائص النظام تماما إذا كانت كثافة الدنيا للإلكترونات معروفة. $\begin{array}{ccc}{V}_{ext}(r)& ⟸& n_0(r)\\ \downarrow & & \uparrow \\ {\mathrm{\Psi }}_i(\{r\})& \to & {\mathrm{\Psi }}_0(\{r\})\end{array}$ التقريبات المستعملة كما هو موضح لاحقا فان نظرية الكثافة الوظيفية نظرية دقيقة تماما تدور في معادلات Khon Sham، ; (وبصرف النظر عن تقريب Born-Oppenheimer) عندما نقيس الكثافة الإلكترونية فإنها تمثل بدقة كثافة نظام مكون من N إلكترون في حالة تفاعل. ومع ذلك، تبقى DFT غير قابل للتطبيق لان كمون التبادل-الارتباط (الذي يحتوي أيضا على تصحيح للطاقة الحركية) ما زال مجهولا.ولذلك من الضروري إيجاد أحسن تقريب لكمون التبادل-الارتباط لكي نستطيع العمل لأن كمون التبادل-الارتباط يمثل قيمة ضعيفة حوالي 20mh/é ولكنه مسئول في الحساب الكمي. يوجد بشكل عام نوعين من التقريب:(تقريب الكثافة المحلية Local Density Approximation LDA)و تقريب الانحدار المعمم Generalized gradient approximations مراجع