Hartree-Fock-művelet: a módszer matematikai háttere és operátorai

A kvantummechanika alapvető célja az atomok és molekulák elektronikus szerkezetének leírása, ami elengedhetetlen a kémiai reakciók, anyagtulajdonságok és spektroszkópiai jelenségek megértéséhez. Azonban a valóságban létező, több elektronból és több atommagból álló rendszerek pontos, analitikus megoldása a Schrödinger-egyenletnek gyakorlatilag lehetetlen. Emiatt a kvantumkémia számos közelítő módszert fejlesztett ki, melyek közül a Hartree-Fock-módszer az egyik legfontosabb és leggyakrabban használt kiindulópont.

Ez a módszer adja a modern kvantumkémiai számítások alapját, és bár vannak korlátai, elengedhetetlen a bonyolultabb, korrelált módszerek megértéséhez. A Hartree-Fock-módszer lényege egy sokelektronos hullámfüggvény közelítése egyetlen, antiszimmetrikus determinánssal, amely egyelektronos spin-orbitálokból épül fel. Ez a megközelítés lehetővé teszi a Schrödinger-egyenlet közelítő megoldását egy iteratív, önkonzisztens mező (SCF) eljárás keretében.

A módszer ereje abban rejlik, hogy egy komplex soktest-problémát egy sor egyelektronos problémára redukál, melyek megoldhatók. Ez a redukció azonban magával vonja az elektronok közötti korreláció egy részének elhanyagolását, ami a Hartree-Fock-módszer legfőbb korlátja. Ennek ellenére a matematikai elegancia és a viszonylag alacsony számítási költség miatt továbbra is alapvető eszköz marad a kvantumkémikusok kezében.

A Born-Oppenheimer közelítés alapkővének szerepe

Mielőtt belemerülnénk a Hartree-Fock-módszer részleteibe, elengedhetetlen megérteni azt az alapvető közelítést, amely lehetővé teszi a molekuláris rendszerek kvantummechanikai kezelését: a Born-Oppenheimer közelítést. Ez a közelítés azon a tényen alapul, hogy az atommagok sokkal nehezebbek, mint az elektronok, így mozgásuk sebessége nagyságrendekkel kisebb. Ez a jelentős tömegkülönbség lehetővé teszi az atommagok és az elektronok mozgásának szétválasztását.

A Born-Oppenheimer közelítés szerint az atommagok helyzete rögzítettnek tekinthető a molekula elektronikus szerkezetének számítása során. Ez azt jelenti, hogy az elektronok mozgását egy rögzített magtérben írjuk le, és a magok közötti taszítást egyszerűen egy konstans tagként kezeljük az energiafüggvényben. Ez a szétválasztás jelentősen leegyszerűsíti a sokrészecskés Schrödinger-egyenletet, mivel az atommagok kinetikus energiáját elhanyagoljuk, és a hullámfüggvényt csak az elektronkoordináták függvényeként tekintjük.

Matematikailag ez azt jelenti, hogy a teljes molekuláris Hamilton-operátor két részre bontható: egy elektronikus részre és egy magi részre. Az elektronikus Schrödinger-egyenletet oldjuk meg először, amelynek megoldása egy adott magkonfigurációhoz tartozó elektronikus energiát és hullámfüggvényt adja meg. Ez az energia aztán a magok közötti potenciális energia felületet alkotja, amelyen a magok mozognak. A Hartree-Fock-módszer kizárólag az elektronikus Schrödinger-egyenlet megoldására fókuszál.

„A Born-Oppenheimer közelítés nélkül a molekuláris kvantummechanika gyakorlatilag megoldhatatlan lenne. Ez a fundamentális lépés teszi lehetővé, hogy az elektronok viselkedését vizsgáljuk, mintha a magok mozdulatlanok lennének.”

A sokelektronos Schrödinger-egyenlet és a variációs elv

A Born-Oppenheimer közelítés alkalmazása után a probléma leegyszerűsödik a sokelektronos Schrödinger-egyenlet megoldására. Ez az egyenlet írja le az N számú elektron viselkedését egy rögzített atommagi térben:

$\hat{H} \Psi = E \Psi$

Ahol $\hat{H}$ az elektronikus Hamilton-operátor, $\Psi$ a sokelektronos hullámfüggvény, és $E$ az elektronikus energia. A Hamilton-operátor a következő tagokból áll:

• Az elektronok kinetikus energiája: $\sum_{i=1}^{N} -\frac{1}{2} \nabla_i^2$ (atomai egységekben).
• Az elektronok és az atommagok közötti vonzás: $\sum_{i=1}^{N} \sum_{A=1}^{M} -\frac{Z_A}{|\mathbf{r}_i – \mathbf{R}_A|}$.
• Az elektronok közötti taszítás: $\sum_{i

Ez az egyenlet N elektron és M mag esetén N elektronkoordinátától függ. A fő nehézséget az elektron-elektron taszítási tag jelenti, amely a részecskék közötti korrelációt írja le, és megakadályozza az egyenlet szétválasztását egyelektronos részekre.

Mivel a sokelektronos Schrödinger-egyenlet egzakt megoldása lehetetlen, közelítő módszerekre van szükség. Itt lép be a képbe a variációs elv, amely a kvantumkémiai számítások egyik alappillére. A variációs elv kimondja, hogy bármely normalizált közelítő hullámfüggvényre ($\Psi_{közelítő}$) számított energia ($E_{közelítő}$) mindig nagyobb vagy egyenlő lesz, mint a rendszer valódi alapállapoti energiája ($E_{valódi}$):

$E_{közelítő} = \frac{\langle \Psi_{közelítő} | \hat{H} | \Psi_{közelítő} \rangle}{\langle \Psi_{közelítő} | \Psi_{közelítő} \rangle} \ge E_{valódi}$

Ez az elv adja a módszer matematikai alapját: a legjobb közelítő hullámfüggvényt úgy találhatjuk meg, hogy minimalizáljuk az energiát. A Hartree-Fock-módszer célja egy olyan hullámfüggvény megtalálása, amely minimalizálja az energiát a megengedett hullámfüggvények egy speciális osztályán belül.

A hullámfüggvény közelítése: a Slater-determináns

A variációs elv alkalmazásához először egy jól definiált közelítő hullámfüggvényre van szükségünk. A Hartree-Fock-módszer kulcsfontosságú lépése a sokelektronos hullámfüggvény közelítése egyetlen Slater-determinánssal. Ez a választás nem véletlen, hanem a kvantummechanika két alapvető elvének figyelembevételéből fakad: a Pauli-elv és a részecskék megkülönböztethetetlensége.

A Pauli-elv kimondja, hogy két fermion (például elektron) nem foglalhatja el ugyanazt a kvantumállapotot. Ennek következménye, hogy a fermionokból álló sokrészecskés hullámfüggvénynek antiszimmetrikusnak kell lennie a részecskék felcserélésére. A Slater-determináns automatikusan teljesíti ezt a követelményt.

A Slater-determináns egy N elektronból álló rendszer hullámfüggvényét N darab egyelektronos spin-orbitál ($\chi_i$) szorzatának antiszimmetrizált formájaként írja le:

$\Psi = \frac{1}{\sqrt{N!}} \begin{vmatrix}
\chi_1(\mathbf{x}_1) & \chi_2(\mathbf{x}_1) & \cdots & \chi_N(\mathbf{x}_1) \\
\chi_1(\mathbf{x}_2) & \chi_2(\mathbf{x}_2) & \cdots & \chi_N(\mathbf{x}_2) \\
\vdots & \vdots & \ddots & \vdots \\
\chi_1(\mathbf{x}_N) & \chi_2(\mathbf{x}_N) & \cdots & \chi_N(\mathbf{x}_N)
\end{vmatrix}$

Ahol $\mathbf{x}_i$ az i-edik elektron térbeli és spin koordinátáját jelöli, $\chi_i(\mathbf{x}_i)$ pedig az i-edik egyelektronos spin-orbitált. Minden spin-orbitál egy térbeli orbitál ($\phi$) és egy spinfüggvény ($\alpha$ vagy $\beta$) szorzataként fejezhető ki.

Ez a determinánsos forma biztosítja, hogy ha két spin-orbitál azonos lenne, a determináns értéke zérus lenne, ami összhangban van a Pauli-elvvel. A $\frac{1}{\sqrt{N!}}$ tényező a normalizációt biztosítja.

A Hartree-Fock-módszer lényege tehát abban rejlik, hogy a legjobb, energiát minimalizáló spin-orbitálokat keressük, feltételezve, hogy a sokelektronos hullámfüggvény egyetlen Slater-determináns formájában írható le. Ez a közelítés az úgynevezett független részecske modell, ahol minden elektron mozgását a többi elektron átlagos, statikus mezőjében írjuk le, elhanyagolva az elektronok közötti pillanatnyi korrelációt.

A Hartree-Fock energiafüggvény levezetése

A Slater-determináns formájában felírt hullámfüggvénnyel a variációs elv alkalmazása azt jelenti, hogy minimalizálnunk kell az energia kifejezést az egyelektronos spin-orbitálok változtatásával. Az energia kifejezése, amelyet minimalizálunk, a következő:

$E = \langle \Psi | \hat{H} | \Psi \rangle$

Ahol $\Psi$ egy normalizált Slater-determináns. A Hamilton-operátor, ahogy már említettük, egyelektronos és kételektronos operátorokból áll. Az egyelektronos operátorok (kinetikus energia és mag-elektron vonzás) additívak, azaz minden elektronra külön hatnak, míg a kételektronos operátorok (elektron-elektron taszítás) páronként hatnak.

Az energia kifejezésének kifejtése egy Slater-determinánsra meglehetősen bonyolult, de a végeredmény egy elegáns formula, amely az egyelektronos spin-orbitálok integráljait tartalmazza:

$E = \sum_{i=1}^{N} \langle \chi_i | \hat{h} | \chi_i \rangle + \sum_{i

Nézzük meg részletesebben a tagokat:

• $\sum_{i=1}^{N} \langle \chi_i | \hat{h} | \chi_i \rangle$: Ez az egyelektronos energiák összege. $\hat{h}$ az egyelektronos Hamilton-operátor, amely tartalmazza az i-edik elektron kinetikus energiáját és az atommagokkal való vonzását.
• $\sum_{iCoulomb-integrálok összege (más néven $J_{ij}$). Ez a tag az i-edik és j-edik elektron közötti elektrosztatikus taszítást írja le, mintha az elektronok töltéseloszlásai elmosódottak lennének, és egymástól függetlenül hatnának.
• $\sum_{icserélődési integrálok összege (más néven $K_{ij}$). Ez a tag tisztán kvantummechanikai eredetű, és a Pauli-elv következménye. Nincs klasszikus analógiája, és az azonos spinű elektronok közötti taszítás csökkenését fejezi ki, ami a hullámfüggvény antiszimmetriájából adódik.

A minimalizálás elvégzéséhez alkalmazzuk a Lagrange-féle multiplikátorok módszerét, figyelembe véve azt a mellékfeltételt, hogy a spin-orbitáloknak ortonormáltnak kell lenniük ($\langle \chi_i | \chi_j \rangle = \delta_{ij}$). A minimalizálás eredményeként kapjuk meg a Hartree-Fock-egyenleteket, amelyek egy sor differenciálegyenletet jelentenek az egyelektronos spin-orbitálokra.

A Fock-operátor és komponensei

A Hartree-Fock-egyenletek a következő formában írhatók fel:

$\hat{f} \chi_i = \epsilon_i \chi_i$

Ahol $\hat{f}$ a Fock-operátor, $\chi_i$ az i-edik spin-orbitál, és $\epsilon_i$ a hozzá tartozó spin-orbitál energia. Ez az egyenlet formálisan hasonlít egy egyelektronos Schrödinger-egyenletre, de a Fock-operátor magában foglalja a többi elektron átlagos mezőjének hatását is, ami teszi az egyenletet önkonzisztenssé.

A Fock-operátor a következő komponensekből épül fel:

$\hat{f} = \hat{h} + \sum_{j=1}^{N} (\hat{J}_j – \hat{K}_j)$

Kinetikus energia operátor

Az első tag az egyelektronos Hamilton-operátor, $\hat{h}$, ami a kinetikus energiát és a mag-elektron vonzást tartalmazza. A kinetikus energia operátor az i-edik elektron mozgási energiáját írja le:

$\hat{T}_i = -\frac{1}{2} \nabla_i^2$

Ez a tag az elektron hullámfüggvényének görbületével arányos, és a kvantummechanika alapvető eleme.

Mag-elektron vonzási operátor

A mag-elektron vonzási operátor az i-edik elektron és az atommagok közötti vonzást írja le. Mivel a magok rögzítettek a Born-Oppenheimer közelítésben, ez a tag potenciális energiát jelent:

$\hat{V}_{ne,i} = \sum_{A=1}^{M} -\frac{Z_A}{|\mathbf{r}_i – \mathbf{R}_A|}$

Ahol $Z_A$ az A-adik mag töltése, $\mathbf{R}_A$ az A-adik mag pozíciója, és $\mathbf{r}_i$ az i-edik elektron pozíciója. Ez a tag biztosítja, hogy az elektronok az atommagok vonzása alatt maradjanak.

Tehát az egyelektronos Hamilton-operátor:

$\hat{h}_i = \hat{T}_i + \hat{V}_{ne,i}$

Coulomb-operátor (J)

A Coulomb-operátor, $\hat{J}_j$, leírja az i-edik elektron és a j-edik elektron töltéseloszlása közötti átlagos elektrosztatikus taszítást. Ez egy „átlagos mező” hatás, ahol a j-edik elektront egy elmosódott töltéseloszlásként kezeljük:

$\hat{J}_j \chi_i(\mathbf{x}_i) = \left( \int \frac{|\chi_j(\mathbf{x}_j)|^2}{|\mathbf{r}_i – \mathbf{r}_j|} d\mathbf{x}_j \right) \chi_i(\mathbf{x}_i)$

A $|\chi_j(\mathbf{x}_j)|^2$ kifejezés a j-edik elektron valószínűségi sűrűségét reprezentálja. A Coulomb-operátor egy lokális operátor, ami azt jelenti, hogy az i-edik elektronra ható potenciál csak az i-edik elektron pozíciójától függ.

Cserélődési operátor (K)

A cserélődési operátor, $\hat{K}_j$, sokkal komplexebb és tisztán kvantummechanikai eredetű. Nincs klasszikus analógiája. Ez az operátor az azonos spinű elektronok közötti antiszimmetria következménye, és csökkenti az azonos spinű elektronok közötti taszítást. Ez a Pauli-elv megnyilvánulása:

$\hat{K}_j \chi_i(\mathbf{x}_i) = \left( \int \frac{\chi_j^*(\mathbf{x}_j) \chi_i(\mathbf{x}_j)}{|\mathbf{r}_i – \mathbf{r}_j|} d\mathbf{x}_j \right) \chi_j(\mathbf{x}_i)$

Látható, hogy a cserélődési operátor nem lokális, azaz az i-edik elektronra ható potenciál nem csak az i-edik elektron pozíciójától, hanem a j-edik elektron pozíciójától is függ, és „felcseréli” az i-edik és j-edik spin-orbitált. Ez a nem-lokalitás teszi a cserélődési operátor számítását sokkal bonyolultabbá, mint a Coulomb-operátorét.

A Fock-operátor tehát az egyelektronos Hamilton-operátorból, valamint az elektronok közötti átlagos taszítást leíró Coulomb- és cserélődési operátorokból tevődik össze. Mivel a Coulomb- és cserélődési operátorok függenek a keresett spin-orbitáloktól ($\chi_j$-ktől), a Hartree-Fock-egyenletek egy integrál-differenciálegyenlet-rendszert alkotnak, amelyet iteratív módon kell megoldani.

A Coulomb- és cserélődési integrálok értelmezése

A Hartree-Fock-módszer alapvető építőkövei a Coulomb- és cserélődési integrálok. Ezek az integrálok adják a matematikai kifejezést az elektronok közötti kölcsönhatásokra a Slater-determináns közelítés keretében. Mélyebb megértésük kulcsfontosságú a módszer fizikai tartalmának felfogásához.

Coulomb-integrál (J_ij)

A Coulomb-integrál $J_{ij}$ az i-edik és j-edik elektron közötti elektrosztatikus taszítás energiáját reprezentálja. Matematikailag a következőképpen definiálható:

$J_{ij} = \langle \chi_i \chi_j | \chi_i \chi_j \rangle = \iint \chi_i^*(\mathbf{x}_1) \chi_j^*(\mathbf{x}_2) \frac{1}{|\mathbf{r}_1 – \mathbf{r}_2|} \chi_i(\mathbf{x}_1) \chi_j(\mathbf{x}_2) d\mathbf{x}_1 d\mathbf{x}_2$

Ez a kifejezés úgy értelmezhető, mint két töltéssűrűség ($\rho_i(\mathbf{r}_1) = |\chi_i(\mathbf{x}_1)|^2$ és $\rho_j(\mathbf{r}_2) = |\chi_j(\mathbf{x}_2)|^2$) közötti klasszikus elektrosztatikus taszítás. Az integrál azt méri, hogy mennyi energiába kerülne, ha a két elektron töltéseloszlása kölcsönhatna egymással. Mindig pozitív értékű, mivel a taszítási energia pozitív.

A Coulomb-integrál független az elektronok spinjétől, csak a térbeli eloszlásuktól függ. Ez azt jelenti, hogy két elektron taszítja egymást, függetlenül attól, hogy spinjük azonos vagy ellentétes.

Cserélődési integrál (K_ij)

A cserélődési integrál $K_{ij}$ a hullámfüggvény antiszimmetriájából fakad, és csak az azonos spinű elektronok között létezik. Matematikailag a következőképpen írható fel:

$K_{ij} = \langle \chi_i \chi_j | \chi_j \chi_i \rangle = \iint \chi_i^*(\mathbf{x}_1) \chi_j^*(\mathbf{x}_2) \frac{1}{|\mathbf{r}_1 – \mathbf{r}_2|} \chi_j(\mathbf{x}_1) \chi_i(\mathbf{x}_2) d\mathbf{x}_1 d\mathbf{x}_2$

A kulcsfontosságú különbség a Coulomb-integrálhoz képest az, hogy a második részecske (itt $\mathbf{x}_2$) hullámfüggvényei felcserélődnek az integrálon belül. Ez a felcserélődés a Pauli-elv következménye, amely szerint két azonos spinű elektron nem lehet ugyanazon a helyen, azaz kerülik egymást. Ez a „kerülés” azt eredményezi, hogy az azonos spinű elektronok közötti taszítás effektíve csökken, ami az energia csökkenésében nyilvánul meg. A cserélődési integrál mindig pozitív, de az energia kifejezésben negatív előjellel szerepel, így az azonos spinű elektronok közötti taszítást csökkenti.

A cserélődési integrál tehát egy kvantummechanikai korrelációt ír le, amelyet Fermi-korrelációnak is neveznek. Ez a korreláció abból adódik, hogy az azonos spinű elektronok hullámfüggvényei kénytelenek nullává válni, amikor a részecskék pozíciói megközelítik egymást. Ez egy „lyukat” hoz létre az elektron sűrűségében a referencia elektron körül, az úgynevezett Fermi-lyukat. Ezt a jelenséget a Hartree-Fock-módszer automatikusan figyelembe veszi.

Fontos megkülönböztetni a Fermi-korrelációt az úgynevezett dinamikus korrelációtól, amelyet a Hartree-Fock-módszer elhanyagol. A dinamikus korreláció az elektronok közötti pillanatnyi, klasszikus értelemben vett taszításra utal, ahol az elektronok pillanatnyi pozíciói befolyásolják egymás mozgását. Ezt a Hartree-Fock-módszer nem képes leírni, mivel az elektronokat csak a többi elektron átlagos mezőjében mozognak.

A Hartree-Fock egyenletek felépítése

A Fock-operátor bevezetésével a Hartree-Fock-egyenletek egy sor egyelektronos egyenletté egyszerűsödnek:

$\hat{f} \chi_i = \epsilon_i \chi_i$

Ez egy eigenérték egyenlet, ahol $\chi_i$ az i-edik spin-orbitál (eigenfüggvény), és $\epsilon_i$ a hozzá tartozó spin-orbitál energia (eigenérték). Ezek az egyenletek egy N x N méretű rendszert alkotnak, ahol N a spin-orbitálok száma (ami egyenlő az elektronok számával egy zárt héjú rendszerben).

Mivel a Fock-operátor a keresett spin-orbitáloktól függ, az egyenletek nem lineárisak, és nem oldhatók meg közvetlenül. Ezért van szükség egy iteratív eljárásra, az úgynevezett önkonzisztens mező (SCF) módszerre.

Az $\epsilon_i$ spin-orbitál energiák fizikai jelentőséggel bírnak. A Koopmans-tétel szerint egy zárt héjú rendszerben a legmagasabb foglalt spin-orbitál energiája közelítőleg megegyezik a legkönnyebben ionizálható elektron ionizációs energiájával. Hasonlóképpen, a legalacsonyabb üres spin-orbitál energiája közelítőleg megegyezik az elektronaffinitással.

A Hartree-Fock energia (E) és a spin-orbitál energiák ($\epsilon_i$) közötti kapcsolat a következő:

$E = \sum_{i=1}^{N} \epsilon_i – \sum_{i

Ez a formula azt mutatja, hogy az elektronok közötti kölcsönhatásokat kétszeresen számoljuk az $\sum \epsilon_i$ összegben (mivel minden $\epsilon_i$ tartalmazza a többi elektron átlagos mezőjét), ezért a Coulomb- és cserélődési integrálokat egyszer le kell vonni. Ez a képlet rávilágít arra is, hogy az egyelektronos energiák önmagukban nem adnak pontos képet a teljes molekuláris energiáról, mivel az elektron-elektron kölcsönhatásokat tartalmazzák.

Az önkonzisztens mező (SCF) eljárás

Mivel a Hartree-Fock-egyenletek nem lineárisak, és a Fock-operátor függ a keresett spin-orbitáloktól, azokat iteratív módon kell megoldani. Ezt az eljárást önkonzisztens mező (SCF) módszernek nevezik. Az SCF eljárás célja, hogy olyan spin-orbitálokat találjon, amelyek önkonzisztensek, azaz a bemeneti orbitálokból számított Fock-operátor ugyanazokat az orbitálokat adja ki megoldásként.

Az SCF eljárás lépései a következők:

1. Kezdő orbitálok becslése: Először szükségünk van egy kezdeti becslésre a spin-orbitálokra. Ezt gyakran a mag-elektron Hamiltonianus (csak kinetikus energia és mag-elektron vonzás) megoldásával, vagy egyszerűbb, empirikus módszerekkel tesszük. Ezekből a kezdeti orbitálokból egy kezdeti sűrűségmátrix is képezhető.
2. Fock-operátor felépítése: Az aktuális spin-orbitálok (vagy sűrűségmátrix) felhasználásával kiszámítjuk a Coulomb- és cserélődési operátorokat, majd felépítjük a teljes Fock-operátort.
3. Hartree-Fock-egyenletek megoldása: Megoldjuk a $\hat{f} \chi_i = \epsilon_i \chi_i$ egyenleteket. Ez egy eigenérték probléma, amely új spin-orbitálokat és spin-orbitál energiákat eredményez.
4. Konvergencia ellenőrzése: Összehasonlítjuk az új spin-orbitálokat (vagy az általuk képzett sűrűségmátrixot, vagy a teljes energiát) az előző iterációból származókkal. Ha a különbség egy előre meghatározott küszöbérték alá esik, az eljárás konvergált, és megkaptuk az önkonzisztens megoldást.
5. Iteráció: Ha az eljárás nem konvergált, az új spin-orbitálokat (vagy sűrűségmátrixot) használjuk a következő iteráció kiindulópontjaként, és visszatérünk a 2. lépésre.

Az SCF eljárás sikeressége nagymértékben függ a kezdeti becsléstől és a konvergencia gyorsítási technikáktól. A konvergencia nem mindig garantált, különösen bonyolult rendszerek vagy nyílt héjú esetek esetén. Az eljárás végén kapott spin-orbitálok adják a rendszer elektronikus szerkezetének Hartree-Fock közelítését.

„Az önkonzisztens mező eljárás a Hartree-Fock-módszer szíve, egy iteratív tánc az elektronok között, amíg mindenki meg nem találja a helyét a közös, átlagos mezőben.”

Bázisfüggvények szerepe és a Roothaan-Hall egyenletek

A Hartree-Fock-egyenletek, $\hat{f} \chi_i = \epsilon_i \chi_i$, differenciálegyenletek formájában vannak, és analitikusan csak a legegyszerűbb rendszerekre (pl. hidrogénatom) oldhatók meg. Molekuláris rendszerek esetében numerikus megoldásra van szükség, amihez a spin-orbitálokat egy véges számú, előre definiált bázisfüggvény lineáris kombinációjaként fejezzük ki.

Ez az úgynevezett lineáris kombináció atomi orbitálokból (LCAO) közelítés. A spin-orbitálokat a következőképpen írjuk fel:

$\chi_i = \sum_{\mu=1}^{K} C_{\mu i} \phi_\mu$

Ahol $\phi_\mu$ a bázisfüggvények, és $C_{\mu i}$ a lineáris kombináció együtthatói. A bázisfüggvények jellemzően Gauss-típusú függvények (Gaussian Type Orbitals, GTOs), amelyek matematikai formája kényelmes az integrálok számításához. A bázisfüggvények száma (K) határozza meg a számítás pontosságát és költségét.

A bázisfüggvények bevezetése a differenciálegyenleteket egy mátrix-eigenérték problémává alakítja át, amelyet Roothaan-Hall egyenleteknek nevezünk (zárt héjú rendszerekre) vagy Pople-Nesbet egyenleteknek (nyílt héjú rendszerekre). A Roothaan-Hall egyenletek a következők:

$\mathbf{F C} = \mathbf{S C} \mathbf{\epsilon}$

Ahol:

• $\mathbf{F}$ a Fock-mátrix, amelynek elemei $F_{\mu\nu} = \langle \phi_\mu | \hat{f} | \phi_\nu \rangle$. Ez a mátrix tartalmazza az egyelektronos és kételektronos integrálokat a bázisfüggvények között.
• $\mathbf{C}$ a koefficiens mátrix, amelynek oszlopai a spin-orbitálok lineáris kombinációjának együtthatóit tartalmazzák ($C_{\mu i}$).
• $\mathbf{S}$ az átfedési mátrix, amelynek elemei $S_{\mu\nu} = \langle \phi_\mu | \phi_\nu \rangle$. Ez a mátrix méri a bázisfüggvények közötti átfedést. Ha a bázisfüggvények ortonormáltak, $\mathbf{S}$ egységmátrix.
• $\mathbf{\epsilon}$ egy diagonális mátrix, amelynek elemei a spin-orbitál energiák ($\epsilon_i$).

A Roothaan-Hall egyenletek megoldása az SCF eljárásba ágyazva történik. Minden iterációban felépítjük a Fock-mátrixot az aktuális koefficiens mátrixból (ami a spin-orbitálokat adja), majd megoldjuk az eigenérték problémát az új koefficiens mátrix és spin-orbitál energiák meghatározásához. Ez az eljárás folytatódik, amíg a koefficiens mátrix vagy az energia konvergál.

A bázisfüggvények kiválasztása kritikus fontosságú. A nagyobb báziskészletek (több bázisfüggvénnyel) pontosabb eredményeket adnak, de jelentősen növelik a számítási költségeket. Különböző típusú báziskészletek léteznek, mint például a minimális bázisok (pl. STO-3G), polarizációs függvényeket tartalmazó bázisok (pl. 6-31G(d,p)), és diffúz függvényeket tartalmazó bázisok (pl. 6-31+G(d,p)), amelyek mindegyike más-más kompromisszumot kínál a pontosság és a számítási költség között.

Korlátozások és az elektronkorreláció hiánya

Bár a Hartree-Fock-módszer forradalmi volt és alapvető betekintést nyújt a molekuláris szerkezetbe, fontos tisztában lenni a korlátaival. A legfőbb korlátozás az elektronkorreláció nem teljes kezelése.

Ahogy korábban említettük, a Hartree-Fock-módszer a sokelektronos hullámfüggvényt egyetlen Slater-determinánssal közelíti. Ez azt jelenti, hogy minden elektron mozgását a többi elektron átlagos, statikus mezőjében írja le. Ez a „független részecske” modell elhanyagolja az elektronok közötti pillanatnyi, dinamikus korrelációt.

Az elektronok valóságban nem mozognak egymástól függetlenül egy átlagos mezőben. Ehelyett pillanatról pillanatra taszítják egymást, és mozgásuk szinkronban van, hogy elkerüljék egymást. Ezt a pillanatnyi kölcsönhatást, amely nem írható le egyetlen Slater-determinánssal, nevezzük elektronkorrelációnak.

A Hartree-Fock-módszer csak a Fermi-korrelációt (az azonos spinű elektronok Pauli-elv által diktált elkerülését) veszi figyelembe a cserélődési operátoron keresztül. Azonban a dinamikus korrelációt (az ellentétes spinű, vagy akár azonos spinű elektronok közötti pillanatnyi elkerülést) figyelmen kívül hagyja. Ez az elhanyagolás azt eredményezi, hogy a Hartree-Fock energia mindig magasabb, mint a valódi energia, és az energiakülönbséget korrelációs energiának nevezzük.

A korrelációs energia hiánya számos problémát okoz:

• Kötési energiák alulbecslése: A Hartree-Fock-módszer általában alulbecsüli a kötési energiákat, mert a molekuláris állapotban nagyobb a korrelációs energia, mint az atomi állapotban.
• Geometriai paraméterek: A kötéshosszak és kötésszögek néha pontatlanok lehetnek.
• Reakcióenergiák: A reakcióenergiák számítása problémás lehet, ha a reaktánsok és termékek korrelációs energiája jelentősen eltér.
• Disszociáció helytelen leírása: A Hartree-Fock-módszer nem képes helyesen leírni a kovalens kötések disszociációját, mivel a disszociált termékek gyakran nyílt héjú jelleget mutatnak, amit egyetlen determináns nem tud pontosan leírni.
• Spektroszkópiai tulajdonságok: Az elektronikus átmenetek energiái, gerjesztett állapotok és polarizálhatóságok pontatlanok lehetnek.

Ezen korlátozások ellenére a Hartree-Fock-módszer továbbra is rendkívül fontos. Egyrészt egy jól definiált, reprodukálható alapállapotot szolgáltat, amelyre a korrelált módszerek épülhetnek. Másrészt számos rendszerre (különösen zárt héjú molekulák alapállapotára a molekula egyensúlyi geometriájának közelében) elfogadható pontosságú eredményeket ad.

Restricted és unrestricted Hartree-Fock (RHF és UHF)

A Hartree-Fock-módszer alkalmazása során fontos megkülönböztetni a különböző megközelítéseket attól függően, hogy a rendszer zárt vagy nyílt héjú-e, és hogyan kezeljük a spin-orbitálokat. A két leggyakoribb megközelítés a Restricted Hartree-Fock (RHF) és az Unrestricted Hartree-Fock (UHF).

Restricted Hartree-Fock (RHF)

Az RHF módszert zárt héjú rendszerekre alkalmazzák, ahol minden térbeli orbitált két elektron foglal el, ellentétes spinnel ($\alpha$ és $\beta$). Ez azt jelenti, hogy a térbeli orbitálok azonosak az $\alpha$ és $\beta$ spinű elektronok számára. Matematikailag:

$\chi_{i\alpha}(\mathbf{x}) = \phi_i(\mathbf{r}) \alpha(\omega)$

$\chi_{i\beta}(\mathbf{x}) = \phi_i(\mathbf{r}) \beta(\omega)$

Ahol $\phi_i$ a térbeli orbitál, $\alpha$ és $\beta$ a spinfüggvények. Az RHF megközelítésben a Fock-operátor ugyanaz az $\alpha$ és $\beta$ spinű elektronok számára. Ez a megközelítés garantálja, hogy a hullámfüggvény szimmetrikus legyen a spinoperátorok (S² és S_z) eigenfüggvényeihez képest, és tiszta spinállapotot eredményez (azaz a spin szennyeződés nulla). Az RHF a legtöbb zárt héjú molekula esetében jó kiindulópontot ad, és viszonylag egyszerűen implementálható.

Unrestricted Hartree-Fock (UHF)

Az UHF módszert nyílt héjú rendszerekre (pl. szabad gyökök, gerjesztett állapotok) alkalmazzák, ahol az $\alpha$ és $\beta$ spinű elektronok számára különböző térbeli orbitálokat engedélyezünk. Matematikailag:

$\chi_{i\alpha}(\mathbf{x}) = \phi_i^\alpha(\mathbf{r}) \alpha(\omega)$

$\chi_{i\beta}(\mathbf{x}) = \phi_i^\beta(\mathbf{r}) \beta(\omega)$

Ahol $\phi_i^\alpha$ és $\phi_i^\beta$ különböző térbeli orbitálok lehetnek. Ez a rugalmasság lehetővé teszi a spinpolarizáció jobb leírását, és gyakran alacsonyabb energiát eredményez, mint az RHF nyílt héjú rendszerek esetén. Az UHF-ben két külön Fock-egyenletet oldunk meg, egyet az $\alpha$ spinű elektronoknak és egyet a $\beta$ spinű elektronoknak, mindkettőt a másik spinű elektronok átlagos mezőjében.

Azonban az UHF-nek van egy jelentős hátránya: az így kapott hullámfüggvény nem feltétlenül tiszta spinállapot. Gyakran „spin szennyezett”, ami azt jelenti, hogy nem csak a kívánt spinállapotot (pl. dublett gyökök esetén) tartalmazza, hanem magasabb spinállapotokat (pl. kvartett, szextett) is. Ez a spin szennyeződés torzíthatja az energiákat és más fizikai tulajdonságokat. Léteznek korrekciós módszerek (pl. spin-projekció) a spin szennyeződés eltávolítására, de ezek további számítási költséggel járnak.

Létezik egy harmadik megközelítés is, a Restricted Open-shell Hartree-Fock (ROHF), amelyet szintén nyílt héjú rendszerekre használnak. Az ROHF megközelítésben a zárt héjú orbitálokat RHF-módon kezeljük (két elektron/orbitál), míg a nyílt héjú (egyszeresen foglalt) orbitálokhoz egyetlen elektront rendelünk. Az ROHF hullámfüggvények tiszta spinállapotot eredményeznek, de az UHF-nél magasabb energiát adhatnak, mivel kevésbé rugalmasak. Az ROHF-et gyakran használják kiindulópontként a korrelált módszerekhez nyílt héjú rendszerek esetén, mivel tiszta spinállapotot biztosít.

Jellemző	Restricted Hartree-Fock (RHF)	Unrestricted Hartree-Fock (UHF)
Alkalmazás	Zárt héjú rendszerek	Nyílt héjú rendszerek
Térbeli orbitálok	Azonos $\alpha$ és $\beta$ spineknek	Különböző $\alpha$ és $\beta$ spineknek
Spin állapot	Tiszta spin állapot	Spin szennyezett lehet
Rugalmasság	Kisebb	Nagyobb
Energia	Magasabb lehet nyílt héjú rendszereknél	Alacsonyabb lehet nyílt héjú rendszereknél

A Hartree-Fock-módszer jelentősége és továbbfejlesztések alapja

A Hartree-Fock-módszer, annak ellenére, hogy számos korláttal rendelkezik, a modern kvantumkémiai számítások sarokköve. Jelentősége túlmutat azon, hogy önmagában képes közelítő megoldásokat szolgáltatni az elektronikus szerkezetre. Sokkal inkább egy kiindulópontként szolgál a pontosabb, korrelált módszerek számára.

A Hartree-Fock-módszer által kapott orbitálok és energiák adják az alapját szinte minden utó-Hartree-Fock (post-Hartree-Fock) módszernek. Ezek a módszerek célja az elektronkorreláció figyelembe vétele, amelyet a Hartree-Fock elhanyagol. Ilyen módszerek a következők:

• Konfiguráció kölcsönhatás (CI): A CI módszerek a Hartree-Fock alapállapoti hullámfüggvényt kiegészítik gerjesztett determinánsokkal (ahol elektronokat emelünk magasabb energiájú, üres orbitálokra), hogy jobban leírják az elektronkorrelációt. A CI-ben a hullámfüggvényt több Slater-determináns lineáris kombinációjaként írjuk fel.
• Möller-Plesset perturbációs elmélet (MPn): Az MPn módszerek a perturbációs elméletet alkalmazzák a Hartree-Fock-módszerre, hogy szisztematikusan beépítsék az elektronkorrelációt. Az MP2 a leggyakrabban használt, és viszonylag olcsón ad korrelációs energiát.
• Kombinált klaszter (CC): A CC módszerek a legpontosabb és legkifinomultabb korrelált módszerek közé tartoznak. Exponenciális operátorok segítségével írják le az elektronkorrelációt, és rendkívül pontos eredményeket szolgáltatnak.

Ezek a fejlettebb módszerek mind a Hartree-Fock orbitálokat használják kiindulópontként, és azok energiáira építenek. A Hartree-Fock-módszer adja a „nulladik rendű” közelítést, amelyre a korrelációs hatások ráépülnek. Ezért a Hartree-Fock-módszer megértése alapvető fontosságú a kvantumkémia mélyebb megértéséhez és a modern számítási módszerek alkalmazásához.

A Hartree-Fock-módszer nemcsak elméleti jelentőséggel bír, hanem gyakorlati alkalmazásai is széles körűek. Segít az alapszerkezetek (pl. molekuláris geometriák, dipólusmomentumok) előrejelzésében, és gyakran elegendő pontosságot biztosít nagy molekulák esetében, ahol a korrelált módszerek számítási költségei megfizethetetlenek lennének. Továbbá, a sűrűségfunkcionál elmélet (DFT) népszerű módszerei is gyakran Hartree-Fock-szerű Kohn-Sham egyenleteket oldanak meg, és a hibrid funkcionálok (pl. B3LYP) részben tartalmazzák a Hartree-Fock cserélődési energiát.

Összességében a Hartree-Fock-módszer egy elegáns és hatékony közelítés a sokelektronos problémára, amely a kvantumkémia alapjait képezi. Matematikai háttere és az operátorok részletes megértése elengedhetetlen ahhoz, hogy a modern kvantumkémiai számításokat ne csak alkalmazzuk, hanem azok mélyebb fizikai tartalmát is felfogjuk.