Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Statistical thermodynamics of liquids using the Monte Carlo method. II. Liquid chloroform

Monte Carlo simulation; free energy simulation; liquid chloroform

DIVULGAÇÃO

Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Vania Elisabeth Barlette

Instituto de Química - Universidade Estadual de Campinas - CP 6154 - 13 083-970 - Campinas - SP

Luiz Carlos Gomide Freitas

Departamento de Química - Universidade Federal de São Carlos - CP 676 - 13565-905 - São Carlos - SP

Recebido em 19/2/97; aceito em 14/4/98

Statistical thermodynamics of liquids using the Monte Carlo method. II. Liquid cloroform. Thermodynamic properties and radial distribution functions for liquid chloroform were calculated using the Monte Carlo method implemented with Metropolis algorithm in the NpT ensemble at 298 K and 1 atm. A five site model was developed to represent the chloroform molecules. A force field composed by Lennard-Jones and Coulomb potential functions was used to calculate the intermolecular energy. The partial charges needed to represent the Coulombic interactions were obtained from quantum chemical ab initio calculations. The Lennard-Jones parameters were adjusted to reproduce experimental values for density and enthalpy of vaporization for pure liquid. All thermodynamic results are in excelent agreement with experimental data. The correlation functions calculated are in good accordance with theoretical results avaliable in the literature. The free energy for solvating one chloroform molecule into its own liquid at 298 K and 1 atm was computed as an additional test of the potential model. The result obtained compares well with the experimental value. The medium effects on cis/trans convertion of a hypotetical solute in water TIP4P and chloroform solvents were also accomplished. The results obtained from this investigation are in agreement with estimates of the continuous theory of solvation.

Keywords: Monte Carlo simulation; free energy simulation; liquid chloroform.

INTRODUÇÃO

O líquido clorofórmio é largamente empregado como solvente para estudos de propriedades físicas e químicas das substâncias¹. Existe na literatura uma grande quantidade de estudos teóricos e experimentais sobre esse líquido². No que concerne à simulação computacional, a dinâmica e estrutura do líquido clorofórmio foi investigada por Evans³, Dietz e Heinzinger^4,5, Böhm e Ahlrichs⁶ e Laaksonen e col.⁷. A técnica de Dinâmica Molecular (MD) foi utilizada nessas investigações, e as moléculas de clorofórmio (CHCl₃) representadas por estruturas rígidas com cinco sítios atômicos de interação. Com esses modelos, foram obtidos dados elucidativos sobre a estrutura do líquido, mas os resultados para propriedades termodinâmicas diferem do comportamento observado experimentalmente. Uma representação mais simplificada foi adotada por Jorgensen e col.⁸ em seus estudos de Monte Carlo (MC): a molécula de CHCl₃ érepresentada por quatro sítios, sendo o grupo CH considerado um único sítio. O potencial para a interação intermolecular foi parametrizado empiricamente. Com esse modelo, foram obtidos resultados termodinâmicos para o líquido em boa concordância com dados experimentais. No entanto, a disponibilidade de um modelo para a molécula de clorofórmio com o átomo de hidrogênio considerado explicitamente é desejável, desde que propriedades termodinâmicas e estruturais do líquido clorofórmio, e de suas misturas binárias, são afetadas pela interação de sítios específicos com este átomo. Um exemplo é a importância da interação entre o grupo carbonila da acetona e o átomo de hidrogênio do clorofórmio para explicar o abaixamento da pressão de vapor nesta mistura binária.

O desenvolvimento de um modelo de cinco sítios para a molécula de clorofórmio é o principal objetivo deste trabalho. Para testar a qualidade do modelo proposto, uma série de testes foi realizada utilizando o método de MC. Dada a importância clorofórmio como solvente em reações químicas, foi estudada uma aplicação relevante para invetigar efeitos deste solvente sobre a estrutura molecular.

MODELO

Há dois métodos bastante conhecidos na literatura para obter potenciais de pares do tipo Lennard-Jones mais Coulomb, e consistem basicamente na determinação dos parâmetros de Lennard-Jones e de cargas. Esses métodos são referenciados como métodos empírico e mecânico-quântico ab initio.

O procedimento empírico para gerar parâmetros para as funções potenciais depende da disponibilidade de dados experimentais. Os potenciais são parametrizados para reproduzir resultados experimentais. Seguindo essa abordagem, Jorgensen e col. têm trabalhado no desenvolvimento de modelos para uma série de compostos orgânicos. São os modelos conhecidos na literatura como modelos OPLS (Optimized Potential for Liquid Simulations), especialmente desenvolvidos para estudos de propriedades termodinâmicas e estruturais de líquidos puros e de suas soluções. O modelo de quatro sítios atômicos de interação para a molécula de clorofórmio faz parte do grupo de potenciais OPLS.

A outra estratégia discutida na literatura, e primeiramente proposta por Beveridge e col.¹⁰, envolve cálculos quânticos ab initio. Nessa abordagem, a função potencial produzida é representativa da superfície de potencial mecânico-quântica, e consiste basicamente em determinar um conjunto de dados energia x distâncias atômicas que compõem um conjunto de pontos do espaço configuracional. Questões importantes como qual região do espaço configuracional amostrar em cálculos quânticos ab initio com moléculas isoladas, e, uma vez determinada esta região, quão bem esta região pode ser amostrada, consistem nos principais critérios para a determinação de potenciais intermoleculares usados em simulação em fase líquida. A região amostrada irá refletir a qualidade dos resultados obtidos em simulação. Em simulação de MC, as regiões do espaço configuracional preferencialmente amostradas são as regiões de mínimo de potencial, as quais correspondem ao equilíbrio termodinâmico. Portanto, tais regiões devem ser bem descritas nos cálculos quânticos. Isso implica em cálculos de alta qualidade, os quais incluem correlação eletrônica e um conjunto de funções de base estendido, especialmente para descrever a parte atrativa do potencial intermolecular. Isso demanda um alto custo computacional, sendo particularmente importante a escolha das geometrias dos complexos a serem estudados. Nessa direção, no início da década de 80, Jorgensen e col. desenvolveram potenciais para os líquidos amônia¹¹ e metanol¹². Os potenciais para a água, desenvolvidos por Clementi e col.^13,14 seguem essa mesma abordagem. Recentemente, nosso grupo desenvolveu nessa linha um potencial para a interação íon cloreto-metanol, com o objetivo de investigar a contribuição das interações íon-solvente no problema da solvatação iônica.

No presente estudo, adotamos o procedimento empírico, como também informações geradas a partir de cálculos mecânico-quânticos ab initio, para desenvolver os parâmetros de potencial de Lennard-Jones e Coulomb. O modelo desenvolvido conjuga esses dois procedimentos às simulações de MC, guardando o compromisso entre o custo computacional e a qualidade dos resultados em simulação.

