JBCS

INTRODUÇÃO

O tópico "ácidos e bases" é apresentado, de forma ampla, durante as disciplinas de química analítica do ensino superior para cursos de Química e Engenharia Química. Contudo, a termodinâmica inerente a esse tema não é abordada nos livros de Química Analítica criando uma lacuna entre a termodinâmica e o equilíbrio químico envolvido nas correspondentes reações.^1-3 De acordo com a literatura pesquisada, não há nenhuma associação das constantes de equilíbrio das reações de dissociação de ácidos ou bases e a relação com a estabilidade das moléculas envolvidas, assim como não se aborda a influência do solvente nas constantes de equilíbrio, como no caso da série de acidez do álcoois que tem completa inversão entre a fase aquosa e a fase gás.⁴

Além disso, nos livros de Química Analítica, não encontrou-se menção sobre propriedades que possam se correlacionar com a acidez ou basicidade como, por exemplo, informações a respeito de possíveis correlações entre acidez/basicidade e a hidrólise de suas correspondentes bases/ácidos conjugados, respectivamente.

Sabe-se que, cada reação, em equilíbrio, tem uma constante de equilíbrio específica e única. Isto ocorre porque cada reação tem uma variação da energia livre de Gibbs (∆G) específica. A relação entre a constante de equilíbrio e a variação da energia livre de Gibbs é: ∆G = -RT lnK. Ou seja, a acidez ou basicidade podem ser calculadas a partir da variação da energia livre de Gibbs da dissociação do ácido ou da base, seja em fase condensada, seja em fase gás. Então, pequenas mudanças nos valores das variações de energia livre de qualquer uma das etapas desse processo podem gerar grandes alterações nos valores K_a ou K_b.⁵ Além disso, uma análise comparativa e qualitativa acerca da acidez ou basicidade de um conjunto de ácidos ou bases da mesma família (mesmo grupo funcional) pode ser feita a partir das estabilidades relativas das espécies químicas envolvidas na reação.⁶

A química computacional é uma poderosa ferramenta que pode auxiliar a compreender melhor algumas das questões levantadas neste trabalho. Com o avanço da química computacional, cada vez mais é necessário e viável explicar fenômenos microscópicos com a ajuda da química quântica.⁷ Então, a área de ensino da Química pode contar com um poderoso recurso, a química computacional, para gerar melhor compreensão acerca de experimentos que compõem um determinado tema ensinado nos livros didáticos.⁸

Neste trabalho, utilizamos os métodos computacionais B3LYP e o CBS-QB3 para calcular valores de variação de energia livre de Gibbs para as reações de dissociação de ácidos e bases e as reações de hidrólise das correspondentes bases ou ácidos conjugados, em fase gás. Com esse trabalho foi estabelecida uma a relação entre estabilidade de cátions metálicos - os ácidos conjugados da reação de dissociação das bases - frente a reação de hidrólise, a correspondente variação da energia livre de Gibbs, e, implicitamente, a correspondente constante de equilíbrio da reação. Também se estabeleceu boa concordância entre as variações da energia livre de Gibbs para a reação de dissociação de ácidos e a reação de hidrólise das correspondentes bases conjugadas. Por fim, comparando-se dados teóricos, em fase gás, com dados experimentais, em fase condensada, mostramos a importância da solvatação no processo de hidrólise de ácidos conjugados e de dissociação de ácidos.

METODOLOGIA

Para esse estudo, foram selecionadas as seguintes moléculas ou espécies químicas: cátions (K⁺; Ca⁺²; Ga⁺³; Ti⁺²; Co⁺²); bases (KOH; Ca(OH)₂; Ga(OH)₃; Ti(H₂O)₆⁺²; Ti(H₂O)₅(OH)⁺¹; Co(H₂O)₆⁺²; Co(H₂O)₅(OH)⁺¹); ácidos (H₂SO₄; H₂SO₃; H₃PO₄; H₃PO₃; HNO₃; HNO₂; HClO₄; HClO); e bases conjugadas (SO₄^-2; SO₃^-2; HSO₄^-; HSO₃^-; PO₄^-3; HPO₄^-2; HPO₃^-2; H₂PO₄^-; H₂PO₃^-; NO₃^-; NO₂^-; ClO₄^-; ClO^-).

