Vaihto- ja korrelaatiofunktionaali

Vaihto- ja korrelaatiofunktionaali (exchange and correlation functional) on tiheysfunktionaaliteoriassa esiintyvä, kokonaisenergiafunktionaaliin liittyvä tuntematon funktionaali, jota käytännön laskuissa joudutaan approksimoimaan. Erilaisia approksimaatiomenetelmiä on esitetty alla.

Toisin kuin aaltofunktiomenetelmässä, jossa Schrödingerin yhtälöön on olemassa useita, teoreettisesti täsmällisiä ratkaisuja, joilla oikeata ratkaisua voidaan systeemaattisesti lähestyä (esimerkiksi konfiguraatiovuorovaikutus tai häiriöteoria menetelmät) ei tiheysfunktionaaliteoriassa ole olemassa systemaattista menetelmää oikean tai paremman funktionaalin löytämiseksi. Funktionaalille on mahdollista johtaa erilaisia reunaehtoja, joita funktionaalin tulisi toteuttaa, mutta paradoksaalisesti useampia reunaehtoja toteuttava funktionaali ei ole välttämättä parempi funktionaali. Osa olemassa olevista funktionaaleista on kehitetty ilman parametreja, etsien funktionaalia joka toteuttaa tiettyjä reunaehtoja, osa taas on luotu parametrisoimalla, jolloin funktionaali kuvaa paremmin tietynlaista systeemiä (esimerkiksi heikosti vuorovaikuttavia systeemiä tai pinnalla tapahtuvia ilmiöitä).

TaustaaMuokkaa

Hohenbergin ja Kohnin teoreemoihin nojaten, monielektronisysteemin perustilan kokonaisenergia voidaan lausua elektronitiheyden funktionaalina:

 

missä   on elektronien kineettinen energia,   on elektronien keskinäinen vuorovaikutus ja   on ulkoiseen potentiaaliin liittyvä potentiaalienergia (esimerkiksi molekyylisysteemissä atomien ydinten muodostama potentiaali). Ulkoinen potentiaali riippuu tutkittavasta systeemistä, mutta kineettinen energia ja elektronien välinen vuorovaikutus ovat universaaleita funktionaaleja. Niiden summaa kutsutaan Hohenberg–Kohn-funktionaaliksi  

 

missä   on Kohn–Sham-menetelmällä määritettävä kineettinen energia,   on elektronitiheyksien välinen coulombinen repulsio (ns. Hartree-repulsio) ja   on vaihto- ja korrelaatiofunktionaali, jota menetelmässä joudutaan approksimoimaan. Funktionaali sisältää osia kineettisestä energista, jota Kohn–Sham-menetelmässä ei voida kuvata (Kohn–Shamin kineettinen energia on aina pienempi tai yhtäsuuri kuin oikea kineettinen energia).

Adiabaattinen kytkentäMuokkaa

Adiabaattinen kytkentä (adiabatic connection) antaa vaihto- ja korrelaatiofunktionaalille eksaktin muodon ja sitä käytetään monesti lähtökohtana XC-funktionaalia approksimatiivista funktionaalia rakennettaessa. Adiabaattinen kytkentä muistuttaa monessa mielessä häiriöteoriaa. Ideana tuoda elektronien välinen vuorovaikutus Hohenberg--Kohn funktionaaliin asteittain, kirjoittamalla funktionaali muodossa

 

missä   kineettisen energian operaattori ja   on elektronien välinen vuorovaikutus. Kun  , funktionaali kuvaa Kohn--Sham referenssitilaa, jossa elektronien välillä ei ole vuorovaikutuksia ja kun   funktionaali vastaa oikeata vuorovaikuttavaa systeemiä. Vuorovaikutus kytketään päälle (eli  ) hitaasti, siten että systeemin perustilan elektronitiheys   säilyy muuttumattomana. Tarkastelemalla vuorovaikuttavan ja Kohn--Sham systeemien funktionaaleja (ja hyödyntämällä integraalin ominaisuuksia) saadaan

 

Kohn--Sham systeemille  , joten yhtälön vasemmalla puolella kineettisen energian termit supistuvat pois ja jäljelle jää

 
 

Vaihto- ja korrelaatioenergia voidaan kirjoittaa myös elektronitiheyden   ja vaihto- ja korrelaatioaukon   välisenä vuorovaikutuksena

 


Paikallistiheysapproksimaatio – LDAMuokkaa