As moléculas de CHCl₃ foram modeladas por estruturas rígidas compostas por cinco sítios de interação. Neste modelo, cada átomo foi representado por um sítio. A Figura 1 mostra o modelo molecular em termos dos seus átomos constituintes. As esferas de Lennard-Jones foram definidas de acordo com os raios s das interações de curto alcance de cada átomo. No presente modelo, o sítio H é representado apenas por uma carga parcial. A geometria adotada para as moléculas de clorofórmio está listada na Tabela 1.

Interações Eletrostáticas

Para fazer a simulação do líquido puro, a parte Coulombica foi representada por cargas parciais obtidas com o método CHELPG¹⁶, implementado no programa GAUSSIAN 92¹⁷. Esse método foi discutido em trabalho anterior, Seção 4. Parte I¹⁸. As cargas resultantes estão listadas na Tabela 1. Os potenciais eletrostáticos quânticos foram calculados em uma malha de pontos sobre camadas envolvendo a molécula CHCl₃. Usamos uma malha com densidade de 1 ponto/Ų em quatro camadas situadas a distâncias de 1,4, 1,6, 1,8 e 2,0 vezes o raio de van der Waals da molécula. Os raios de van der Waals usados para os átomo de C, H e Cl foram, respectivamente, 1,50 Å, 1,45 Å e 2,00 Å. Com esse procedimento para a seleção de pontos, 5024 pontos foram selecionados com um desvio quadrático médio RMS = 3,2 kcal mol^-1, e um desvio quadrático médio relativo RRMS = 46%. Os pontos foram escolhidos em regiões fora da superfície de van der Waals da molécula, uma vez que são essas as regiões acessíveis ao solvente, e portanto, particularmente importantes para descrever interações intermoleculares. O potencial eletrostático quântico foi calculado para cada um dos pontos das amostras a partir da função de onda calculada e da geometria experimental dada acima. Esses resultados foram utilizados como dados de entrada na subrotina de mínimos quadrados, com a condição de reproduzir a carga total da molécula, mas não o seu momento de dipolo. Utilizamos uma função de onda correlacionada ao nível de teoria de perturbação de Møller-Plesset de segunda ordem¹⁹ (MP2) com o conjunto de funções de base MC-311G(d) para os átomos de cloro²⁰ e 6-311G(d,p) para os átomos de carbono de hidrogênio²¹.

É interessante observar que as cargas CHELPG resultantes fornecem um momento de dipolo de 1,32 D para a molécula de CHCl₃, cerca de 20% acima do resultado experimental em fase gasosa² (1,07 D). Esse valor superestimado para o momento de dipolo é desejável para representar efeitos de polarização existentes no estado líquido^23-25. Como os cálculos dos potenciais moleculares quânticos são realizados em fase gasosa, eles não carregam informações sobre efeitos de polarização devido ao solvente. Uma melhoria substancial neste modelo seria usar cargas atômicas variáveis que fossem capazes de reproduzir a polarização da molécula devido a interação com o solvente.

Parâmetros de Lennard-Jones

O potencial de Lennard-Jones para o presente modelo é caracterizado por quatro parâmetros independentes: s_C e e_C para o átomo de carbono, e s_Cl e e_Cl para os átomos de cloro. Não usamos parâmetros de Lennard-Jones para o átomo de hidrogênio neste modelo. O parâmetro s_H para o átomo de H é cerca de 2 Ų⁵. Como adotamos s_C = 3,8 Å (Tabela 1) e o comprimento de ligação entre os átomos de C e H é 1,1 Å^26-28, o sítio H encontra-se interno ao raio de esfera dura do potencial de Lennard-Jones centrado sobre o sítio C.

Os parâmetros de Lennard-Jones do potencial OPLS para o clorofórmio de quatro sítios, desenvolvido por Jorgensen e col.⁸, foram usados como valores iniciais no processo de otimização. Os parâmetros do modelo OPLS para o clorofórmio estão listados na Tabela 1. A otimização consistiu de simulações interativas de MC. Os parâmtros de Lennard-Jones foram ajustados de modo a reproduzir resultados experimentais^29,30 de densidade, r, e entalpia de vaporização molar à pressão constante, DH_V, para o líquido à temperatura de 298 K e à pressão de 1 atm. Os parâmetros cruzados s_C-Cl e e_C-Cl estão relacionados aos parâmetros C - C e Cl - Cl por regras de cruzamento geométricas, conforme equação 31, Parte I¹⁸. Encontramos esses valores como sendo s_C-Cl e e_C-Cl = 0,15 kcal mol^-1. A Tabela 1 mostra os parâmetros obtidos neste modelo.

SIMULAÇÕES DE MONTE CARLO

Procedimento Computacional

As propriedades para o líquido puro foram calculadas usando a metodologia de simulações de MC com amostragem de Metroplolis e condições periódicas de contorno. O ensemble escolhido para o presente estudo foi o isobárico e isotérmico. A esolha do ensemble NpT foi motivada pela sua correspondência com as condições experimentais usuais de temperatura e pressão constantes. Os cálculos foram realizados usando o programa DIADORIM³¹. Uma caixa cúbica de simulação contendo N = 256 moléculas foi utilizada. A geometria das moléculas foi mantida rígida durante a simulação. O volume L³ inicial da caixa de simulação foi calculado usando a densidade experimental do líquido puro. Após especificado o comprimento L da caixa, as N moléculas foram distribuidas ao longo de uma rede cúbica de face centrada. Com essa definição de configuração inicial, cada nova configuração foi gerada transladando e rodando aleatoriamente uma molécula também escolhida aleatoriamente. A translação foi feita para o centro-de-massa nos três eixos cartezianos, conforme equações 15, 16, e 17, Parte I¹⁸. A rotação foi feita para a molécula inteira em torno de um dos eixos, conforme equações 18, e 19, Parte I¹⁸.

As translações e rotações foram permitidas dentro de certos limites máximos definidos pelo parâmetros Dr_max = ±0,23 Å e DQ_max = 0,15°, a fim de que a razão entre as configurações aceitas e as configurações geradas ficasse entre 40-50%. Como no ensemble NpT o volume pode variar, novas configurações também foram geradas por movimentos de volume, de acordo com a equação 24, Parte I¹⁸. As coordenadas do centro-de-massa das N moléculas foram escalonadas, como descrito na seção 2.3, Parte I¹⁸. Tentativas para mover o volume foram feitas a cada 10³ configurações dentro dos limites definidos por DV_max = ±680ų, de forma que fosse mantida a mesma razão acima para a aceitação de novas configurações.

