**קירוב היקל** (Hückel method) או **שיטת אורביטלים מולקולריים של היקל** (HMO), על שמו של אריך היקל (Erich Hückel) היא שיטת צירוף ליניארי של אורביטלים אטומיים לאורביטלים מולקולריים (LCAO MO) הבסיסית ביותר. שיטה זו מאפשרת חישוב אנרגיות של אורביטלים מולקולריים של קשר פאי במערכות פחממנים מצומדות, כגון אתילן (Ethylene), בנזן (Benzene) ובוטאדיאן (1,3-Butadiene). שיטה זו מהווה בסיס תאורטי לכלל היקל.
בהמשך, השיטה הורחבה למערכות מצומדות הטרואטומיות, למשל פירידין (Pyridine), פירול (Pyrrole) ופוראן (Furan)
.
  

אמנם, בכימיה הקוונטית קיימות שיטות חישוב מתקדמות הרבה יותר, אך חשיבות מודל היקל היא דידאקטית ותיאורו מופיע בספרי כימיה קוונטית כאחת מאבני הדרך לשיטות מתקדמות יותר.

מבוא ומאפייני השיטה
-------------------


הנחת היסוד של השיטות סמי-אמפיריות המוקדמות הייתה כי אלקטרוני ה-\pi וה-\sigma  מופרדים לחלוטין. זהו רק מודל מקורב ועל מגבלותיו יורחב בהמשך. ההצדקה לקירוב זה מבוססת על העובדה כי אורביטלי ה-\sigma  וה-\pi הם בעלי סימטריה שונה. הטיפול ב-{\displaystyle n\_{\pi }} האלקטרוני ה-\pi של המערכת מתבצע באמצעות ההמילטוניאן הבא (הכתוב ביחידות אטומיות):




H
^

π

=

∑

i
=
1

n

π



H
^

π

c
o
r
e

(
i
)
+

∑

i
=
1

n

π


∑

j
>
i

1

r

i
j



{\displaystyle {\hat {H}}\_{\pi }=\sum \_{i=1}^{n\_{\pi }}{\hat {H}}\_{\pi }^{\mathrm {core} }(i)+\sum \_{i=1}^{n\_{\pi }}\sum \_{j>i}{\frac {1}{r\_{ij}}}}

{\displaystyle {\hat {H}}_{\pi }=\sum _{i=1}^{n_{\pi }}{\hat {H}}_{\pi }^{\mathrm {core} }(i)+\sum _{i=1}^{n_{\pi }}\sum _{j>i}{\frac {1}{r_{ij}}}}
כאשר




H
^

π

c
o
r
e

(
i
)
=
−

1
2

∇

i

2

+
V
(
i
)

{\displaystyle {\hat {H}}\_{\pi }^{\mathrm {core} }(i)=-{\frac {1}{2}}\nabla \_{i}^{2}+V(i)}

{\displaystyle {\hat {H}}_{\pi }^{\mathrm {core} }(i)=-{\frac {1}{2}}\nabla _{i}^{2}+V(i)}
ו-{\displaystyle V(i)} הוא פוטנציאל אפקטיבי של אלקטרון ה-i הנובע מגרעינים ואלקטרוני ה-\sigma .
  

מכאן, בשיטת היקל מזניחים את הדחיה הבין אלקטרונית ({\displaystyle 1/r\_{ij}}) בהמילטוניאן עצמו, כך שההמילטוניאן של מערכת ה-\pi הוא סכום המילטוניאנים חד אלקטרוניים:




H
^

π

=

∑

i
=
1

n

π

=

H
^


e
f
f

(
i
)

{\displaystyle {\hat {H}}\_{\pi }=\sum \_{i=1}^{n\_{\pi }}={\hat {H}}^{\mathrm {eff} }(i)}

{\displaystyle {\hat {H}}_{\pi }=\sum _{i=1}^{n_{\pi }}={\hat {H}}^{\mathrm {eff} }(i)}
כיוון שההמילטוניאן הוא סכום המילטוניאנים החד-אלקטרוניים, ניתן להפריד בין משתנים. כלומר, אורביטלים מולקולריים מקיימים




