תורת הקירוב המקומית של הצפיפות

מתוך המכלול, האנציקלופדיה היהודית
קפיצה לניווט קפיצה לחיפוש

קירוב הצפיפות המקומית (באנגלית: Local Density Approximation ובראשי תיבות LDA) היא שיטה חישובית בכימיה קוונטית, פיזיקת המצב המוצק ומדע החומרים שפותחה על ידי הפיזיקאים וולטר קוהן ולו ז'ו שאם (אנ') בשנות ה-60 של המאה ה-20.

רקע כללי

שיטות תורת הקירוב המקומית הן קבוצה של קירובים פונקציונליים לאנרגיית ההחלפה-קורלציה (Exchange-Correlation ובראשי תיבות XC) בתורה הכללית של פונקציונל הצפיפות. קירובים אלו תלויים אך ורק בצפיפות המטען המרחבית, המסומנת הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho(\bold{r})} .

ישנן גישות שונות לסוג כזה של קירובים, אך הנפוצה ביותר זו הנגזרת ממודל הג'ליום (Jellium), בו מניחים שצפיפות האלקטרונים אחידה בכל המערכת. קירובים אלו, המבוססים על מודל הג'ליום, מאפשרים לחשב פונקציונליים כמו החלפה וקורלציה בצורה מדויקת, דבר שאינו אפשרי באמצעות תורת פונקציוניל הצפיפות הכללית.

בצורתה הכללית, עבור מערכות עם ספין לא מקוטב, ניתן לכתוב את אנרגיית ההחלפה-קורלציה בעזרת תורת הקירוב המקומית של הצפיפות באופן הבא:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle E^{\text{LDA}}_{\text{xc}}[\rho] = \int \rho(\bold{r})\epsilon_{\text{xc}}(\rho (\bold{r})) d\bold{r}}

כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho(\bold{r})} היא צפיפות המטען בנקודה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bold{r}} במרחב ו- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon_{\text{xc}}} היא אנרגיית ההחלפה-קורלציה, המהווה פונקציה של צפיפות המטען.

ניתן להכליל קירוב זה בצורה ישירה על ידי הכללת אינטראקציות ספיניות, אזי הוא נקרא Local Spin-Density Approximation:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\text{xc}}^{\text{LSDA}}[\rho_\uparrow:\rho_\downarrow]= \int \rho(\bold{r})\epsilon_{\text{xc}}(\rho_\uparrow,\rho_\downarrow)d\bold{r}}

אם לוקחים בחשבון גם את גרדיאנט הצפיפות, מקבלים את הקירוב המקומי הנוסף שנקרא Generalized Gradient Approximations:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\text{xc}}^{\text{GGA}}[\rho_\uparrow:\rho_\downarrow]= \int \rho(\bold{r})\epsilon_{\text{xc}}(\rho_\uparrow,\rho_\downarrow,\nabla\rho_\uparrow,\nabla\rho_\downarrow)d\bold{r}}

אנרגיית ההחלפה-קורלציה ניתנת לכתיבה בצורה הבאה:

הפענוח נכשל (שגיאת המרה. השרת ("https://wikimedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle E_{\text{xc}}^{\text{LDA}}=E_{\text{x}}^{\text{LDA}}+E_{\text{c}}^{\text{LDA}}}

יישומים

תורת הקירוב המקומית של הצפיפות, ובמיוחד הגרסה שמתייחסת לגרדיאנט הצפיפות (GGA), מצויה בשימוש רב בקרב פיזיקאים בתחומי המצב המוצק על מנת לבחון אינטראקציות אלקטרוניות ומגנטיות בחומרים מוליכים למחצה, ובכלל זה תחמוצות מוליכות למחצה. בנוסף, משתמשים בתורה זו יחד עם תוכנות סימולציה בכדי לחזות את רמת פרמי ומבנה הפסים בתחמוצות מוליכות למחצה מסוימות.[1][2][3]

תורת הקירוב המקומית של הצפיפות עובדת בצורה מדויקת למדי עבור מגוון רחב של חומרים. בפרט, היא שימושית עבור חישובים של משוואות מצב, קבועים אלסטיים, מבנים אלקטרוניים ותכונות פיזיקליות נוספות של החומר.[4]

גז אלקטרונים הומוגני

הקירוב עבור הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon_{\text{xc}}} תלוי רק בצפיפות וניתן לפתח אותו בדרכים רבות. הגישה המוצלחת ביותר מבוססת על מודל גז האלקטרונים ההומוגני, שנבנה על ידי הצבת הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle N} אלקטרונים בעלי אינטראקציה בתוך נפח הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle V} , עם מטען רקע חיובי השומר על המערכת נייטרלית. לאחר מכן נלקחים הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle N} ו-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle V} לאינסוף, אך באופן כזה ששומר על הצפיפות (הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho=\frac{N}{V}} ) סופית. זהו קירוב שימושי, שכן האנרגיה הכוללת מורכבת רק מהאנרגיה הקינטית, אנרגיית האינטראקציה האלקטרוסטטית ואנרגיית החלפה שנובעת מהאינטראקציה בין האלקטרונים. פונקציית הגל ניתנת לביטוי במונחים של גלים מישוריים. בפרט, עבור צפיפות קבועה, צפיפות אנרגיית ההחלפה מתכונתית ל-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho^{1/3}} .

פונקציונל ההחלפה

כשדנים באנרגיית ההחלפה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\text{x}}^{\text{LDA}}} במודל הג'ליום, ישנן גם משואות אנליטיות פשוטות יחסית, שניתן לכתוב כך:[5][6]

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\text{x}}^{\text{LDA}}[\rho] = -\frac{3}{4}\left ( \frac{3}{\pi} \right) ^{\frac{1}{3}}\int \rho(\bold{r})^{\frac{4}{3}}d\bold{r}}

פונקציונל הקורלציה

לאנרגיית הקורלציה יש פתרון מדויק שמתקיים רק בגבולות של צפיפות גבוהה או נמוכה מאוד. בגבול של צפיפות גבוהה, אפשר לכתוב את אנרגיית הקורלציה כך:[5]

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon_{C}=A\ln(r_{s})+B+r_{s}(C\ln(r_{s})+D)}

ובגבול של צפיפות נמוכה:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon_C = \frac {1}{2}\left(\frac {g_0}{r_{s}} + \frac {g_0}{r_{s}^{\frac{3}{2}}} + \cdots\right)}

כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle r_{s}} , המכונה פרמטר ויגנר-זייץ, הוא פרמטר חסר-ממדים שמוגדר כרדיוס של כדור שמקיף אלקטרון בודד, מחולק ברדיוס בוהר.[7] פרמטר ויגנר-זייץ קשור לצפיפות עצמה:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{4}{3}\pi r_{\rm s}^{3} = \frac{1}{\rho}}

כמו כן, הוצע ביטוי אנליטי המשמש עבור כל טווח הצפיפויות, והמבוסס על תורת ההפרעות הרב-גופית. אנרגיות הקורלציה המחושבות עולות בקנה אחד עם התוצאות של סימולציות מונטה קרלו קוונטיות עד כדי שגיאה של 2 מילי-הרטרי. סימולציות מונטה קרלו מדויקות עבור האנרגיה של גז אלקטרונים הומוגני, מבוצעות עבור מספר ערכי ביניים של הצפיפות ומספקות בתורן ערכים מדויקים של צפיפות אנרגיית הקורלציה.[8]

קיטוב ספין

ההרחבה של פונקציונל הצפיפות גם למערכות מקוטבות-ספין היא מיידית עבור אנרגיית ההחלפה, שבה ה-scaling ע"פ הספין הוא ידוע, אולם עבור אנרגיית הקורלציה יש להשתמש בקירובים נוספים. בתורת פונקציונל הצפיפות, מערכת מקוטבת-ספין עושה שימוש בשתי צפיפויות ספין, הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{\alpha}} ו- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{\beta}} , כאשר הצפיפות הכוללת היא סכום שלהן. הצורה המתמטית של ה- Local Spin-Density Approximation תהיה במקרה זה:[9]

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\rm xc}^{\mathrm{LSDA}}[\rho_{\alpha},\rho_{\beta}] = \int\mathrm{d}\mathbf{r}\ \rho(\mathbf{r})\epsilon_{\rm xc}(\rho_{\alpha},\rho_{\beta})\ }