Usamos a metodologia de raio de corte esférico para calcular a energia configuracional. A energia de interação clorofórmio-clorofórmio foi calculada sempre que a distância r_ij entre dois sítios de duas moléculas de clorofórmio dentro de um raio de corte definido por r_c = 15Å. A energia intermolecular para distâncias além do raio de corte foi corrigida somente para as interações de Lennard-Jones. Entretanto, é necessário ressaltar que a correção de longo alcance para o potencial de Lennard-Jones inclui também uma correção para a interação dipolo-dipolo: segundo Keesom³², a componente radial da interação dipolo-dipolo, obtida após a média angular para a orientação entre os dipolos, comporta-se assintoticamente como 1/r⁶. Logo, a contribuição da energia de interação entre dipolos além do raio de corte pode ser incluida na correção de longo alcance para o termo de Lennard-Jones de potência correspondente na componente radial. Com essa correção de longo alcance, a energia total para o presente sistema diminui de 1 - 2%. Podemos acrescentar também que o clorofórmio é um líquido com baixa constante dielétrica (e_r = 4,8), o que está associado a uma fraca correlação dipolo-dipolo.

Esse comportamento justifica as premissas introduzidas por Keesom para efetuar a média angular sobre a orientação dipolo-dipolo, e que a energia de interação comporta-se como uma função que decai radialmente com 1/r⁶. Logo, forças eletrostáticas de longo alcance têm pouca contribuição para as propriedades termodinâmicas e estruturais desse líquido. Como um resultado desses efeitos, podemos citar a grande volatilidade do líquido clorofórmio.

Na fase de equilíbrio 0,3 x 10⁷ configurações foram geradas. Os resultados obtidos nessa fase foram descartados. A seguir, 0,1 x 10⁸ configurações extras foram geradas para o cálculo das médias, completanto um total de 0,13 x 10⁸ configurações.

RESULTADOS E DISCUSSÕES

Termodinâmica

As propriedades termodinâmicas calculadas estão listadas na Tabela 2. Resultados teóricos e experimentais existentes na literatura também são mostrados nessa Tabela. Todas as propriedades calculadas, segundo o modelo proposto neste trabalho, concordam muito bem com os resultados experimentais e teóricos disponíveis.

A entalpia de vaporização foi calculada a partir da energia configuracional média para o líquido, E, usando-se a aproximação,

DH_V » -E + RT, (1)

onde RT=0,592 kcal mol^-1 a T=298K. Deve ser observado que a energia -E calculada inclui correções para as interações de Leonard-Jones além do raio de corte. O valor "experimental" de E na Tabela 2 foi estimado usando-se a relação acima com o valor experimental³⁰ para DH_V.

Os valores médios de densidade, r, e de entalpia de vaporização molar à pressão constante, DH_V, calculados mo presente modelo são aproximadamente 1% menores do que os correspondentes valores de MC obtidos por Jorgensen e col.⁸ com o modelo de quatro sítios para as moléculas de clorofórmio.

Comparados às observações experimentais, os resultados obtidos com o presente modelo para r e DH_V são maior e menor dessa mesma quantidade, respectivamente. O valor simulado para o volume molar, V_m, é também aproximadamente 1% maior do que o valor experimental.

O valor para a capacidade calorífica molar à pressão constante, C_p, que aparece na Tabela 2 já contém as contribuições vibracionais, rotacionais e translacionais obtidas para o monômeros em fase gasosa, a saber, 7,32 cal mol^-1K^-1, 2,98 cal mol^-1K^-1 e 4,97 cal mol^-1K^-1, respectivamente. Esses cálculos foram realizados usando-se o programa MOPAC com Hamiltoniano AM1³³. A contribuição média de para C_p configuracional obtida através do programa DIADORIM³¹ foi de 12,57 cal mol^-1K^-1 e adicionada aos valores acima para compor o valor apresentado na Tabela 2. O valor assim obtido para C_p está em bom acordo com o resultado experimental (Tabela 2). Os valores para capacidade calorífica molar, C_p, coeficiente de expansão isobárica, a_p e compressibilidade isotérmica, k_T, são calculados através das flutuações de entalpia e volume. As incertezas associadas a essas quantidades são maiores do que as incertezas associadas a propriedades como densidade e entalpia e requerem simulações mais longas^34-36. Isso pode ser observado nas Figuras 2 e 3 para C_p e r. As médias foram efetuadas no intervalo de 0,3 x 10⁷ a 0,13 x 10⁸ configurações. Pode-se observar que r já tem seu valor convergido a partir do início desse intervalo. O valor de C_p alcança a convergência por volta de 0,6 x 10⁷ configurações, ou seja, após 0,3 x 10⁸ configurações terem sido geradas na fase de médias. Comportamentos semelhantes foram observados para a_p e k_T neste trabalho e não estão mostrados aqui. Uma simulação preliminar de MC sob as condições do ensemble NpTa 298 K e 1 atm foi feita usando os parâmetros de potencial para o clorofórmio puro reportado por Laaksonen e col.⁷. Esses autores realizaram simulações de MD à volume constante com um modelo de cinco sítios para o clorofórmio. Os resultados de MC obtidos neste trabalho para r e DH_V no ensemble NpT usando o modelo proposto por esses autores foram, respectivamente, 1,493±0,002 g cm^-3 e 9,09±0,01 kcal mol^-1. Esse último resultado é aproximadamente 20% maior do que o valor experimental. Devido a essa discrepância, neste trabalho é desenvolvido um modelo de cinco sítios para o clorofórmio mais apropriado para estudos de propriedades termodinâmicas. Vale salientar que o modelo de cinco sítios utilizado por Laaksonen e col. é essencialmente o mesmo reportado por Evans³, e Dietz e Heizinger⁴, exceção feita para as cargas. Entretanto, a contribuição das interações eletrostáticas para a energia intermolecular calculada com o presente modelo é apenas cerca de 4%. Este resultado está em bom acordo com os valores de contribuição Coulombica reportados por Laaksonen e col. (cerca de 6%), e Böhm e Ahlrichs⁶ (cerca de 5%).

Estrutura

Seis funções de distribuição radial de pares foram construídas a partir do potencial proposto: g_C-C, g_C-H, g_H-H, g_C-Cl, g_H-Cl, g_Cl-Cl. Essas funções são apresentadas nas Figuras 4 e 5. As curvas para os números de vizinhos como uma função da distância de separação entre os sítios são apresentadas na Figura 6.

Informações sobre a disposição média das moléculas de clorofórmio podem ser obtidas através da função g_C-C (linha sólida na Figura 4), uma vez que o átomo de carbono ocupa aproximadamente o centro da molécula. As posições dos dois picos que aparecem nessa função correspondem aos raios médios das camadas de moléculas que se formam em torno de uma molécula no líquido. Observamos que o primeiro pico é alargado, com um máximo em r = 5,3 Å. O segundo pico, um pouco mais alargado, aparece com um máximo aproximadamente em r = 10 Å. O raio médio da segunda camada sendo um pouco menor do que duas vezes o raio médio da primeira camada de moléculas. Esses resultados foram os mesmos encontrados por Dietz e Heinzinger⁴, e por Laaksonen e col.⁷. As larguras dos picos indicam como variam as distâncias médias entre um átomo de carbono e os átomos de carbono das moléculas nas camadas vizinhas, devido ao empacotamento das moléculas. A área sob os picos fornece o número médio de moléculas nas camadas, ou números de coordenação. Portanto, a integração de r²g_C-C para o primeiro pico, de r = 0 até o primeiro mínimo correspondente em r = 7,6 Å fornece, em média, aproximadamente 13 moléculas de clorofórmio na primeira camada, como pode ser visto pela Figura 6. Esse resultado concorda com o resultado de 12,6 moléculas obtido por Dietz e Heinzinger⁴ com simulações MD.

