Termocromismo em soluções de alcóxidos de Vanádio(IV): uma abordagem pela modelagem molecular

Thermochromic behavior of vanadium(IV) alkoxides in solution: a molecular modeling aproach

thermochromism; vanadium(IV) alkoxides; molecular modeling

ARTIGO

Termocromismo em soluções de alcóxidos de Vanádio(IV): uma abordagem pela modelagem molecular

Thermochromic behavior of vanadium(IV) alkoxides in solution: a molecular modeling aproach

Alexandre C. de Freitas; Kátia C. M. Westrup; Giovana G. Nunes; Denis J. Gulin; Roberto L. A. Haiduke; Jaísa F. Soares; Eduardo L. de Sá^*

Departamento de Química, Setor de Ciências Exatas, Universidade Federal do Paraná, CP 19081, 81530-900 Curitiba - PR, Brasil

ABSTRACT

The thermochromic behavior exhibited by vanadium(IV) alkoxides, [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] and [V₂(μ-ONep)₂(ONep)₆ ], OPrⁱ = isopropoxide and ONep = neopentoxide, was studied by molecular modeling using DFT, TDDFT and INDO/S methods. The vibrational and electronic spectra calculated for [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] were very similar to the experimental data registered for crystalline samples of the complex and for its solutions at low temperature (< 210 K), while spectra recorded at high temperature (> 315 K) were compatible with those calculated for the monomeric form, [V(OPrⁱ)₄]. These results consistently point to a monomer/dimer equilibrium as an explanation for the solution thermochromism of {V(OPrⁱ)₄}_n. In spite of the structural similarity between [V₂(μ-ONep)₂(ONep)₆ ] and [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] in the solid state, the thermochromic behavior of the former could not be explained by the same model, and the possibility of tetranuclear aggregation at low temperatures was also investigated.

Keywords: thermochromism; vanadium(IV) alkoxides; molecular modeling.

INTRODUÇÃO

Quando a variação de temperatura provoca mudanças estruturais em nível microscópico, quer sejam elas drásticas como alterações conformacionais, ou sutis como na variação do estado de spin, temos o fenômeno de termocromismo. Assim como em todo equilíbrio químico, processos reversíveis que envolvem a agregação e a dissociação de moléculas - que podem explicar o comportamento termocrômico de um sistema - são dependentes da temperatura. Obviamente, para que estes processos impliquem em termocromismo, as formas agregadas e dissociadas devem exibir diferenças marcantes nos espectros de absorção ou emissão de fótons.

Como a temperatura não é uma variável que aparece na Equação de Schrödinger, não há uma correlação direta entre o uso da mecânica quântica e da modelagem molecular como ferramentas de interpretação do comportamento termocrômico de um sistema. Quando a realização de medidas espectroscópicas com temperatura variável não é possível, a proposição da natureza das espécies químicas envolvidas e do mecanismo pelo qual o termocromismo opera deve ser embasada no senso químico do investigador. Apesar de não ser um método direto ou que conduza necessariamente a conclusões definitivas, a abordagem do termocromismo pela modelagem molecular permite que o investigador proponha várias estruturas e mecanismos, teste-os através do confronto com resultados experimentais ou mesmo sugira outras explicações para o fenômeno.

Os alcóxidos homolépticos de metais de transição têm sido empregados como precursores de óxidos metálicos de grande aplicação tecnológica.¹ Este uso pode ser explicado pela facilidade com que óxidos de alta pureza e de estequiometria e estrutura definidas podem ser obtidos a partir destes alcóxidos, principalmente através do processo sol-gel ou pelos métodos de deposição química ou física de vapor (CVD e PVD).^1,2 As características dos produtos advindos destes processos são amplamente dependentes das estruturas moleculares dos precursores em solução.

Acredita-se que muitos alcóxidos de metais de transição (bloco d) sejam protagonistas de equilíbrios de agregação e dissociação entre as formas monomérica e dimérica,³ comportamento que pode ser relacionado com a facilidade de formação de ligações M-(μ-OR)-MOR, com alcóxidos em pontes μ₂ ou μ₃. Já na química de alcóxidos homolépticos de metais do bloco f, a oligomerização é um fenômeno comum, pois o grande raio iônico do centro de coordenação permite que o número de ligantes ao seu redor seja elevado. Apesar destes precedentes, o sistema em estudo no presente trabalho é o primeiro em que o fenômeno da oligomerização é acompanhado de termocromismo.

