parte vi teoria de funcional densidade dft
Download
Skip this Video
Download Presentation
Parte VI  Teoria de Funcional Densidade (DFT)

Loading in 2 Seconds...

play fullscreen
1 / 46

Parte VI  Teoria de Funcional Densidade (DFT) - PowerPoint PPT Presentation


  • 76 Views
  • Uploaded on

Parte VI  Teoria de Funcional Densidade (DFT). Joaquim Delphino Da Motta Neto Departamento de Química, Cx. Postal 19081 Centro Politécnico, Universidade Federal do Paraná (UFPR) Curitiba, PR 81531-990, Brasil.

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about ' Parte VI  Teoria de Funcional Densidade (DFT)' - tasya


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
parte vi teoria de funcional densidade dft

Parte VI  Teoria de Funcional Densidade (DFT)

Joaquim Delphino Da Motta Neto

Departamento de Química, Cx. Postal 19081

Centro Politécnico, Universidade Federal do Paraná (UFPR)

Curitiba, PR 81531-990, Brasil

slide2

Já temos uma certa familiaridade com teoria Hartree-Fock e suas correções (CI, MBPT, CC)...

Todos estes métodos são baseados no determinante de Slater — uma complicada função de muitas variáveis...

XXIX ENEQui, Curso 9, Aula #6

resumo
Resumo
  • Racional
  • Modelo de Thomas & Fermi (1928)
  • Teoremas de Hohenberg & Kohn (1965)
  • Equações de Kohn & Sham (1965)
  • Aproximação da densidade local (LDA)
  • Aproximação de gradiente generalizado
  • Funcional de Becke-Lee-Yang-Perdew

XXIX ENEQui, Curso 9, Aula #6

racional
Racional

Em cálculos de estrutura eletrônica, os núcleos são fixos (aproximação de Born & Oppenheimer ), gerando um potencial externo estático V no qual os elétrons se movem. Um estado eletrônico estacionário é então descrito por uma função de onda que satisfaz a equação de Schrödinger multi-eletrônica,

XXIX ENEQui, Curso 9, Aula #6

slide5

A diferença real entre o problema de uma partícula e o complicadíssimo problema de muitas partículas está estreitamente relacionada com o termo de interação U.

Como vimos anteriormente, existem muitos métodos sofisticados para resolver a equação de Schrödinger. O problema com tais métodos é o enorme custo computacional, que torna virtualmente impossível aplicar tais métodos a sistemas mais complexos tais como proteínas.

XXIX ENEQui, Curso 9, Aula #6

slide6

DFT oferece uma alternativa bem atraente, sendo muito mais versátil por que oferece uma maneira de mapear sistematicamente o problema de muitos corpos (com U) em um problema de um corpo (sem U).

A variável essencial é a densidade dada por

XXIX ENEQui, Curso 9, Aula #6

slide7

P. Hohenberg e W. Kohn provaram em 1964 [1] que a relação acima pode ser revertida, ou seja, para uma densidade de estado fundamental n0 é possível calcular a função de onda

correspondente. Em outras palavras, 0 é um funcional unívoco de n0.

XXIX ENEQui, Curso 9, Aula #6

slide8

Funcional da densidade é uma idéia antiga. Os físicos tentaram explicar o gás de elétrons

em termos parecidos com

os que usamos hoje, e muitos

conceitos permanecem...

XXIX ENEQui, Curso 9, Aula #6

llewelyn h thomas 1903 1992
Llewelyn H. Thomas (1903-1992)

Obteve seu Ph.D. em 1927. Ficou famoso pela precessão de Thomas. Entrou para o Watson Lab (Columbia) em 1945. Inventou a memória de “core” em 1946, dois anos antes de An Wang.

XXIX ENEQui, Curso 9, Aula #6

enrico fermi 1901 1954
Enrico Fermi (1901-1954)

Em 1942 conseguiu levar a cabo a primeira reação nuclear controlada. Foi recrutado para o Projeto Manhattan.

XXIX ENEQui, Curso 9, Aula #6

modelo de thomas fermi 1928
Modelo de Thomas & Fermi (1928)

A idéia é usar a teoria do gás uniforme de elétrons não-interagentes, e aplicar este funcional a átomos e moléculas para descrever sólidos.

XXIX ENEQui, Curso 9, Aula #6

slide12

Em 1928 L. Thomas e E. Fermi postularam independentemente que os elétrons estão distribuídos uniformemente no espaço de fase, com dois elétrons para cada elemento h3 de volume.

Assim, para cada elemento d3r podemos encher uma esfera do espaço de momento até o momento de Fermi pf , totalizando

XXIX ENEQui, Curso 9, Aula #6

slide13

Igualando o número de elétrons no espaço de coordenadas ao número de elétrons no espaço de fase, fica

