Índices de reatividade química a partir da teoria do funcional de densidade: formalismo e perspectivas

Chemical reactivity indexes from density functional theory: formalism and perspectives

density functional theory; reactivity indexes; hardness; softness

Revisão

ÍNDICES DE REATIVIDADE QUÍMICA A PARTIR DA TEORIA DO FUNCIONAL DE DENSIDADE: FORMALISMO E PERSPECTIVAS

Hélio Anderson Duarte*

Departamento de Química - ICEx, Universidade Federal de Minas Gerais, Av. Antônio Carlos, 6627, Pampulha, 31270-901 Belo Horizonte - MG

* e-mail: [e-mail]

Recebido em 30/5/00; aceito em 6/11/00

CHEMICAL REACTIVITY INDEXES FROM DENSITY FUNCTIONAL THEORY: FORMALISM AND PERSPECTIVES. The fundaments of the modern Density Functional Theory (DFT), its basic theorems, principles and methodology are presented. This review also discuss important and widely used concepts in chemistry but that had not been precisely defined until the development of the DFT. These concepts were proposed and used from an empirical base, but now their precise definition are well established in the DFT formalism. Concepts such as chemical potential (electronegativity), hardness, softness and Fukui function are presented and their consequences to the understanding of chemical reactivity are discussed.

Keywords: density functional theory; reactivity indexes; hardness; softness.

INTRODUÇÃO

A Teoria de Drude¹ de metais foi provavelmente a primeira tentativa de se usar a densidade eletrônica como variável básica na descrição de um sistema eletrônico. Isto aconteceu apenas três anos depois da descoberta do elétron por Thompson em 1897. Drude aplicou a teoria cinética dos gases a um metal, considerado como um gás homogêneo de elétrons, para desenvolver a sua teoria sob condução térmica e elétrica. Da mesma forma que na teoria cinética dos gases, as interações elétron-elétron e elétron-íon eram desprezadas. Mais tarde, Sommerfeld usou a distribuição quântica de Fermi-Dirac no lugar da distribuição clássica de Maxwell-Boltzman para a distribuição da velocidade eletrônica no modelo de Drude. O modelo de Drude-Sommerfeld continha várias simplificações e, consequentemente, não foi possível fazer uma descrição precisa de um sólido a partir deste modelo, embora ele tenha sido considerado como um grande avanço naquela época.

O próximo passo foi o modelo de Thomas-Fermi², proporcionando uma sensível melhora em relação ao modelo de Drude-Sommerfeld. Neste modelo, partiu-se do pressuposto de que os elétrons estão distribuídos uniformemente no espaço. Usando argumentos estatísticos para aproximar a distribuição dos elétrons, Thomas e Fermi desenvolveram o funcional de energia de Thomas-Fermi (TF) (Eq. 1 com Cx = 0). O modelo de Thomas-Fermi-Dirac (TFD) inclui, em adição, a energia de troca para um gás de elétrons desenvolvido por Dirac³. O funcional de energia (E) de TFD é dado por

onde e .

Os quatro termos da direita da Eq. (1) correspondem à energia cinética, ao potencial externoⁱ, ao potencial de Coulomb e à energia de troca, respectivamente. r e r correspondem a densidade eletrônica e as coordenadas, respectivamente.

O modelo TFD pode ser resolvido apenas numericamente e a solução da equação variacional

onde m é o potencial químico e N o número de partículas, é essencialmente independente do número de partículas. Porém, estes modelos de funcionais clássicos (TF e TFD) não podem competir em precisão com métodos auto-consistentes como, por exemplo, o método de Hartree-Fock. O modelo de TFD falha ao reproduzir a estrutura quântica de camadas dos átomos. Além disso, o método não prediz ligação molecular e os resultados para átomos não são muito precisos. A TFD é supersimplificada e não permite predições quantitativas^4,5.

Foi em 1964, com a publicação de dois teoremas por Hohenberg e Kohn , que o uso da densidade eletrônica como variável básica foi rigorosamente legitimado. Eles forneceram os fundamentos da teoria do funcional de densidade moderna (DFT, do Inglês Density Functional Theory) e mostraram que os modelos baseados no funcional de energia TFD devem ser vistos como uma aproximação de uma teoria exata. Os dois teoremas⁶ mostram que existe um funcional de energia exato da densidade eletrônica E[r] e um princípio variacional exato para este funcional, semelhante a Eq. (2). E, em 1965, Kohn e Sham⁷ propuseram uma forma de contornar o problema de se encontrar o funcional de energia exato (o método KS). Desde a publicação destes dois artigos, a DFT tem atraído cada vez mais a atenção da comunidade científica e, atualmente, ela é largamente usada para se estudar sistemas cada vez mais complexos.