באשר לאנרגיית ההחלפה, התוצאה הסופית (ולא רק עבור קירובי צפיפות מקומית) ידועה במונחי הפונקציונל שלא כולל קיטוב ספין:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\rm x}[\rho_{\alpha},\rho_{\beta}] = \frac{1}{2}\bigg( E_{\rm x}[2\rho_{\alpha}] + E_{\rm x}[2\rho_{\beta}] \bigg)\ }

התלות הספינית של צפיפות אנרגיית הקורלציה ניתנת לכתיבה על ידי הגדרת קיטוב-הספין היחסי:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \zeta(\mathbf{r}) = \frac{\rho_{\alpha}(\mathbf{r})-\rho_{\beta}(\mathbf{r})}{\rho_{\alpha}(\mathbf{r})+\rho_{\beta}(\mathbf{r})}\ }

ניתן לראות כי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \zeta = 0\,} מייצג מצב דיאמגנטי של אי-קיטוב-ספין, שבו הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha\,} ו-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta\,} הן צפיפויות ספין זהות, בעוד הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \zeta = \pm 1} מייצג מצב פרומגנטי שבו אחת מצפיפויות הספין איננה. צפיפות אנרגיית הקורלציה הספינית עבור ערכים נתונים של הצפיפות הכוללת והקיטוב היחסי, המסומנת כ-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon_{c}(\rho,\zeta)} , מחושבת כך שיהיה ניתן לבצע אינטרפולציה בין ערכי קצה. צורות מסוימות של התלות הזו פותחו במקביל לפונקציונלי הקורלציה תחת קירוב הצפיפות המקומית.[10]

פוטנציאל קורלציה-החלפה

פוטנציאל הקורלציה-החלפה מייצג את אנרגיית הקורלציה-החלפה בקירוב הצפיפות המקומית, והוא נתון ע"י:[5]

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle v_{\rm xc}^{\mathrm{LDA}}(\mathbf{r}) = \frac{\delta E^{\mathrm{LDA}}}{\delta\rho(\mathbf{r})} = \epsilon_{\rm xc}(\rho(\mathbf{r})) + \rho(\mathbf{r})\frac{\partial \epsilon_{\rm xc}(\rho(\mathbf{r}))}{\partial\rho(\mathbf{r})}\ }

במערכות סופיות, הפוטנציאל הזה דועך אסימפטוטית באופן מעריכי. זוהי תוצאה שגויה; פוטנציאל ההחלפה-קורלציה האמיתי דועך לאט בהרבה, באופן קולומבי (כלומר היפרבולי). הדעיכה המהירה המלאכותית מביעה עצמה במספר של אורביטלי קוהן-שאם שהפוטנציאל יכול לקשור (משמע, האורביטלים שהאנרגיה שלהם שלילית). פוטנציאל זה לא יכול לתמוך בסדרת רידברג, והמצבים שהוא כן קושרים גבוהים מדי באנרגיה. הדבר מביא לכך שאנרגיית אורביטל ה-HOMO גבוהה מדי, ולכן כל חיזוי לגבי אנרגיית היינון כנובע ממשפט קופמנס, אינו מדויק. מעבר לכך, הקירוב מספר תיאור לא מדויק לצורונים עשירים באלקטרונים כגון אניונים, שבהם לעיתים קרובות לא ניתן לקשור אלקטרון נוסף; הקירוב גורס כי צורונים אלה יהיו בלתי יציבים.[11] במקרה של קיטוב-ספין, פוטנציאל ההחלפה-קורלציה דורש אינדקסי ספין. אולם, אם רק מחשיבים את אנרגיית ההחלפה, ניתן לקבל פוטנציאל המלוכסן באינדקסי הספין:[12]

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle v_{\rm xc, \alpha \beta}^{\mathrm{LDA}}(\mathbf{r}) = \frac{\delta E^{\mathrm{LDA}}}{\delta\rho_{\alpha \beta}(\mathbf{r})} = \frac{1}{2}\delta_{\alpha\beta}\frac{\delta E^{\mathrm{LDA}}[2\rho_{\alpha}]}{\delta \rho_{\alpha}} = - \delta_{\alpha\beta}\Big(\frac{3}{\pi}\Big)^{1/3}2^{1/3}\rho_{\alpha}^{1/3}}

