פרק אחד-עשר: קשר כימי ומולקולות — מכניקה קוונטית בפועל

הבעיה: משוואת שרדינגר לא ניתנת לפתרון מדויק

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

עבור מולקולת מים ($H_{2}O$) — עשרה אלקטרונים. עבור חלבון ממוצע — עשרות אלפי אלקטרונים. הפתרון המדויק אינו רק קשה — הוא בלתי-אפשרי מעשית לכל מולקולה אמיתית.

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

---

הפרדת בורן-אופנהיימר: ראשית כל, הקפיאו את הגרעינים

הצעד הראשון והחשוב ביותר הוא קירוב בורן-אופנהיימר (Born-Oppenheimer approximation, 1927).

הרעיון: פרוטון כבד מאלקטרון פי כ-1836. לכן — גרעינים נעים לאט בהרבה מאלקטרונים. מבחינתם של האלקטרונים, הגרעינים כמעט עומדים במקומם. מבחינתם של הגרעינים, האלקטרונים “כבר מצאו את מצבם” ברגע שהגרעינים זזו קצת.

זה מאפשר להפריד את הבעיה לשני שלבים:

1. בעיה אלקטרונית: עבור כל סידור קבוע של גרעינים, פותרים את משוואת שרדינגר לאלקטרונים בלבד ומקבלים את האנרגיה האלקטרונית.
2. בעיה גרעינית: הגרעינים נעים על משטח פוטנציאל (potential energy surface, PES) שנקבע על-ידי האנרגיה האלקטרונית בכל נקודה.

משטח הפוטנציאל הוא מושג מרכזי בכימיה: הוא “הנוף” שעליו מתרחשות תגובות כימיות. שקע מינימלי — מבנה יציב. עמק בין שני שקעים — מצב מעבר (transition state). גובה העמק מעל המינימום — אנרגיית השפעול. כל זה, מבחינה מתמטית, הוא פונקציה של קואורדינטות הגרעינים בלבד — ובורן-אופנהיימר הוא מה שמאפשר לחשב אותו.

---

הרטרי-פוק: כל אלקטרון בשדה הממוצע של האחרים

לאחר שהפרדנו את הגרעינים, נשארת עדיין הבעיה האלקטרונית. אפילו היא לא ניתנת לפתרון מדויק עבור יותר מאלקטרון אחד, בגלל האינטראקציות בין האלקטרונים.

שיטת הרטרי-פוק (Hartree-Fock, HF) היא הקירוב הפשוט ביותר שעדיין לוקח ברצינות את עיקרון האיסור של פאולי. הרעיון:

במקום לטפל בכל האינטראקציות בין אלקטרונים במדויק, כל אלקטרון נע בשדה ממוצע הנוצר על-ידי כל האלקטרונים האחרים. הבעיה ה”רב-גופית” (many-body) הופכת לאוסף של בעיות חד-גופיות (one-body) — כל אחת פותרת תנועה של אלקטרון יחיד בפוטנציאל יעיל.

כדי לקיים את עיקרון פאולי (שני אלקטרונים לא יכולים להיות באותו מצב), פונקציית הגל הרב-אלקטרונית כתובה כדטרמיננטת סלייטר (Slater determinant) — מכפלה מסוימת של פונקציות חד-אלקטרוניות שמבטיחה אנטי-סימטריה.

הרטרי-פוק נותן תוצאות סבירות לגיאומטריות מולקולות ותדירויות רטט, אך מחטיא גורם חשוב: קורלציית אלקטרונים. אלקטרונים אינם באמת עצמאיים — כל אחד “יודע” בזמן אמת היכן האחרים. שדה ממוצע הוא, בדיוק, ממוצע — ומחמיץ את הדינמיקה של הזוגות. שגיאת הקורלציה הזו קריטית עבור קשרים כימיים, מצבי עירור, ותכונות מגנטיות.

---

DFT: הגאונות של הצפיפות

ב-1964 פרסמו וולטר קון (Walter Kohn) ופייר הוהנברג (Pierre Hohenberg) תוצאה מפתיעה:

כל התכונות של מערכת אלקטרונים בשיווי משקל נקבעות באופן יחיד על-ידי צפיפות האלקטרונים $n(\mathbf{r})$ — פונקציה של שלושה משתנים בלבד.