H
^


e
f
f

(
i
)

ϕ

i

=

e

i

ϕ

i

{\displaystyle {\hat {H}}^{\mathrm {eff} }(i)\phi \_{i}=e\_{i}\phi \_{i}}

{\displaystyle {\hat {H}}^{\mathrm {eff} }(i)\phi _{i}=e_{i}\phi _{i}}
והאנרגיה של מערכת ה-\pi ניתנת על ידי


E

π

=

∑

i
=
1

n

π


e

i

{\displaystyle E\_{\pi }=\sum \_{i=1}^{n\_{\pi }}e\_{i}}

{\displaystyle E_{\pi }=\sum _{i=1}^{n_{\pi }}e_{i}}
יש לשים לב, כי שיטת היקל לא מציינת במפורש את צורתו של {\displaystyle {\hat {H}}^{\mathrm {eff} }(i)}, במקום זאת היא מספקת כללים לבנייתו, אשר מתוארים בהמשך.
  

קירוב נוסף הוא פריסה של {\displaystyle \phi \_{i}} באמצעות בסיס אטומי מינימלי (מכאן סיווג השיטה LCAO MO) בצורה הבאה:


ϕ

i

=

∑

r
=
1

n

C



c

r
i

f

r

{\displaystyle \phi \_{i}=\sum \_{r=1}^{n\_{\mathbb {C} }}c\_{ri}f\_{r}}

{\displaystyle \phi _{i}=\sum _{r=1}^{n_{\mathbb {C} }}c_{ri}f_{r}}
כאשר f\_{r} הוא אורביטל אטומי מסוג {\displaystyle 2p\_{z}} ו-{\displaystyle n\_{C}} הוא מספר הפחמנים במולקולה. המקדמים מתקבלים בתהליך מינימיזציה על אינטגרל הוריאציה.
  

נסכם את מאפייני השיטה:

* בדרך כלל השיטה מוגבלת למערכות מצומדות.
* כפי שתואר לעיל, השיטה מטפלת רק באלקטרוני ה-\pi.
* השיטה לא לוקחת בחשבון את הגאומטריה של המולקולה אלא רק את יחסי שכן קרוב/רחוק כפי שיתואר בהמשך.
* השיטה מוגבלת למערכות מישוריות.

פיתוח שיטת היקל
---------------

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



E

0

≤

⟨
Ψ

|

H
^

|

Ψ
⟩

⟨
Ψ

|

Ψ
⟩


{\displaystyle E\_{0}\leq {\frac {\langle \Psi |{\hat {H}}|\Psi \rangle }{\langle \Psi |\Psi \rangle }}}

{\displaystyle E_{0}\leq {\frac {\langle \Psi |{\hat {H}}|\Psi \rangle }{\langle \Psi |\Psi \rangle }}}כאשר {\displaystyle \Psi } היא פונקציית הווריאציה. כלומר, ערך תצפית סביב המילטוניאן של פונקציית וריאציה כלשהי יהיה תמיד חסם עליון לאנרגיית מצב יסוד.

בשיטת היקל נציב את:



ϕ

π
i

=

∑

r
=
1

n

C



c

r

f

r

{\displaystyle \phi \_{\pi i}=\sum \_{r=1}^{n\_{\mathbb {C} }}c\_{r}f\_{r}}

{\displaystyle \phi _{\pi i}=\sum _{r=1}^{n_{\mathbb {C} }}c_{r}f_{r}}בתור פונקציית הווריאציה ונחשב את אנרגיית הווריאציה:



ε

i

=

⟨

∑

r
=
1

n

C


c

r
i

f

r

|


H
^


e
f
f

|

∑

s
=
1

n

C


c

s
i

f

s

⟩

⟨


∑

r
=
1

n

C


c

r
i

f

r

|

∑

s
=
1

n

C


c

s
i

f

s

⟩

=

∑

r
=
1

n

c



∑

s
=
1

n

C


c

r

∗

c

s

H

r
s