Interessante e cuidadosa comparação entre os métodos computacionais usados para fazer cálculos ab initio convencionais e DFT é apresentada no artigo de Morgon e Custódio⁸ publicado neste periódico. O desenvolvimento da metodologia computacional para fazer cálculos DFT leva, invariavelmente, as equações matemáticas semelhantes às equações do método de Hartree-Fock-Roothaan. Porém, apesar da semelhança, os resultados de cálculos DFT devem ser interpretados à luz da Teoria do Funcional de Densidade. Diferentemente de outros métodos, tais como métodos semi-empíricos, que tentam se aproximar do método Hartre-Fock tanto quanto possível, a DFT está relacionada com a solução exata do problema de muitos elétrons.

Esta revisão apresenta os fundamentos da teoria moderna do funcional de densidade, seus teoremas básicos, princípios e metodologia. Conceitos importantes na descrição teórica de reações químicas, tais como potencial químico, conceitos de maciez e dureza (hardness e softness) e outros índices de reatividade química redefinidos, de uma maneira exata a partir da Teoria do Funcional de Densidade são discutidos.

TEOREMAS DE HOHEMBERG-KOHN

Primeiro Teorema

No artigo de 1964, Hohenberg e Kohn⁶ (HK) consideraram um sistema de N-elétrons descrito pelo Hamiltoniano não relativístico:

onde é a energia cinética, o operador de repulsão elétron-elétron o qual inclui a repulsão Coulombiana e todos os termos não clássicos (troca e correlação) e é o potencial externo com respeito aos elétrons, normalmente devido às cargas dos núcleos Za:

A densidade eletrônica é definida pela equação seguinte:

onde Y(r₁,r₂,...r_n) é a solução do estado fundamental do Hamiltoniano.

Originalmente, supos-se que o estado fundamental não é degenerado. A energia total do sistema é então dada por

ou, se o potencial externo é separado em um funcional trivial da densidade eletrônica, por

e aplica-se universalmente à todos os sistemas eletrônicos. O número de elétrons, N, e o potencial externo, , no qual os elétrons se movem, definem completamente o sistema de muitos elétrons. O primeiro teorema de HK estabelece que o potencial externo, , é um funcional único de r(r), além de uma constante aditiva.

Segundo Teorema

O segundo teorema estabelece que qualquer aproximação da densidade eletrônica, , de tal modo que e , a energia total .

A forma do operador de energia cinética e do operador da repulsão eletrônica é a mesma para qualquer sistema eletrônico com um dado número de partículas e potencial externo. Pode-se, então, definir um funcional universal,

No entanto, é preciso observar que define seu próprio e, consequentemente, o Hamiltoniano e . A função , por sua vez, pode ser usada como uma função tentativa para o sistema com o potencial externo u(r). De acordo com o princípio variacional, temos

E_u se refere ao funcional da energia total com relação ao potencial externo u(r).

Os dois teoremas de HK mostram como se pode determinar o estado fundamental de um sistema, com um dado potencial externo, usando a densidade eletrônica tridimensional como variável básica ao invés de usar a função de onda de N-elétrons, que é muito mais complexa.

Na prática, nós queremos usar o princípio variacional para encontrar r(r) em um dado potencial externo u(r), e não o oposto. O primeiro teorema de HK estabelece que se r(r) é conhecido, pode-se determinar u(r). Porém, dado uma densidade eletrônica tentativa, , como podemos garantir que é "n-representável", ou seja, que corresponda à um potencial fisicamente realizável? Uma outra questão surge em relação a "N-representabilidade" da densidade eletrônica. A densidade eletrônica é "N-representável" se ela pode ser obtida a partir de uma função de onda antissimétrica. Gilbert⁹ mostrou que qualquer densidade eletrônica razoável é "N-representável", se as seguintes condições forem satisfeitas:

Felizmente, a Teoria do Funcional de Densidade pode ser formulada de tal forma que a condição da "n-representabilidade" pode ser contornada através da formulação de Levy¹⁰ da procura restrita (Levy Constrained Search). Seguindo Parr e Yang⁴,

ou seja, a procura restrita consiste em testar todas as densidades tentativas, e para cada densidade encontrar as funções de onda Y(r₁,r₂,...r_n) que minimizem F[r] e que gerem aquela densidade. Substituindo a Eq. (11) na Eq. (9) temos

A minimização interna está restrita a todas as funções de onda, Y, que fornece r, enquanto a minimização externa é feita sobre todos os possíveis r.