Características sobre orientações preferenciais entre as moléculas no líquido podem ser inferidas, examinado-se a função g_H-Cl. O pico pouco intenso que aparece nessa função (linha pontilhada na Figura 5) mostra que a distância intermolecular média mais próxima no líquido ocorre em r = 3,2 Å entre os átomos de hidrogênio de uma molécula de clorofórmio e os átomos de cloro da molécula vizinha, como também encontrado por Laaksonen e col⁷. A área sob essa curva até o primeiro mínimo, em torno de 3,7 Å, mostra que, em média, dois átomos de cloro são primeiros vizinhos ao átomo de hidrogênio, como pode ser observado pela Figura 6. Um cálculo mecânico-quântico ab initio ao nível MP2/MC-311G(d,p) mostra que o vetor momento de dipolo para a molécula de clorofórmio, na geometria experimental, está preferencialmente orientado na direção da ligação CH. Esses resultados indicam duas características importantes desse líquido. Primeiramente, o alinhamento parcial dos dipolos entre moléculas vizinhas. Como segundo resultado importante, consideramos que, em média, o átomo de hidrogênio de uma molécula está localizado entre dois átomos de cloro de uma mesma molécula vizinha. A Figura 7 ilustra esse comportamento. Nessa Figura, mostramos em destaque a molécula de clorofórmio mais próxima do centro da caixa de simulação coordenando uma outra molécula de clorofórmio, mais próxima a ela. Essa Figura foi retirada de uma configuração composta das 256 moléculas de clorofórmio da caixa de simulação em um instante arbitrário do líquido. As conclusões acima concordam com os resultados obtidos por Bertagnolli³⁷ a partir da análise de funções de distribuição radial de pares. Esse autor encontrou que os eixos dos dipolos moleculares estão inclinados entre 0^o e 50^o e deslocados de modo que o átomo de hidrogênio está direcionado para dois átomos de cloro. Bertagnolli e col. investigaram a estrutura do líquido usando técnicas combinadas de difração de raios-X e espalhamento de nêutrons^27,28,37-40. Como pode ser observado pelo crescimento lento dos números de coordenação para as funções g_C-H e g_H-H, arranjos entre os átomos H - H e C - H são quase inexistentes a curtas distâncias. A função g_C-H apresenta um ombro aproximadamente em 4,3 Å, e um pico alargado aproximadamente em 6 Å. A integração da curva r²g_C-H em torno desse pico, até a região de mínimo correspondente, mostra, em média, 14 moléculas vizinhas. Esse resultado está em acordo com o número de primeiros vizinhos encontrado como 13 para a correlação C - C.

A função g_Cl-Cl apresenta um pico em 3,8 Å com um mínimo aproximadamente em 5 Å, e um segundo pico em 6,2 Å. Esses resultados para as posições dos máximos estão em bom acordo com os correspondentes valores de 3,75 Å e 6,2 Å reportados por Dietz e Heinzinger. Esses resultados estão também em bom acordo com os correspondentes valores de 3,79 Å e 4,75 Å para as posições dos primeiros máximo e mínimo reportados para a correlação Cl - Cl no líquido CCl₄ em simulação de MD⁴¹.

De um modo geral, os resultados obtidos com o presente modelo para a estrutura do líquido CHCl₃ são similares aos resultados de MD existentes na literatura. Pequenas diferenças são observadas para as funções g_C-H e g_H-Cl obtidas neste trabalho (Figuras 4 e 5) e aquelas obtidas por MD. Na posição radial 4,4 Å para a função g_C-H encontramos um ombro, enquanto os resultados de MD indicam um pico definido, ainda que pequeno. Na posição radial para a função g_H-Cl, correpondente à distância mais curta de aproximação entre duas moléculas de CHCl₃ no líquido, encontramos um pico menos acentuado comparado àqueles obtidos por MD. Ainda, a distância média de separação entre os sítios H e Cl encontrada por Dietz e Heinzinger é levemente menor do que 3 Å. O valor de densidade adotado por Dietz e Heinzinger foi de 1,484 g cm^-3 para uma faixa de temperatura de 295-299 K. Pelas observações acima, o modelo desenvolvido neste trabalho parece indicar um líquido menos estruturado quando comparado aos resultados existentes na literatura para os modelos de cinco sítios. Vale lembrar que dos modelos existentes na literatura, somente o modelo apresentado por Laaksonsen e col. foi testado neste trabalho através de simulações de MC no ensemble NpT, e que entre os modelos existentes há diferença apenas nas cargas dos sítios. Laaksonsen e col. relizaram simulações de MD no ensemble NpT com uma densidade de 1,40 g cm^-3 para uma temperatura de 288 ± 15 K. Esse valor de densidade é aproximadamente 7% menor do que o resultado encontrado neste trabalho (1,493 g cm^-3, apresentado na seção anterior) em simulações de MC no ensemble NpT com a mesma função potencial reportada por esses autores e para a temperatura de 298 K. Uma vez que flutuações de volume não são permitidas com o uso do ensemble N_VT, o sistema não pode relaxar até encontrar a densidade apropriada para a temperatura em questão. Com isso, são esperadas algumas diferenças nas funções distribuição radial obtidas através das simulações N_VT e NpT.

SIMULAÇÕES DE ENERGIA LIVRE

Energia Livre de Solvatação

O termo solvatação está associado à interação entre um soluto e um solvente. Seguindo Ben-Naim⁴², definimos a solvatação de uma molécula de soluto por um solvente como o processo de transferir a molécula de soluto da fase de gás ideal para o interior de um líquido à temperatura T e à pressão p constantes. Como um teste adicional para o modelo proposto, nesta seção calculamos a energia livre de Gibbs de solvatação de uma molécula de clorofórmio pelo seu próprio líquido a temperatura 298 K e pressão de 1 atm.

Procedimento Computacional

A teoria de perturbação termodinâmica com a amostragem preferencial foi implementada no programa DIADORIM³¹ e utilizada no presente estudo. Essa metodologia tem sido correntemente empregada em estudos que envolvem energia livre^23,43-45. A molécula de clorofórmio mais próxima do centro da caixa cúbica de simulação, previamente equilibrada e contendo 256 moléculas de clorofórmio, foi considerada como o soluto para compor a configuração inicial. A estratégia adotada para calcular a diferença de energia livre de solvatação foi de aniquilar a molécula de soluto na solução. Adotamos para isso um procedimento denominado electrostatic decoupling^25,46. Esse procedimento divide a perturbação em dois estágios. Primeiramente, as cargas sobre os sítios são removidas, mantendo-se os parâmetros de Lennard-Jones fixos em seus valores iniciais. Nessa etapa, que chamamos de conversão de A ® A', o conjunto de parâmetros a ser perturbado é o conjunto de cargas {q_j} para a molécula constituida de j átomos. O conjunto de cargas inicial está listado na Tabela 1, juntamente com os parâmetros de Lennard-Jones. Posteriormente, o processo é completado com o desaparecimento dos parâmetros de Lennard-Jones. Denominamos essa etapa de conversão de A' ® B.

