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

Technische Samenvatting: Hydratatievrije Energieën van Lineaire Alkanen: Systematische Afwijkingen in Gangbare Watermodellen en Hun Correctie

Probleemstelling
Gangbare moleculaire krachtvelden overschatten systematisch de hydratatievrije energieën ($\Delta G_{hyd}$) van hydrofobe opgeloste stoffen, wat leidt tot een overdreven hydrofoob effect in simulaties. Deze discrepantie manifesteert zich als onderschatte kritische micelconcentraties voor surfactanten en overschatte vrije energieën voor adsorptie aan de olie-water-interface. Hoewel alkanen dienen als een fundamenteel modelsysteem voor het begrijpen van niet-polaire solvatatie, falen standaardparameteriseringen—specifiek de combinatie van het TraPPE-UA krachtveld voor alkanen en diverse watermodellen met standaard Lorentz–Berthelot mengregels—om experimentele oplosbaarheidsdata nauwkeurig te reproduceren. Eerdere studies hebben deze afwijkingen opgemerkt, maar een systematische correctie die toepasbaar is op een breed scala aan alkaanketenlengtes (van methaan tot eicosaan) en meerdere watermodellen, ontbrak.

Methodologie
De auteurs berekenden hydratatievrije energieën voor lineaire alkanen variërend van methaan ($C_1$) tot eicosaan ($C_{20}H_{42}$) met behulp van Vrije Energie Perturbatie (FEP) binnen het isotherm-isobaar ensemble. De studie gebruikte het TraPPE-UA united-atom krachtveld voor alkanen, gecombineerd met vier verschillende watermodellen:

• Drie-site modellen: SPC/E en OPC3.
• Vier-site modellen: TIP4P/2005 en OPC.

Simulaties werden uitgevoerd met LAMMPS, gebruikmakend van een soft-core Lennard-Jones potentiaal om singulariteiten tijdens alchemische transformaties te voorkomen. De koppelingsparameter $\lambda$ werd gevarieerd in 40 vensters. De studie gebruikte standaard Lorentz–Berthelot mengregels voor initiële vergelijkingen en paste lange-afstand staartcorrecties toe voor niet-verschoven potentialen.

Om de systematische overschatting te corrigeren, gebruikten de auteurs een reparameterisatiestrategie gebaseerd op de decompositie van de hydratatievrije energie:
$$\Delta G_{hyd} = \Delta G_{cavity} + \Delta H_{att}$$
waarbij $\Delta G_{cavity}$ de vrije energie is om een opgeloste-stof-grootte holte te creëren en $\Delta H_{att}$ de aantrekkende opgeloste-stof-oplosmiddel-interacties vertegenwoordigt (benaderd als de Lennard-Jones interactie-energie, $\Delta U_{att}$). Door gebruik te maken van gevestigde waarden voor $\Delta G_{cavity}$, leidden de auteurs de doelwit experimentele $\Delta H_{att}$ af en pasten ze de Lennard-Jones putdiepte-parameter ($\epsilon$) voor alkaan-water aan om deze doelen te matchen bij 300 K. De afstand-parameter ($\sigma$) werd constant gehouden, aangezien holte-vrije energieën zeer gevoelig zijn voor $\sigma$.

Daarnaast evalueerde de studie:

1. Het General Amber Force Field (GAFF) gecombineerd met TIP4P/2005 en SPC/E watermodellen.
2. Het eerder voorgestelde HH-alkaan model (Ashbaugh et al.) uitgebreid tot $C_{20}$.
3. De impact van het toepassen van een verschoven Lennard-Jones potentiaal (waarbij de potentiaal op de cutoff wordt genulld) versus niet-verschoven potentialen met staartcorrecties.

Belangrijkste Resultaten