e
f
f


∑

r
=
1

n

C



∑

s
=
1

n

C


c

r

∗

c

s

S

r
s



≡

A
B

{\displaystyle \varepsilon \_{i}={\frac {\left\langle \displaystyle \sum \_{r=1}^{n\_{C}}c\_{ri}f\_{r}\right|{\hat {H}}^{\mathrm {eff} }\left|\displaystyle \sum \_{s=1}^{n\_{C}}c\_{si}f\_{s}\right\rangle }{\left\langle \left.\displaystyle \sum \_{r=1}^{n\_{C}}c\_{ri}f\_{r}\right|\displaystyle \sum \_{s=1}^{n\_{C}}c\_{si}f\_{s}\right\rangle }}={\frac {\displaystyle \sum \_{r=1}^{n\_{c}}\displaystyle \sum \_{s=1}^{n\_{C}}c\_{r}^{\*}c\_{s}H\_{rs}^{\mathrm {eff} }}{\displaystyle \sum \_{r=1}^{n\_{C}}\displaystyle \sum \_{s=1}^{n\_{C}}c\_{r}^{\*}c\_{s}S\_{rs}}}\equiv {\frac {A}{B}}}

{\displaystyle \varepsilon _{i}={\frac {\left\langle \displaystyle \sum _{r=1}^{n_{C}}c_{ri}f_{r}\right|{\hat {H}}^{\mathrm {eff} }\left|\displaystyle \sum _{s=1}^{n_{C}}c_{si}f_{s}\right\rangle }{\left\langle \left.\displaystyle \sum _{r=1}^{n_{C}}c_{ri}f_{r}\right|\displaystyle \sum _{s=1}^{n_{C}}c_{si}f_{s}\right\rangle }}={\frac {\displaystyle \sum _{r=1}^{n_{c}}\displaystyle \sum _{s=1}^{n_{C}}c_{r}^{*}c_{s}H_{rs}^{\mathrm {eff} }}{\displaystyle \sum _{r=1}^{n_{C}}\displaystyle \sum _{s=1}^{n_{C}}c_{r}^{*}c_{s}S_{rs}}}\equiv {\frac {A}{B}}}{\displaystyle \phi \_{\pi i}} נדרשת להיות מנורמלת. כעת, ניתן לבצע מינימיזציה של {\displaystyle \varepsilon }, למשל, לפי המקדמים {\displaystyle \left\{c\_{r}^{\*}\right\}}, כך שעבור כל {\displaystyle r=1,2,\dots ,n\_{C}} מתקבל:



∂

ε

i

∂

c

r

∗


=

∑

s
=
1

n

C


c

s
i

(

H

r
s

e
f
f

−
ε

S

r
s

)

B

=
0

{\displaystyle {\frac {\partial \varepsilon \_{i}}{\partial c\_{r}^{\*}}}={\frac {\displaystyle \sum \_{s=1}^{n\_{C}}c\_{si}(H\_{rs}^{\mathrm {eff} }-\varepsilon S\_{rs})}{B}}=0}

{\displaystyle {\frac {\partial \varepsilon _{i}}{\partial c_{r}^{*}}}={\frac {\displaystyle \sum _{s=1}^{n_{C}}c_{si}(H_{rs}^{\mathrm {eff} }-\varepsilon S_{rs})}{B}}=0}וזה מוביל למערכת של {\displaystyle n\_{C}} משואות סקולריות מהצורה:



∑

s
=
1

n

C


c

s
i

(

H

r
s

e
f
f

−

ε

i

S

r
s

)

=
0

for

r
=
1
,
2
,
…
,

n

C

.

{\displaystyle \sum \_{s=1}^{n\_{C}}c\_{si}\left(H\_{rs}^{\mathrm {eff} }-\varepsilon \_{i}S\_{rs}\right)=0\quad {\text{for}}\quad r=1,2,\dots ,n\_{C}.}

