Hydration Free Energies of Linear Alkanes: Systematic Deviations in Common Water Models and Their Correction

Resumo Técnico: Energias Livres de Hidratação de Alcanos Lineares: Desvios Sistemáticos em Modelos de Água Comuns e Sua Correção

Declaração do Problema
Campos de força moleculares comuns superestimam sistematicamente as energias livres de hidratação ($\Delta G_{hyd}$) de solutos hidrofóbicos, levando a um efeito hidrofóbico exagerado em simulações. Essa discrepância manifesta-se como concentrações críticas de micelização subestimadas para surfactantes e energias livres de adsorção na interface óleo-água superestimadas. Embora os alcanos sirvam como um sistema modelo fundamental para compreender a solvatação não polar, as parametrizações padrão — especificamente a combinação do campo de força TraPPE-UA para alcanos com vários modelos de água e as regras de mistura padrão de Lorentz–Berthelot — falham em reproduzir com precisão os dados experimentais de solubilidade. Estudos anteriores notaram esses desvios, mas uma correção sistemática aplicável a uma ampla gama de comprimentos de cadeia de alcano (de metano a eicosano) e múltiplos modelos de água tem sido inexistente.

Metodologia
Os autores calcularam as energias livres de hidratação para alcanos lineares variando de metano ($C_1$) a eicosano ($C_{20}H_{42}$) usando Perturbação de Energia Livre (FEP) dentro do ensemble isotérmico-isobárico. O estudo empregou o campo de força united-atom TraPPE-UA para alcanos combinado com quatro modelos de água distintos:

• Modelos de três sítios: SPC/E e OPC3.
• Modelos de quatro sítios: TIP4P/2005 e OPC.

As simulações foram conduzidas usando LAMMPS com um potencial de Lennard-Jones de núcleo macio para evitar singularidades durante transformações alquímicas. O parâmetro de acoplamento $\lambda$ foi variado em 40 janelas. O estudo utilizou as regras de mistura padrão de Lorentz–Berthelot para comparações iniciais e aplicou correções de cauda de longo alcance para potenciais não deslocados.

Para corrigir a superestimação sistemática, os autores utilizaram uma estratégia de reparametrização baseada na decomposição da energia livre de hidratação:
$$\Delta G_{hyd} = \Delta G_{cavity} + \Delta H_{att}$$
onde $\Delta G_{cavity}$ é a energia livre de criação de uma cavidade do tamanho do soluto e $\Delta H_{att}$ representa as interações atrativas soluto-solvente (aproximadas como a energia de interação de Lennard-Jones, $\Delta U_{att}$). Ao utilizar valores estabelecidos para $\Delta G_{cavity}$, os autores derivaram o $\Delta H_{att}$ experimental alvo e ajustaram o parâmetro de profundidade do poço de Lennard-Jones entre alcano e água ($\epsilon$) para corresponder a esses alvos a 300 K. O parâmetro de distância ($\sigma$) foi mantido constante, pois as energias livres de cavidade são altamente sensíveis a $\sigma$.

Adicionalmente, o estudo avaliou:

1. O Campo de Força Geral Amber (GAFF) combinado com os modelos de água TIP4P/2005 e SPC/E.
2. O modelo HH-alkane proposto anteriormente (Ashbaugh et al.) estendido até $C_{20}$.
3. O impacto de aplicar um potencial de Lennard-Jones deslocado (onde o potencial é zerado no corte) versus potenciais não deslocados com correções de cauda.

Principais Resultados

• Superestimação Sistemática: Todos os modelos de água testados (SPC/E, OPC3, TIP4P/2005, OPC) combinados com TraPPE-UA superestimaram $\Delta G_{hyd}$ em relação aos dados experimentais e às estimativas de contribuição de grupo. O desvio foi mais pronunciado para modelos de três sítios e aumentou com o comprimento da cadeia de alcano.
• Sucesso da Reparametrização: O ajuste da profundidade do poço de Lennard-Jones entre alcano e água ($\epsilon$) em aproximadamente 5% (aumentando-a em relação ao valor de Lorentz–Berthelot) alinhou com sucesso os resultados da simulação aos valores experimentais a 300 K.
  • Os modelos reparametrizados mantiveram excelente concordância com os dados experimentais em uma faixa de temperatura de 290–350 K.
  • Os modelos TIP4P/2005, SPC/E e OPC3 reparametrizados mostraram desvios quadráticos médios (RMSDs) de 0,07, 0,10 e 0,06 $kT$, respectivamente.
• Desempenho do GAFF: O campo de força GAFF combinado com os modelos de água TIP4P/2005 ou SPC/E produziu desvios menores em relação aos valores experimentais comparado à combinação original TraPPE-UA/TIP4P/2005, embora ainda exibisse tendências sistemáticas (superestimação para alcanos pequenos com TIP4P/2005, subestimação para alcanos grandes).
• Modelo HH-Alcano: A extensão do modelo HH-alkane (originalmente validado apenas até o butano) para eicosano confirmou sua capacidade de reproduzir dados experimentais/de contribuição de grupo até $C_{20}$.
• Deslocamento de Potencial: O estudo demonstrou que deslocar o potencial de Lennard-Jones pelo seu valor no corte esférico introduz desvios sistemáticos em $\Delta G_{hyd}$. Embora adicionar correções de cauda ao potencial deslocado possa recuperar matematicamente os resultados não deslocados, o ato de deslocar em si exacerba erros se não for rigorosamente corrigido, destacando a sensibilidade das energias livres alquímicas a esquemas de truncamento de potencial.

Significado e Alegações
O artigo afirma que a superestimação sistemática da hidratação hidrofóbica em campos de força comuns é um fenômeno robusto que pode ser corrigido através de um ajuste de uma única etapa no parâmetro de profundidade do poço de Lennard-Jones entre alcano e água. Os autores enfatizam que:

1. Transferibilidade da Correção: O aumento de ~5% em $\epsilon$ é consistente entre diferentes modelos de água (tanto de 3 sítios quanto de 4 sítios), sugerindo um mecanismo de correção generalizável para solutos não polares.
2. Limitações das Regras de Mistura: Os resultados destacam que campos de força parametrizados em termodinâmica de componentes puros frequentemente performam mal em misturas quando regras de mistura heurísticas (como Lorentz–Berthelot) são aplicadas sem calibração específica de interações cruzadas.
3. Insight Metodológico: O estudo fornece um caminho prático para melhorar as previsões de energia livre de hidratação ao calibrar parâmetros de interação diretamente contra energias livres de formação de cavidade e dados experimentais de solubilidade, em vez de depender exclusivamente de regras de mistura padrão.
4. Deslocamentos de Potencial: O trabalho serve como um alerta quanto ao uso de potenciais deslocados em cálculos de energia livre, demonstrando que tais deslocamentos podem introduzir erros sistemáticos significativos na termodinâmica de hidratação.

Os autores concluem que energias livres de hidratação precisas para espécies não polares como alcanos podem ser alcançadas calibrando diretamente os parâmetros de interação usando energias livres de formação de cavidade, melhorando assim a confiabilidade de simulações moleculares para auto-organização hidrofóbica e fenômenos de solvatação.