Parte I  o Problema da Correlação Eletrônica

1. Parte I  o Problema da Correlação Eletrônica 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

2. Neste curso vamos examinar os diferentes métodos para obter funções de onda acuradas a partir dos conhecidos determinantes de Slater... XXIX ENEQui, Curso 9, Aula #2

3. Resumo • Motivação • O Problema da Correlação • Determinantes excitados • Interação de configurações • Função de onda exata • Método MCSCF, CAS-SCF • Exemplo de aplicação XXIX ENEQui, Curso 9, Aula #2

4. Que é Química Quântica?... XXIX ENEQui, Curso 9, Aula #2

5. Encontrar soluções aproximadas para a equação de Schrödinger eletrônica é a principal ocupação dos quânticos. O primeiro passo é conseguir uma função de ordem zero razoável, o que é usualmente feito usando-se a aproximação Hartree-Fock. XXIX ENEQui, Curso 9, Aula #2

6. A mais simples função de onda anti-simétrica usada para descrever o estado fundamental de sistemas com N elétrons é o determinante de Slater, XXIX ENEQui, Curso 9, Aula #2

7. De acordo com o princípio variacional, a melhor função de onda desta forma funcional é aquela que dá a mais baixa energia possível A flexibilidade variacional está apenas na escolha dos spin-orbitais. XXIX ENEQui, Curso 9, Aula #2

8. Minimizando E0 com respeito à escolha de spin-orbitais, podemos deduzir uma equação (a equação Hartree-Fock) que determina os spin-orbitais ótimos. Em vários textos (ex. Szabo, capítulo 3) mostra-se que esta equação é uma equação do tipo de autovalor com a forma XXIX ENEQui, Curso 9, Aula #2

9. Método SCF Como a equação Hartree-Fock é não-linear, ela deve ser resolvida iterativamente. O procedimento para resolver a equação Hartree-Fock é chamado “método do campo autoconsistente” (SCF). A idéia básica é simples. Fazendo um chute inicial dos spin-orbitais, pode-se calcular o campo médio experimentado por cada elétron e resolver a equação de autovalor para um novo conjunto de spin-orbitais. XXIX ENEQui, Curso 9, Aula #2

10. Usando estes novos spin-orbitais, podemos obter novos campos e repetir o procedimento até que os campos não mudem mais. Os N spin-orbitais com energias mais baixas são chamados “ocupados” ou “buracos”. O determinante de Slater formado por estes spin-orbitais é chamado “função de onda do estado fundamental Hartree-Fock”, e é a melhor aproximação para o estado fundamental do sistema. XXIX ENEQui, Curso 9, Aula #2

11. XXIX ENEQui, Curso 9, Aula #2

12. Que propriedades interessantes podem ser calculadas? XXIX ENEQui, Curso 9, Aula #2

13. Valores esperados de operadores Os autovalores de F’ são as energias orbitais i (que como vimos podem ser associadas com os potenciais de ionização). Como já vimos bastante, o valor esperado da energia é dado por XXIX ENEQui, Curso 9, Aula #2

14. Com a definição do operador de Fock, temos e portanto podemos escrever a energia como XXIX ENEQui, Curso 9, Aula #2

15. Este é um resultado conveniente; se substituímos a expansão nas funções de base, obtemos uma fórmula fechada para a energia que é prontamente calculada de quantidades disponíveis em qualquer estágio do procedimento autoconsistente: XXIX ENEQui, Curso 9, Aula #2

16. Se a energia eletrônica E0 é calculada usando a mesma matriz ℙ que foi usada para formar F, então E0 será um limite superior para a verdadeira energia em qualquer estágio da iteração. Se somamos a repulsão nuclear à energia eletrônica E0 , obtemos a energia total XXIX ENEQui, Curso 9, Aula #2

17. Operadores monoeletrônicos A maior parte das propriedades interessantes (momento de dipolo, momento de quadrupolo, gradiente de campo no núcleo, suscetibilidade magnética etc.) são descritas por somas de operadores monoeletrônicos. XXIX ENEQui, Curso 9, Aula #2

18. Das regras de elementos de matriz que vimos anteriormente, os valores esperados para tais operadores terão a forma de forma que, além da matriz densidade ℙ, precisamos calcular apenas o conjunto de integrais de um elétron para termos os valores esperados dos operadores de um elétron. XXIX ENEQui, Curso 9, Aula #2

19. Momento de dipolo (  ) A definição clássica do momento de dipolo de uma coleção de cargas qi com vetores posição ri é A definição correspondente do operador quântico é XXIX ENEQui, Curso 9, Aula #2

20. O primeiro termo é a contribuição dos elétrons e o segundo termo é a contribuição (clássica) dos núcleos. Portanto, o operador dipolo eletrônico é claramente uma soma de operadores monoele- trônicos. Assim, usando a equação vista anteriormente, Esta é uma equação vetorial com componentes nas direções x, y e z. XXIX ENEQui, Curso 9, Aula #2

21. Análise de população A densidade de carga representa a probabilidade de encontrar um elétron em várias regiões do espaço, e é usualmente representada por mapas de contorno de planos cortando a molécula. XXIX ENEQui, Curso 9, Aula #2

22. Não há uma definição única do número de elétrons associados com um átomo ou núcleo dentro da molécula. E agora ?... XXIX ENEQui, Curso 9, Aula #2

23. Robert S. Mulliken (1896-1986) Em 1921, recebeu o Ph.D. na Univ. de Chicago por seu trabalho com W.D. Harkens na separação de isótopos. Em 1922 foi para Harvard. Em 1927 foi para a Europa e discutiu a teoria de Hund. Logo conseguiu associar números quânticos aos elétrons em moléculas. XXIX ENEQui, Curso 9, Aula #2