Resolvendo para o momento de Fermi e substituindo, a energia cinética do gás de elétrons se torna um funcional da densidade:

XXIX ENEQui, Curso 9, Aula #6

slide14

Logo, neste método o funcional de energia cinética é

A interação elétron-elétron é dividida em dois termos:

A que se referem estas partes?...

XXIX ENEQui, Curso 9, Aula #6

j repuls o cl ssica
J Repulsão clássica

EX funcional de troca

XXIX ENEQui, Curso 9, Aula #6

slide16

O sucesso inicial do modelo fez com que diversos pesquisadores se ocupassem intensamente do gás de elétrons...

XXIX ENEQui, Curso 9, Aula #6

norman h march
Norman H. March

XXIX ENEQui, Curso 9, Aula #6

paul a m dirac 1902 1984
Paul A.M. Dirac (1902-1984)

Em 1925 publicou seus cinco artigos sobre Mecânica Quântica.

Em 1928 estabeleceu a relação entre Mecânica Quântica e Relatividade com seu famoso artigo do spin ½.

XXIX ENEQui, Curso 9, Aula #6

john c slater 1908 1976
John C. Slater (1908-1976)

Filho de uma família de acadêmicos, recebeu seu doutorado em Harvard em 1923. Foi o chefe do Depto. de Física do MIT desde 1930.

XXIX ENEQui, Curso 9, Aula #6

slide20

Após um hiato de várias décadas, o problema foi revivido com vigor por dois pesquisadores...

XXIX ENEQui, Curso 9, Aula #6

pierre c hohenberg 1934
Pierre C. Hohenberg (1934-...)

Ficou famoso por seu trabalho em DFT em Yale. Atualmente é provost na NY Univ.

XXIX ENEQui, Curso 9, Aula #6

walter kohn 1923
Walter Kohn (1923-...)

Recebeu seu Ph.D. em Harvard em 1948. Nos anos 60 lançou as bases do método DFT. Desde 1991 está na UC em Santa Barbara.

XXIX ENEQui, Curso 9, Aula #6

teoremas de hohenberg kohn 1965
Teoremas de Hohenberg-Kohn (1965)

1) O potencial externo V sentido pelos elétrons é um funcional único da densidade eletrônica.

2) A energia do estado fundamental E[n] é mínima para a densidade exata.

P. Hohenberg and W. Kohn, Phys. Rev. 1964, 136, B864

XXIX ENEQui, Curso 9, Aula #6

equa es de kohn sham 1965
Equações de Kohn-Sham (1965)

O problema variacional de minimizar o funcional de energia E[n] pode ser resolvido aplicando-se o método dos multiplicadores indeterminados de Lagrange. Usa-se o fato de que o funcional na equação pode ser escrito como um funcional de uma densidade hipotética de um sistema não-interagente.

XXIX ENEQui, Curso 9, Aula #6

slide25

Ts denota a energia cinética não-interagente, e Vs é um potencial efetivo externo no qual as partículas se movem. Obviamente

se Vs for escolhido como sendo

XXIX ENEQui, Curso 9, Aula #6

slide26

Logo, pode-se resolver as equações de Kohn-Sham deste sistema não-interagente auxiliar

e determinar os orbitais i que reproduzem a densidade do sistema original.

XXIX ENEQui, Curso 9, Aula #6

slide27

O potencial de uma particular efetivo Vs pode ser escrito em mais detalhe como

onde VXC é o famigerado potencial de “troca e correlação” que inclui todas as interações de muitas partículas.

W. Kohn and L. J. Sham, Phys. Rev. 1965, 140, A1133

XXIX ENEQui, Curso 9, Aula #6

slide28

O maior problema com DFT é que os funcionais exatos de troca e correlação não são conhecidos, exceto para o gás de elétrons.

Entretanto, existem várias aproximações que permitem o cálculo de certas quantidades físicas com bastante acurácia...

XXIX ENEQui, Curso 9, Aula #6

aproxima o da densidade local
Aproximação da densidade local

Em Física a aproximação mais usada é a aproximação da densidade local (local density approximation, LDA), na qual o funcional depende apenas da densidade na coordenada onde o funcional é avaliado.

XXIX ENEQui, Curso 9, Aula #6

slide30

Por sua vez, a aproximação da densidade de spin local ( local spin-density approximation, LSDA) é apenas uma generalização da LDA para incluir o spin eletrônico:

XXIX ENEQui, Curso 9, Aula #6

slide31

Nas últimas décadas, fórmulas altamente acuradas para a energia de troca-correlação

foram construídas a partir

de simulações sobre o modelo do gás de elétrons livres...

XXIX ENEQui, Curso 9, Aula #6

aproxima o gga 1988
Aproximação GGA (1988)