A formulação de Levy também elimina a limitação dos teoremas de HK, segundo os quais, o estado fundamental tem que ser não degenerado. Além do mais, nenhuma referência a "n-representabilidade" da densidade eletrônica é necessária, uma vez que r(r) é originada a partir de uma função de onda antissimétrica. A questão da "N-representabilidade" da densidade eletrônica é uma condição que tem que ser satisfeita, porém é uma condição mais fácil de ser atendida como foi demonstrada por Gilbert⁹, qualquer densidade eletrônica razoável é N-representável se as Eqs (10a), (10b) e (10c) forem satisfeitas.

EQUAÇÕES DE KOHN-SHAM

De acordo com os teoremas de HK, o funcional energia é escrito como segue:

O estado fundamental satisfaz o princípio estacionário como requerido pelo princípio variacional (Eq. (9)),

O multiplicador de Lagrange, m, é o potencial químico dos elétrons e será discutido mais tarde. A função universal F(r) é determinada pelo procedimento de minimização da Eq. (11).

Uma importante idéia introduzida por Kohn e Sham⁷ é usar o conceito de um sistema de referência de partículas independentes, ou seja, que não interagem. Eles rescreveram a Eq. (13) tornando explícita a repulsão elétron-elétron de Coulomb, e definindo uma nova função universal G(r),

onde,

e T_s[r] é o funcional de energia cinética de um sistema de elétrons que não interagem tendo a mesma densidade eletrônica do sistema que interage. E_xc inclui o termo de interação elétron-elétron não clássica (troca e correlação) mais a parte residual da energia cinética, T[r]-T_s[r], onde T[r] é a energia cinética exata para o sistema de elétrons que interagem.

Torna-se importante ressaltar que as diferenças entre as definições da correlação eletrônica tradicionais (métodos ab initio convencional) e da DFT são geralmente ignoradas e funcionais de correlação são usualmente julgados pelo seu desempenho em reproduzir a energia de correlação tradicional, o qual é geralmente o único conhecido com um razoável grau de precisão¹¹. Embora o formalismo DFT seja exato, a expressão para o funcional de troca-correlação não é conhecido em sua forma exata.

Orbital de um elétron pode, então, ser utilizado permitindo calcular precisamente a energia cinética via um procedimento auto-consistente:

Os orbitais são auto-funções do Hamiltoniano KS efetivo de um elétron,

onde, o potencial efetivo, é definido como segue,

onde,

A densidade eletrônica, r(r) , é dada por

Portanto, os orbitais Kohn-Sham (KS), y_i, são obtidos a partir da equação de Schrödinger de um elétron

As Equações (19), (20) e (21) representam o esquema Kohn-Sham auto-consistente (KS-SCF). Uma vez que o potencial efetivo, u_eff(r) depende da densidade eletrônica, r(r) , as equações de Kohn-Sham são resolvidas através de um procedimento auto-consistente.

No esquema Kohn-Sham, a densidade exata do estado fundamental de um sistema de elétrons que interagem é gerada a partir da solução de um problema auxiliar do sistema de elétrons que não-interagem definido de forma a fornecer a mesma densidade do estado fundamental. Este aspecto do método KS é discutido em detalhe nas referências 4 e 5. As equações KS, de maneira semelhante às equações de Hartree-Fock, fornecem equações de um elétron que descrevem sistemas de muitos elétrons. Além disso, as equações Kohn-Sham, em princípio, são exatas uma vez que incorporam totalmente os efeitos da correlação eletrônica (troca e correlação) e solução dessas equações formalmente equivalem a resolver exatamente o problema variacional da DFT. Em outras palavras, as equações KS vão além da aproximação de Hartree-Fock e ainda, oferecem uma grande vantagem: o potencial de troca-correlação é local. Na prática, isto significa que as integrais de quatro centros do método de Hartree-Fock usadas para calcular o termo de troca (que é um operador não local) não precisam ser realizadas, reduzindo-se desta forma o custo computacional.

A energia total é determinada usando a Eq. (15) ou pela fórmula

O esquema KS permite calcular a densidade eletrônica do estado fundamental. Todas as outras propriedades do sistema podem ser calculadas desde que os funcionais da densidade eletrônica apropriados sejam conhecidos.

A parte mais difícil para resolver as equações KS é determinar o funcional de troca-correlação (XC)u_xc(r) , uma vez que a exata forma analítica do potencial de troca-correlação não é conhecida. Muitos trabalhos têm sido devotados ao desenvolvimento e teste das diferentes aproximações dos funcionais de troca-correlação^12,13. A aproximação da densidade local baseada no modelo de gás homogêneo de elétrons^14,15,a aproximação do gradiente generalizado^16-21, o funcional LAP (Laplaciano), que usa a densidade de energia cinética e o Laplaciano da densidade eletrônica, são alguns dos funcionais de XC mais usados no cálculo de estrutura eletrônica. Os funcionais híbridos que incluem parte da energia de troca exata (termo de troca do método de Hartree-Fock) têm sido largamente utilizados principalmente no estudo de sistemas orgânicos. No entanto, vale ressaltar que estes funcionais não são puros uma vez que o potencial efetivo u_eff(r) deixa de ser um operador local.