{\displaystyle \sum _{s=1}^{n_{C}}c_{si}\left(H_{rs}^{\mathrm {eff} }-\varepsilon _{i}S_{rs}\right)=0\quad {\text{for}}\quad r=1,2,\dots ,n_{C}.}במשואות אלה, אנרגיה {\displaystyle \varepsilon } והמקדמים {\displaystyle \left\{c\_{si}\right\}} הם הנעלמים. ביחס לווקטור המקדמים, זוהי מערכת משוואות הומוגנית וקיים פתרון לא טריויאלי אם ורק אם הדטרמיננטה מתאפסת, כלומר כאשר:


det

(

H

r
s

e
f
f

−

ε

i

S

r
s

)

=
0.

{\displaystyle \det \left(H\_{rs}^{\mathrm {eff} }-\varepsilon \_{i}S\_{rs}\right)=0.}

{\displaystyle \det \left(H_{rs}^{\mathrm {eff} }-\varepsilon _{i}S_{rs}\right)=0.}קיימים {\displaystyle n\_{C}} ערכי {\displaystyle \varepsilon }, אשר מקיימים את המשוואה (שורשי המשוואה). הערך הקטן ביותר מקבוצת הפתרונות {\displaystyle \left\{\varepsilon \_{i}\right\}} זהו הקירוב הטוב ביותר לאנרגיית מצב היסוד במסגרת שיטת היקל. שאר הערכים של האנרגיות הם הערכה לאנרגיות מצבים המעוררים.

ניתן לקבוע את {\displaystyle \phi \_{\pi i}} על ידי קביעת וקטור המקדמים המתאימים {\displaystyle \left\{c\_{si}\right\}\_{s=1}^{n\_{C}}}.

רכיבי מטריצה המייצגת המילטוניאן ואינטגרלי החפיפה
------------------------------------------------

כפי שצוין קודם, שיטת היקל לא מציינת במפורש את צורתו של {\displaystyle {\hat {H}}^{\mathrm {eff} }(i)}. פתרון הבעיה המנוסחת לעיל כרוך בבניית מטריצה מייצגת של המילטוניאן בבסיס של אורביטלים אטומיים. כלומר יש לחשב את רכיבי המטריצה {\displaystyle H\_{rs}^{\mathrm {eff} }}. בנוסף, יש לחשב את אינטגרלי החפיפה {\displaystyle S\_{rs}}.

שיטת היקל מספקת כללים פשוטים לקביעת הרכיבים לעיל.

* מודל היקל מניח כי {\displaystyle H\_{rr}^{\mathrm {eff} }} הוא זהה לכל הפחמנים במולקולה. נשים לב, כי זהו ערך תצפית לאנרגיה של אורביטל אטומי. ערך זה ניתן לקירוב על ידי אנרגיית יוניזציה והוא זה שיקבע את סקאלת האנרגיה במערכת. נתון זה נמדד אמפירית, מכאן האופי הסמי-אמפירי של השיטה.
* ערכיהם של האלמנטים מחוץ לאלכסון הראשי נקבעים באופן הבא: 

H

r
s

e
f
f

=