Para os estudos realizados nessa seção, e na seção "Efeitos do solvente..." a seguir, cada simulação teve um período de equilíbrio em que 0,4 x 10⁷ configurações foram geradas e descartadas. As médias foram obtidas com 0,6 x 10⁷ configurações.

Resultados e Discussões

A diferença de energia livre total, DG_total, obtida com o procedimento descrito acima, com sinal oposto, pode ser diretamente comparada com o resultado experimental para a energia livre de solvatação do clorofórmio em seu próprio líquido. Para a conversão de A ® A', foram necessários somente dois passos, resultando em dois incrementos de energia livre. Para a conversão de A' ® B, vinte simulações foram necessárias, resultando em vinte incrementos de energia livre. A mudança total de energia foi obtida somando-se cada mudança incremental de energia. Os resultados parciais para cada etapa, bem como DG_total, estão mostrados na Tabela 3.

Observamos que a mudança na energia livre associada às forças Coulombicas (transformação A ® A') é pequena comparada à mudança na energia livre associada às forças de curto alcance (transformação A' ® B), as quais contribuem com 92% da mudança total de energia no processo. Esse resultado está de acordo com o resultado já obtido da contribuição Coulombica para a energia intermolecular. Observamos, também, que o número de perturbações necessárias para desparecer com a molécula no interior do líquido (transformação A' ® B) foi superior àquela necessária para "descarregar" a mesma (transformação A ® A'). As perturbações nos parâmetros de Lennard-Jones, que são os parâmetros associados ao desaparecimento da molécula de soluto, envolvem pequenas variações nas distâncias devido aos rearranjos das moléculas de solvente nas imediações do soluto. Uma vez que são as interações de curto alcance que variam rapidamente com as distâncias intermoleculares, e não as interações Coulombicas, evidentemente que o sistema deve ser lentamente perturbado ao longo da transformação A' ® B.

A escolha de l foi parcialmente indicada pela incerteza estastística calculada para um incremento de energia livre. A escolha Dl = 0,05 resultou em uma incerteza de cerca de 0,1 kcal mol^-1para cada valor individual de DG e uma incerteza estimada de 0,51 kcal mol^-1 para o processo total, englobando as duas etapas. Como pode ser visto na Tabela 3, o valor calculado para a energia livre de solvatação, DG_total, foi de - 4,32 kcal mol^-1. Este resultado está em bom acordo com o valor experimental⁴⁷ de - 4,182 kcal mol^-1.

O procedimento adotado de mudar primeiramente a parte eletrostática do potencial de interação soluto-solvente, mantendo-se inalterados os parâmetros de Lennard-Jones do soluto, foi motivado por alguns resultados inicialmente obtidos nas simulações. Como o átomo de hidrogênio da molécula de soluto foi representado unicamente por uma carga parcial positiva, este átomo deixa de experimentar repulsão de van der Waals, devido à aproximação com os átomos de cloro do solvente, à medida que o parâmetro s do átomo de carbono do soluto é reduzido no processo de aniquilação. Nessa situação, as mudanças de energia livre resultam em valores excessivamente grandes.

Na conversão A' ® B, ou seja, a partir da remoção das cargas do soluto até sua aniquilação na solução, a amostragem adotada foi aquela de passo único, e não de passo duplo. Na amostragem de passo único, os estados intermediários são sempre posteriores ao estado de referência. Na amostragem de passo duplo, os estados intermediários são os estados imediatamente anterior e posterior em relação a este. As simulações dos estados intermediários anterior e posterior em relação ao estado de referência são realizadas utilizando-se as configurações representativas do estado de referência. Significa dizer que o solvente não relaxa quando o soluto sofre as modificações características dos estados anterior e posterior. Com relação aos valores de s, temos que s_i-1 > s_i > s_i+1. Logo, no estado anterior, o soluto é mais volumoso que no estado de referência. Os parâmetros s sobre os átomos do soluto são maiores nesse estado. Uma vez que o solvente está equilibrado para a situação em que o soluto é menos volumoso (caracterizado por s_i) que no estado anterior (caracterizado por s_i-1), é esperada superposição artificial dos átomos de soluto e átomos de solvente na fase anterior, resultando em energias configuracionais altamente repulsivas. Uma amostragam de passo único foi, então, usada para contornar esse problema, embora tenha elevado o custo computacional.

Efeitos do Solvente sobre o Equilíbrio Conformacional

O fato de um meio líquido exercer influência sobre a conformação de um soluto dissolvido é bem conhecido na literatura^1,48-50. O enrolamento de proteínas, e os efeitos de hidratação sobre as estruturas do DNA, são exemplos dessas manifestações⁵¹. Para levar em conta os efeitos do solvente sobre a estabilidade conformacional em estudos teóricos via simulação computacional, as quantidades de interesse são a coordenada conformacional, F, e o potencial de torção média, W(F), que é definido como W(F) = V^gas(F) + w^liq(F). O termo V^gas(F) é a energia intermolecular em fase gasosa, e o termo w^liq(F) descreve os efeitos induzidos pelo solvente. Podemos dizer que w^liq(F) é o trabalho induzido pelo solvente, necessário para levar a coordenada conformacional de um estado de referência F_i para um estado F_j. A energia induzida pelo solvente sobre o soluto é um resultado das interações soluto-solvente e solvente-solvente. Essa contribuição é genericamente denominada de efeitos do solvente, e pode ser medida pela variação de energia livre devido às interações soluto-solvente e solvente-solvente, DDG_solv.

Como solventes foram utilizados o meio aquoso e um meio de constante dielétrica baixa, representado por moléculas de clorofórmio. O modelo TIP4P para as moléculas de água⁵², e o modelo de cinco sítios para o clorofórmio proposto neste trabalho, foram adotados neste estudo. Os parâmetros de potencial e geometria para a água estão listados na Tabela 4. Como o principal objetivo deste estudo é investigar a contribuição do solvente para a energia livre de rotação da molécula hipotética nesses dois meios, a contribuição intramolecular, V^gas(F), foi desprezada.

Procedimento Computacional

