next up previous
Next: Aproximação de Hückel Up: Algumas Aplicações Previous: Energia Eletrônica de Hartree-Fock

Matriz Densidade de Ordem-N

Introduz-se nesta seção o funcional da densidade de ordem N [18] a partir do produto interno entre os elementos tex2html_wrap_inline2239 tex2html_wrap_inline2295 H e tex2html_wrap_inline2981 tex2html_wrap_inline2295 H tex2html_wrap_inline2987 . Mostra-se, como caso particular desta formulação para N=2, que pode-se obter uma relação entre o índice de ligação (ou de Wiberg) tex2html_wrap_inline2991 para os átomos a e b e uma correlação entre as flutuações das cargas de a e b [19].

Do requerimento de antissimetria de um sistema de partículas resultaram simplificações do formalismo da matriz densidade, como ve-se a seguir. Definindo o orbital de spin tex2html_wrap_inline2873 a partir os geradores de um espaço multilinear alternado e usando suas propriedades antissimétricas, procura-se reinterpretar o conceito de funcional da densidade de ordem N. Assume-se para isso o espaço vetorial de Hilbert H, munido de um produto interno e considera-se a forma 2N-linear tex2html_wrap_inline3009 definida, na aproximação de Hartree-Fock, por

displaymath3011

  equation622

Para orbitais de spin tex2html_wrap_inline2239 mono eletrônicos, a matriz tex2html_wrap_inline3015 tex2html_wrap_inline3017 é representada por,

  equation627

sendo tex2html_wrap_inline3019 um produto interno (veja seção 1). Assim, defini-se a matriz densidade de primeira ordem tex2html_wrap_inline3015 tex2html_wrap_inline3023 , a qual contém a contribuição da parte de spin pela expressão

  equation633

onde tex2html_wrap_inline3025 a contribuição spinorial referente ao orbital tex2html_wrap_inline2873 e tex2html_wrap_inline3029 tex2html_wrap_inline3023 é

  equation637

tex2html_wrap_inline3033 e tex2html_wrap_inline3035 são respectivamente os coeficientes covariantes e contravariantes da expansão LCAO. Como foi mostrado no capítulo , tex2html_wrap_inline3015 tex2html_wrap_inline3039 é uma forma multilinear alternada nas linhas e nas colunas da matriz tex2html_wrap_inline3015 tex2html_wrap_inline3017 . Isto quer dizer que tex2html_wrap_inline3015 é alternada nos tex2html_wrap_inline2239 e tex2html_wrap_inline2981 separadamente, como mostra a relação abaixo

  equation646

Em particular observa-se que, para N=2, a positividade do determinante na equação (0.95), para os geradores tex2html_wrap_inline2239 e tex2html_wrap_inline3055 H se exprime pela seguinte inequação:

  equation657

ou

  equation660

sendo a igualdade válida somente se os geradores tex2html_wrap_inline2239 são linearmente dependentes. Esta equação é conhecida como desigualdade e Cauchy-Schwarz. O número de desigualdades diferentes para N partículas é dado por tex2html_wrap_inline3063

Para um sistema de 2 partículas o determinante da matriz tex2html_wrap_inline3065 (0.95) pode ser escrito em função da matriz tex2html_wrap_inline3067 , da seguinte forma:

displaymath3069

displaymath3071

  equation669

Somando sobre todos os orbitais tex2html_wrap_inline3073 , tex2html_wrap_inline3075 e todos os spins tex2html_wrap_inline3025 e tex2html_wrap_inline3079 , tem-se

displaymath3081

  equation681

ou

displaymath3083

  equation686

onde tex2html_wrap_inline3085 é a carga do átomo a [20] e tex2html_wrap_inline2991 é o índice de ligação[21] entre os átomos a e b,

  equation692

e

  equation701

O lado direito da equação (0.103) expressa a correlação entre as flutuações das cargas dos átomos a e b [19]. Para bases ortogonais, tex2html_wrap_inline2991 se reduz ao índice de Wiberg [22]. Um caso bastante interessante é quando o espaço vetorial em consideração for decomponível, i.e. tex2html_wrap_inline3101 e cada um dos elementos de uma base tex2html_wrap_inline3103 H for ortogonal a cada um dos elementos da base tex2html_wrap_inline3107 H; assim,

displaymath3111

  equation716

No caso de sistemas moleculares, cuja base no espaço vetorial obedece a condição acima, tem-se uma simplificação nos cálculos, já que os termos de interação entre os elementos tex2html_wrap_inline3113 e tex2html_wrap_inline3115 são nulos por construção. Este fato é devido a ortogonalidade entre os dois subespaços.


next up previous
Next: Aproximação de Hückel Up: Algumas Aplicações Previous: Energia Eletrônica de Hartree-Fock

Kleber Mundim
Sun Jul 13 18:14:36 CDT 1997