• Systematische Overschatting: Alle geteste watermodellen (SPC/E, OPC3, TIP4P/2005, OPC) gecombineerd met TraPPE-UA overschatten $\Delta G_{hyd}$ ten opzichte van experimentele data en schattingen op basis van groep-bijdragen. De afwijking was duidelijker voor drie-site modellen en nam toe met de alkaan-ketenlengte.
• Succes van Reparameterisatie: Het aanpassen van de alkaan-water Lennard-Jones putdiepte ($\epsilon$) met ongeveer 5% (verhogend ten opzichte van de Lorentz–Berthelot waarde) bracht de simulatieresultaten succesvol in overeenstemming met experimentele waarden bij 300 K.
  • De gereparameteriseerde modellen behielden uitstekende overeenstemming met experimentele data over een temperatuurbereik van 290–350 K.
  • De gereparameteriseerde TIP4P/2005, SPC/E en OPC3 modellen vertoonden root mean squared deviations (RMSD's) van respectievelijk 0,07, 0,10 en 0,06 $kT$.
• GAFF Prestaties: Het GAFF krachtveld gecombineerd met TIP4P/2005 of SPC/E water leverde kleinere afwijkingen op ten opzichte van experimentele waarden dan de originele TraPPE-UA/TIP4P/2005 combinatie, hoewel het nog steeds systematische trends vertoonde (overschatting voor kleine alkanen met TIP4P/2005, onderschatting voor grote alkanen).
• HH-Alkaan Model: Het uitbreiden van het HH-alkaan model (oorspronkelijk alleen gevalideerd tot butaan) tot eicosaan bevestigde het vermogen om experimentele/groep-bijdrage data te reproduceren tot $C_{20}$.
• Potentiaal Verschuivingen: De studie toonde aan dat het verschuiven van de Lennard-Jones potentiaal met zijn waarde bij de sferische cutoff systematische afwijkingen in $\Delta G_{hyd}$ introduceert. Hoewel het toevoegen van staartcorrecties aan de verschoven potentiaal wiskundig de niet-verschoven resultaten kan herstellen, verergert de verschuiving zelf de fouten als deze niet strikt wordt gecorrigeerd, wat de gevoeligheid van alchemische vrije energieën voor potentiaal-truncatieschema's benadrukt.

Betekenis en Beweringen
Het artikel beweert dat de systematische overschatting van hydrofobe hydratatie in gangbare krachtvelden een robuust fenomeen is dat kan worden gecorrigeerd door een een-staps aanpassing van de alkaan-water Lennard-Jones putdiepte-parameter. De auteurs benadrukken dat:

1. Transferabiliteit van Correctie: De ~5% verhoging van $\epsilon$ consistent is over verschillende watermodellen (zowel 3-site als 4-site), wat wijst op een generaliseerbaar correctiemechanisme voor niet-polaire opgeloste stoffen.
2. Beperkingen van Mengregels: De resultaten onderstrepen dat krachtvelden geparameteriseerd op thermodynamica van zuivere componenten vaak slecht presteren in mengsels wanneer heuristische mengregels (zoals Lorentz–Berthelot) worden toegepast zonder specifieke kruis-interactie kalibratie.
3. Methodologisch Inzicht: De studie biedt een praktische weg om voorspellingen van hydratatievrije energieën te verbeteren door interactieparameters direct te kalibreren tegen vrije energieën voor holte-vorming en experimentele oplosbaarheidsdata, in plaats van uitsluitend te vertrouwen op standaard mengregels.
4. Potentiaal Verschuivingen: Het werk dient als een waarschuwing voor het gebruik van verschoven potentialen in vrije energie-berekeningen, en toont aan dat dergelijke verschuivingen significante systematische fouten kunnen introduceren in hydratatiethermodynamica.

De auteurs concluderen dat nauwkeurige hydratatievrije energieën voor niet-polaire soorten zoals alkanen kunnen worden bereikt door interactieparameters direct te kalibreren met behulp van vrije energieën voor holte-vorming, waardoor de betrouwbaarheid van moleculaire simulaties voor hydrofobe zelfassemblage en solvatatieverschijnselen wordt verbeterd.