Utilizamos a metodologia de teoria de perturbação termodinâmica com amostragem preferencial implementada no programa DIADORIM para investigar quantitativamente a influência do solvente sobre a barreira de energia livre para a conversão cis/trans de uma molécula hipotética com apenas um grau de liberdade conformacional. Preparamos, incialmente, duas soluções diluídas para as simulações. Para isso, usamos dois sistemas: o sistema I foi composto de 255 moléculas de água TIP4P, encerradas em uma célula cúbica a 298 K e pressão de 1 atm, onde o volume da célula pode variar (ensemble NpT); e o sistema II foi composto de 256 moléculas de clorofórmio sob as mesmas condições de temperatura e pressão que o sistema I. As configurações iniciais para as soluções diluídas foram obtidas a partir dos sistemas I e II já preparados (equilibrados) acima. No sistema I, três moléculas de água mais próximas do centro da caixa foram substituidas pela molécula de soluto hipotética; e no sistema II, duas moléculas de clorofórmio mais próximas do centro da caixa de simulação foram recolocadas pelo soluto. Imaginamos as soluções suficientemente diluídas de modo que apenas uma molécula de soluto é considerada. A geometria adotada para o soluto hipotético foi baseada em uma cadeia de quatro sítios com comprimento de ligação de 1,45 Å, e ângulo de ligação de 120^o. Os parâmetros de potencial OPLS para o metano³⁴ foram adotados para cada sítio da cadeia hipotética. Esses parâmetros estão listados na Tabela 4. Para cada sítio foram adotadas cargas parciais de 0,5 |e|. Cargas com sinais negativos foram usadas nas extremidades da cadeia para manter coerência com a distribuição de carga em uma cadeia finita, e cargas positivas para os sítios centrais. Uma representação esquemática do soluto hipotético nas suas conformações cis/trans é apresentada na Figura 8.

As configurações seguintes foram geradas movendo o volume, as moléculas de solvente e a molécula de soluto. Tentativas para mover o soluto, tanto na solução soluto-água quanto na solução soluto-clorofórmio, foram feitas a cada 100 configurações. Os limites máximos permitidos para a translação e rotação do soluto nas duas soluções foram, respectivamente, ±0,12 Å e ±12^o. Para a solução soluto-água, foram tentados movimentos de volume a cada 10³ configurações, e os parâmetros DV_max, Dr_max e DQ_max assumiram, respectivamente, os valores ±130 ų, ±0,15 Å e ±15^o. Para a solução soluto-clorofórmio, o volume variou segundo os limites ±580 ų, e as moléculas de clorofórmio ficaram sujeitas às mesmas condições para rotação e translação que no líquido clorofórmio puro.

O procedimento de corte esférico foi utilizado para ambas as soluções. Na solução soluto-água, as interações água-água e soluto-água foram calculadas até a distância r_c = 8,5 Å. Para as interações soluto-clorofórmio e clorofórmio-clorofórmio na solução soluto-clorofórmio, o raio de corte adotado foi o mesmo que para a interação clorofórmio-clorofórmio no líquido clorofórmio puro (15 Å). As simulações foram realizadas quando o soluto variou sua conformação de cis para trans. Ou, equivalentemente, podemos dizer que na coordenada conformacional F, o soluto variou de um ângulo diedro F = 0^o para F = 180^o. Para converter o soluto da forma de cis para a forma trans, usamos o parâmetro de acoplamento l variando de l = 0 (para F = 0) a l = 1 (para F = 180^o). A conversão foi realizada em incrementos de DF = 9^o a partir da conformação cis, nos sentidos contrário e a favor dos ponteiros do relógio, simultaneamente. Esse procedimento foi necessário para minimizar as perturbações nos solventes devido à rotação da molécula de soluto. Como medida de economia de tempo computacional, adotamos a amostragem de passo duplo. Com a amostragem dupla, as simulações foram baseadas nos sistemas soluto-água e soluto-clorofórmio de referência, determinados pela coordenada conformacional F_i do soluto. A partir desses sistemas de referência, a conformação do soluto foi perturbada em ambas as direções, anterior e posterior, ou de F_i para F_j = F_i + DF e de F_i para F_j = F_i - DF. Computacionalmente, três solutos foram monitorados simultaneamente: o soluto de referência com a coordenada conformacional F_i, e os dois solutos perturbados com F_j = F_i ± DF. Essa estratégia permitiu cobrir duas vezes DF ou 18^o, em uma única simulação. Com isso, as simulações foram equilibradas para F_i = 9^o, 27^o, 45^o, 63^o, 81^o, 99^o, 117^o, 135^o, 153^o e 171^o, perfazendo um total de 10 simulações e 20 incrementos de energia livre relativa, DDG_solv, para cada um dos sistemas, soluto-água e soluto-clorofórmio. Esses resultados estão apresentados na Tabela 3.

RESULTADOS E DISCUSSÕES

Dentre os líquidos orgânicos correntemente estudados, o clorofórmio constitui um meio modelo para processos físico-químicos que ocorrem em um meio líquido de baixa constante dielétrica (e_r = 4,8), ao contrário da água que é comumente representativa de um meio de alta constante dielétrica (e_r = 78,3). Os solventes água e clorofórmio são considerados meios modelo para estudos das mudanças conformacionais assumidas pelas moléculas de soluto nos meios de baixa constante dielétrica (lipídico) e aquoso⁸.

A abordagem tradicional para efeitos de solvente sobre o equilíbrio conformacional é baseada na teoria de campo de reação de Onsager para dipolos moleculares em fase condensada^1,48. No modelo de Onsager, a molécula de soluto é tratada como um dipolo no interior de uma cavidade esférica embebida em um solvente modelado por um dielétrico contínuo, sem estrutura. O potencial de reação produzido pelo meio, em resposta à polarização induzida pelo dipolo, depende da constante dielétrica do meio, e_r, e é usado como uma estimativa para a energia de solvatação via equação de Onsager,

Onde µ é o momento de dipolo do soluto, e r o raio da cavidade. Para duas espécies em equilíbrio, A e B, consideradas de mesmo tamanho, a diferença de energia de solvatação é dada por,

Qualitativamente, a expressão acima diz que o confôrmero com momento de dipolo maior será preferencialmente estabilizado em um meio de constante dielétrica mais alta.

Os resultados obtidos via simulação para as diferenças de energia livre de solvatação cumulativas, relativas a 180^o, DG_solv, estão apresentados na Tabela 5, e na Figura 9. Nessa Figura, c representa clorofórmio, e w representa água.