O SIGNIFICADO DAS ENERGIAS DOS ORBITAIS KS E A EXTENSÃO DA DFT AO NÚMERO DE OCUPAÇÃO NÃO-INTEIRO

O teorema de Janak²² e a conseqüente extensão da DFT ao número de ocupação não-inteiro é apresentada, seguindo a formulação de Mermin²³ e Perdew et al.²⁴.

Primeiramente, a definição de T_s(r) é generalizada, onde todos os possíveis números de ocupação n_i e os orbitais ortonormais y_i são procurados:

o subscrito J indica Janak. A densidade r(r) é, agora, definido como segue

Portanto, é possível expressar a energia total em termos de n_i e y_i,

com

Seguindo Janak²², usando a regra da cadeia para derivadas de funcional, a expressão da energia total pode ser diferenciada com respeito ao número de ocupação n_i para fornecer,

Esta equação estabelece uma regra aufbau, o que significa que o funcional de energia, Eq. (26), é minimizado quando se coloca n_i = 1 para os N orbitais y_i mais baixos em energia e n_i = 0 para os outros orbitais. Ocupação fracionária dos orbitais é possível apenas dentro de um conjunto de orbitais degenerados ao nível de Fermi. Porém, para sistemas de elétrons que não interagem entre si e que são n-representáveis, esta regra não se aplica²⁵. A Eq. 28 também se aplica para sistemas com número não inteiro de elétrons, ou seja, número fracionário . De acordo com Perdew et al.²⁴, a Eq. (28) fornece um significado físico do orbital KS ocupado de maior energia,e_max, ou seja,

onde I, A e N_Z são o potencial de ionização, a afinidade eletrônica e a soma dos números atômicos dos átomos no sistema com N elétrons, respectivamente.

A extensão da DFT para sistemas com um número de elétrons fracionário derivado por Mermin²³ e Janak²² e sua conseqüência é discutida em detalhe por Russier²⁶ e Perdew²⁷. Recentemente, Duarte e Salahub²⁸ mostraram em detalhes estas extensões e como elas podem ajudar a aperfeiçoar o modelo de cluster para estudar adsorção de moléculas em superfícies metálicas.

POTENCIAL QUÍMICO E A DEFINIÇÃO DE ELETRONEGATIVIDADE

É uma vantagem muito grande o número de elétrons N e a densidade eletrônica desempenharem um papel central na Teoria do Funcional de Densidade. É muito mais difícil relacionar funções de onda de muitos elétrons com a estrutura e reatividade molecular que quantidades observáveis mais simples tais como r(r) e o número de elétrons N. A DFT permite definir também importantes conceitos químicos relacionados com a estrutura e a reatividade química de moléculas. São apresentados aqui os conceitos mais importantes que resultam da DFT e que constituem área de intensa atividade de pesquisa^4,29,30.

Na seção 3, ao estabelecer-se o princípio estacionário na DFT (Eq. (14)), o multiplicador de Lagrange m foi relacionado ao potencial químico,

Portanto, m=m(N,u) mede a tendência dos elétrons de escaparem de um sistema em equilíbrio³¹. Se aproximar a definição do potencial químico mostrado na equação acima pela diferença finita (veja Figura 1), chega-se a

onde I e A correspondem ao Potencial de Ionização vertical e à Afinidade Eletrônica vertical, para que seja levado em conta o fato de que o potencial externo é mantido constante (veja Eq. (30)). A Figura 1 mostra graficamente a definição de I e A, e como estas grandezas podem ser usadas para estimar a primeira e segunda derivada da energia em relação ao número de elétrons.

Esta é a mesma definição de eletronegatividade (c_M) introduzida por Mulliken^32,33, em 1934. O conceito de eletronegatividade é importante na Teoria do Funcional de Densidade, uma vez que ele é o potencial químico e, por isso, desempenha um papel importante no princípio variacional, Eq. (14). Esta expressão foi identificada primeiramente por Parr et al.³⁴.

Sanderson propôs em 1954³⁵ o princípio da igualação da eletronegatividade: "quando dois ou mais átomos, inicialmente possuindo eletronegatividades diferentes, ligam-se para formar um composto, suas eletronegatividades tornam-se ajustadas a um valor intermediário dentro do composto"³⁶, ou seja, o potencial químico em um sistema em equilíbrio é igual em toda a parte. Quando duas moléculas se aproximam, o potencial químico deve ser igualado para os dois sistemas em equilíbrio. O potencial químico (eletronegatividade) é uma propriedade de estado do sistema e pode ser calculada via DFT, Eq.(31). A transferência de elétrons entre duas moléculas é dirigida pelas diferenças do potencial químico. O formalismo DFT leva rigorosamente ao conceito de eletronegatividade e ao princípio da igualação da eletronegatividade, mesmo embora todos estes conceitos tenham sido formulados sobre uma base empírica no passado.