O alcóxido de vanádio(IV) formado em presença de íons isopropóxido é conhecido desde 1962 quando foi sintetizado por Bradley, que inicialmente acreditava se tratar de [V(OPrⁱ)₄], onde OPrⁱ = isoproprila.⁴ Entretanto, a determinação da massa molar do produto por crioscopia em benzeno mostrou que se tratava de uma estrutura de agregação com um grau de associação de 1,17. Hoje sabe-se que, em fase cristalina, se tem um composto dimérico de fórmula [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆].⁵ Recentemente, utilizando uma outra rota sintética, [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] foi obtido em forma de monocristais,⁶ que, dissolvidos em tolueno ou hexano, apresentaram um comportamento termocrômico marcante, com colorações que variavam do amarelo-ouro (abaixo de 210 K) ao azul-cobalto (a partir de 315 K), passando pelo verde nas temperaturas intermediárias.⁷ Um comportamento termocrômico semelhante foi observado para um outro alcóxido de vanádio(IV), o [V₂(μ-ONep)₂(ONep)₆ ], cuja estrutura cristalina é análoga à do [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆].⁸ Entretanto, em especial no sistema com neopentóxido, não há evidências espectroscópicas inequívocas sobre a natureza das espécies moleculares que mudam a composição do sistema em resposta à variação da temperatura.

Este trabalho propôs-se a sugerir a natureza das espécies químicas responsáveis pelo comportamento termocrômico tão incisivo destes alcóxidos de vanádio(IV), [V₂(μ-OR)₂(OR)₆ ], R = isopropila ou neopentila. Para que este objetivo fosse alcançado, as ferramentas de modelagem molecular, com ênfase no emprego da Teoria do Funcional de Densidade (DFT, Density Functional Theory), da Teoria do Funcional de Densidade Dependente do Tempo (TDDFT, Time Dependent Density Functional Theory) e da parametrização semi-empírica INDO/S, foram utilizadas.

PARTE EXPERIMENTAL

Todas as operações relativas à síntese e à caracterização dos complexos foram realizadas sob atmosfera de N₂ (pureza 99,999%, White Martins) ou sob vácuo de 10^-3-10^-4 Torr, utilizando técnicas de Schlenk ou glove-box. Os solventes foram purificados por métodos padronizados,⁹ secados exaustivamente e destilados antes do uso. Cristais dos alcóxidos de vanádio [V₂(μ-ONep)₂(ONep)₆ ] e [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] foram obtidos em rendimentos de 40 e 60%, respectivamente, por uma rota sintética simples e reprodutível.⁸ Os espectros de absorção na região do infravermelho (400 a 4000 cm^-1) com transformada de Fourier (IV-TF) foram obtidos em um instrumento BIO-RAD, Excalibur Series (FTS 3500 GX). Os espectros foram registrados em soluções com hexano ou em emulsões das amostras pulverizadas em óleo mineral (Nujol). Os espectros eletrônicos (UV-Vis, 190 a 900 nm) foram obtidos em um instrumento Shimadzu UV-2401 PC, em soluções com hexano, tolueno ou a partir de emulsões em óleo mineral.

MÉTODOS COMPUTACIONAIS

Os complexos de vanádio(IV) com grupos isopropóxido e neopentóxido, em fase sólida cristalina, apresentam-se sob a forma de dímeros cujas fórmulas moleculares são, respectivamente, [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] e [V₂(μ-ONep)₂(ONep)₆ ]. As estruturas destas espécies foram determinadas por difratometria de raios X e estão mostradas nas Figuras 1 e 2.^6,8 Estes resultados cristalográficos serviram de ponto de partida para a otimização das geometrias das formas diméricas em solução com tolueno. A inclusão do solvente foi feita pelo emprego do método PCM (Polarizable Continuum Model).¹⁰ As geometrias das formas monoméricas foram geradas a partir das estruturas cristalográficas das formas diméricas, através da eliminação de átomos escolhidos com base em critérios de simetria molecular. Os cálculos foram feitos inicialmente utilizando-se a aproximação HF/6-31G¹¹ disponível no pacote de programas GAMESS,¹² e depois no nível B3LYP/LANL2DZ,¹³ disponível no pacote Gaussian03,¹⁴ ambos disponibilizados no Laboratório Central de Processamento de Alto Desempenho da Universidade Federal do Paraná (LCPAD/UFPR). A partir das estruturas otimizadas, os espectros vibracionais foram calculados utilizando a metodologia DFT com a mesma função de base empregada na otimização da geometria. As estruturas moleculares constituíram-se como mínimos locais, o que pôde ser caracterizado pela ausência de frequências vibracionais negativas nos espectros calculados. Os espectros eletrônicos foram obtidos utilizando-se as metodologias TDDFT e ZINDO,^15,16 implementadas no programa Gaussian03, empregando-se 30 estados excitados como espaço ativo.