Esses resultados mostram que para os dois solventes, o confôrmero gauche (F próximo de 60^o) é melhor estabilizado por solvatação. A energia livre de solvatação decresce do confôrmero trans para gauche, e volta a crescer até o confôrmero cis. Esse comportamento é verificado para os dois solventes. O confôrmero gauche tem DG_solv mais negativo em relação ao confôrmero trans por 6,3 kcal mol^-1 em clorofórmio, e por 11,4 kcal mol^-1 em água TIP4P. O confôrmero gauche é também melhor solvatado em relação ao confôrmero cis por 4,1 kcal mol^-1 em clorofórmio, e por 5,7 kcal mol^-1 em água TIP4P. Observa-se desses dados que a solvatação do confôrmero gauche é mais forte em água TIP4P do que em clorofórmio. Os resultados para as diferenças de energias de interação soluto-água e soluto-clorofórmio estão também apresentados na Figura 9. Observa-se que a diferença de energia de interação soluto-solvente entre a conformação gauche relativa à conformação cis do soluto é de 7,1 kcal mol^-1 em água TIP4P, e -5,4 kcal mol^-1 em clorofórmio. Esses resultados indicam que uma interação preferencial das moléculas de solvente com o confôrmero gauche é responsável pela estabilização da molécula de soluto. Essa conclusão é reforçada pelos resultados obtidos para as energias de interação solvente-solvente apresentados na Tabela 6. Nessa Tabela, observa-se que a interação água-água na solução soluto-água, bem como a interação clorofórmio-clorofórmio na solução soluto-clorofórmio, não mudam significativamente quando o soluto muda sua conformação de cis para trans.

O efeito dos dois meios sobre a barreira de energia para a rotação cis/trans do soluto foi medido pela energia livre de solvatação para a conversão cis/trans, DG_solv(trans®cis ). Como pode ser observado pela Tabela 5 e pela Figura 9, o resultado obtido para essa conversão em água TIP4P foi de -5,6 kcal mol^-1, e em clorofórmio foi de -2,2 kcal mol^-1. Portanto, o valor obtido para a barreira de energia livre para a rotação interna da molécula da forma cis para a forma trans do soluto é maior em água TIP4P do que em clorofórmio, que significa que a conversão cis/trans ocorre com menor custo energético em clorofórmio do que em água, implicando em maior estabilização da conformação cis pela água. Esse resultado está de acordo com a teoria de contínuo Onsager, que diz que a solvatação do confôrmero de maior momento de dipolo, o confôrmero cis, é mais favorecida em um meio de maior constante dielétrica, no caso, a água. É importante observar na Figura 9 que ao longo da interconversão cis/trans, a variação obtida para a energia mecânica é maior que a calculada para a energia livre. A variação de energia livre pode ser escrita na forma DG = DH - TDS, a qual pode ser aproximada para DG = DE - TDS, onde DE é a variação da energia mecânica. Desta equação, e com os dados obtidos pela Figura 9, a variação de entropia pode ser obtida. Obtivemos TDS = 7,36 kcal mol^-1 para a interconversão no meio aquoso, e TDS = 2,88 kcal mol^-1 para a interconversão em clorofórmio. Esses dados mostram o aumento de entropia ao longo do processo de interconversão. Essa perda de ordem deve-se ao fato das moléculas de solvente estarem mais organizadas devido à interação com o momento de dipolo do soluto na conformação cis. É também interessante observar que a variação de entropia é maior no meio aquoso.

CONCLUSÕES

Um modelo de cinco sítios para representar a interação entre moléculas de clorofórmio foi desenvolvido. Esse modelo utiliza um potencial efetivo de pares da forma Lennard-Jones mais Coulomb para calcular a interação intermolecular. As cargas parciais necessárias para representar as interações eletrostáticas foram derivadas do potencial eletrostático molecular quântico calculado para a molécula de clorofórmio. É importante ressaltar que os sítios de carga foram definidos em posições coincidentes com os núcleos dos átomos constituintes da molécula de clorofórmio. Uma função de onda correlacionada ao nível de teoria de pertubação de Møller-Plesset de segunda ordem, com o conjunto de funções de base MC-311G(d,p). foi utilizada para gerar o potencial eletrostático. Os parâmetros de Lennard-Jones foram otimizados por meio de simulações interativas de Monte Carlo para reproduzir dados termodinâmicos experimentais de densidade e entalpia molar de vaporização para o líquido puro. Essas simulações foram realizadas no ensemble isobárico-isotérmico a 298 K e 1 atm. Utilizamos regras de cruzamento para obter os parâmetros de Lennard-Jones para a interação entre átomos de carbono e cloro. A partir desse modelo, obtivemos resultados termodinâmicos para o líquido puro em muito boa concordância com resultados teóricos e experimentais disponíveis na literatura. Os resultados estruturais obtidos indicam pouca organização espacial entre as moléculas de clorofórmio no líquido, e globalmente estão de acordo com resultados da literatura.

Como um teste adicional para o modelo proposto neste trabalho, calculamos a energia livre de solvatação de uma molécula de clorofórmio em seu próprio líquido, mantendo-se a temperatura e a pressão constantes. A metodologia de teoria de perturbação termodinâmica, dicutida em trabalho anterior, foi utilizada nesses estudos. O resultado obtido concorda muito bem com o resultado experimental existente na literatura. Para obter-se uma boa estimativa para a energia livre de solvatação, apontamos como principais dificuldades a determinação de quais, e quantos, estados de referência e estados intermediários devem ser utilizados para descrever a transformação, segundo as equações 40 e 41 apresentadas na Parte I¹⁸.

O líquido clorofórmio é utilizado como um meio protótipo para representar a dissolução de moléculas em meio baixa constante dielétrica. Dada a importância da investigação de efeitos do meio sobre a estrutura molecular, fizemos um estudo para prever mudanças conformacionais de solutos orgânicos nos meios aquoso e de baixa constante dielétrica. O modelo de água TIP4P, desenvolvido por Jorgensen e col., foi utilizado para representar o meio aquoso. Utilizando a teoria de perturbação termodinâmica, investigamos a influência do meio sobre a barreira de energia para uma rotação cis/trans de uma molécula hipotética. Resultados qualitativos esperados pela teoria de contínuo de Onsager para a solvatação de dipolos moleculares em um meio definido por uma constante dielétrica foram confirmados a partir desses estudos. Os dados obtidos mostram também a reorganização das moléculas do solvente devido à influência do momento de dipolo do soluto. Os valores obtidos para a variação de entropia mostram que essa reorganização molecular é mais intensa no solvente aquoso.

AGRADECIMENTOS

À FAPESP e ao CNPq pelo apoio financeiro. Este trabalho teve suporte computacional parcial do CENAPAD-SP.

e-mail: [e-mail]

www.lqt.dq.ufscar.br