Todas essas moléculas foram otimizadas utilizando o nível de teoria B3LYP/6-311G++(2d,2p) da teoria do funcional da densidade (DFT).^9-14 Usamos o pacote GAUSSIAN09 para realização dos cálculos quânticos.¹⁵ A partir da geometria otimizada, foi realizado novo cálculo de otimização, além do cálculo de frequência e propriedades termodinâmicas das correspondentes moléculas usando o método CBS-QB3,¹⁶ cujos cálculos de propriedades termodinâmicas são comparáveis aos valores experimentais, além de ser usado para cálculos de metais de transição na forma iônica e cálculos de pKa.^17,18 É importante ressaltar que os resultados obtidos foram para fase gás e nenhuma relação pode ser estabelecida com dados experimentais obtidos em fase condensada devido à influência do efeito do solvente nos resultados obtidos.

RESULTADOS E DISCUSSÃO

Validação do método

Na Figura 1 são apresentados valores de comprimentos de ligação entre o metal central (Fe⁺², Co⁺², Cu²⁺) e o átomo de oxigênio dos grupos aqua dos complexos M(H₂O)₆⁺² obtidos por meio da difração de raios-X, por cálculos computacionais em que a otimização foi feita inicialmente pelo nível de teoria B3LYP/6-311G++(2d,2p), seguida da otimização pelo método CBS-QB3, ou seja, B3LYP/6-311G++(2d,2p)//CBS-QB3. Para reforçar a comparação, foram acrescidos também dados teóricos calculados por outros métodos.

Figura 1. Geometria otimizada dos complexos Co(H₂O)₆ ⁺², Cu(H₂O)₆ ⁺² e Fe(H₂O)_{6 +2} com distâncias interatômicas, em Angstroms, obtidos pelos métodos B3LYP/6-311G++(2d,2p)//CBS-QB3, por outros métodos teóricos^20,21 (entre parênteses) e obtidos experimentalmente (entre colchetes)¹⁹

Pode-se observar a eficácia do método B3LYP/6-311G++(2d,2p)//CBS-QB3 para cálculos envolvendo metais de transição na forma iônica ao compararmos dados geométricos obtidos teoricamente com aqueles experimentais. A diferença entre os valores obtidos experimentalmente e aqueles obtidos a partir da otimização via B3LYP/6-311G++(2d,2p)//CBS-QB3 não ultrapassa 0,15 Å (Figura 1).^19,20 Além disso, os valores calculados teoricamente nesse trabalho também concordam com obtidos teoricamente por outros métodos de cálculo.^20,21

Relação entre variação da energia livre de Gibbs de hidrólise de cátions metálicos e seus correspondentes raios atômicos

Na Tabela 1, são apresentadas as reações de hidrólise dos cátions, Ti⁺², Fe⁺², Co⁺², K⁺, Ca⁺² e Ga⁺³, metais do quarto período da tabela periódica, e suas correspondentes variações de energia livre de Gibbs (∆G_hid) e de entalpia (∆H_T) calculadas teoricamente, em fase gás, e variação de entalpia obtida experimentalmente, em fase condensada (∆H_EXP).²² Seus raios atômicos do tipo iônico também são apresentados na Tabela 1.²³ De acordo com a literatura, a hidrólise de cátions metálicos de elementos de transição, como titânio, ferro e cobalto, leva à formação de complexos hexacoordenados por ligantes de molécula de água do tipo [M(H₂O)₆]⁺² de geometria bipirâmide tetraédrica.^24,25

À medida que avançamos da esquerda para direita no quarto período da tabela periódica, o raio atômico diminui, e a estabilidade dos correspondentes cátions metálicos também diminui, porque quanto mais localizada a carga, mais instável é a espécie química. Consequentemente, as reações de hidrólise desses cátions ficam mais exergônicas e exotérmicas quando o tamanho do correspondente cátion diminui (ver Tabela 1). Podemos esperar que a basicidade das correspondentes bases siga tendência inversa, como será mostrado posteriormente.