זהו משפט הוהנברג-קון הראשון. המשמעות: במקום לחפש פונקציית גל של $N$ אלקטרונים — שתלויה ב-$3N$ קואורדינטות ועשויה להיות פונקציה של מאות-אלפי משתנים — מספיק לדעת את $n(\mathbf{r})$: כמה אלקטרונים יש, בממוצע, בכל נקודה במרחב. שלושה משתנים, לא $3N$.

המשפט השני של הוהנברג-קון אומר שהאנרגיה הכוללת היא פונקציונל של הצפיפות — $E[n]$ — ושצפיפות מצב היסוד (ground state) היא זו שממזערת אנרגיה זו. זה עיקרון וריאציוני: ה”צפיפות הנכונה” היא זו שנותנת את האנרגיה הנמוכה ביותר.

משוואות קון-שם

ב-1965 המשיכו קון ולו-ג’או שם (Lu Jeu Sham) ותרגמו את הרעיון לשיטה חישובית מעשית: משוואות קון-שם (Kohn-Sham equations).

הרעיון: מדמיינים מערכת היפותטית של אלקטרונים ללא אינטראקציות, שהצפיפות שלה זהה לצפיפות של המערכת האמיתית. עבור מערכת כזו, הבעיה פתירה בדיוק. את כל מה שלא יודעים — כל הפיסיקה של האינטראקציות האמיתיות — דוחסים לתוך פונקציה אחת שנקראת פונקציונל החילוף-קורלציה (exchange-correlation functional), $E_{xc}[n]$.

בפועל: $E_{xc}$ אינו ידוע במדויק — הוא הוא “הלא-ידוע” היחיד בכל שיטת DFT. כל ה”שיפורים” ב-DFT לאורך עשרות השנים האחרונות הם, בעיקרם, שיפורים בקירוב של $E_{xc}$.

סולם ג’ייקוב — הפונקציונלים

בקהילת DFT, מקובל לדבר על “סולם ג’ייקוב” (Jacob’s ladder) של פונקציונלים — כל שלב גבוה יותר, מדויק יותר, ויקר יותר חישובית:

• LDA (Local Density Approximation): $E_{xc}$ תלוי רק בצפיפות המקומית בכל נקודה. פשוט, מהיר, לא מדויק מספיק לכימיה אורגנית.
• GGA (Generalized Gradient Approximation): מוסיפים גם את גרדיינט הצפיפות — שיפור משמעותי. הפונקציונל PBE הוא ה-GGA הנפוץ ביותר; אם ראיתם VASP מחשב משהו, כנראה שהשתמשתם ב-PBE.
• Hybrid functionals: מערבים חלק מהחילוף של הרטרי-פוק (מדויק) עם GGA. פונקציונל B3LYP הוא הנפוץ בכימיה אורגנית; HSE נפוץ בחישובי פסים. יקרים יותר, מדויקים יותר.
• מעל לכך: meta-GGA, RPA — שיטות מדויקות עוד יותר, אך יקרות מאוד חישובית, שמורות לבעיות ספציפיות.

בפועל: PBE לגאומטריה, HSE06 לפערי פסים, B3LYP לכימיה אורגנית — זה ה”תרגום לשפת השדה”. אם תגיעו ל-VASP או Quantum ESPRESSO, תדעו עכשיו על מה מדובר כשיגידו לכם “הריצה PBE”.

---

מה DFT יודע לעשות — ומה לא

עובד היטב:

• גאומטריית מולקולות וגבישים: מרחקי קשר, זוויות, פרמטרי תא — בדיוק של ~1% ביחס לניסוי
• אנרגיות תגובה ואנרגיות קשר: לרוב בדיוק של כמה kJ/mol
• ספקטרוסקופיה IR וראמן: תדירויות רטט, דרך חישוב פונונים
• מבנה פסי אנרגיה: לחומרים מוליכים ומוליכי-למחצה עם פונקציונלים מתאימים
• תכונות מגנטיות: מומנטים מגנטיים, סדר מגנטי

עובד פחות טוב (דורש שיטות מתקדמות יותר):

• פערי פסים של מבודדים: LDA/GGA לרוב מזעירים אותם משמעותית; דרוש HSE או GW
• מצבים מעוררים: DFT מחשב מצב יסוד; עבור ספקטרוסקופיה UV-Vis דרוש TD-DFT
• כוחות ון-דר-ולס ודיספרסיה: DFT רגיל מחמיץ אותם; דרושות תיקוני DFT-D
• מערכות מתואמות חזק (strongly correlated, כגון אוקסידים של מתכות מעבר): DFT+U או שיטות מעבר