Vejamos o exemplo da formação do composto NaF a partir dos seus elementos, no estado neutro em fase gasosa (Na e F). A Figura 2 mostra a formação deste composto exemplificando o princípio da igualação da eletronegatividade e o papel do potencial químico na formação da ligação química. A Figura 2a mostra a variação da energia total de cada fragmento em relação a variação da carga eletrônica e a Figura 2b corresponde a 1a. derivada da energia total com relação a carga. Inicialmente, os dois elementos estão infinitamente separados, cada qual com a carga q igual a zero. O Flúor apresenta potencial químico em torno de 9,3 e o Sódio em torno de 3,1. Em seguida, carga eletrônica flui do flúor (com maior m) para o sódio (com menor m). Em torno de q = 0,35 (m = 4,8) os potenciais químicos são igualados e o fluxo de carga é interrompido. A definição do potencial químico em cada fragmento desta molécula é dado pela Eq. (32). Igualado o potencial químico, o sistema se encontra em equilíbrio e o composto NaF é formado. Vale notar que pela análise de Mulliken, o cálculo DFT do composto NaF (distância de ligação calculada para 1,96 Å) prevê carga de 0,30 sobre o átomo de sódio e -0,30 sobre o átomo de flúor.

Obviamente, esta é uma visão simplificada da formação da ligação química. Alguém poderia se perguntar como dois átomos de oxigênio (que apresentam o mesmo potencial químico) podem reagir para formar a molécula de oxigênio. A explicação é dada pelo princípio da dureza máxima sugerida por Pearson em 1987³⁷ e que será discutida na próxima seção.

O POTENCIAL QUÍMICO E SUAS DERIVADAS: DUREZA, MACIEZ E FUNÇÃO DE FUKUI

A seguinte equação diferencial descreve um sistema eletrônico quando este muda de um estado fundamental para o outro:

deduzida a partir da generalização da teoria de perturbação de 1ª ordem para incluir mudança no número total de elétrons do sistema⁴.

A energia total é um funcional do número de elétrons N e do potencial externo u(r) , ou seja, E = E[N,u]. Logo, pode-se escrever a equação diferencial total da energia E que relaciona a mudança do sistema de um estado fundamental para um outro:

Lembrando que a partir da densidade eletrônica r(r) determina-se o número de elétrons N, podemos também escrever a expressão da energia total como um funcional da densidade eletrônica: E=E[r,u]. Logo a Eq. (34) pode ser rescrita como segueⁱⁱ:

Porém, torna-se necessário mencionar que embora a Eq. (34) seja exata, a Eq. (35) é apenas uma aproximação, pois estamos assumindo que r(r) e u(r) são independentes. No entanto, os teoremas de Hohenberg-Kohn demonstram que u(r) completamente determina a densidade eletrônica do estado fundamental.

Comparando as Eqs. (33) e (34), obtem-se as seguintes igualdades:

O potencial químico m, por sua vez, é um funcional do número de elétrons N e do potencial externo u(r) , isto é, m=m[N,u]. Diferenciando esta expressão temos:

A Eq. (38) poderia ser derivada levando-se em conta que m é um funcional de r(r) e, consequentemente, uma equação mais complexa seria obtida.

Importantes conceitos químicos surgem a partir das Eqs. (33),(34),(36) e (37):

Dureza (Hardness)³⁸

Se utilizarmos a aproximação da diferença finita (veja Figura 1), chegamos a equação:

Esta definição é a mesma que a sugerida por Pearson em 1963³⁹.

Recentemente, vários trabalhos foram publicados mostrando a aplicação destes conceitos para a compreensão do processo de isomerização em compostos como peróxido de hidrogênio, tioperóxido de hidrogênio e persulfeto de hidrogênio. Toro-Labbé et al.⁴⁰ mostraram como o potencial químico e a dureza variam com a ângulo de torção da molécula. Embora eles tenham usado em seus cálculos a aproximação da diferença finita, ou seja, h=1/2(E_L-E_H), onde E_L e E_H correspondem a energia do LUMO e HOMO, respectivamente, os resultados se comportam como o esperado. A dureza apresenta um máximo onde se encontra os mínimos de energia na superfície de energia potencial. O princípio da dureza máxima foi também exemplificada no estudo do rearranjo intramolecular de uma série de moléculas tais como HCN, HClO, HONS, H₂SO, F₂S₂, H₃PO, CH₂SH₂⁴¹. Em todos estes sistemas a dureza é máxima nos pontos correspondentes aos mínimos na superfície de energia potencial.