A aproximação de gradiente generalizado (Generalized gradient approximation, GGA) também é local, mas leva em conta o gradiente da densidade na mesma coordenada:

XXIX ENEQui, Curso 9, Aula #6

slide33

Muito bons resultados para geometrias moleculares e energias de estado fundamental são obtidos usando funcionais dentro da GGA.

Vários melhoramentos têm sido alcançados, desenvolvendo melhores representações dos funcionais...

XXIX ENEQui, Curso 9, Aula #6

axel d becke 1953
Axel D. Becke (1953-...)

Em 2000 recebeu a medalha de Schrödinger da WATOC.

Em 2006 foi eleito

fellow da Royal Society of London.

XXIX ENEQui, Curso 9, Aula #6

funcional b3lyp 1993
Funcional B3LYP (1993)

a0 = 0,20 aX = 0,72 aC = 0,81

A. D. Becke, J. Chem. Phys. 98 (1993) 5648

C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 1988, 37, 785

XXIX ENEQui, Curso 9, Aula #6

aplica es
Aplicações

O grande interesse em cálculos DFT vem do fato que ele torna possível alcançar otimizações de geometria para sistemas maiores que 50 átomos.

Todos os programas comerciais de uso geral (Gaussian, GAMESS, C4 etc.) têm a opção de otimizar geometrias com bases Gaussianas de qualidade double- ou superior.

Consequentemente, muitas aplicações em biologia têm usado este protocolo.

XXIX ENEQui, Curso 9, Aula #6

slide37

Exemplos

XXIX ENEQui, Curso 9, Aula #6

problema c ncer
Problema: Câncer
  • Câncer ainda é uma das doenças com maiores índices de mortalidade no planeta.
  • Existem vários tipos diferentes de tumores que atacam diferentes pontos do organismo.
  • Pouco se sabe sobre os mecanismos que disparam o quadro patológico.
  • Em geral protocolos que envolvem quimioterapiasão severos aos pacientes.

Seminário do DQUI, Julho 2002

terapia fotodin mica pdt
Terapia Fotodinâmica (PDT).
  • Protocolos baseados nas propriedades fotofísicas de sistemas orgânicos já foram aprovados e estão em uso nos EUA e Japão.
  • Sensitizadores geralmente são derivados de porfirinas e/ou ftalocianinas.

M. Ochsner, Arzneim.-Forsch. Drug Design 47(II), 1185-1194 (1997).

R. Bonnett, Chem. Soc. Reviews 24, 19-33 (1995).

D. Wöhrle & A. Hirth, Russ.Chem. Bull. 47(5), 807-16 (1998).

Seminário do DQUI, Julho 2002

mecanismos atribu dos
Mecanismos atribuídos:
  • Tipo I: o sensitizador excitado (S*) reage diretamente com um substrato fornecendo radicais que podem dar origem a radicais livres.
  • Tipo II: transferência de energia do sensitizador excitado produzindo 1DgO2 .

M. Ochsner, Arzneim.-Forsch. Drug Design 47(II), 1185-1194 (1997).

R. Bonnett, Chem. Soc. Reviews 24, 19-33 (1995).

Seminário do DQUI, Julho 2002

problemas dos fotosensitizadores conhecidos ex fotofrina
Problemas dos fotosensitizadores conhecidos (ex.: fotofrina)
  • Baixa intensidade de absorção na região de interesse, vermelho (l > 650 nm).
  • Pouca solubilidade em água, e portanto pouca penetrabilidade no tecido alvo.
  • Tendência para acumular na região da aplicação, levando a danos no tecido (toxicidade no escuro).
  • Falta de homogeneidade química.

Seminário do DQUI, Julho 2002

como podemos contribuir
Como podemos contribuir ???
  • Cálculos INDO/S para racionalizar os espectros de absorção/emissão de compostos conhecidos.
  • Estabelecer correlações com dados expe- rimentais (propriedades fotofísicas).
  • Relações estrutura-atividade (QSAR) sobre os dados de atividade biológica (principal- mente necrose dos tumores).

Seminário do DQUI, Julho 2002

slide44

Inicialmente: cálculos semi-empíricos em sistemas conhecidos (TPP, THPP, THPC, THBPC, AlPc, ZnNPc e outros).1) otimização de geometria com AM12) espectro de absorção com INDO/S

Ref.: M.J.S. Dewar & E.G. Zoebisch JACS 107(13), 3902-3909 (1985).

Ref.: J.E. Ridley & M.C. Zerner, Theoret. Chim. Acta 32, 111 (1973).

Seminário do DQUI, Julho 2002

alvos iniciais porfirinas
Alvos iniciais: porfirinas

Seminário do DQUI, Julho 2002

resultados para tpp geometria am1
Resultados para TPP: geometria AM1...

Seminário do DQUI, Julho 2002

ad