---

קשר כימי — תמונה מהירה דרך MO

לא ניכנס עמוק לתורת MO (זה נושא לקורסים אחרים), אבל כמה עקרונות שווה להכיר:

כאשר שני אטומים מתקרבים, האורביטלים שלהם מתאחדים: כל זוג AOs (atomic orbitals) מסוג מתאים יוצר זוג MOs — קשר (bonding, $\sigma$, $\pi$) ואנטי-קשר (antibonding, ${\sigma }^{\ast }$, ${\pi }^{\ast }$).

• MO קשר: צפיפות אלקטרונים בין הגרעינים — מושך שניהם, מייצב את הקשר
• MO אנטי-קשר: צפיפות אלקטרונים מחוץ לאזור הבין-גרעיני — מדחה, מערער קשר

סדר קשר = (אלקטרונים בקשר − אלקטרונים באנטי-קשר) / 2

כמה תוצאות מפורסמות שנגזרות ישירות:

• $N_2$: סדר קשר 3 (קשר משולש) — אחד מהקשרים החזקים בכימיה, לכן $N_2$ כל כך אינרטי
• $O_2$: שני האלקטרונים האחרונים ממלאים שני MOs אנטי-קשר מנוונים, כל אחד עם ספין שונה — לכן $O_2$ פארמגנטי, עובדה שלא מתחייבת ממבנה לואיס, אך יוצאת אוטומטית מתורת MO
• $F_2$: קשר בודד חלש, כי שני ${\pi }^{\ast }$ כמעט מבטלים שני $\pi$ קשרים

מהאטום למוצק: הגשר לפרק 12

בפרק 10 ראינו שכאשר $10^{23}$ אטומים נמצאים זה לצד זה, רמות האנרגיה הבדידות מתפצלות לפסים. עכשיו אפשר לנסח זאת בשפת MO:

כשמולקולה דו-אטומית נוצרת — שתי AOs → שתי MOs. כשרשת של $N$ אטומים נוצרת — $N$ AOs → $N$ MOs כל-כך צפופות שהן פס. רוחב הפס תלוי בעוצמת החפיפה בין AOs שכנות — ובמולקולות חזקות-קשר, הפס רחב.

פס הולכה = פסי MOs ריקים (אנטי-קשר) שאלקטרונים יכולים “לזרום” בהם. פס ערכיות = פסי MOs מלאים. פסת אסור (band gap) = הפרש האנרגיה בין השניים.

• מוליך: פס הולכה חלקית מלא, אין פסת אסור
• מבודד: פסת אסור רחבה
• מוליך-למחצה: פסת אסור צרה

זה הכל. הפרק הבא ילמד תורת פסים בפרטים — אבל הרעיון הכימי שמאחוריה הוא בדיוק תורת MO, מוחלת על מוצק.

---

מה למדנו בפרק זה

הפרק הזה היה, במידה רבה, גשר — בין הפיזיקה הבסיסית שבנינו (פרקים 1–10) לבין הכימיה החישובית המעשית שתפגשו בחיי המחקר.

• קירוב בורן-אופנהיימר מאפשר להפריד תנועת גרעינים מאלקטרונים — ובכך מגדיר את משטח הפוטנציאל, שעליו תגובות כימיות “מטפסות ויורדות”.
• הרטרי-פוק הוא הקירוב הפשוט ביותר: כל אלקטרון בשדה ממוצע. עובד, אך מחמיץ קורלציה.
• DFT ומשפטי הוהנברג-קון: כל הפיזיקה גלומה בצפיפות $n(\mathbf{r})$ — שלושה משתנים, לא $3N$. משוואות קון-שם הפכו זאת לשיטה חישובית. הפונקציונל $E_{xc}$ הוא הבלתי-ידוע היחיד; LDA, GGA (PBE), ופונקציונלים היברידיים (B3LYP, HSE) הם קירובים ברמות שונות.
• תורת MO בקצרה: AOs → MOs קשר ואנטי-קשר; סדר הקשר מוסבר, פארמגנטיות $O_2$ מוסברת.
• מהמולקולה למוצק: $N$ AOs → $N$ MOs צפופות = פס אנרגיה. זו ראשית תורת הפסים, שתתפתח בפרק 12.