Maciez (Softness), que é o inverso da Dureza⁴²:

Estes conceitos propostos inicialmente sobre uma base empírica estão, a partir do formalismo DFT, precisamente definidos.

O princípio enunciado por Pearson³⁹: "ácidos macios preferem bases macias, ácidos duros preferem bases duras", também chamado de princípio HSAB (do inglês: Hard Soft Acid Base principle) também está contemplado dentro da DFT. Porém não é uma tarefa fácil estabelecer teoricamente este princípio. Recentemente, Parr e co-autores⁴³, forneceram duas demonstrações deste princípio com a restrição de que os dois reagentes tenham a mesma eletronegatividade, ao invés do potencial externo constante. A caracterização do tautomerismo keto-enol dos derivados do acetil foram analisados recentemente em termos do potencial químico e dureza⁴⁴. O princípio HSAB foi invocado mostrando o bom acordo com os dados termodinâmicos. De acordo com o princípio HSAB, a forma keto é favorecida em relação a forma enol.

Índice de Reatividade de Fukui (ou Função de Fukui) :

Este índice de reatividade mede como uma perturbação externa afeta o potencial químico de um sistema. A função de Fukui, f(r), envolve a densidade eletrônica do átomo ou molécula na sua região de valência. Na realidade, a função de Fukui pode ser usada para medir a reatividade em relação a um ataque nucleofílico, eletrofílico ou mesmo frente a um reagente radicalar⁴⁵. Três definições numéricas foram propostas para as funções de Fukui por Yang e Parr⁴⁶. Estas três definições são devidas ao fato de que a Eq. 42 apresenta descontinuidade para N igual a um número inteiro.

Em outras palavras, se usarmos a aproximação da diferença finita nas definições apresentadas nas Eqs 43 e 44, tem-se que:

A terceira função, f⁰(r), é estimada como uma média das funções f⁺(r) e f^-(r). Em contraste com os conceitos de dureza e maciez que são quantidades globais, a função de Fukui é uma quantidade local e reflete as propriedades dos diferentes sítios dentro de uma molécula. Usando da aproximação dada pelas Eqs 46 e 47, a Figura 3 mostra os orbitais HOMO e LUMO do cloreto de t-butila indicando onde um ataque nucleofílico ou eletrofílico pode acontecer.

A energia total de um sistema pode ser expandida em termos destas grandezas, como se segue abaixo⁴⁷:

Outras grandezas importantes para a química e física surgem a partir da DFT tais como os conceitos de maciez local (Eq. 49) e sua inversa, chamada de dureza local.

A Eq. 49 mostra que a maciez pode ser determinada em diferentes regiões da molécula e, deste modo, avaliando a reatividade nos diferentes sítios de um sistema molecular.

Os conceitos relacionados ao potencial químico e suas derivadas foram discutidas em detalhes na caracterização de reações químicas por Toro-Labbé⁴⁸. Estes conceitos foram aplicados no estudo da transferência de prótons em complexos do tipo CHX-XH ... CHX-XH (X=O, S)⁴⁹ e para estimar energias de ligação de complexos ligados por ligações de hidrogênio⁵⁰. O sucesso na aplicação dos conceitos de potencial químico, dureza e maciez para se caracterizar reações químicas é notório, apesar de na grande maioria dos trabalhos utilizarem aproximações da diferença finita para avaliar estas propriedades. Em relação as propriedades locais tais como a dureza e maciez locais além da função de Fukui ainda são muito incipientes, embora sejam de óbvio interesse.

CONSIDERAÇÕES FINAIS

Sabemos há muito tempo que a química teórica já chegou a um ponto em seu desenvolvimento no qual cálculos teóricos tem, em muitas situações, o status de experimento, dada a sua precisão. A química teórica tem sido utilizada intensamente no auxílio da interpretação de resultados experimentais, na avaliação da viabilidade de determinadas rotas sintéticas, na compreensão de fenômenos nos quais as técnicas experimentais não podem ainda ser diretamente utilizadas, no estudo de caminhos de reações, de barreiras de reação, da influência de ligantes em compostos de coordenação, na espectroscopia do infra-vermelho (IR) e no VIS-UV, etc⁵¹. No entanto, sabemos que ainda é um grande desafio interpretar e transmitir as informações, geradas pelos cálculos teóricos, numa linguagem capaz de ser facilmente compreendida por um químico experimental. Em relação a este último ponto, a Teoria do Funcional de Densidade apresenta uma grande vantagem. A descrição da energia eletrônica e das propriedades moleculares a partir de uma quantidade observável, a densidade eletrônica, permite uma comunicação mais objetiva entre teóricos e experimentais. Vale ressaltar que a Teoria do Funcional de Densidade, em seu formalismo, levanta o véu sobre importantes aspectos da química que ainda eram vistos à luz do empirismo. Conceitos como potencial químico, m, maciez, dureza e índice de reatividade formalmente definidos a partir da DFT abrem uma nova perspectiva na compreensão dos fenômenos químicos. Tais quantidades podem ser combinadas em uma equação fundamental:

O potencial químico e seus derivados são atualmente objeto de intensa pesquisa (veja, por exemplo, as referências 29,30,52,53) A DFT fornece uma descrição da estrutura dos átomos, moléculas, clusters, superfícies e sólidos. De acordo com Kohnⁱⁱⁱ et al.²⁹, e com eles estamos em pleno acordo, DFT é uma linguagem conveniente e universal para a teoria de estrutura eletrônica, a qual ajuda substancialmente a unificar a química orgânica, a inorgânica, química de superfície e a ciência dos materiais. Podemos ainda acrescentar que as dificuldades em implementar o cálculo de índices de reatividade (maciez, dureza, funções de Fukui, etc.) de forma conveniente a um custo computacional acessível estão sendo devidamente tratados. Deste modo, acreditamos que a DFT ocupará um espaço ainda maior, tornando-se também uma ferramenta de grande utilidade no estudo de sistemas bioquímicos, interações droga-aceptor e mecanismos de atividade biológica.

AGRADECIMENTOS

Gostaria de agradecer igualmente aos Profs. Eucler B. Paniago e Wagner B. De Almeida pelas contribuições e sugestões a este trabalho. Agradecemos também ao Centro Nacional de Processamento de Alto Desempenho - CENAPAD-MG/CO pelo acesso aos seus recursos computacionais e ao Conselho Nacional de Desenvolvimento Científico - CNPq, pelo apoio financeiro.

REFERÊNCIAS

1. Drude, P.; Annalen der Physik 1900, 1, 566. Drude, P.; Annalen der Physik 1900, 3, 369.

2. Fermi, E.; Rend. Accad. Lincei 1927, 6, 602. Thomas, L. H.; Proc. Camb. Phil. Soc. 1927, 23, 542.

3. Dirac, P. A. M.; Proc. Camb. Phil. Soc. 1930, 26, 376.

4. Parr, R. G.; Yang, W.; Density Functional Theory of Atoms and Molecules, Oxford University Press, Oxford, 1989.

5. Dreizler, R, M,; Gross, E. K. U.; Density Functional Theory. An Approach to the Quantum Many-Body Problem, Springer-Verlag, Berlin, 1990.

6. Hohenberg, P.; Kohn, W.; Phys. Rev. 1964, 136, B864.

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

8. Morgon, N. H.; Custódio, R.; Quim. Nova 1995, 18, 44.

9. Gilbert, T. L.; Phys. Rev. B 1975, 12, 2111.

10. Levy, M.; Proc. Natl. Acad. Sci. USA 1979, 76, 6062. Levy, M.; Phys. Rev. A. 1982, 26, 1200. Levy, M; Perdew, J. P.; Phys. Rev. A 1985, 32, 2010.

11. Baerends, E. J.; Gritsenko, O. V.; J. Phys. Chem. A 1997, 30, 5383.

12. Parr, R. G.; Yang, W., Annu. Rev. Phys. Chem. 1995, 46, 701.

13. Perdew, J. P.; Burke, K.; Intern. J. Quantum Chem. 1996, 57, 309.

14. Vosko, S. H.; Wilk, L.; Nusair, M.; Can. J. Phys. 1980, 58, 1200.

15. Slater J. C.; Phys. Rev. 1951, 81, 385.

16. Becke, A. D.; J. Chem. Phys. 1986, 84, 4524.

17. Perdew, J. P.; Phys. Rev. Lett. 1985, 55, 1665.

18. Perdew, J. P.; Yue, W.; Phys. Rev. B 1986, 33, 8800.

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

20. Becke, A. D.; Phys. Rev. A 1988, 38, 3098.

21. Perdew, J. P.; Phys. Rev. B 1986, 33, 8822. Ibid. 1986, 34, 7406.

22. Janak, J. F.; Phys. Rev. B 1978, 18, 7165.

23. Mermin, N. D.; Phys. Rev. 1965, 137, A1441.

24. Perdew, J. P.; Parr, R. G.; Levy, M.; Balduz, J. L. Jr.; Phys. Rev. Lett. 1982, 49, 1691.

25. Perdew, J. P.; Zunger, A.; Phys. Rev. B 1981, 23, 5048.

26. Russier, V.; Phys. Rev. B 1992, 45, 8894.