Homogeeniselle elektronikaasulle on mahdollista ratkaista vaihto- ja korrelaatiofunktionaali hyvin tarkasti. Paikallistiheysapproksimaatiossa (local density approximation, LDA) tuntematon vaihto- ja korrelaatiofunktionaali korvataan homogeenisen elektronikaasun funktionaalilla. Tasainen elektronikaasu ei luonnollisestikaan kuvaa mitään realistista systeemiä, mutta approksimaatio toimii kohtalaisen hyvin systeemeissä, joissa muutokset elektronitiheydessä ovat kohtalaisen pieniä, kuten esimerkiksi kiinteä metalli. Sen sijaan esimerkiksi kemian kannalta oleellisissa systeemissä, lähinnä molekyyleissä, elektronitiheydessä voi esiintyä rajujakin muutoksia ja LDA:lla määritettävät energiat ovat yleensä kaukana oikeasta.

Paikallistiheysapproksimaatiossa vaihto- ja korrelaatioenergia kirjoitetaan vaihto- ja korrelaatio-osien summana. Homogeenisen elektronikaasun vaihto-osa voidaan ratkaista eksaktisti ja korrelaatio-osaa eri elektronitiheyksissä on arvioitu laskennallisesti [Monte Carlo -simulaatio|Monte Carlo] -simulaatioiden avulla [1].

 

Ilmiselvistä puutteistaan huolimatta LDA soveltuu melko hyvin kiinteän olomuodon mallintamiseen ja antaa muun muassa kohtalaisen tarkkoja molekyylirakenteita. LDA epäonnistuu kuitenkin energetiikassa ja esimerkiksi sidosenergiat voivat poiketa todellisista arvoista useita elektronivoltteja (1 eV ≈ 100 kJ mol−1). Moderneissa funktionaaleissa LDA:n puutteita korjataan huomioimalla elektronitiheyden epähomogeenisuus elektronitiheyden gradientilla tai korkeammilla derivaatoilla. Vaikka teoriassa oikea funktionaali ei riipukaan elektronien spinistä, kirjoitetaan elektronitiheys yleensä kahden spintiheyden summana,

 

Gradienttikorjaus – GGAMuokkaa

Luonteva parannus LDA-malliin, on huomioida elektronitiheyden epähomogeenisuus tekemällä funktionaalista riippuvainen ympäristöstään. Huomioimalla funktionaalissa elektronitiheyden gradientti saadaan niin kutsuttu yleistetty gradienttiapproksimaatio (generalized gradient approximation, GGA). Jos funktionaalissa huomioidaan korkeampia derivaattoja, kutsutaan funktionaalia meta-GGA funktionaaliksi. Kuten LDA-funktionaalitkin, GGA-funktionaalit kirjoitetaan yleensä vaihto- ja korrelaatio-osien summana. Gradientti huomioidaan niin kutsutaan redusoidun gradientin   avulla, joka periaatteessa kuvaa elektronitiheyden poikkeavuutta homogeenisuudestä. Gradienttikorjattu funktionaali rakennetaan yleensä paikallistiheysapproksimaation päälle.

 

Osa gradienttikorjatuista funktionaaleista on rakennettu ab initio eli nojaten teoriaan, osa on parametrisoitu laskuihin tai sovitettu dataan. Yleinen ei-empiirinen yleisfunktionaali on Perdew–Burke–Ernzerhofin PBE-funktionaali [2], jossa esimerkiksi vaihtoenergian lauseke on muotoa

 

missä funktion   parametrit on valittu siten, että   (Lieb–Oxford-raja).

LuettavaaMuokkaa

  • W. Koch, M. C. Holthausen: A chemist's guide to density functional theory. Wiley-Vch Weinheim, 2001. ISBN 978-3-527-30372-4.
  • C. J. Cramer: Essentials of computational chemistry: theories and models. John Wiley & Sons, Ltd, 2004. ISBN 978-0-470-09182-1.
  • D. S. Sholl, J. A. Steckel: Density functional theory, a practical introduction. John Wiley & Sons, Inc, 2009. ISBN 978-0-470-37317-0.

LähteetMuokkaa

  1. D. M., Ceperley, B. J. Alder: Ground State of the Electron Gas by a Stochastic Method. Physical Review Letters, 1980, nro 45, s. 566-569. Artikkelin verkkoversio Viitattu 4.3.2014.
  2. J. Perdew, K. Burke, M. Ernzerhof: Generalized gradient approximation made simple. Physical Review Letters, 1996, nro 77, s. 3865-3868. Artikkelin verkkoversio Viitattu 12.4.2014.