Hartree–Fock-teoria

Hartree–Fock-teoria on approksimatiivinen menetelmä kvanttimekaanisen monielektronisysteemin Schrödingerin yhtälön ratkaisemiseksi. Lähtökohtana teoriassa on kirjoittaa monielektronisysteemin aaltofunktio Slaterin determinantin avulla. Yritefunktio johtaa niin kutsuttuihin Hartree–Fock-yhtälöihin, jotka kuvaavat elektronien liikettä toisten elektronien muodostamassa keskimääräisessä elektrostaattisessa kentässä. Hartree–Fock-approksimaatio on monissa käytännön sovelluksissa huono ja antaa lähinnä kvalitatiivisia tuloksia. Menetelmää on kuitenkin mahdollista parantaa lukuisin eri keinoin, jolloin eksaktia ratkaisua Schrödingerin yhtälöön on mahdollista lähestyä systemaattisesti.
Kemiassa yleisesti käytetty orbitaalimalli, jossa elektronit sijaitsevat eri muotoisilla ja lokaaleilla orbitaaleilla, perustuu Hartree–Fock-approksimaatioon. Yleisesti ottaen Hartree–Fock-teoria ennustaa hyvin molekyylien rakenteet, mutta epäonnistuu sidosenergioiden ja värähtelytaajuuksien määrittämisessä. Teorian vajavaisuudet johtuvat elektronikorrelaation puuttumisesta, joka on mahdollista huomioida erilaisilla Hartree–Fock-johdannaisilla.
Slaterin determinantti.
Slaterin determinanttia käytetään fermionisysteemin antisymmetrisen aaltofunktion kirjoittamiseen. Kvanttimekaanisessa systeemissä ei voida tehdä erottelua kahden hiukkasen välillä vaan monihiukkassysteemiä kuvaa yksi todennäköisyysjakauma. Tämän vuoksi monihiukkassysteemin aaltofunktion tulee olla verrannollinen aaltofunktioon, jossa kahden hiukkasen koordinaatteja on vaihdettu eli formula_1. Mikäli muutos tehdään uudelleen, saadaan formula_2 eli formula_3. Aaltofunktio voi olla siis joko "symmetrinen" (+) tai "antisymmetrinen" (-) koordinaatteja vaihdettaessa. Bosonit määritellään symmetrisiksi ja fermionit antisymmetrisiksi hiukkasiksi. Antisymmetrisyyden vuoksi voidaan monifermionisysteemin aaltofunktio kirjoittaa determinanttimuodossa, sillä determinantti muuttaa merkkiään riviä tai saraketta vaihdettaessa aivan kuten antisymmetrisen aaltofunktion kuuluukin. Lisäksi, jos kaksi riviä ovat samat, on determinantti nolla eli fermionien kieltosääntö pätee.
formula_4
Esimerkiksi kahden elektronin systeemille Slaterin determinantti on
formula_5
Yllä formula_6 on spinorbitaali, joka sisältää elektronin paikkakoordinaatin lisäksi myös elektronin spinin:
formula_7
Periaatteessa Slaterin determinantti on separoituva yrite, koska yksielektroniorbitaalit kirjoitetaan tulomuodossa. Antisymmetrisyyden vuoksi elektronien spinit ovat kuitenkin korreloituneet, joten elektronien välillä on olemassa jo alkujaan heikko korrelaatio. Koska systeemiä kuvaavaa aaltofunktiota ei tiedetä, joudutaan yritettä approksimoimaan eri keinoin. Yritefunktiossa esiintyviä spinorbitaaleja approksimoidaan kantafunktioiden avulla, jotka riippuvat joistain vapaista parametreista. Näitä parametreja muuttumalla, voidaan systeemin perustilan aaltofunktio löytää. Kvanttimekaniikan variaatioperiaate asettaa perustilan aaltofunktion systeemin kokonaisenergian ehdottomaksi alarajaksi: yritefunktio, joka ei ole oikea perustilan aaltofunktio, antaa aina perustilan energiaksi oikeaa energiaa suuremman arvon. Paras yritefunktio (käytettyjen kantafunktioiden rajoissa) löydetään siis etsimällä kokonaisenergiafunktionaalin minimi.
Hartree–Fock-approksimaatio.
Oikea monielektronisysteemin aaltofunktio voidaan kirjoittaa spinorbitaalien tulojen lineaarikombinaationa. Tämä lineaarikombinaatio voidaan kirjoittaa myös Slaterin determinanttien avulla.
formula_8
Hartree–Fock-teoriassa aaltofunktiota — joka koostuu Slaterin determinanttien lineaarikombinaatiosta — approksimoidaan ainoastaan yhdellä Slaterin determinantilla eli yllä esitetty sarja katkaistaan. 
Tarkoituksena Hartree–Fock-menetelmässä on löytää perustilaa vastaava aaltofunktio yhden determinantin approksimaation rajoissa. Oikean aaltofunktion muotoa ei tunneta, mutta se voidaan löytää soveltamalla variaatioperiaatetta. Elektronisysteemin kokonaisenergia saadaan Hamiltonin operaattorin odotusarvosta. Approksimoimalla aaltofunktiota Slaterin determinantilla kokonaisenergia on
formula_9
missä nimittäjään on kirjoitettu integraali normittamista varten. Jatkossa aaltofunktion oletetaan olevan normitettu, eli :formula_10. formula_11 on niin kutsuttu yksielektronioperaattori, joka sisältää Hamiltonin operaattorin termit, jotka riippuvat ainoastaan yhdestä elektronista (kineettinen energia sekä vuorovaikutus atomiydinten kanssa). formula_12 kaksielektronioperaattori, joka sisältää elektronien väliset vuorovaikutukset.
formula_13
formula_14. Kokonaisenergiafunktionaalin minimi löydetään varioimalla spinorbitaaleja eli tekemällä pieni muutos orbitaliin, formula_15. Orbitaaleja varioitaessa vaaditaan, että spinorbitaalit säilyvät ortogonaalisina eli formula_16. Tämä ehto toteutetaan Lagrangen kertoimien avulla:
formula_17
Muutos energiassa on
formula_18
Koska suoritettavat variaatiot eri orbitaaleihin ovat toisistaan riippumattomia, täytyy energiafunktionaalin differentiaalin tulla nollaksi jokaisella orbitaalilla formula_19. Tämä johtaa niin kutsuttuihin Hartree–Fock-yhtälöihin:
formula_20
Yllä formula_21 on niin kutsuttu Fock-operaattori. formula_11 on yksielektronioperaattori joka on määritelty yllä. formula_23 on Coulomb-operaattori, joka huomioi elektronien välisen sähköstaattisen vuorovaikutuksen. formula_24 on vaihto-operaattori (exchange operator), jolle ei ole olemassa klassista vastinetta. Operaattori on seuraus fermionisysteemin aaltofunktion antisymmetrisyydestä. Huomioi, että kun juokseva indeksi formula_25 Coulomb- ja vaihto-operaattorit kumoavat toisensa, toisin sanoen elektroni ei vuorovaikuta itsensä kanssa. Tämä on intuitiivinen ja triviaalin oloinen tulos, mutta ongelma muun muassa tiheysfunktionaaliteoriassa.
Hartree–Fock-yhtälöt ovat periaatteessa Schrödingerin yhtälöt kirjoitettuna elektronille, joka liikkuu toisten elektronien muodostamassa keskimääräisessä kentässä. Fock-operaattori kirjoitetaan jokaiselle orbitaalille erikseen, mutta operaattori riippuu kuitenkin kaikista muista orbitaaleista (Coulomb- ja vaihto-operaattorit). Jotta perustilan energia saataisiin ratkaistua, joudutaan orbitaaleista ensin luomaan järkevä alkuarvaus, jota sitten iteratiivisella menetelmällä optimoidaan. Tätä kutsutaan itseytyvän kentän menetelmäksi (self-consistent field, SCF). Alkuarvaus luodaan yleensä lineaarikombinaationa atomiorbitaaleista (linear combination of atomic orbitals, LCAO). 
Slaterin determinantin separoituva luonne johtaa vaillinaiseen elektroni–elektroni-korrelaatioon, mikä voidaan havaita muun muassa Coulomb-operaattorissa. Muiden elektronien muodostama varauspilvi saadaan suoraan laskemalla elektronien elektronitiheydet yhteen, aivan kuten riippumattomien satunnaistapahtumien kohdalla. Hartree–Fock-mallista puuttuva korrelaatioenergia on suureena hieman hankala määritellä fysikaalisesti intuitiivella tavalla. Yleinen, ehkä hieman kiistelty, määritelmä on määrittää korrelaatioenergia erotuksena oikean, eksaktin energian ja Hartree–Fock energian välillä.
formula_26
Täsmällisen elektronikorrelaation puuttuminen ajaa elektronit lähemmäs toisiaan, minkä vuoksi Coulombinen repulsio systeemissä kasvaa todellista suuremmaksi. Tämän vuoksi Hartree–Fock kokonaisenergia on aina oikeaa perustilan energiaa suurempi. Korrelaatioenergia on yleensä hyvin pieni suhteessa systeemiin ja kattaa suurimmillaankin vain muutaman prosentin kokonaisenergiasta, yleensä alle prosentin. Kemiallisissa systeemeissä tämä pieni erotus on kuitenkin keskeinen.