27. Perdew, J. P., In.: Density Functional Methods in Physics, R. M. Dreizler and J. da Providência, Eds; NATO ASI Series-B Bol. 123, 1985.

28. Duarte, H. A.; Salahub, D. R.; J. Chem. Phys. 1998, 108, 743.

29. Kohn, W; Becke, A. D.; Parr, R. G.; J. Phys. Chem. 1996, 100, 12974.

30. Chermette, H; J. Comp. Chem. 1999, 20, 129.

31. Parr, R. G.; Donnelly, R. A.; Levy, M. and Palke, W. E.; J. Chem. Phys. 1978, 68, 3801.

32. Mulliken, R. S.; J. Chem. Phys. 1934, 2, 782.

33. Gyftopoulos, E. P.; Hatsopoulos, G. N.; Proc. Natl. Acad. Sci. USA 1965, 60, 786.

34. Parr, R. G.; Donnelly, R. A.; Levy, M.; Palke, W. E.; J. Chem. Phys. 1978, 68, 3801.

35. Sanderson, R. T.; J. Chem. Educ. 1952, 29, 539. Ibid. 1954, 31, 238.

36. Sanderson, R. T.; Chemical Bounds and Bound Energy; Academic: New York, 1971.

37. Pearson, R. G.; J. Chem. Educ. 1987, 64, 561.

38. Parr, R. G.; Pearson, R. G.; J. Am. Chem. Soc. 1983, 105, 105.

39. Pearson, R. G.; J. Am. Chem. Soc. 1963, 85, 3533.

40. Gutiérrez-Oliva, S.; Letelier, J. R.; Toro-Labbé, A.; Mol. Phys. 1999, 96, 61.

41. Solá, M.; Toro-Labbé, A.; J. Phys. Chem. A 1999, 103, 8847.

42. Yang, W.; Parr, R. G.; Proc. Natl. Acad. Sci. USA 1985, 82, 6723.

43. Chattaraj, P. K.; Lee, H.; Parr, R. G.; J. Am. Chem. Soc. 1991, 113, 1855.

44. Péreza, P.; Toro-Labbé, A.; J. Phys. Chem. A 2000, 104, 1557.

45. Lee, C., Yang, W., Parr, R. G., J. Mol. Structure (Theochem) 1987, 163, 305.

46. Perdew, J. P.; Parr, R. G. ; Levy, M.; Balduz, J. L.; Phys. Rev. Lett. 1982, 49, 1691.

47. Liu, S; Parr, R. G.; J. Chem. Phys. 1997, 106, 5578, 5586.

48. Toro-Labbé, A.; J. Phys. Chem. A. 1999, 103, 4398.

49. Jaque, P.; Toro-Labbé, A.; J. Phys. Chem. A 2000, 104, 995.

50. Gutiérrez-Oliva, S.; Jaque, P.; Toro-Labbé, A.; J. Phys. Chem. A 2000, 104, 8955.

51. Veja, por exemplo, alguns dos trabalhos desenvolvidos em nosso laboratório: Haslett, T. L., Fedrigo, S.; Bosnick, K.A.; Bier, K. D.; Moskovits, M.; Duarte, H. A.; Salahub, D.; J. Am. Chem. Soc. 2000, 122, 6039. Ladeira, A. C. Q.; Alves, M. C. M.; Duarte, H. A.; Ciminelli, V. S. T.; Geochim. Cosmochim. Acta 2001, 65, 1211. Oliveira, M. A.; Pernaut, J-M.; Duarte, H. A.; De Almeida, W. B.; J. Phys. Chem. A 2000, 104, 8256. Duarte, H. A.; Salahub, D. R.; Haslett, T.; Moskovits, M.; Inorg. Chem. 1999, 38, 3895. Duarte, H. A.; Salahub, D. R.; Top. Catal. 1999, 9, 123. Duarte, H. A.; Carvalho, S.; Paniago, E. B.; De Almeida; W. B.; J. Inorg. Biochem. 1998, 72, 71.

52. De Proft, F.; Liu S.; Parr, R. G.; J. Chem. Phys. 1997, 107, 3000.

53. Mineva, T.; Russo, N.; Sicilia, E.; Toscano, M.; J. Chem. Soc. Far. Trans. 1997, 93, 3309.

ⁱⁱ O diferencial de um funcional F[f(x)] é dado por

ⁱⁱⁱ Agraciado, juntamento com o J. Pople, com o prêmio Nobel de química, 1998.

ⁱ

Potencial externo normalmente consiste na posição e na carga dos núcleos dos átomos que formam a molécula.

Datas de Publicação

• Publicação nesta coleção
  31 Ago 2001
• Data do Fascículo
  Ago 2001

Histórico

• Aceito
  06 Nov 2000
• Recebido
  30 Maio 2000