הערות שוליים

  1. ↑ Segall, M.D.; Lindan, P.J (2002). "First-principles simulation: ideas, illustrations and the CASTEP code". Journal of Physics: Condensed Matter. 14 (11): 2717. Bibcode:2002JPCM...14.2717S. doi:10.1088/0953-8984/14/11/301. S2CID 250828366.
  2. ↑ Zhao, G. L.; Bagayoko, D.; Williams, T. D. (1999-07-15). "Local-density-approximation prediction of electronic properties of GaN, Si, C, and RuO2". Physical Review B. 60 (3): 1563–1572. Bibcode:1999PhRvB..60.1563Z. doi:10.1103/physrevb.60.1563. ISSN 0163-1829.
  3. ↑ Bagayoko, Diola (בדצמבר 2014). "Understanding density functional theory (DFT) and completing it in practice". AIP Advances. 4 (12): 127104. Bibcode:2014AIPA....4l7104B. doi:10.1063/1.4903408. ISSN 2158-3226. {{cite journal}}: (עזרה)
  4. ↑ Assadi, M.H.N; et al. (2013). "Theoretical study on copper's energetics and magnetism in TiO2 polymorphs". Journal of Applied Physics. 113 (23): 233913–233913–5. arXiv:1304.1854. Bibcode:2013JAP...113w3913A. doi:10.1063/1.4811539. S2CID 94599250.
  5. ^ 5.0 5.1 5.2 Parr, Robert G; Yang, Weitao (1994). Density-Functional Theory of Atoms and Molecules. Oxford: Oxford University Press. ISBN 978-0-19-509276-9.
  6. ↑ Dirac, P. A. M. (1930). "Note on exchange phenomena in the Thomas-Fermi atom". Proc. Camb. Phil. Soc. 26 (3): 376–385. Bibcode:1930PCPS...26..376D. doi:10.1017/S0305004100016108.
  7. ↑ Murray Gell-Mann and Keith A. Brueckner (1957). "Correlation Energy of an Electron Gas at High Density" (PDF). Phys. Rev. 106 (2): 364–368. Bibcode:1957PhRv..106..364G. doi:10.1103/PhysRev.106.364. S2CID 120701027.
  8. ↑ D. M. Ceperley and B. J. Alder (1980). "Ground State of the Electron Gas by a Stochastic Method". Phys. Rev. Lett. 45 (7): 566–569. Bibcode:1980PhRvL..45..566C. doi:10.1103/PhysRevLett.45.566. S2CID 55620379.
  9. ↑ Oliver, G. L.; Perdew, J. P. (1979). "Spin-density gradient expansion for the kinetic energy". Phys. Rev. A. 20 (2): 397–403. Bibcode:1979PhRvA..20..397O. doi:10.1103/PhysRevA.20.397.
  10. ↑ von Barth, U.; Hedin, L. (1972). "A local exchange-correlation potential for the spin polarized case". J. Phys. C: Solid State Phys. 5 (13): 1629–1642. Bibcode:1972JPhC....5.1629V. doi:10.1088/0022-3719/5/13/012.
  11. ↑ Fiolhais, Carlos; Nogueira, Fernando; Marques Miguel (2003). A Primer in Density Functional Theory. Springer. p. 60. ISBN 978-3-540-03083-6.
  12. ↑ Giustino, Feliciano (2014). Materials Modelling Using Density Functional Theory: Properties and Predictions. Oxford University Press. p. 229.
הערך באדיבות ויקיפדיה העברית, קרדיט,
רשימת התורמים
רישיון cc-by-sa 3.0

תורת הקירוב המקומית של הצפיפות36410447Q898241

אוחזר מתוך "https://www.hamichlol.org.il/w/index.php?title=תורת_הקירוב_המקומית_של_הצפיפות&oldid=3279240"