Hartree-Fock-approksimationen (HF) (også Hartree-Fock-metoden eller den selvkonsistente felt-metode (engelsk: self-consistent field method, SCF) er en variationsmetode inden for kvantekemi. Den bruges til at løse den tidsuafhængige Schrödinger-ligning for et system af flere elektroner omkring atomkerner, da en analytisk løsning ikke findes. Approksimation behandler hver elektron som om, den bevæger sig i et middelfelt skabt af alle de andre elektroner.[1]
For et system af
elektroner i fx et molekyle antages det, at den samlede bølgefunktion
kan approksimeres som en Slater-determinant, hvor hver elektron beskrives af såkaldte spinorbitaler, der er elementer i Slater-determinanten:

Det kan vises, at der for hver spinorbital er en såkaldt Hartree-Fock-ligning:

hvor
er Fock-operatoren for elektron 1. Den er givet ved den isolerede elektrons Hamilton-operator
plus et effektivt potentiale
fra de andre elektroner:

mens
er en energi associeret med spinorbitalen.
Hartree-Fock-metoden går ud på at løse dette ligningssystem for derved at finde en approksimativ løsning til
.[2][1]
I de følgende to afsnit vil problemet først få en mere deltaljeres behandling, hvorefter Hartree-Fock-ligningen vil blive udledt.
Et generelt system af atomkerner og elektroner beskrives tidsuafhængigt med Schrödinger-ligningen:

hvor
er Hamilton-operatoren,
er bølgefunktionen, og
er energien. Hamilton-operatoren er summen af alle bidrag til den kinetiske og potentielle energi:

Første led
er elektronernes kinetiske energi,
er atomkernernes kinetiske energi,
er elektron-elektron-frastødning,
er atomkerne-atomkerne-frastødning, mens
er elektron-atomkerne-tiltrækning. Vha. Born-Oppenheimer-approksimationen kan en simplere Hamilton-operator opskrives, der kun inkluderer termer, der afhænger af elektronerne:

Skrevet fuldt ud repræsenteres alle interaktionerne med Coulomb-potentialer:

For overskuelighedens skyldes bruges fra nu af desuden atomare enheder, hvor mange af faktorerne i ligningen sættes til 1:

Dette kan gøres endnu mere kompakt, idet de første to led kan betragtes som en sum af Hamilton-operatoren
for hver elektron uden vekselvirkning:

Schrödinger-ligningen for elektronerne er altså:
Det er denne ligning, der kan løses med Hartree-Fock-metoden.[2][1]
Grunden, til at en Slater-determinant bruges til at approksimere bølgefunktionen, er, at den er en funktion af samtlige elektroner, og at den er antisymmetrisk, hvilket vil sige, at den får et negativt fortegn, hvis to elektroner ombyttes. Dette er nødvendigt, da elektroner er fermioner. I princippet kan bølgefunktionen skrives eksakt som en uendelig sum af Slater-determinanter, men dette er ikke praktisk muligt. I stedet antages det, at bølgefunktionen kan beskrives tilstrækkeligt med én enkelt Slater-determinant. Denne approksimation er central for Hartree-Fock-metoden. Andre metoder såsom konfigurationsvekselvirkningsmetoden (CI) bruger flere Slater-determinanter.
Energien
i grundtilstanden derfor givet ved:
![{\displaystyle E_{0}[\Psi _{0}]=\left\langle \Psi _{0}\right|{\hat {H}}\left|\Psi _{0}\right\rangle =\left\langle \chi _{1}\chi _{2}\dots \chi _{N}\right|{\hat {H}}\left|\chi _{1}\chi _{2}\dots \chi _{N}\right\rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/8ebf4158b7c2071ac131723baebbc5d7a0ffc113)
hvor
er et funktional af bølgefunktionen. For at finde de optimale spinorbitaler, skal energien jf. variationsprincippet minimeres. Det kræves yderligere, at spinorbitalerne er ortonormale, hvilket vil sige, at det indre produkt skal være lig med Kroneckers delta:

Variationen af dette er

da Kroneckers delta er en konstant og derfor har en variation på nul. Vha. Lagrange-multiplikatorer kan denne betingelse tilføjes for samtlige kombinationer af spinorbitaler. Funktionalet
, der skal minimeres, er derfor:
![{\displaystyle L[{\chi _{i}}]=\left\langle \chi _{1}\chi _{2}\dots \chi _{N}\right|{\hat {H}}\left|\chi _{1}\chi _{2}\dots \chi _{N}\right\rangle -\sum _{a}\sum _{b}\varepsilon _{ab}\left(\left\langle \chi _{a}|\chi _{b}\right\rangle -\delta _{ab}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/65bc9494d325d70ca165c7e5e4c2a18193dd429b)
Variationen skal være nul, når spinorbitalerne er optimerede:
![{\displaystyle {\begin{aligned}0&=\delta L[{\chi _{i}}]=\delta E_{0}-\sum _{a}\sum _{b}\varepsilon _{ab}\left(\left\langle \delta \chi _{a}|\chi _{b}\right\rangle +\left\langle \chi _{a}|\delta \chi _{b}\right\rangle \right)\\0&=\delta L[{\chi _{i}}]=\delta E_{0}-\sum _{a}\sum _{b}\varepsilon _{ab}\left\langle \delta \chi _{a}|\chi _{b}\right\rangle -\sum _{a}\sum _{b}\varepsilon _{ab}\left\langle \chi _{a}|\delta \chi _{b}\right\rangle \\0&=\delta L[{\chi _{i}}]=\delta E_{0}-\sum _{a}\sum _{b}\varepsilon _{ab}\left\langle \delta \chi _{a}|\chi _{b}\right\rangle -\sum _{a}\sum _{b}\varepsilon _{ba}\left\langle \chi _{b}|\delta \chi _{a}\right\rangle \\0&=\delta L[{\chi _{i}}]=\delta E_{0}-\sum _{a}\sum _{b}\varepsilon _{ab}\left\langle \delta \chi _{a}|\chi _{b}\right\rangle -\sum _{a}\sum _{b}\left(\varepsilon _{ab}\left\langle \delta \chi _{a}|\chi _{b}\right\rangle \right)^{*}\\0&=\delta L[{\chi _{i}}]=\delta E_{0}-2\sum _{a}\sum _{b}\varepsilon _{ab}\left\langle \delta \chi _{a}|\chi _{b}\right\rangle \end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4c1c89df4109c97dfc3af152a5224cbb0b989647)
Det er her udnyttet, at summernes indekser er arbitrære, og at både integralerne og
er Hermitiske. Det andet led er allerede evalueret, men for at finde variationen af
skal størrelsen først udtrykkes direkte vha. spinorbitalerne.
For at finde grundtilstandsenergien anvendes det, at Slater-determinanten er en sum af Hartree-produkter, der gennemgår samtlige permutationer af elektroner i de forskellige orbitaler, og at leddene skifter fortegn for hver permutation. Dette kan skrives som:

hvor
er permutationsoperatoren, og summen er over antallet af samtlige permutationer. For
spinorbitaler, er der
mulige permutationer. Dette kan alternativt ses som en operator
, der virker på et Hartree-produkt for at gøre det antisymmetrisk:

hvor
derfor kaldes for antisymmetriseringsoperatoren:

Med denne notation kan
skrives som:

Vha. regnereglerne for
kan dette omskrives:

Udtrykket for Hamiltonoperatoren indsættes:

Det første led repræsenterer energien
for de isolerede elektroner, men det andet led er energien
pga. vekselvirkningen. Summerne kan simplificeres ved at betragte de enkelte led. Det første led for permutationen
for elektron
skrives ud som integraler:

Da
kun afhænger af elektron
, kan udtrykket skrives som et produkt af integraler:

Integralet med Hamilton-operatoren er den isolerede elektrons energi
. Spinorbitalerne er ortonormale, så de andre integraler givet 1:

For alle andre permutterede led vil der være mindst et integrale af to forskellige spinorbitaler, hvilket giver 0. Derfor er det første led udtrykket for
blot en sum af
:

Uden vekselvirkningerne er energien altså blot summen af de enkelte elektroners energi, hvilket også er forventeligt. Elektronernes koordinater
er blevet integrerede ud - resultatet ville være det samme, så længe orbitaler er ens - så
angiver nu, hvilke orbitaler der integreres.
Vekselvirkningsenergien kan evalueres på en lignende måde. For elektronerne
og
ved
gælder:

Alle integraler, der ikke har noget med
og
at gøre, giver 1:

Energibidraget for to elektroner er her altså forventningsværdien af den inverse afstand - pga. atomare enheder - imellem dem. Dette udtrykker blot den elektrostatiske frastødning
, som også kendes fra klassisk fysik.
Men det er ikke hele historien. Det tilsvarende led for
- en ombytning af elektron
og
giver:

Dette minder om det forrige bidrag, men elektronerne på højre side af afstanden er nu i forskellige orbitaler. Dobbeltintegralet er derfor ikke blot en gennemsnitlig afstand. Dette bidrag er ombytningsenergien
og findes ikke i den klassisk fysik.
Alle andre permutationer giver derimod integraler af forskellige spinorbitaler og giver derfor nul. Dvs. at

hvor
og
refererer til orbitalerne. Det ses, at
og
er ens, når
er ens:

Dvs. at summen ikke behøver at undgå dette led:

Faktoren
er tilføjet for at opveje, at summerne nu tæller alle interaktioner to gange.
Den samlede energi
er derfor:

Der altså tre forskellige interaktioner, der bidrager til energien.[2][1]
Variationen af
er nu en sum af variationer:

For
er variationen:

Som forklaret tidligere er det arbitrært, hvilken elektron der integreres over. For nemhedens skyld kaldes den her blot for elektron 1.
Da Hamilton-operatoren er Hermitisk, er de to led ens:

For
integreres der over elektron 1 og elektron 2:

Her er den Hermitiske egenskab igen udnyttet. For
findes et tilsvarende udtryk:

Her er det yderligere udnyttet, at det er arbitrært, om elektronen kaldes 1 eller 2. Samlet bliver variationen derfor:

Det er her igen udnyttet, at integrationsvariablerne og summernes indekser er arbitrære.
Variationen inkl. ortonormalitetsbetingelsen er derfor:
![{\displaystyle {\begin{aligned}0=&\delta L[{\chi _{i}}]=2\sum _{a=1}^{N}\int \delta \chi _{a}^{*}(1){\hat {h}}(1)\chi _{a}(1)\mathrm {d} \tau _{1}+2\sum _{a=1}^{N}\sum _{b=1}^{N}\int \int \delta \chi _{a}^{*}(1)\chi _{b}^{*}(2){\frac {1}{r_{12}}}\chi _{a}(1)\chi _{b}(2)\mathrm {d} \tau _{1}\mathrm {d} \tau _{2}-2\sum _{a=1}^{N}\sum _{b=1}^{N}\int \int \delta \chi _{a}^{*}(1)\chi _{b}^{*}(2){\frac {1}{r_{12}}}\chi _{a}(2)\chi _{b}(1)\mathrm {d} \tau _{1}\mathrm {d} \tau _{2}-2\sum _{a=1}^{N}\sum _{b=1}^{N}\varepsilon _{ab}\int \delta \chi _{a}^{*}(1)\chi _{b}(1)\mathrm {d} \tau _{1}\\0=&\sum _{a=1}^{N}\left(\int \delta \chi _{a}^{*}(1){\hat {h}}(1)\chi _{a}(1)\mathrm {d} \tau _{1}+\sum _{b=1}^{N}\int \int \delta \chi _{a}^{*}(1)\chi _{b}^{*}(2){\frac {1}{r_{12}}}\chi _{a}(1)\chi _{b}(2)\mathrm {d} \tau _{1}\mathrm {d} \tau _{2}-\sum _{b=1}^{N}\int \int \delta \chi _{a}^{*}(1)\chi _{b}^{*}(2){\frac {1}{r_{12}}}\chi _{a}(2)\chi _{b}(1)\mathrm {d} \tau _{1}\mathrm {d} \tau _{2}-\sum _{b=1}^{N}\varepsilon _{ab}\int \delta \chi _{a}^{*}(1)\chi _{b}(1)\mathrm {d} \tau _{1}\right)\\0=&\sum _{a=1}^{N}\int \delta \chi _{a}^{*}(1)\left({\hat {h}}(1)\chi _{a}(1)+\sum _{b=1}^{N}\int \chi _{b}^{*}(2){\frac {1}{r_{12}}}\chi _{b}(2)\mathrm {d} \tau _{2}\chi _{a}(1)-\sum _{b=1}^{N}\int \chi _{b}^{*}(2){\frac {1}{r_{12}}}\chi _{a}(2)\mathrm {d} \tau _{2}\chi _{b}(1)-\sum _{b=1}^{N}\varepsilon _{ab}\chi _{b}(1)\right)\mathrm {d} \tau _{1}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a5ecac7ba9e1f6b1e317b596c2e93b0d2548880b)
Integralerne i midten kan omskrives til operatorer, der virker på
:
Pr. definition ændrer
altså spinorbitalen fra
til
. Dermed bliver udtrykket:
![{\displaystyle {\begin{aligned}0=&\sum _{a=1}^{N}\int \delta \chi _{a}^{*}(1)\left({\hat {h}}(1)\chi _{a}(1)+\sum _{b=1}^{N}{\hat {J}}_{b}(1)\chi _{a}(1)-\sum _{b=1}^{N}{\hat {K}}_{b}(1)\chi _{a}(1)-\sum _{b=1}^{N}\varepsilon _{ab}\chi _{b}(1)\right)\mathrm {d} \tau _{1}\\0=&\sum _{a=1}^{N}\int \delta \chi _{a}^{*}(1)\left(\left({\hat {h}}(1)+\sum _{b=1}^{N}\left[{\hat {J}}_{b}(1)-{\hat {K}}_{b}(1)\right]\right)\chi _{a}(1)-\sum _{b=1}^{N}\varepsilon _{ab}\chi _{b}(1)\right)\mathrm {d} \tau _{1}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9d26561691b3f78b52db3678b02c8b5d06e60ec6)
Det samlede udtryk skal være nul for alle variationer af
. For at det er muligt, må udtrykket i parentesen være nul:
![{\displaystyle {\begin{aligned}\left({\hat {h}}(1)+\sum _{b=1}^{N}\left[{\hat {J}}_{b}(1)-{\hat {K}}_{b}(1)\right]\right)\chi _{a}(1)-\sum _{b=1}^{N}\varepsilon _{ab}\chi _{b}(1)&=0\\\left({\hat {h}}(1)+\sum _{b=1}^{N}\left[{\hat {J}}_{b}(1)-{\hat {K}}_{b}(1)\right]\right)\chi _{a}(1)&=\sum _{b=1}^{N}\varepsilon _{ab}\chi _{b}(1)\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd9fbe641a212ac6f23008cc51ed000cd7521ff6)
På højre side kan summen ophæves ved at rotere spinorbitalerne.[2][1]
Spinorbitalerne i Hartree-Fock-approksimationen skal altså selv overholde en egenværdi-ligning kaldet Hartree-Fock-ligningen:
![{\displaystyle \left({\hat {h}}(1)+\sum _{b=1}^{N}\left[{\hat {J}}_{b}(1)-{\hat {K}}_{b}(1)\right]\right)\chi _{a}(1)=\varepsilon _{a}\chi _{a}(1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10b7b0b6d9d18db86e5741b93b7fd808eea062ec)
Denne ligning minder om Schrödinger-ligningen for en enkelt elektron, hvor elektronens bølgefunktion er
. Operatoren er Fock-operatoren
, der består af den isolerede elektrons Hamilton-operator plus korrektioner, der tager højde for elektron-elektron-interaktionerne:
Den samlede Slater-determinant skal altså bestå af spinorbitaler, der er egenfunktioner til Fock-operatoren. Interaktionsbidraget kan også ses som et effektivt potentiale kaldet Hartree-Fock-potentialet:
![{\displaystyle {\hat {v}}^{\text{HF}}(1)=\sum _{b=1}^{N}\left[{\hat {J}}_{b}(1)-{\hat {K}}_{b}(1)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/52a96f7fb251dd2e8ba406c4167b40191eaf674a)
Selvom
virker på spinorbitalerne, ses det dog også, er operatoren afhænger af spinorbitalerne pga. definitionerne på
og
. Hartree-Fock-ligningerne må derfor nødvendigvis løses vha. en iterativ metode såsom en fikspunktsiteration. For at optimere spinorbitalerne skrives de normalt som en lineær kombination af såkaldte basisfunktioner, der er centrerede omkring hvert enkelt atom. Jo større basissættet er, jo bedre bliver spinorbitalerne, men jo mere omfattende bliver udregningen også.[2][1]
Hvordan, Hartree-Fock-ligningerne løses, afhænger af, hvilke spinorbitaler benyttes. Da elektronerne er fermioner, kan de ikke eksisterer i ens spinorbitaler, men to elektroner kan godt eksistere i den samme rumlige orbital
, hvis de blot har forskellige spin
og
. Inden for den såkaldte begrænsede Hartree-Fock-metode antages det, at
- og
- elektronerne netop har ens rumlige orbitaler, mens det i den ubegrænsede Hartree-Fock-metode tillades, at de rumlige orbitaler er forskellige. Den begrænsede metode er best egnet til systemer med et lige antal elektroner (closed shell), mens den ubegrænsede metode er bedre til at håndtere et ulige antal elektroner (open shell).[3]
- ^ a b c d e f Szabo, Attila; Ostlund, Neil S. (1996). "Chapter 3. The Hartree-Fock Approximation". Modern Quantum Chemistry (revideret 1. udgave). Dover Publications. s. 108-230. ISBN 0486691861.
- ^ a b c d e Rauk, Arvi (2001). "Appendix A: Derivation of Hartree-Fock Theory". Orbital Interaction Theory of Organic Chemistry (e-bog) (2. udgave). s. 218-226. ISBN 0-471-22041-8. Hentet 31. januar 2021.
- ^ Szabo, Attila; Ostlund, Neil S. (1996). "3.4 Restricted Closed-Shell Hartree-Fock: The Roothaan Equations". Modern Quantum Chemistry (revideret 1. udgave). Dover Publications. s. 131-132. ISBN 0486691861.