RESULTADOS E DISCUSSÃO

Os espectros eletrônicos e vibracionais calculados não se alteraram sensivelmente com a inclusão do solvente utilizando-se o método PCM. Este resultado é congruente com aquele obtido por Nemykin e colaboradores no estudo dos fatores que influenciam as energias de excitação vertical de ftalocianinas.¹⁷ Desta forma, nas discussões que se seguem não serão feitas distinções entre resultados obtidos considerando-se ou não a presença do solvente.

Sistema [V₂(μ-OPrⁱ)₂ (OPrⁱ)₆]/[V(OPrⁱ)₄]

Para o composto binuclear, [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆], a geometria de partida foi a estrutura determinada por difratometria de raios X, com multiplicidade de spin (2S+1) igual a 3. Após convergência, a geometria otimizada deste dímero manteve uma precisão, em relação aos dados experimentais, de cinco casas decimais para comprimentos de ligação e de três casas para ângulos de ligação.

Para a espécie monomérica, [V(OPrⁱ)₄], uma estrutura gerada a partir de comprimentos e ângulos de ligação padrão, com multiplicidade de spin do estado fundamental igual a 2, foi o ponto inicial do cálculo. Semelhantemente ao descrito por Haaland et al. para a estrutura molecular do complexo [V(OBu^t)₄] (Bu^t = terc-butila) determinada em fase gasosa por difração de elétrons,¹⁸ o poliedro de coordenação calculado para o [V(OPrⁱ)₄] apresenta-se ligeiramente achatado em relação a um tetraedro perfeito, adotando a simetria S₄. Os ângulos médios de ligação (O-V-O) calculados para o [V(OPrⁱ)₄] são 112,2(16)º e 108,4(2)º, comparados com 115,1(19)º e 106,7(11)º no [V(OBu^t)₄], o que indica uma distorção menor no complexo com isopropóxido. Esta distorção pode ser atribuída a uma combinação dos efeitos (i) da doação π do alcóxido para os orbitais d_x²_-y² e/ou d_z² do vanádio(IV) e (ii) da repulsão espacial entre os ligantes volumosos.¹⁸

O espectro eletrônico experimental de uma solução de [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] em tolueno apresenta uma mudança nítida e reversível com a variação da temperatura (Figura 3). Nota-se claramente a presença de uma banda em 674 nm em temperaturas altas, que diminui em intensidade com o resfriamento do sistema; concomitantemente, observa-se o aparecimento de uma banda em 470 nm. O ponto isosbéstico em 530 nm sugere a presença de somente duas espécies químicas em equilíbrio na solução (formas azul e amarela); a cor verde resultaria da coexistência destas formas em temperaturas intermediárias. A espécie formada a altas temperaturas seria o complexo mononuclear [V(OPrⁱ)₄], a partir do qual se cristaliza o [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆]. A espécie predominante em baixas temperaturas seria uma molécula binuclear, como a da estrutura cristalina. Ocorreria assim um equilíbrio de agregação monômero-dímero. A Tabela 1 mostra os valores dos sinais experimentais e dos calculados para estas duas espécies em solução.

A análise da Tabela 1 mostra que a proposta mecanística feita no parágrafo anterior é condizente com os resultados experimentais. No caso da espécie monomérica, são previstos três sinais, com um máximo em 678 nm calculado por TDDFT (B3LYP/LANL2DZ). Experimentalmente, há somente um sinal bastante alargado com máximo em 674 nm, que pode perfeitamente estar ocultando os sinais de menor intensidade previstos pelo cálculo TDDFT. Para a espécie supostamente presente à baixa temperatura, há a previsão da banda observada experimentalmente em 470 nm por ambos os métodos teóricos usados. É também comum aos dois métodos teóricos a previsão, para a espécie de baixa temperatura, de uma banda na região de menor energia, entre 700 e 800 nm, também detectada experimentalmente (Figura 3). Assim, os resultados INDO/S e TDDFT reforçam-se mutuamente e apontam para uma confirmação da proposta de agregação.