• 1. Reichardt, C.; Solvents and Solvent Effects in Organic Chemistry, Verlag Chemie, Weinheim, 2^nd Edition 1988.
• 2. Evans, M. W.; Ed.; Dynamical Processes in Condensed Matter, Willey, New York 1985.
• 3. Evans, M. W.; J. Mol. Liquids 1983, 25, 211.
• 4. Dietz, W.; Heinzinger, K.; Ber. Bunsenges. Phys. Chem. 1985, 89, 968.
• 5. Dietz, W.; Heinzinger, K.; Ber. Bunsenges. Phys. Chem. 1984, 88, 543.
• 6. Böhm, H. J.; Ahlrichs, R.; Mol. Phys. 1985, 54, 1261.
• 7. Kovaks, H.; Kowalewski, J.; Laaksonen, A.; J. Chem. Phys. 1990, 94, 7378.
• 8. Jorgensen, W. L.; Briggs, J. M.; Contreras, M. L.; J. Chem. Phys. 1990, 94, 1683.
• 9. Politi, J. R. S., Dissertaçăo de Mestrado, Departamento de Química, Universidade Federal de Săo Carlos, S. Carlos, 1994.
• 10. Swaminathan, S.; Whitehead, R. J.; Guth, E.; Beveridge, D. L.; J. Am. Chem. Soc. 1977, 99, 7817.
• 11. Jorgensen, W. L.; Ibrahim, M.; J. Am. Chem. Soc. 1980, 102, 3309.
• 12. Jorgensen, W. L.; J. Am. Chem. Soc. 1980, 102, 543.
• 13. Lie, G. C.; Clementi, E.; Yoshimine, M.; J. Chem. Phys. 1976, 64, 2314.
• 14. Lie, G. C.; Clementi, E.; Phys. Rev. A 1986, 33, 2679.
• 15. Silva, L. B., Dissertaçăo de Mestrado, Instituto de Química, Universidade Estadual de Campinas, Campinas, 1996
• 16. Breneman, C. M.; Wiberg, K. B.; J. Comp. Chem. 1990, 11, 361.
• 17. Frisch, M. J.; Trucks, G. W.; Head-Gordon, M.; Gill, P. M. W.; Wong, M. W.; Foresman, J. B.; Johnson, B. G.; Schlegel, H. B.; Robb, M. A.; Repogle, E. S.; Gomperts, R.; Andres, J. L.; Ragharachari, K.; Binkley, J. S.; Gonzalez, C.; Martin, R. L.; Fox, D. J.; Defrees, D. J.; Baker, J.; Stewart, J. J. P.; Pople, J. A. 1992, Gaussian 92, Revision C, Gaussian, Inc.; Pittsburg, PA.
• 18. Barlette, V. E.; Freitas, L. C. G.; Quim. Nova 1999, 22, 254.
• 19. Hehre, W. J.; Radom, L.; Schleyer, P. V.; Pople, J. A.; Ab Initio Molecular Orbital Theory, John Wiley & Sons, New York 1986.
• 20. McLean, A. D.; Chandler, G. S.; J. Chem. Phys. 1980, 72, 5639.
• 21. Krishnan, R.; Binkley, J. S.; Seeger, R.; Pople, J. A.; J. Chem. Phys. 1980, 72, 650.
• 22. McClellan, A. L.; Tables of Experimental Dipole Moments, Freeman, London 1963.
• 23. Kollman, P.; Chem. Rev. 1993, 93, 2395.
• 24. Merz Jr., K. M.; J. Comp. Chem. 1992, 13, 749.
• 25. Carlson, H. A.; Nguyen, T. B.; Orozco, M.; Jorgensen, W. L.; J. Comp. Chem. 1993, 14, 1240.
• 26. Jen, M.; Lide Jr.; D. R.; J. Chem. Phys. 1962, 36, 2525.
• 27. Bertagnolli, H.; Leicht, D. O.; Zeidler, M. D.; Mol. Phys. 1978, 35, 193.
• 28. Bertagnolli, H.; Leicht, D. O.; Zeidler, M. D.; Chieux, P.; Mol. Phys. 1978, 36, 1769.
• 29. Riddick, J. A.; Bunger, W. B.; Organic Solvents, Wiley, New York 1970.
• 30. Majer, V.; Sváb, L.; Svoboda, V.; J. Chem. Thermodynamics 1980, 12, 843.
• 31. Freitas, L. C. G.; Diadorim Program, Version 2.0, Departamento de Química, Universidade Federal de Săo Carlos, 1995.
• 32. Israelachvili, J. N.; Intermolecular and Surface Forces, Academic Press, London 1985.
• 33. Stewart, J. J. P., Mopac 6.0, United States Air Force Academy, F. J. Seiler Research Laboratory, 1990.
• 34. Jorgensen, W. L.; Madura, J. D.; Swenson, C. J.; J. Am. Chem. Soc. 1984, 106, 6638.
• 35. Jorgensen, W. L.; J. Chem. Phys. 1986, 90, 1276.
• 36. Jorgensen, W. L.; Briggs, J. M.; Mol Phys 1988, 63, 547.
• 37. Bertagnolli, H.; Ber. Bunsenges. Phys. Chem. 1981, 85, 644.
• 38. Bertagnolli, H.; Leicht, D. O.; Zeidler, M. D.; Chieux, P.; Mol. Phys. 1978, 35, 199.
• 39. Bertagnolli, H.; Chieux, P.; Ber. Bunsenges. Phys. Chem. 1980, 84, 1225.
• 40. Bertagnolli, H.; Chieux, P.; Mol. Phys. 1984, 51, 617.
• 41. McDonald, I. R.; Bounds, D. G.; Klein, M. L.; Mol. Phys. 1982, 45, 521.
• 42. Ben-Naim, A.; Solvation Thermodynamics, Plenum Press, New York 1987.
• 43. Gunsteren, W. F. V.; Weiner, P. K.; Eds.; Computer Simulation of Biomolecular Systems: Theoretical and Experimental Applications, ESCOM, Leiden 1989.
• 44. Beveridge, D. L.; DiCapua, F. M.; Ann. Rev. Biophys. Biophys. Chem. 1989, 18, 431.
• 45. Reynolds, C. A.; King, P. M. 1989, in Richards, W. G.; Ed.; Computer-Aided Molecular Design, Cap. 4, p. 51, IBC Technical Services Ltd, Great Britain.
• 46. Bash, P. A.; Singh, U. C.; Langridge, R.; Kollman, P. A.; Science 1987, 236, 564.
• 47. Ben-Naim, A; Marcus, Y.; J. Chem. Phys. 1984, 81, 2016.
• 48. Abraham, R. J.; Bretschneider, E. 1974, in Orville-Thomas, W. J.; Ed., Internal Rotation in Molecules, Cap. 13, Wiley, New York.
• 49. Jorgensen, W. L.; J. Phys. Chem. 1983, 87, 5304.
• 50. Jorgensen, W. L.; Buckner, J. K.; J. Phys. Chem. 1987, 91, 6083.
• 51. Vázquez. M.; Némethy, G.; Scheraga, H. A.; Chem. Rev. 1994, 94, 2183.
• 52. Jorgensen, W. L.; Chandrasenkhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L.; J. Phys. Chem. 1983, 79, 926.
• 53. Marcus, Y.; Ion Solvation, Wiley, New York 1985.
• 54. Harmony, M. D.; Laurie, V. W.; Kuczkowsky, R. L.; Schwendeman, R. H.; Ramsay, D. A.; Phys. Chem. Ref. Data 1978, 8

Datas de Publicação

• Publicação nesta coleção
  08 Fev 2000
• Data do Fascículo
  Jul 1999

Histórico

• Aceito
  14 Abr 1998
• Recebido
  19 Fev 1997