Pode-se observar que os valores teóricos de variação da entalpia para a reação de hidrólise dos cátions seguem a mesma tendência dos valores de ∆G_T. Contudo, a tendência dessas grandezas em fase gás é contrária àquela do ∆H_EXP em fase condensada. A inversão das tendências de entalpia em fase gás e em fase condensada para a reação de hidrólise mostra que o efeito do solvente desempenha importante papel no valor experimental das propriedades termodinâmicas e constante de equilíbrio das reações de hidrólise. Como dito anteriormente, quanto menor o cátion metálico, mais localizado a sua carga atômica, o que leva a maior instabilidade em fase gás. Mas, em fase condensada, a diminuição do raio atômico e maior localização da carga atômica leva à maior interação com as moléculas de água e à maior estabilidade do sistema cátion e camada de solvatação, o que explica a tendência inversa de valores de entalpia comparando-se a fase gás e a fase condensada da reação de hidrólise dos cátions estudados. Assim, com esses resultados, estabelecemos a relação entre estabilidade, entalpia/energia livre de Gibbs, e, implicitamente, a constante de equilíbrio, envolvendo as reações de hidrólise de cátions metálicos.

A Figura 2 mostra uma relação exponencial entre ∆G_hid e raio atômico dos correspondentes metais envolvidos na hidrólise. Pelo que pode ser observado na tendência exponencial da Figura 2, com o aumento do raio atômico do metal, mais estável o correspondente cátion e, consequentemente, mais endergônica (ou menos favorável) é a reação de hidrólise do correspondente cátion metálico.

Figura 2. Variaçao da energia livre de Gibbs da reaçao de hidrólise de Ti⁺², Fe⁺², Co⁺², K⁺, Ca⁺² e Ga⁺³ versus seus correspondentes raios atômicos

Relação da variação da energia livre de Gibbs com as reações de dissociação de bases e os raios atômicos dos metais envolvidos na base

A Tabela 2 mostra os dados da variação da energia livre de Gibbs, em fase gás, para a reação de dissociação da base (∆G_B) envolvendo os metais dos elementos representativos do quarto período da tabela periódica.

Devido à complexidade das estruturas e reações de equilíbrio envolvidas para os cátions metálicos de transição hidratados, nós optamos por não calcular as possíveis reações de dissociação dessas espécies químicas.^26-28

A análise dos dados da Tabela 2 nos mostra que também há uma relação entre variação da energia livre de Gibbs da dissociação da base e o correspondente raio atômico do metal do quarto período da tabela periódica envolvido na base. Nesse caso, há uma relação inversa entre tamanho do raio atômico e ∆G_B. Quanto menor o cátion e maior sua carga, mais instável ele será, o que, a partir da equação ∆G = -RT lnK, leva a uma maior variação da energia livre de dissociação da correspondente base. Essa relação entre tamanho do cátion e sua estabilidade está de acordo com o fato de que quanto mais deslocalizada a carga de uma espécie química (o que ocorre com o maior tamanho do cátion, por exemplo), maior será sua estabilidade. Sabendo-se que as reações de dissociação das bases mostradas na Tabela 2 são reações endergônicas, e que isso leva a uma relação inversa entre a constante de equilíbrio e energia livre de Gibbs, consequentemente, quanto menor o cátion (ácido conjugado da base), maior é a variação de energia livre de Gibbs da reação de dissociação da base e menor é a correspondente constante de basicidade. Novamente, estabelecemos a relação entre estabilidade, propriedade termodinâmica e equilíbrio químico.

Relação da variação da energia livre de Gibbs com as reações de dissociação de ácidos e as hidrólises de suas correspondentes bases conjugadas

Foi analisada a relação entre energia livre de Gibbs entre duas reações: a dissociação dos ácidos e a reação de hidrólise da correspondente base conjugada, cujos valores são mostrados nas Tabelas 3 e 4, respectivamente. As equações químicas gerais estudadas são mostradas abaixo:

Os dados das Tabelas 3 e 4 foram plotados na Figura 3, mostrando uma relação inversa entre as magnitudes das duas energias livres de Gibbs (de dissociação do ácido e da hidrólise da correspondente base conjugada).

Figura 3. Variaçao da energia livre de Gibbs da reaçao de dissociaçao de ácidos (∆G_A) versus variaçao da energia livre de Gibbs de hidrólises das correspondentes bases conjugadas (∆G_hid)

O índice de correlação no gráfico plotado da Figura 3 foi igual a um, indicando que, para o caso dos ácidos, quanto mais exergônica (ou menos endergônica) a reação de hidrólise, mais endergônica será a reação de dissociação do ácido. Ou seja, existe então uma relação inversa e perfeita entre os valores de diferença de energia livre de Gibbs da reação de dissociação do ácido e aqueles da correspondente reação de hidrólise da base conjugada. Então, pode-se observar que quanto mais estável for a base conjugada do ácido frente à reação de hidrólise (o que se reflete numa maior endergonicidade da correspondente reação), maior será a constante de acidez do correspondente ácido (que se reflete na menor endergonicidade da correspondente reação).

Em última análise, nós observamos que não há relação linear entre os valores de variação de energia livre de Gibbs da dissociação dos ácidos estudados em fase gás obtidos teoricamente e os correspondentes valores experimentais em fase condensada. O fato de não haver relação linear entre os dados experimental e teórico indica que a solvatação das moléculas envolvidas na reação de dissociação de ácidos (obtida experimentalmente) tem importante influência na variação de energia livre de Gibbs para essa reação em fase condensada.

A acidez e basicidade são propriedades termodinâmicas em que as constantes de acidez e basicidade dependem da variação da energia livre de Gibbs das reações de dissociação de ácidos ou bases. Em fase gás, os fatores que afetam essas propriedades são apenas intrínsecos (dependendo da energia eletrônica da molécula, da entalpia e da entropia), enquanto que em fase condensada, além de fatores intrínsecos, existem fatores extrínsecos que dependem das interações entre solvente e soluto, em que o soluto é o ácido/base conjugada ou base/ácido conjugado. Os fatores extrínsecos podem alterar total ou parcialmente o comportamento da acidez ou basicidade. A falta de relação linear entre os dados experimental (fase condensada) e teórico (fase gás) dos ácidos estudados indica que o fator extrínseco altera parcialmente o comportamento da acidez dos ácidos inorgânicos em comparação com o comportamento em fase gás.

Pelos resultados apresentados, é recomendável que o aluno entenda primeiramente o comportamento da acidez e basicidade em fase gás. E, somente então, numa análise seguinte, seja ensejada a compreensão do comportamento da acidez e basicidade em fase aquosa por parte do aluno, pois assim ele poderá perceber o quanto o fator extrínseco altera o comportamento da acidez e basicidade frente aos fatores intrínsecos.

CONCLUSÕES

Com esse trabalho, mostramos como a Físico-Química e a Química Computacional podem ajudar a ampliar a compreensão dos fatores associados às reações de dissociação entre ácidos e base. Dentre as relevantes contribuições, apontamos que:

1 - Existe uma relação direta entre ∆G_hid e raio atômico dos correspondentes cátions metálicos envolvidos na hidrólise, indicando que quanto maior o raio atômico, menos espontânea será a reação de hidrólise do correspondente íon;

2 - Quanto mais favorável termodinamicamente for a reação de hidrólise de uma espécie química, menor será a constante de acidez do correspondente ácido. Isso ocorre por causa da relação inversa entre constante de acidez e energia livre de Gibbs para reações endergônicas, em que o aumento da endergonicidade da reação de dissociação do ácido leva a uma menor constante de acidez;

3 - O efeito do solvente desempenha papel importante nas reações de dissociação de ácidos em fase condensada.

Acreditamos então que esse trabalho poderá contribuir para o maior entendimento do tópico ácido-base por parte de um aluno que cursa disciplinas de Química Analítica especificamente.

AGRADECIMENTOS