A comparação dos espectros experimentais de absorção no infravermelho, obtidos tanto para o sólido cristalino quanto para a solução em baixa, média e alta temperaturas (Figura 4), com os espectros calculados para as espécies [V(OPrⁱ)₄] e [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] (Tabela 2) também confirma a hipótese proposta como explicação do efeito termocrômico. A diminuição gradativa da intensidade dos sinais espectroscópicos em 720 (este mascarado pelo sinal intenso, em 724 cm^-1, do hexano), 640, 613, 583 e 512 cm^-1 - característicos da espécie dimérica - com o aumento da temperatura e sem a substituição por sinais característicos da espécie monomérica é indicativa de que o comportamento termocrômico possa ser entendido com base no equilíbrio de dimerização, no qual a espécie dimérica predomina no estado sólido e em baixas temperaturas, transformando-se na espécie monomérica em temperaturas superiores. Ambas as metodologias HF e DFT conseguem prever com razoável acurácia a posição dos sinais na espectroscopia vibracional, com erros geralmente inferiores a 5%.

Para o sistema [V(OPrⁱ)₄] / [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆], o valor experimental do ΔH de dimerização, -54,6 kJ mol^-1, foi calculado a partir das constantes de equilíbrio determinadas em função da temperatura.⁶ Este valor também é compatível com um processo de agregação molecular em que a formação das pontes μ-OPrⁱ contribui para o abaixamento da energia livre do sistema.

Sistema [V₂(μ-ONep)₂(ONep)₆]/[V(ONep) ₄]

De acordo com evidências oriundas da espectroscopia paramagnética eletrônica (RPE), a multiplicidade de spin (2S+1) da forma dimérica [V₂(μ-ONep)₂(ONep)₆ ] corresponde a um estado fundamental tripleto. Esta informação, juntamente com a geometria molecular determinada por difratometria de raios X, formaram o ponto de partida para os cálculos relativos à espécie binuclear. Para o monômero, [V(ONep)₄], uma estrutura gerada a partir de distâncias e ângulos de ligação médios relatados para outros compostos de vanádio e uma multiplicidade de spin do estado fundamental igual a 2 foram as premissas adotadas no início do cálculo.

Também neste caso, como no complexo com isopropóxido, a geometria experimental da espécie dimérica foi mantida, nos cálculos, com uma precisão de cinco casas decimais para distâncias de ligação, e de três casas para ângulos de ligação. Já o [V(ONep)₄] otimizado apresentou uma geometria tetraédrica praticamente perfeita, com um ângulo de ligação O-V-O médio de 109,1(8)º e uma distância de ligação V-O média de 1,7824(3) Å. A ausência de distorção, neste caso, pode ser explicada pela presença, no neopentóxido, de um grupo "espaçador" -CH₂- que não existe no isopropóxido e que afasta o metal da porção mais volumosa da cadeia carbônica do ligante. Em consequência disso, a repulsão intramolecular diminui e o centro de vanádio(IV) assume uma geometria de coordenação mais regular.

Soluções de [V₂(μ-ONep)₂(ONep)₆ ] (Figura 2) em tolueno apresentam, da mesma forma que soluções do complexo análogo [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆], uma variação espectral significativa com a mudança de temperatura (compare as Figuras 3 e ⁵). Neste caso são observados também dois padrões nos espectros: em altas temperaturas, há uma banda com máximo em 665 nm, que diminui de intensidade com o abaixamento da temperatura, ao passo que uma banda alargada em torno de 478 nm se faz mais pronunciada, sendo acompanhada por um crescimento na intensidade de absorção em torno de 730 nm. Devido à semelhança estrutural e dos espectros eletrônicos, o comportamento termocrômico observado para este neopentóxido de vanádio(IV) também foi atribuído, inicialmente, a um equilíbrio de agregação monômero/dímero.

As espécies [V(ONep)₄] e [V₂(μ-ONep)₂(ONep)₆ ] tiveram os seus espectros eletrônicos simulados. Para o [V(ONep)₄], os resultados foram semelhantes aos já apresentados para o [V(OPrⁱ)₄]. Assim, a transição experimental que ocorre com máximo em 665 nm foi calculada pelos métodos TDDFT e INDO/S como ocorrendo em 616 e 681 nm, respectivamente. Para a espécie [V₂(μ-ONep)₂(ONep)₆ ], por sua vez, as transições foram calculadas pelo método TDDFT em 397, 611 e 778 nm. Para ambas as espécies há também uma absorção na região do infravermelho, com força de oscilador pequena, que foi atribuída a uma transição eletrônica do orbital HOMO para o LUMO + 1, ambos baseados em orbitais atômicos d do metal. As demais transições entre orbitais de fronteira não apresentaram força de oscilador significativa.