24. População de Mulliken Ás vezes é útil fazer uma “análise de população”. Observe que divide o número total de elétrons em dois elétrons por orbital molecular. XXIX ENEQui, Curso 9, Aula #2

25. Substituição da expansão na base fornece e assim é possível interpretar o elemento diagonal ( ℙ S ) como o número de elétrons “ocupando” o -ésimo orbital de base. XXIX ENEQui, Curso 9, Aula #2

26. Para que serve isso?... XXIX ENEQui, Curso 9, Aula #2

27. Gilbert Newton Lewis (1875-1946) Obteve seu Ph.D. em Harvard em 1899 trabalhando com amálgamas. Em 1908 publicou seu primeiro artigo sobre relatividade, onde deduziu a relação massa-energia de uma forma diferente da de Einstein. Em 1916 formulou a idéia de que uma ligação covalente consiste de um par de elétrons compartilhado entre átomos (o famoso artigo sobre a regra do octeto). Em 1919 realizou o famoso experimento em que descobriu o O4. Em 1923 lançou sua teoria de ácidos e bases. XXIX ENEQui, Curso 9, Aula #2

28. Lembrete: notação de Lewis A molécula é desenhada com os átomos em suas posições aproximadas, e os elétrons de valência (pertencentes à camada mais externa do átomo) são representados como pontos ao redor do átomo. Regra do octeto (1916) Diferentes átomos compartilham elétrons sempre no sentido de completar oito elétrons em sua camada mais externa. O hidrogênio “quer” completar dois elétrons. XXIX ENEQui, Curso 9, Aula #2

29. A regra é muito bem-sucedida em explicar qualitativamente a estrutura e as propriedades de muitos compostos químicos diferentes. XXIX ENEQui, Curso 9, Aula #2

30. Às vezes pode acontecer de, para um determinado arranjo espacial dos átomos, a estrutura eletrônica poder ser descrita por mais de uma estrutura de Lewis válida... XXIX ENEQui, Curso 9, Aula #2

31. Neste caso, aparece a “ressonância” entre mais de uma estrutura de Lewis. Por exemplo, veja o caso clássico do ozônio: A teoria VB mostra claramente por que este composto (que é um gás) deve ser angular e não linear. XXIX ENEQui, Curso 9, Aula #2

32. Químicos gostariam de manter o quadro fácil da teoria de valência de Lewis, logo uma população que suporte esta visão é bem vinda... XXIX ENEQui, Curso 9, Aula #2

33. XXIX ENEQui, Curso 9, Aula #2

34. XXIX ENEQui, Curso 9, Aula #2

35. Logicamente, essa definição não é única... XXIX ENEQui, Curso 9, Aula #2

36. Per-Olov Löwdin (1916-2000) Participou dos maiores desenvolvimentos da Mecânica Quântica. Estabeleceu diversos conceitos usados até hoje. Fundou o QTP e foi editor do IJQC por anos. XXIX ENEQui, Curso 9, Aula #2

37. População de Löwdin A definição não é única. Como tr AB = tr BA, podemos escrever para qualquer . XXIX ENEQui, Curso 9, Aula #2

38. Particularmente, com  = 1/2, temos XXIX ENEQui, Curso 9, Aula #2

39. Podemos mostrar que ℙ’ é a matriz densidade em termos de uma base simetricamente ortogonalizada: Os elementos diagonais de ℙ’ são usados para uma análise de população de Löwdin. XXIX ENEQui, Curso 9, Aula #2

40. Ok. Só nos resta então escolher uma base. Que tipos de bases poliatômicas existem? XXIX ENEQui, Curso 9, Aula #2

41. Bases poliatômicas Várias sugestões de bases já foram colocadas para cálculos moleculares, mas apenas duas são de uso comum. A função de Slater 1s norma- lizada tem a forma XXIX ENEQui, Curso 9, Aula #2

42. Outro tipo de base é a função do tipo Gaussiano centrada em RA onde  é o “expoente orbital Gaussiano” caracterizando a função de base. Este tipo de base é mais usado em cálculos ab initio. XXIX ENEQui, Curso 9, Aula #2

43. Diferenças entre STOs e GTOs As maiores diferenças entre estas funções ocorrem em r = 0 e r muito grande. Por inspeção apenas verifica-se que em r = 0 (a chamada “condição de cusp”) as funções de Slater têm uma derivada finita e as funções Gaussianas têm derivada nula. Em grandes valores de r, a função Gaussiana decai muito mais rapidamente que a função de Slater. XXIX ENEQui, Curso 9, Aula #2

44. STOs vs. GTOs XXIX ENEQui, Curso 9, Aula #2

45. Para cálculos de estrutura eletrônica molecular, é preferível usar uma base de funções de Slater pois elas se parecem mais com os orbitais hidrogenóides. Por exemplo, é possível mostrar que em grandes valores de r os orbitais moleculares decaem como e-ar e não e-ar2 . XXIX ENEQui, Curso 9, Aula #2

46. Orbitais hidrogenóides Por exemplo, no estado fundamental ( 2S ) do hidrogênio temos um elétron ocupando um orbital 1s, o qual é descrito pela função esfericamente simétrica: XXIX ENEQui, Curso 9, Aula #2

47. XXIX ENEQui, Curso 9, Aula #2

48. XXIX ENEQui, Curso 9, Aula #2

49. XXIX ENEQui, Curso 9, Aula #2

50. XXIX ENEQui, Curso 9, Aula #2