Agradecemos à CAPES, CNPq e FAPERN, projeto PPP 005/2011, pelo o auxílio financeiro para realização dessa pesquisa.

REFERÊNCIAS

1. Skoog, D. A.; West, D. M.; Holler, F. J.; Crouch, S. R.; Fundamentos da Química Analítica, trad. 8ª ed., Thomson: São Paulo, 2003.

2. Harris, D. C.; Análise Química Quantitativa, trad. 8ª ed., LTC: Rio de Janeiro, 2012.

3. Christian, G. D.; Analytical Chemistry, 6ª ed., Wiley: New York, 2003.

4. Boand, G.; Houriet, R.; Gaumann, T.; J. Am. Chem. Soc. 1983, 105, 2203.

5. Bell, R. P.; The Proton in Chemistry, 1ª ed., Chapman and Hall: Londres, 1973.

6. Anslyn, E. V.; Dougherty, D. A.; Modern Physical Organic Chemistry, 1ª ed., University Science Books: Sausalito, 2006.

7. Cramer, C. J.; Essentials of Computational Chemistry: Theories and Models, 2ª ed., Wiley: Chichester, 2004.

8. Leach, A. R.; Molecular Modelling: Principles and Applications, 2ª ed., Prentice Hall: Harlow, 2001.

9. Becke, A. D.; J. Chem. Phys. 1993, 98, 1372.

10. Becke, A. D.; J. Chem. Phys. 1993, 98, 5648.

11. Dunning, T. H.; J. Chem. Phys. 1989, 90, 1007.

12. Lee, C.; Yang, W.; Parr, R. G.; Phys. Rev. B: Condens. Matter 1988, 37, 785.

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

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

15. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, N. J.; Klene, M.; Knox, J. E.; 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.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J.; Gaussian 09w; Revision C.01; Gaussian, Inc., Wallingford CT, 2009.

16. Montgomery, J. A.; Frisch, M. J.; Ochterski, J. W.; Petersson, G. A.; J. Chem. Phys. 1999, 110, 2822.

17. Charmet, A. P.; Quartarone, G.; Ronchin, L.; Tortato, C.; Vavasori, A,; J. Phys. Chem. A 2013, 117, 6846.

18. Altinay, G.; Citir, M.; Metz, R. B.; J. Phys. Chem. A 2010, 114, 5104.

19. Ohtaki, H.; Yamaguchi, T.; Maeda, M.; Bull. Chem. Soc. Jpn. 1976, 49, 701.

20. Semrouni, D.; Isley, W. C.; Clavaguéra, C.; Dognon, J. P.; Cramer, C. J.; Gagliardi, L.; J. Chem. Theory Comput. 2013, 9, 3062.

21. van Duin, A. C. T.; Bryantsev, V. S.; Diallo, M. S.; Goddard, W. A.; Rahaman, O.; Doren, D. J.; Raymand, D.; Hermansson, K.; J. Phys. Chem. A 2010, 114, 9507.

22. Smith, D. W.; J. Chem. Educ. 1977, 54, 540.

23. Shannon, R. D.; Acta Crystallogr. Sect. A 1976, 32, 751.

24. Lincoln, S. F.; Richens, D. T.; Sykes, A. G. Em Comprehensive Coordination Chemistry II; McCleverty, J. A.; Meyer, T. J., eds.; Pergamon: Oxford, 2003, vol.1, section II, cap. 1.25.

25. Ogden, M. I.; Beer, P. D. Em Encyclopedia of Inorganic Chemistry; John Wiley & Sons: Chichester, 2006.

26. Burgess, J.; Metal Ions in Solution, Ellis Horwood: Chichester, 1978.

27. Richens, D. T.; The Chemistry of Aqua Ions: Synthesis, Structure and Reactivity: A Tour Through the Periodic Table of the Elements, 1ª ed., Wiley: New York, 1997.

28. Baes, C. F.; Mesmer, R. E.; The Hydrolysis of Cations, John Wiley & Sons: New York, 1976.

29. Weast, R. C.; Handbook of Chemistry and Physics, 59ª ed., CRC Press: Boca Raton, 1978.