A análise dos resultados revela claramente que o espectro eletrônico calculado para a espécie [V₂(μ-ONep)₂(ONep)₆ ] difere do espectro experimental obtido em baixas temperaturas: a absorção que deveria ser observada em 478 nm não aparece como resultado dos cálculos. Além disto, os sinais calculados para [V₂(μ-ONep)₂(ONep)₆ ] são aqueles que experimentalmente decrescem com a diminuição da temperatura, quando um comportamento contrário seria o esperado. A hipótese de que a forma dimérica não participa do mecanismo que origina o termocromismo é também sugerida pelo fato de que os espectros vibracionais registrados em solução (Figura 6) não são similares ao espectro do [V₂(μ-ONep)₂(ONep)₆ ] sólido (cristalino), pelo menos nas temperaturas investigadas experimentalmente. Neste contexto, uma explicação alternativa para o efeito termocrômico residiria num equilíbrio entre a forma monomérica e outra de maior grau de agregação, diferente do dímero. O menor grau de agregação superior ao dímero, compatível com a fórmula mínima dos alcóxidos moleculares de vanádio(IV) é o do tetrâmero [V₄(μ-ONep)₆(ONep)₁₀ ], representado tentativamente na Figura 7. Na estrutura do tetrâmero pode-se observar que há centros metálicos em dois ambientes químicos diferentes, ligados a cinco ou seis ligantes neopentóxidos em geometria bipiramidal (base triangular) ou octaédrica distorcidas, respectivamente. A estrutura tetramérica mostrada na Figura 7 foi obtida através da duplicação da estrutura dimérica otimizada, preservando-se ao máximo as distâncias e ângulos de ligação.

Para a espécie tetramérica, as transições foram calculadas pelo método INDO/S, considerando um quinteto (2S+1=5) como o estado fundamental, sendo que os principais sinais calculados se localizam em 428 (extremamente intenso), 499, 523, 775 e 783 nm (intensidades baixas). A análise destes resultados sugere que, se assumirmos que a banda experimental alargada em torno de 730 nm (^{Figura 5}) decorre da sobreposição dos sinais calculados em 775 e 783 nm, e que o mesmo vale para os sinais em 499 e 523 nm em relação ao ombro experimental em 478 nm, a hipótese de uma espécie tetramérica predominante em baixas temperaturas poderá ser considerada consistente.

A análise dos espectros de absorção no infravermelho (Figura 6) permite mostrar que há boa concordância entre a curva experimental para o [V₂(μ-ONep)₂(ONep)₆ ] cristalino e a curva calculada para a espécie dimérica. Porém, a concordância entre os pares (curvas simuladas para as espécies monomérica e tetramérica e as suas possíveis congêneres experimentais, curvas a +45 ºC e a -25 ºC, respectivamente) foi insuficiente para que se pudesse afirmar conclusivamente que o termocromismo neste sistema opera segundo um mecanismo de oligomerização entre a forma monomérica e a tetramérica. Em resumo, com base nas diferenças entre os espectros experimentais e calculados de absorção no infravermelho, é possível afirmar que as espécies ou formas de agregação presentes em solução, em baixas temperaturas, não são as mesmas do estado cristalino.

CONCLUSÕES

O comportamento termocrômico das espécies [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆] e [V₂(μ-ONep)₂(ONep)₆] em solução foi estudado por espectroscopia eletrônica, espectroscopia vibracional e técnicas de modelagem molecular. Os resultados indicam claramente que, nas soluções de [V₂(μ-OPrⁱ)₂(OPrⁱ) ₆], o termocromismo se origina de um equilíbrio de dimerização, enquanto nas soluções de [V₂(μ-ONep)₂(ONep)₆ ] os indícios são para um equilíbrio de tetramerização. Entretanto, a precária concordância entre os espectros de absorção no infravermelho calculados para as espécies [V(ONep)₄] e [V₄(μ-ONep)₆(ONep)₁₀ ] e os espectros experimentais registrados em alta e baixa temperaturas, respectivamente, não permite afirmar com segurança que a forma tetramérica seja majoritária quando o sistema é submetido ao resfriamento.

AGRADECIMENTOS

