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

Résumé technique : Énergies libres d'hydratation des alcanes linéaires : Déviations systématiques dans les modèles d'eau courants et leur correction

Énoncé du problème
Les champs de force moléculaires courants surestiment systématiquement les énergies libres d'hydratation ($\Delta G_{hyd}$) des solutés hydrophobes, conduisant à un effet hydrophobe exagéré dans les simulations. Cette discordance se manifeste par des concentrations micellaires critiques sous-estimées pour les tensioactifs et des énergies libres d'adsorption à l'interface huile–eau surestimées. Bien que les alcanes servent de système modèle fondamental pour comprendre la solvatation non polaire, les paramétrisations standard — spécifiquement la combinaison du champ de force TraPPE-UA pour les alcanes et de divers modèles d'eau avec les règles de mélange standard de Lorentz–Berthelot — échouent à reproduire avec précision les données de solubilité expérimentales. Des études antérieures ont noté ces déviations, mais une correction systématique applicable sur une large gamme de longueurs de chaînes d'alcanes (du méthane à l'éicosane) et de plusieurs modèles d'eau a fait défaut.

Méthodologie
Les auteurs ont calculé les énergies libres d'hydratation pour des alcanes linéaires allant du méthane ($C_1$) à l'éicosane ($C_{20}H_{42}$) en utilisant la perturbation d'énergie libre (FEP) dans l'ensemble isotherme–isobare. L'étude a employé le champ de force united-atom TraPPE-UA pour les alcanes combiné à quatre modèles d'eau distincts :

• Modèles à trois sites : SPC/E et OPC3.
• Modèles à quatre sites : TIP4P/2005 et OPC.

Les simulations ont été réalisées avec LAMMPS en utilisant un potentiel de Lennard-Jones à noyau doux pour éviter les singularités lors des transformations alchimiques. Le paramètre de couplage $\lambda$ a été varié sur 40 fenêtres. L'étude a utilisé les règles de mélange standard de Lorentz–Berthelot pour les comparaisons initiales et a appliqué des corrections de queue à longue portée pour les potentiels non décalés.

Pour corriger la surestimation systématique, les auteurs ont utilisé une stratégie de reparamétrisation basée sur la décomposition de l'énergie libre d'hydratation :
$$\Delta G_{hyd} = \Delta G_{cavity} + \Delta H_{att}$$
où $\Delta G_{cavity}$ est l'énergie libre de création d'une cavité de la taille du soluté et $\Delta H_{att}$ représente les interactions attractives soluté–solvant (approximées comme l'énergie d'interaction de Lennard-Jones, $\Delta U_{att}$). En utilisant des valeurs établies pour $\Delta G_{cavity}$, les auteurs ont dérivé le $\Delta H_{att}$ expérimental cible et ajusté le paramètre de profondeur de puits de Lennard-Jones entre l'alcane et l'eau ($\epsilon$) pour correspondre à ces cibles à 300 K. Le paramètre de distance ($\sigma$) a été maintenu constant, car les énergies libres de cavité sont hautement sensibles à $\sigma$.

De plus, l'étude a évalué :

1. Le champ de force General Amber Force Field (GAFF) combiné avec les modèles d'eau TIP4P/2005 et SPC/E.
2. Le modèle HH-alcane précédemment proposé (Ashbaugh et al.) étendu jusqu'à $C_{20}$.
3. L'impact de l'application d'un potentiel de Lennard-Jones décalé (où le potentiel est ramené à zéro à la coupure) par rapport aux potentiels non décalés avec corrections de queue.

Résultats clés

• Surestimation systématique : Tous les modèles d'eau testés (SPC/E, OPC3, TIP4P/2005, OPC) combinés à TraPPE-UA ont surestimé $\Delta G_{hyd}$ par rapport aux données expérimentales et aux estimations par contribution de groupe. La déviation était plus prononcée pour les modèles à trois sites et augmentait avec la longueur de la chaîne d'alcane.
• Succès de la reparamétrisation : L'ajustement de la profondeur de puits de Lennard-Jones entre l'alcane et l'eau ($\epsilon$) d'environ 5 % (en l'augmentant par rapport à la valeur de Lorentz–Berthelot) a aligné avec succès les résultats de simulation sur les valeurs expérimentales à 300 K.
  • Les modèles reparamétrés ont maintenu un excellent accord avec les données expérimentales sur une plage de température de 290–350 K.
  • Les modèles TIP4P/2005, SPC/E et OPC3 reparamétrés ont montré des écarts quadratiques moyens (RMSD) de 0,07, 0,10 et 0,06 $kT$, respectivement.
• Performance du GAFF : Le champ de force GAFF combiné avec les modèles d'eau TIP4P/2005 ou SPC/E a produit des écarts plus faibles par rapport aux valeurs expérimentales que la combinaison originale TraPPE-UA/TIP4P/2005, bien qu'il ait encore présenté des tendances systématiques (surestimation pour les petits alcanes avec TIP4P/2005, sous-estimation pour les grands alcanes).
• Modèle HH-Alcane : L'extension du modèle HH-alcane (initialement validé uniquement jusqu'au butane) à l'éicosane a confirmé sa capacité à reproduire les données expérimentales/de contribution de groupe jusqu'à $C_{20}$.
• Décalage des potentiels : L'étude a démontré que le décalage du potentiel de Lennard-Jones par sa valeur à la coupure sphérique introduit des déviations systématiques dans $\Delta G_{hyd}$. Bien que l'ajout de corrections de queue au potentiel décalé puisse mathématiquement récupérer les résultats non décalés, l'acte de décalage lui-même exacerbe les erreurs si elles ne sont pas rigoureusement corrigées, soulignant la sensibilité des énergies libres alchimiques aux schémas de troncature des potentiels.

Signification et revendications
L'article revendique que la surestimation systématique de l'hydratation hydrophobe dans les champs de force courants est un phénomène robuste qui peut être corrigé par un ajustement en une seule étape du paramètre de profondeur de puits de Lennard-Jones entre l'alcane et l'eau. Les auteurs soulignent que :

1. Transférabilité de la correction : L'augmentation d'environ 5 % de $\epsilon$ est cohérente à travers différents modèles d'eau (à la fois 3 et 4 sites), suggérant un mécanisme de correction généralisable pour les solutés non polaires.
2. Limites des règles de mélange : Les résultats mettent en évidence que les champs de force paramétrés sur la thermodynamique des composants purs fonctionnent souvent mal dans les mélanges lorsque des règles de mélange heuristiques (comme Lorentz–Berthelot) sont appliquées sans calibration spécifique des interactions croisées.
3. Insight méthodologique : L'étude fournit une voie pratique pour améliorer les prédictions d'énergie libre d'hydratation en calibrant directement les paramètres d'interaction contre les énergies libres de formation de cavité et les données de solubilité expérimentales, plutôt que de s'appuyer uniquement sur des règles de mélange standard.
4. Décalages de potentiels : Ce travail sert d'avertissement concernant l'utilisation de potentiels décalés dans les calculs d'énergie libre, démontrant que de tels décalages peuvent introduire des erreurs systématiques significatives dans la thermodynamique de l'hydratation.

Les auteurs concluent que des énergies libres d'hydratation précises pour les espèces non polaires comme les alcanes peuvent être obtenues en calibrant directement les paramètres d'interaction à l'aide d'énergies libres de formation de cavité, améliorant ainsi la fiabilité des simulations moléculaires pour l'auto-assemblage hydrophobe et les phénomènes de solvatation.