Parte vi teoria de funcional densidade dft
Download
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


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 suas correções (CI, MBPT, CC)...

  • 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 suas correções (CI, MBPT, CC)...

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


A suas correções (CI, MBPT, CC)...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


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


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


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) tentaram explicar o gás de elétrons

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) tentaram explicar o gás de elétrons

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) tentaram explicar o gás de elétrons

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


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


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


Logo, neste método o funcional de energia cinética é ao número de elétrons no espaço de fase, fica

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 ao número de elétrons no espaço de fase, fica  Repulsão clássica

EX funcional de troca

XXIX ENEQui, Curso 9, Aula #6


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 pesquisadores se ocupassem intensamente do gás de elétrons...

XXIX ENEQui, Curso 9, Aula #6


Paul a m dirac 1902 1984
Paul A.M. Dirac (1902-1984) pesquisadores se ocupassem intensamente do gás de elétrons...

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) pesquisadores se ocupassem intensamente do gás de elétrons...

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


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-...) com vigor por dois pesquisadores...

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-...) com vigor por dois pesquisadores...

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) com vigor por dois pesquisadores...

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) com vigor por dois pesquisadores...

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


Ts com vigor por dois pesquisadores...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


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


O potencial de uma particular efetivo sistema não-interagente auxiliarVs 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


O maior problema com DFT é que os funcionais sistema não-interagente auxiliarexatos 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 sistema não-interagente auxiliar

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


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

XXIX ENEQui, Curso 9, Aula #6


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) energia de troca-correlação

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


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-...) de estado fundamental são obtidos usando funcionais dentro da GGA.

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) de estado fundamental são obtidos usando funcionais dentro da GGA.

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 de estado fundamental são obtidos usando funcionais dentro da GGA.

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


Exemplos de estado fundamental são obtidos usando funcionais dentro da GGA.

XXIX ENEQui, Curso 9, Aula #6


Problema c ncer
Problema: Câncer de estado fundamental são obtidos usando funcionais dentro da GGA.

  • 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). de estado fundamental são obtidos usando funcionais dentro da GGA.

  • 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: de estado fundamental são obtidos usando funcionais dentro da GGA.

  • 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


Seminário do DQUI, Julho 2002 de estado fundamental são obtidos usando funcionais dentro da GGA.


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


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 conhecidos (TPP, THPP, THPC, THBPC, AlPc, ZnNPc e outros).

Seminário do DQUI, Julho 2002


Resultados para tpp geometria am1
Resultados para TPP: geometria AM1... conhecidos (TPP, THPP, THPC, THBPC, AlPc, ZnNPc e outros).

Seminário do DQUI, Julho 2002


ad