12.07.2015 Views

Teoria do Funcional da Densidade - Departamento de Física - UFMG

Teoria do Funcional da Densidade - Departamento de Física - UFMG

Teoria do Funcional da Densidade - Departamento de Física - UFMG

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Teoria</strong> <strong>do</strong> <strong>Funcional</strong> <strong>da</strong>Densi<strong>da</strong><strong>de</strong>Angélica MeloUniversi<strong>da</strong><strong>de</strong> Fe<strong>de</strong>ral <strong>de</strong> Minas Gerais<strong>Departamento</strong> <strong>de</strong> FísicaFDezembro 2009angelica@fisica.ufmg.br


Introdução:Mecânica quântica• Na mecânica quântica, a equação fun<strong>da</strong>mental que <strong>de</strong>screve o esta<strong>do</strong> <strong>de</strong> umsistema dinâmico é a equação <strong>de</strong> Schroedinger:on<strong>de</strong> (r, t) é a função <strong>de</strong> on<strong>da</strong> total <strong>do</strong> sistema e H éo hamiltoniano.• Para potenciais in<strong>de</strong>pen<strong>de</strong>ntes <strong>do</strong> tempo, po<strong>de</strong>mos separar as variáveis espaciaise temporal na solução:• Isso nos leva à equação <strong>de</strong> Schroedinger in<strong>de</strong>pen<strong>de</strong>nte <strong>do</strong> tempo:on<strong>de</strong> E é a energia total <strong>do</strong> sistema.


Introdução:O hamiltoniano <strong>de</strong> um sistema <strong>de</strong> muitos corpos• Para um sistema interagente <strong>de</strong> N elétrons e M núcleos:on<strong>de</strong>


Introdução:A aproximação <strong>de</strong> Born-Oppeinheimer• Aproximação Born-Oppeinheimer :Hipótese:•Núcleos são muito mais pesa<strong>do</strong>s <strong>do</strong> que os elétrons Núcleos fixosConseqüências:• Despreza-se a energia cinética <strong>do</strong>s núcleos.• Faz a interação núcleo-núcleo uma constante.CONSTANTE• As posições <strong>do</strong>s núcleos <strong>de</strong>finem um potencial externo para o movimento <strong>do</strong>selétrons.• O potencial elétron-núcleo passa a <strong>de</strong>pen<strong>de</strong>r parametricamente <strong>da</strong>s coor<strong>de</strong>na<strong>da</strong>snucleares.