A. C. de Freitas, K. C. M. Westrup e G. G. Nunes agradecem ao CNPq pelas bolsas de estudos (mestrado, iniciação científica e pós-doutorado). E. L. de Sá e J. F. Soares agradecem ao CNPq (projetos 200413/2005-8 e 476715/2008-3), à FINEP (CT-INFRA) e ao Laboratório Central de Processamento de Alto Desempenho (LCPAD/UFPR) pelo apoio financeiro e pela infra-estrutura para a realização deste trabalho.

REFERÊNCIAS

1. Kessler, V. G.; Chem. Commun. 2003, 1213.

2. Hubert-Pfalzgraf, L. G.; Inorg. Chem. 2003, 6, 102.

3. Gulliver, E. A.; Garvey, J. W.; Wark, T. A.; Hampden-Smith, M. J.; Datye, A.; J. Am. Ceram. Soc. 2005, 74, 1091; Barnhart, D. M.; Clark, D. L.; Gordon, J. C.; Huffman, J. C.; Watkin, J. G.; Zwick, B. D.; Inorg. Chem. 1995, 34, 5416.

4. Bradley, D. C.; Mehta, M. L.; Can. J. Chem. 1962, 40, 1183.

5. Kempe, R.; Spannenberg, A.; Z. Kristallogr. - New Cryst. Struct. 1997, 212, 485.

6. Nunes, G. G.; Friedermann, G. R.; dos Santos, J. L. B.; Herbst, M. H.; Vugman, N. V.; Hitchcock, P. B.; Leigh, G. J.; Sá, E. L.; da Cunha, C. J.; Soares, J. F.; Inorg. Chem. Commun. 2005, 8, 83.

7. Nunes, G. G.; Tese de Doutorado, Universidade Federal do Paraná, Brasil, 2004.

8. Reis, D. M.; Westrup, K. C. M.; Nunes, G. G.; Barison, A.; Ribeiro, R. R.; Sá, E. L. de; Hughes, D. L.; Soares, J. F.; J. Braz. Chem. Soc. 2009, 20, 613; Westrup, K. C. M.; Reis, D. M.; Nunes, G. G.; Barison, A.; Ribeiro, R. R.; Sá, E. L.; Hitchcock, P. B.; Soares, J. F.; Cambridge Structural Database, 2008, deposition number CCDC 705491.

9. Perrin, D. D.; Armarego, W. L. F.; Purification of Laboratory Chemicals, 3^rd ed., Butterworth Heenemann: Oxford, 1997.

10. Tomasi, J.; Persico, M.; Chem. Rev.1994, 94, 2027.

11. Ditchfield, R.; Hehre, W. J.; Pople, J. A.; J. Chem. Phys. 1971, 54, 724.

12. Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.; Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S.; Windus, T. L.; Dupuis, M.; Montgomery, J. A.; J. Comput. Chem. 1993, 14, 1347.

13. Hay, P. J.; Wadt, W. R.; J. Chem. Phys. 1985, 82, 270.

14. Gaussian 03, Revision D.01, Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, Jr., J. A.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A.; Gaussian, Inc., Wallingford CT, 2004.

15. Stratmann, R. E.; Scuseria, G. E.; Frisch, M. J.; J. Chem. Phys. 1998, 109, 8218; Bauernschmitt, R.; Ahlrichs, R.; Chem. Phys. Lett. 1996, 256, 454; Casida, M. E.; Jamorski, C.; Casida, K. C.; Salahub, D. R.; J. Chem. Phys. 1998, 108, 4439.

16. Ridley, J. E.; Zerner, M. C.; Theor. Chim. Acta 1973, 32, 111; Ridley, J. E.; Zerner, M. C.; Theor. Chim. Acta 1976, 42, 223; Bacon, A. D.; Zerner, M. C.; Theor. Chim. Acta 1979, 53, 21.

17. Nemykin, V. N.; Hadt, R. G.; Belosludov, R. V.; Mizuseki, H.; Kawazoe, Y.; J. Phys. Chem. A 2007, 111, 12901.

18. Haaland, A.; Rypdal, K.; Volden, H. V.; Andersen, R. A.; J. Chem. Soc., Dalton Trans. 1992, 891.

Recebido em 2/12/08; aceito em 12/7/09; publicado na web em 27/11/09

Datas de Publicação

• Publicação nesta coleção
  05 Fev 2010
• Data do Fascículo
  2010

Histórico

• Aceito
  12 Jul 2009
• Recebido
  02 Dez 2008