{

β

C

r

a
n
d

C

s

b
o
n
d
e
d


0

o
t
h
e
r
w
i
s
e



{\displaystyle H\_{rs}^{\mathrm {eff} }={\begin{cases}\beta &\mathrm {C} \_{r}\;\mathrm {and} \;\mathrm {C} \_{s}\;\mathrm {bonded} \\0&\mathrm {otherwise} \end{cases}}}

{\displaystyle H_{rs}^{\mathrm {eff} }={\begin{cases}\beta &\mathrm {C} _{r}\;\mathrm {and} \;\mathrm {C} _{s}\;\mathrm {bonded} \\0&\mathrm {otherwise} \end{cases}}}כאן ההנחה היא כי החפיפה בין האורביטלים של שכנים רחוקים (פחמנים אשר לא קשורים זה לזה) קטנה, לכן ניתן להזניח את רכיבי המטריצה {\displaystyle H\_{rs}^{\mathrm {eff} }} כאשר מדובר באורביטלים של פחמנים לא שכנים.
* בדומה לרכיבי מטריצה המייצגת המילטוניאן, רכיבי מטריצת החפיפה הם:

S

r
s

=

δ

r
s

{\displaystyle S\_{rs}=\delta \_{rs}}

{\displaystyle S_{rs}=\delta _{rs}}כאשר {\displaystyle \delta \_{rs}} היא הדלתא של קרונקר. במילים אחרות, מטריצת החפיפה היא מטריצת היחידה. הפירוש הפיזיקאלי של הנחה זו היא כי אורביטלים אטומיים על פחמנים שונים אינם חופפים כלל. ברוב המקרים הנחה זו אינה מתקיימת, כמובן. מסירים הנחה זו בשיטת היקל המורחבת (Extended Hückel method).

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

תוצאות חישוב בשיטת היקל - בוטאדיאן
----------------------------------

ניתן לראות בתמונה את האורביטלים המאוכלסים והריקים של המולקולה.אורביטלים מולקולרים של בוטאדיאן
טיפול בשיטת היקל במולקולת בוטאדיאן: פונקציית הגל של אלקטרוני ה-\pi היא צירוף ליניארי של ארבעה אורביטלים אטומיים \phi  ומקדמי הצירוף הם c\_{i}.

{\displaystyle \ \Psi =c\_{1}\phi \_{1}+c\_{2}\phi \_{2}+c\_{3}\phi \_{3}+c\_{4}\phi \_{4}}

המשוואה הסקולרית:


[

α
−
E

β

0

0

β

α
−
E

β

0

0

β

α
−
E

β

0

0

β

α
−
E

]

×

[

c

1


c

2


c

3


c

4


]

=
0

{\displaystyle {\begin{bmatrix}\alpha -E&\beta &0&0\\\beta &\alpha -E&\beta &0\\0&\beta &\alpha -E&\beta \\0&0&\beta &\alpha -E\\\end{bmatrix}}\times {\begin{bmatrix}c\_{1}\\c\_{2}\\c\_{3}\\c\_{4}\\\end{bmatrix}}=0}

{\displaystyle {\begin{bmatrix}\alpha -E&\beta &0&0\\\beta &\alpha -E&\beta &0\\0&\beta &\alpha -E&\beta \\0&0&\beta &\alpha -E\\\end{bmatrix}}\times {\begin{bmatrix}c_{1}\\c_{2}\\c_{3}\\c_{4}\\\end{bmatrix}}=0}אשר מובילה ל-


(
α
−
E
)
(
α
+
β
−
E
)
−

β

2

=
0

{\displaystyle (\alpha -E)(\alpha +\beta -E)-\beta ^{2}=0\,}

{\displaystyle (\alpha -E)(\alpha +\beta -E)-\beta ^{2}=0\,}כך שאנרגיות האורביטליות הן:


E
±
=
α
+

1
±

5

2

β

{\displaystyle E\pm =\alpha +{\frac {1\pm {\sqrt {5}}}{2}}\beta }

{\displaystyle E\pm =\alpha +{\frac {1\pm {\sqrt {5}}}{2}}\beta }

ראו גם
------

* כימיה קוונטית
* מטריצה כייצוג של העתקה ליניארית

לקריאה נוספת
------------

* Ira N.Levine, Quantum Chemistry, Prentice-Hall, 5, 2000

קישורים חיצוניים
----------------

1. N. Goudard, Y. Carissan, D. Hagebaum-Reignier, S. Humbel (2008). "HuLiS : Java Applet – Simple Hückel Theory and Mesomery – program logiciel software" (בצרפתית). נבדק ב-19 באוגוסט 2010. `{{cite web}}`:  (עזרה)תחזוקה - ציטוט: multiple names: authors list (link)

הערות שוליים
------------

1. ↑  בהקשר של מולקולות אורגניות, האטום שאינו פחמן קרוי הטרואטום (Heteroatom).
2. ↑  Andrew Streitwieser, **Molecular Orbital Theory for Organic Chemist**, Wiley, New York, 1961
3. ↑  Ira N.Levine, 16, **Quantum Chemistry**, 5, Prentice-Hall, 2000