Introdução:O hamiltoniano eletrônico• Assim ficamos com o hamiltoniano eletrônico:on<strong>de</strong> o conjunto <strong>de</strong> coor<strong>de</strong>na<strong>da</strong>s {R A } tem que ser <strong>da</strong><strong>do</strong>.• A equação <strong>de</strong> Schroedinger a ser resolvi<strong>da</strong> passa a ser:on<strong>de</strong> e = e (r 1 ... r N ) é a equação <strong>de</strong> on<strong>da</strong> eletrônica.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>• O que é a teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>?Essa teoria é uma alternativa <strong>de</strong> solução para o problema <strong>do</strong>hamiltoniano <strong>de</strong> muitos corpos interagentes.• A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong> promove a <strong>de</strong>nsi<strong>da</strong><strong>de</strong>eletrônica n(r) à variável chave, na qual o cálculo <strong>de</strong> to<strong>do</strong>s osoutros observáveis po<strong>de</strong>m ser basea<strong>do</strong>s.• O que nos garante que a <strong>de</strong>nsi<strong>da</strong><strong>de</strong> eletrônica po<strong>de</strong> exerceresse papel <strong>de</strong> variável fun<strong>da</strong>mental são <strong>do</strong>is teoremaspropostos por Hohenberg e Kohn.Eles formam a base teórica <strong>do</strong> DFT e são enuncia<strong>do</strong>s a seguir.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:Teoremas <strong>de</strong> Hohemberg-Kohn:Teoremas HKTeorema 1: O potencial externo v ext (r) senti<strong>do</strong> pelos elétrons é umfuncional único <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong> eletrônica <strong>do</strong> esta<strong>do</strong> fun<strong>da</strong>mental n 0(r).Teorema 2: A energia <strong>do</strong> esta<strong>do</strong> fun<strong>da</strong>mental E 0[n] é mínima para a<strong>de</strong>nsi<strong>da</strong><strong>de</strong> n 0(r) exata.• Assim, o primeiro teorema garante que po<strong>de</strong>mos escrever a energiacomo um funcional único <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong> eletrônica <strong>do</strong> esta<strong>do</strong>fun<strong>da</strong>mental.E = E[n o]• E o segun<strong>do</strong> garante que a energia <strong>do</strong> esta<strong>do</strong> fun<strong>da</strong>mental tem aproprie<strong>da</strong><strong>de</strong> <strong>do</strong> princípio variacional.E[n 0< E[n]


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:O <strong>Funcional</strong> Energia• Po<strong>de</strong>mos escrever a energia como um funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>eletrônica:E[n] = T[n] + U[n] +V ext[n],on<strong>de</strong>•T[n] é a energia cinética <strong>do</strong> sistema.•U[n] é o termo <strong>de</strong> interação entre elétrons.•V ext[n] é a energia <strong>de</strong>vi<strong>do</strong> ao potencial externo gera<strong>do</strong>pelos núcleos.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:formalismo partícula in<strong>de</strong>pen<strong>de</strong>nte• Kohn e Sham propuseram escrever o funcional energia <strong>de</strong>ntro <strong>do</strong> formalismo<strong>de</strong> partícula in<strong>de</strong>pen<strong>de</strong>nte:E KS [n] = T 0 [n] + U H [n] +V ext [n] +V xc [n]on<strong>de</strong>• T 0 [n] é a energia cinética <strong>de</strong> um sistema não interagente <strong>de</strong>elétrons com <strong>de</strong>nsi<strong>da</strong><strong>de</strong> n.•U H [n] é o termo <strong>de</strong> interação coulombiana média entre elétrons.• V xc [n] é o termo que contém a correção <strong>da</strong> energia cinética ( T -T 0 ) e <strong>da</strong> energia <strong>de</strong> Hartree (U – U H ), e é chama<strong>do</strong> <strong>de</strong> energia <strong>de</strong>troca-correlação.• Assim temos um formalismo <strong>de</strong> partícula in<strong>de</strong>pen<strong>de</strong>nte que inclui os efeitos <strong>de</strong>interação <strong>de</strong> muitos corpos, sen<strong>do</strong>, portanto, formalmente exato.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:As equações <strong>de</strong> Kohn-Sham• Explicitan<strong>do</strong> os termos <strong>do</strong> funcional energia temos:• Minimizan<strong>do</strong> esse funcional em relação a <strong>de</strong>nsi<strong>da</strong><strong>de</strong> com o vínculo <strong>do</strong> número<strong>de</strong> partículas constante, obtemos as equações que <strong>de</strong>screvem o esta<strong>do</strong> <strong>do</strong>sistema, as equações <strong>de</strong> Kohn-Sham:on<strong>de</strong>• i são os autovalores <strong>de</strong> Kohn-Sham e i as autofunções <strong>de</strong> Kohn-Sham.• Então minimizar o funcional energia em termos <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>eletrônica é equivalente a resolvermos uma equação <strong>do</strong> tipoScrhodinger com um potencial efetivo V ef.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:As equações <strong>de</strong> Kohn-Sham• Nesse sistema fictício <strong>de</strong> partículas in<strong>de</strong>pen<strong>de</strong>ntes, a <strong>de</strong>nsi<strong>da</strong><strong>de</strong> eletrônica n(r)é escrita em termos <strong>do</strong>s orbitais <strong>de</strong> Kohn-Sham como:• E note que o potencial efetivo <strong>de</strong> Kohn-Sham é um funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>,V ef = V ef [n]:• Assim, as equações <strong>de</strong> Konh-Sham <strong>de</strong>vem ser resolvi<strong>da</strong>s <strong>de</strong> formaauto-consistente, pois o potencial efetivo <strong>de</strong>pen<strong>de</strong> <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>eletrônica que <strong>de</strong>pen<strong>de</strong> <strong>da</strong>s autofunções <strong>de</strong> Kohn-Sham:


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:<strong>Teoria</strong> Formalmente Exata• A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong> é uma teoria formalmente exata.• Até agora não fizemos nenhuma aproximação sem ter <strong>de</strong>sconta<strong>do</strong>seus erros <strong>de</strong>ntro <strong>do</strong> formalismo.• A per<strong>da</strong> <strong>da</strong> correlação e <strong>da</strong> anti-simetria <strong>da</strong>s funções <strong>de</strong> on<strong>da</strong> <strong>de</strong>vi<strong>do</strong>ao mapeamento por um sistema <strong>de</strong> partículas in<strong>de</strong>pen<strong>de</strong>ntes, éconsi<strong>de</strong>ra<strong>da</strong> no termo <strong>da</strong> energia <strong>de</strong> troca-correlação como umfuncional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>.• A <strong>de</strong>pendência <strong>de</strong>sse termo com a <strong>de</strong>nsi<strong>da</strong><strong>de</strong> não é conheci<strong>do</strong>exatamente até agora.• Por isso <strong>de</strong>vemos trata-lo por meio <strong>de</strong> aproximações.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:Termo <strong>de</strong> troca-correlação• A aproximação mais simples é a chama<strong>da</strong> Aproximação <strong>da</strong> Densi<strong>da</strong><strong>de</strong> Local,LDA, <strong>do</strong> inglês Local Density Approximation.A aproximação LDA consiste em expressar a energia <strong>de</strong> troca-correlação porelétron em um ponto r como aquela em um gás homogêneo que tenha a mesma<strong>de</strong>nsi<strong>da</strong><strong>de</strong> n(r) em r:Essa aproximação funciona bem para sistemas cujas <strong>de</strong>nsi<strong>da</strong><strong>de</strong>seletrônicas não variam rapi<strong>da</strong>mente <strong>de</strong>ntro <strong>de</strong> um região pequena.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:Termo <strong>de</strong> troca-correlação• Para os casos on<strong>de</strong> a <strong>de</strong>nsi<strong>da</strong><strong>de</strong> eletrônica varia no espaço <strong>de</strong> formamenos suave, inclui-se a <strong>de</strong>pendência <strong>da</strong> primeira <strong>de</strong>riva<strong>da</strong> espacial<strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>:Essa é a chama<strong>da</strong> Aproximação <strong>do</strong> Gradiente Generaliza<strong>do</strong>, GGA, <strong>do</strong>inglês Generalized Gradient Approximation. Nela, a energia <strong>de</strong> trocacorrelaçãopor elétron é substituí<strong>da</strong> por uma função local <strong>da</strong><strong>de</strong>nsi<strong>da</strong><strong>de</strong> eletrônica e <strong>do</strong> gradiente <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>.Com o termo <strong>do</strong> gradiente <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong> eletrônica presente nofuncional <strong>de</strong> troca-correlação espera-se que uma melhor <strong>de</strong>scrição<strong>do</strong>s sistemas não homogêneos seja obti<strong>da</strong>.


A teoria <strong>do</strong> funcional <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>• Como teoria a DFT é formalmente exata <strong>de</strong>ntro <strong>do</strong> formalismo <strong>de</strong>Kohn-Sham. Mas na prática <strong>de</strong>vemos fazer aproximações.• O grau <strong>de</strong> precisão <strong>do</strong>s cálculos <strong>de</strong>pen<strong>de</strong>rá <strong>da</strong> quali<strong>da</strong><strong>de</strong> <strong>da</strong>saproximações feitas nos funcionais <strong>da</strong> <strong>de</strong>nsi<strong>da</strong><strong>de</strong>. E sofistica<strong>do</strong>sméto<strong>do</strong>s computacionais nos levam hoje a obter um alto grau <strong>de</strong>confiabili<strong>da</strong><strong>de</strong> <strong>do</strong>s resulta<strong>do</strong>s:Por exemplo:•Comprimentos <strong>de</strong> ligação em moléculas: erro médio menor que0.001 nm.•Constantes <strong>de</strong> re<strong>de</strong> <strong>de</strong> sóli<strong>do</strong>s: erro médio menor que 0.005 nm.•Energia total <strong>de</strong> moléculas: erro médio menor que 0.2 eV. (Paracomparação: uma molécula <strong>de</strong> água tem energia total igual a 2081.1eV)


Referências• P. Hohenberg e W. Kohn, Phys. Rev. 136, B864 (1964)• W. Kohn e L. J. Sham, Phys. Rev. 140, A1133 (1965)• Klaus Capelle, A Bird’s-Eye View of Density-Functional Theory(2006)• Richard M. Martin, Electronic Structure - Basic Theory and PracticalMethods (2004)• J. D. M. Vianna, A. Fazzio e S. Canuto, <strong>Teoria</strong> Quântica <strong>de</strong> Moléculase Sóli<strong>do</strong>s - Simulação Computacional (2004)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!