Lähdekoodi sivulle Tiheysfunktionaaliteoria

'''Tiheysfunktionaaliteoria''' (density functional theory, '''DFT''') on [[laskennallinen fysiikka|laskennallisessa fysiikassa]] ja [[laskennallinen kemia|kemiassa]] käytetty [[kvanttimekaniikka|kvanttimekaaninen]] menetelmä mikroskooppisten systeemien elektronirakenteen tutkimiseen. Menetelmää käytetään erityisesti [[kiinteän olomuodon fysiikka|kiinteän olomuodon fysiikassa]] sekä [[kvanttikemia]]ssa. Perustan teorialle loivat [[Pierre Hohenberg]] ja [[Walter Kohn]] <ref name='hohenberg'>{{Lehtiviite | Tekijä =  P. Hohenberg, W. Kohn | Otsikko = Inhomogeneous electron gas| Julkaisu = Physical review| Ajankohta = 1964 | www-teksti = http://journals.aps.org/pr/abstract/10.1103/PhysRev.136.B864 | Viitattu = 3.3.2014 | Lopetusmerkki = }}</ref> vuonna 1964 osoittaessaan, että elektronisysteemin perustilan elektronitiheys määrittää systeemin tilan yksiselitteisesti. Tällöin systeemin [[sisäenergia]] on mahdollista kirjoittaa [[aaltofunktio]]n sijasta elektronitiheysfunktion ''funktionaalina'', jolloin laskennalliset kustannukset kutistuvat merkittävästi. Kiinteän olomuodon fysiikassa tiheysfunktionaaliteoria on aina ollut suositumpi kuin aaltofunktiopohjaiset menetelmät ja kiinteän olomuodon tutkimisessa tiheysfunktionaaliteoriaa on käytetty jo 1970-luvulta asti. Kvanttikemistien suosioon menetelmä kuitenkin nousi vasta 1980-luvun lopulla uusien, kehittyneempien funktionaalien myötä. Vuonna 1998 [[Nobelin kemianpalkinto|kemian Nobelin palkinto]] myönnettiin Kohnille teorian keksimisestä ja [[John Pople]]lle laskennallisten menetelmien kehittämisestä.
<ref name='review'>{{Lehtiviite | Tekijä =  K. Burke | Otsikko = Perspective on density functional theory | Julkaisu = Journal of Chemical Physics| Ajankohta = 2012 | www-teksti = http://scitation.aip.org/content/aip/journal/jcp/136/15/10.1063/1.4704546 | Viitattu = 5.3.2014 }}</ref>

Tiheysfunktionaaliteoria on teoriassa eksakti, mutta käytännössä kokonaisenergialausekkeessa esiintyvän [[vaihto- ja korrelaatiofunktionaali]]n oikeaa muotoa ei tunneta, joten funktionaalia joudutaan approksimoimaan. Tämän seurauksena myös tulokset ovat riippuvaisia käytetystä approksimaatiosta. Approksimaatiosta seuraavia puutteita tiheysfunktionaaliteoriassa ovat muun muassa:
* Yksikään approksimatiivinen funktionaali ei sisällä heikkoa [[Dispersiovoima|dispersiovuorovaikutusta]], mikä on ongelma esimerkiksi hiilivetyjä mallinnettaessa. Dispersiota voidaan approksimoida esimerkiksi empiirisillä korjauksilla.<ref name='review'/><ref name='grimme'>{{Lehtiviite | Tekijä =  S. Grimme | Otsikko = Density functional theory with London dispersion corrections | Julkaisu = Wiley Interdisciplinary Reviews: Computational Molecular Science | Ajankohta = 2011 | www-teksti = http://onlinelibrary.wiley.com/doi/10.1002/wcms.30/abstract? | Viitattu = 5.3.2014 | Lopetusmerkki = }}</ref><ref name='grimme2'>{{Lehtiviite | Tekijä =  S. Grimme ''et al''. | Otsikko = A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu | Julkaisu = The Journal of Chemical Physics | Ajankohta = 2010 | www-teksti = http://scitation.aip.org/content/aip/journal/jcp/132/15/10.1063/1.3382344 | Viitattu = 5.3.2014}}</ref> Dispersion arvioimista suoraan elektronitiheydestä käsin on myös esitetty <ref name='langrethlundqvist'>{{Lehtiviite | Tekijä =  Y. Andersson, D. C. Langreth, B. I. Lundqvist | Otsikko = van der Waals Interactions in Density-Functional Theory | Julkaisu = Physical review letters | Ajankohta = 1996 | www-teksti = http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.76.102 | Viitattu = 5.3.2014}}</ref>.
* Perinteinen tiheysfunktionaaliteoria ei sisällä virittyneitä tiloja, koska virittyneen tilan elektronitiheys ei yksiselitteisesti määrää potentiaalia. Virittyneitä tiloja on kuitenkin mahdollista käsitellä ajasta riippuvan tiheysfunktionaaliteorian (time-dependent density functional theory, TD-DFT) avulla.<ref name='review'/> 
* Tiheysfunktionaaliteoria on tunnetusti surkea [[Rajaorbitaalit|HOMO–LUMO-välin]] eli ylimmän miehitetyn ja alimman miehittämättömän orbitaalin energiaeron arvioimisessa.

== Hohenberg–Kohn-teoreemat ==

Hohenberg–Kohn-teoreemat antavat teoreettisen oikeutuksen menetelmän käytölle. Ensimmäinen [[teoreema]] osoittaa, että systeemin perustilan elektronitiheys yksiselitteisesti määrää systeemin tilan. Toinen teoreema luo pohjan perustilan elektronitiheyden löytämiselle: oikeaa perustilan elektronitiheyttä voidaan lähestyä minimoimalla perustilan energiaa. Toisin sanoen oikea elektronitiheys antaa alimman mahdollisen perustilan energian. Tiheysfunktionaaliteoria on laskennallisesti huomattavasti kevyempi kuin monet korreloituneet [[Hartree-Fock-teoria|Hartree–Fock-johdannaiset]], kun ''N'':n elektronin systeemin ''3N''-ulotteisen aaltofunktion sijasta ratkaistaan vain 3-ulotteista elektronitiheyttä.

Monielektronisysteemeissä [[atomiydin|ydinten]] ja [[elektroni]]en liike yleensä erotetaan toisistaan (Born–Oppenheimer eli adiabaattinen approksimaatio), jolloin elektronit liikkuvat ydinten muodostamassa "staattisessa" potentiaalissa. Stationääritilaa <math>\Psi (\vec r_1 , ..., \vec r_N) = \Psi </math> voidaan kuvata ajasta riippumattoman Schrödingerin yhtälön avulla:

:<math> \hat H \Psi = (\hat T + \hat V_{ee} + \hat V_{ext}) \Psi = \left( - \sum _i ^N \frac{1}{2} \nabla ^2 + \sum _i ^N \sum _{j > i} \frac{1}{|\vec r_i - \vec r_j|} - \sum _i ^N \sum _{k} \frac{Z_k}{|\vec r_i - \vec R_k|} \right) \Psi = E \Psi </math>

Yllä <math> \hat H </math> on systeemiä kuvaava [[Hamiltonin operaattori]] ja <math> E </math> systeemin kokonaisenergia. <math> \vec R_k </math> on atomiytimen ''k'' paikkavektori ja <math> Z_k </math> ytimen varaus. <math> \hat T </math> on elektronien [[kineettinen energia]], <math>\hat V_{ee}</math> elektronien välinen coulombinen vuorovaikutus ja <math>\hat V_{ext}</math> elektronien ja atomiydinten vuorovaikutus eli atomien virittämä ulkoinen potentiaali. Kineettisen energian <math> \hat T </math> ja elektronirepulsion <math> \hat V_{ee} </math> operaattorit ovat samoja kaikissa ''N''-elektronisysteemeissä, mutta ulkoinen potentiaali, jossa elektronit sijaitsevat vaihtuu systeemistä toiseen. Yllä on käytetty niin kutsuttuja atomiyksiköitä, jolloin vakioita ei tarvitse pyörittää erikseen mukana. Atomiyksiköissä energian yksikkö on Hartree, ''E<sub>h</sub>'' = 27,21138386 eV = 4,35974394·10<sup>−18</sup> J.<ref name='crc'>{{Kirjaviite | Tekijä =  W. M. Haynes, D. R. Lide, T. J. Bruno | Nimeke = CRC Handbook of Chemistry and Physics 2012-2013 | Vuosi = 2012 | Julkaisija = CRC Press | Tunniste = | Isbn = 978-1-4200-6679-1| Viitattu = 4.3.2014 }}</ref>

Ensimmäisen Hohenberg–Kohn teoreeman todistus on suoraviivainen ''reductio ad absurdum'': oletetaan, että kahden toisistaan eroavan ulkoisen potentiaalin, <math>\hat V_{ext} </math> ja <math>\hat V_{ext} '</math>, perustilan elektronitiheys <math>\rho(\vec r) = |\Psi (\vec r)|^2 =  |\Psi' (\vec r)|^2</math> olisi sama. Jos potentiaalit poikkeavat toisistaan enemmän kuin pelkän vakion verran, myös potentiaaleja vastaavat eri aaltofunktiot, <math>\Psi</math> ja <math>\Psi '</math>. Variaatioperiaatteen mukaan oikea aaltofunktio antaa alhaisimman perustilan energian. Jos siis operoidaan aaltofunktiota <math>\Psi'</math> toisen systeemin Hamiltonin operaattorilla <math>\hat H</math> saadaan perustilan energia, joka on oikeaa perustilan arvoa suurempi:

::<math>E_0 < \left\langle \Psi' | \hat H | \Psi ' \right\rangle = \left\langle \Psi' | \hat H' | \Psi ' \right\rangle + \left\langle \Psi' | \hat H - \hat H' | \Psi ' \right\rangle = E_0 ' + \left\langle \Psi' | (\hat T +\hat V_{ee} + \hat V_{ext}) - (\hat T +\hat V_{ee} + \hat V_{ext} ') | \Psi ' \right\rangle</math>

:<math> \Rightarrow E_0  < E_0 ' + \int \rho (\vec r) \left( \hat V_{ext} - \hat V_{ext}' \right) d \vec r </math>

Jos kuitenkin aaltofunktiota <math>\Psi</math> operoidaan toisella Hamiltonin operaattorilla <math>\hat H'</math> saadaan tulokseksi

:<math> E_0'  < E_0 - \int \rho (\vec r) \left( \hat V_{ext} - \hat V_{ext}' \right) d \vec r </math>

Yhdistämällä kaksi edellistä yhtälöä, saadaan tulokseksi <math> E_0 + E_0' < E_0' + E_0 \Leftrightarrow 0 < 0</math> mikä on ristiriita, joten alun oletus on virheellinen: perustila riippuu yksiselitteisesti elektronitiheydestä.0

== Kohn–Sham-menetelmä ==

Vaikka Hohenbergin ja Kohnin teoreemat osoittavat, että elektronitiheys on perustilaan liittyvä yksiselitteinen [[suure]], eivät teoreemat kuitenkaan auta elektronitiheyden löytämisessä. Käytännön menetelmän perustilan löytymiseksi esittivät Kohn ja Sham vuonna 1965 <ref name='kohnsham'>{{Lehtiviite | Tekijä =  W. Kohn, L. J. Sham | Otsikko = Self-Consistent Equations Including Exchange and Correlation Effects | Julkaisu = Physical review| Ajankohta = 1965 | www-teksti = http://journals.aps.org/pr/abstract/10.1103/PhysRev.140.A1133 | Viitattu = 3.3.2014 | Lopetusmerkki = }}</ref>. Kohnin ja Shamin menetelmässä perustilan energia ratkaistaan kahden rinnakkaisen systeemin avulla, joista toinen laskee systeemin kokonaisenergian elektronitiheyden avulla samalla kun toinen määrittää elektronien kineettisen energian Kohn–Sham-orbitaalien avulla. Jälkimmäinen systeemi koostuu [[Hartree-Fock-teoria|Hartree–Fock-menetelmän]] kaltaisista ei-vuorovaikuttavista yksielektroniorbitaaleista, jotka sijaitsevat ulkoisessa potentiaalissa. Orbitaalien avulla kineettisen energian osuus on mahdollista laskea melko tarkasti.

Hohenberg–Kohnin teoreemojen perusteella on perustilan energia elektronitiheyden ''funktionaali'':

:<math> E_0 [\rho _0] = T [\rho _0] + E_{ee} [\rho _0] + E_{ext} [\rho _0]</math>

Yllä <math> T [\rho ] </math> on elektronien kineettinen energia,  <math> E_{ee} [\rho ]</math> on elektronien välisen vuorovaikutuksen [[potentiaalienergia]] ja <math> E_{ext} [\rho ]</math> ulkoiseen potentiaaliin liittyvä energia, <math> E_{ext} [\rho ] = \int \rho (\vec r)  V _{ext} d\vec r </math>. Kineettisen energian ja elektroni–elektroni-termin energia on sama kuin kineettisen ja elektroni–elektroni-operaattorin [[odotusarvo]]:

:<math> \left\langle \Psi | \hat T + \hat V_{ee} | \Psi  \right\rangle = T[\rho] + E_{ee}[\rho] = F_{HK} [\rho] </math>

missä <math> F_{HK} [\rho] </math> on niin kutsuttu Hohenberg–Kohn funktionaali. Koska kineettinen energia eikä elektronien välinen repulsio ovat systeemistä riippumattomia funktionaaleja, on Hohenberg–Kohn funktionaali riippumaton systeemistä. Vaikka tiheysfunktionaaliteoria on itsessään eksakti teoria, on ongelmana juuri Hohenberg–Kohn funktionaali, jolle ei ole olemassa lauseketta vaan funktionaalia joudutaan approksimoimaan.

Funktionaalin elektroni–elektroni-osa voidaan joka tapauksessa jakaa kahteen selvään komponenttiin: varaustiheyksien väliseen coulombiseen vuorovaikutukseen sekä ei-klassiseen osaan, joka sisältää elektronien vaihto- ja coulomb-korrelaation sekä varaustiheyden itseisvuorovaikutuksen:

:<math> F_{HK} [\rho] = T[\rho] + E_{ee}[\rho] = T[\rho] + \frac{1}{2} \int \int \frac{\rho(\vec r_1) \rho(\vec r_2)}{|\vec r_1 - \vec r_2|} d\vec r_1 d\vec r_2 + E_{ncl}[\rho] = T[\rho] + J[\rho]+ E_{ncl}[\rho] </math>

Ainoa suure, joka siis voidaan suoraan laskea elektronitiheydestä on kahden varaustiheyden välinen repulsio. Ei-klassiselle vaihto- ja korrelaatio-osalle ja kineettiselle osalle ei tunneta funktionaaleja. Kohn ja Sham esittivät <ref name='kohnsham'/>, että kineettistä energiaa voitaisiin approksimoida laskemalla se aivan kuten [[Hartree-Fock-teoria|Hartree–Fock-menetelmässä]] eli kineettisen energian operaattorin odotusarvona, <math>\left\langle \Psi | \hat T | \Psi  \right\rangle = \left\langle \Psi | - \frac{1}{2} \nabla ^2 | \Psi  \right\rangle </math>. Tätä varten täytyy luoda ei-vuorovaikuttava referenssisysteemi, joka koostuu yksielektroniorbitaaleista eli niin kutsutuista Kohn–Sham-orbitaaleista <math> \phi ^{KS} </math>. Hartree–Fock menetelmän tapaan referenssisysteemiä kuvaa yksi Slaterin determinantti ja jokaiselle orbitaalille kirjoitetaan Fock-operaattorin kaltainen Kohn–Sham-operaattori:

:<math> \hat f ^{KS} = -\frac{1}{2} \nabla ^2 + V_S(\vec r)</math>

Referenssisysteemin orbitaalien neliön tulee vastata oikean systeemin elektronitiheyttä, eli <math> \sum _i ^N |\phi ^{KS} _i |^2 = \rho (\vec r) </math>. Operaattorin potentiaalitermi <math> V_S(\vec r) </math> muodostaa yhteyden oikean sekä referenssisysteemin välille ja on muotoa

:<math> V_S(\vec r) = \int \frac{\rho(\vec r_2)}{|\vec r_1 - \vec r_2|}d\vec r_2 - \sum _k \frac{Z_k}{|\vec r_1 - \vec R_k|} + V_{XC} (\vec r_1) </math>

Ensimmäinen termi kuvaa orbitaalin vuorovaikutusta muiden orbitaalien kanssa, toinen termi kuvaa atomiydinten muodostamaa ulkoista potentiaalia ja <math> V_{XC} </math> on niin kutsuttu vaihto- ja korrelaatiopotentiaali. Koska jokaisen orbitaalin Kohn–Sham-operaattori riippuu muista orbitaaleista, täytyy ominaisarvotehtävä ratkaista iteratiivisesti itseytyvän kentän (self-consistent field, SCF) avulla. Potentiaalitermi <math> V_{XC} </math> saadaan ''[[Vaihto- ja korrelaatiofunktionaali|vaihto- ja korrelaatioenergian]]'' funktionaaliderivaatasta:

:<math> V_{XC} = \frac{\delta E_{XC}}{\delta \rho}</math>

missä <math> E_{XC} [\rho] = T _C [\rho] + E_{ncl} [\rho] </math>. Toisin sanoen vaihto- ja korrelaatioenergia sisältää uupuvan osan kineettisestä energiasta <math> T _C [\rho]</math>, kun kineettistä energiaa approksimoidaan Kohn–Sham-orbitaaleilla sekä elektronien vuorovaikutusten ei-klassiset osat <math> E_{ncl} [\rho] </math>. Tätä uupuvaa osaa kokonaisenergiafunktionaalista joudutaan approksimoimaan ja saatavat tulokset ovat hyvin riippuvaisia käytetystä funktionaalista, joten funktionaalin valinta on tehtävä huolella ja oltava perusteltavissa. Vaikka oikealle funktionaalille ei ole löydetty lauseketta, on sille johdettavissa ominaisuuksia joita vasten kehitettyjä funktionaaleja voidaan testata. Osa kehitetyistä funktionaaleista pohjautuu ainoastaan teoreettisiin argumentteihin (esimerkiksi Perdew–Burke–Ernzerhofin funktionaali PBE), toiset ovat parametrisoituja esimerkiksi kokeellisiin tuloksiin tai tarkkoihin kvanttimekaanisiin laskuihin.

==Tiheysfunktionaaliteorian ongelmakohtia==

[[Tiedosto:XC-kuvaaja.png|thumb|Systeemin kokonaisenergia elektronien lukumäärän funktiona. Keltaisella piirretty kuvaaja kuvaa eksaktia energiaa ja vihreällä approksimatiivista funktionaalia. Approksimatiivisella funktionaalilla funktio ei ole lineaarinen kokonaislukuelektronien välillä ja konveksin profiilin vuoksi kokonaisenergia voi olla alhaisempi murto-osa elektroneilla.]]

Huolimatta suuresta suosiostaan laskennallisten fyysikkojen ja kemistien keskuudessa, tiheysfunktionaaliteoriassa esiintyy silti lukuisia perustavanlaatuisia ongelmakohtia, jotka tuovat systemaattista virhettä tuloksiin. Virheet eivät sisälly varsinaiseen teoriaan vaan ovat seurausta vaihto- ja korrelaatiofunktionaalin approksimatiivisesta käsittelystä, joka johtaa muun muassa elektronien ylidelokalisoitumiseen sekä staattisen korrelaation virheelliseen kuvaamiseen.<ref name='insights'>{{Lehtiviite | Tekijä =  A. J. Cohen, P. Mori-Sánchez, and W. Yang | Otsikko = Insights into current limitations of density functional theory | Julkaisu = Science| Ajankohta = 2008 | www-teksti = http://www.sciencemag.org/content/321/5890/792.short | Viitattu = 24.7.2014 | Lopetusmerkki = }}</ref>

Approksimatiivisten funktionaalien korrelaatiovirheen vuoksi nykyisellä tiheysfunktionaaliteorialla ei voida luotettavasti kuvata degeneerisiä tai miltei degeneerisiä tiloja, kemiallisen sidoksen katkeamista tai systeemejä, joissa elektronikorrelaatio on erityisen voimakasta. Elektronien ylidelokalisoitumisella on vastaavasti lukuisia seurauksia: reaktioiden liian alhaiset aktivaatioenergiat sekä materiaalien [[Valenssivyö|vyövälin]], ionisten yhdisteiden dissosiaatioenergian ja varauksensiirtoreaktioiden reaktioenergian arvioiminen väärin.

Oikean systeemin kokonaisenergia on elektronien lukumäärän suhteen jatkuva funktio, mutta funktion derivaatta on epäjatkuva kokonaislukumäärien kohdalla johtuen elektronien diskreetistä luonteesta. Kokonaisenergia "murto-osaelektroneilla" saadaan siis interpoloimalla kahden kokonaislukuarvon väliltä. Approksimativiisilla funktionaaleilla kokonaisenergiassa ei epäjatkuvuuskohtia esiinny vaan funktio on konveksi ja systeemin energia on tämän vuoksi oikeaa energiaa alhaisempi.

==Luettavaa==

* {{Kirjaviite | Tekijä = W. Koch, M. C. Holthausen | Nimeke = A chemist's guide to density functional theory | Vuosi = 2001 | Julkaisija = Wiley-Vch Weinheim | Tunniste = | Isbn = 978-3-527-30372-4 | Viitattu = 4.3.2014 }}
* {{Kirjaviite | Tekijä = C. J. Cramer | Nimeke = Essentials of computational chemistry: theories and models | Vuosi = 2004 | Julkaisija = John Wiley & Sons, Ltd | Tunniste = | Isbn = 978-0-470-09182-1 | Viitattu = 4.3.2014 }}
* {{Kirjaviite | Tekijä = D. Cook | Nimeke = Handbook of computational quantum chemistry | Vuosi = 2005 | Julkaisija = Dover publications, Inc | Tunniste = | Isbn = 0-486-44307-8 | Viitattu = 4.3.2014 }}

==  Lähteet  ==

{{Viitteet}}

[[Luokka:Laskennallinen tiede]]