Accurate Solvation Properties in supercritical CO$_2$ with Molecular Density Functional Theory

Technische Samenvatting: Nauwkeurige solvatie-eigenschappen in superkritisch CO2 met Moleculaire Dichtheidsfunctionaaltheorie

Probleemstelling
Superkritisch koolstofdioxide (scCO2) is een uiterst efficiënt, milieuvriendelijk oplosmiddel voor industriële toepassingen variërend van extractie tot de synthese van nanodeeltjes. Het optimaliseren van deze processen vereist nauwkeurige voorspellingen van solvatie-eigenschappen, specif seguitoerde solvatievrije energieën en solvatiestructuren, over variërende thermodynamische condities. Hoewel Molecular Dynamics (MD) simulaties een hoge voorspellende nauwkeurigheid en microscopisch inzicht bieden, lijden zij onder prohibitieve computationele kosten, met name bij het berekenen van thermodynamische eigenschappen via geavanceerde samplingstrategieën (bijv. Widom-insertie, thermodynamische integratie). Daartegenover staan bestaande parametrische methoden (vergelijking-van-toestand, dichtheid-gebaseerd, of solvatiemodellen) die lagere computationele kosten hebben, maar lijden onder beperkte transfereerbaarheid, omdat ze vaak her-parametrisatie vereisen voor nieuwe soluten of systeemwijzigingen. Er is behoefte aan een methode die de nauwkeurigheid van MD combineert met de computationele efficiëntie die vereist is voor grootschalige parametrische studies.

Methodologie
Deze studie maakt gebruik van Moleculaire Dichtheidsfunctionaaltheorie (MDFT), een klassieke cDFT-formalismen aangepast voor moleculaire vloeistoffen. De benadering minimaliseert een solvatiefunctionaal, $\Delta\Omega_{solv}[\rho]$, om de evenwichtsmoleculaire dichtheid $\rho(\mathbf{r}, \omega)$ en de solvatievrije energie te bepalen. De functionaal wordt opgedeeld in drie componenten:

1. Ideale Bijdrage ($\Delta\Omega_{id}$): Houdt rekening met de entropie van de solventdichtheid.
2. Externe Bijdrage ($F_{ext}$): Representeert solute-solvent interacties afgeleid van standaard krachtvelden (EPM2 voor CO2, OPLS-AA voor moleculaire soluten, SPC/E voor water).
3. Excess Bijdrage ($\Delta\Omega_{exc}$): Representeert solvent-solvent interacties. Deze term wordt benaderd met de Homogene Referentievloeistof (HRF) benadering, die ervan uitgaft dat dichtheidsfluctuaties nabij de solute lijken op die van de homogene bulk. De HRF-term is een kwadratische expansie afhankelijk van de directe paarcorrelatiefuncties (DCF) van de zuivere bulk-solvent, verkregen uit voorafgaande MD-simulaties en integraal vergelijkingstheorie.

Om de beperkingen van de HRF-benadering (die hogere-orde fluctuaties verwaarloost) aan te pakken, onderzochten de auteurs drie niveaus van benadering voor de bridge functionaal ($\Delta\Omega_{bridge}$), die het verschil vertegenwoordigt tussen de exacte excess functionaal en de HRF-term:

• Geen Bridge: $\Delta\Omega_{bridge} = 0$.
• Hard-Sphere Bridge (HSB): Benadert de bridge met een hard-sphere excess functionaal (FMT) gekalibreerd op de bulk vrije energie.
• Coarse-Grained Bridge (CGB): Modelleert de bridge als een functionaal van een coarse-grained dichtheidsveld, afgeleid van een polynomiale expansie van de bulk vergelijking-van-toestand.

De MDFT-berekeningen werden uitgevoerd in een kubische box ($L=30$ Å) met een ruimtelijk rooster van $0.3$ Å en een maximale hoekindex van $n_{max}=4$. Deze resultaten werden gebenchmarkt tegen MD-simulaties met 4000 CO2-moleculen en de Bennett acceptance ratio methode voor vrije energieberekeningen.

Belangrijkste Bijdragen

• Toepassing van Formalisme: De studie past succesvol een moleculaire cDFT-formalismen toe, afhankelijk van zowel ruimtelijke posities als moleculaire oriëntaties, op scCO2.
• Parametrisatie van de Bridge Functionaal: De auteurs hebben twee verschillende bridge functionalen geconstrueerd en geparametriseerd (HSB en CGB) voor scCO2. Ze toonden aan dat de CGB-functionaal, met een coarse-graining lengte van $\sigma_g = 1.45$ Å, de homogene bridge bijdrage uit de vergelijking-van-toestand nauwkeurig reproduceert.
• Computationele Efficiëntie: De studie kwantificeert het computationele voordeel van MDFT, waarbij wordt opgemerkt dat een enkele berekening minder dan één minuut vereist op een laptop (ca. 0.01 CPU·uren), vergeleken met ongeveer 640 CPU·uren voor de equivalente MD-simulatie.

Resultaten
Het MDFT-raamwerk werd getest tegen MD-benchmarks voor twee klassen van soluten: sferische (argon-achtige) soluten en negen moleculaire soluten (inclusclusief water, ethanol, tolueen en CO2 zelf) bij een nabij-kritisch staatpunt ($n_0 = 0.8 n_c$, $T = 1.05 T_c$).

• Structurele Eigenschappen: Voor alle soluten vertoonden de door MDFT voorspelde radiale distributiefuncties (RDF's) een vrijwel perfecte overlap met de MD-resultaten. Deze overeenkomst bleef behouden, ongeacht de gebruikte bridge functionaal benadering (HRF, HSB of CGB).
• Solvatievrije Energieën:
  • Voor sferische soluten overschatte de HRF-benadering (geen bridge) de solvatievrije energieën systematisch voor grote soluten ($\sigma > 7$ Å). Echter, zowel de HSB als de CGB correcties leverden resultaten die bijna ononderscheidbaar waren van de MD-benchmarks over een breed scala aan groottes en interactie-energieën.
  • Voor moleculaire soluten reproduceerden alle drie de MDFT-varianten (HRF, HSB, CGB) de MD solvatievrije energieën met een gemiddelde relatieve fout van ongeveer 1%.
  • Opvallend genoeg bleek voor moleculaire soluten in scCO2 de bijdrage van de bridge functionaal verwaarloosbaar te zijn. De HRF-benadering alleen was voldoende om de thermodynamica te vatten, wat contrasteert met eerdere MDFT-studies naar vloeibaar water waar bridge-termen essentieel waren. Dit wordt toegeschreven aan de afwezigheid van een sterk directioneel waterstofbrugnetwerk in scCO2.
• Nauwkeurigheid vs. Kosten: De MDFT-benadering bereikte een bijna kwantitatieve nauwkeurigheid in zowel structurele als thermodynamische eigenschappen, terwijl de computationele inspanning met een factor van ongeveer 60.000 werd verminderd ten opzichte van MD.

Betekenis en Claims
Het artikel claimt dat MDFT een robuust en computationeel efficiënt alternatief biedt voor MD voor het modelleren van solvatie in scCO2. De primaire betekenis ligt in het aantonen dat de Homogene Referentievloeistof (HRF) benadering voldoende is om de MD benchmark resultaten voor scCO2 solvatie-thermodynamica te herstellen, mits de bulk solvent correlaties nauwkeurig bekend zijn. Deze bevinding suggere houdt in dat de complexe bridge correcties die vaak nodig zijn voor andere vloeistoffen (zoals water) onnodig kunnen zijn voor scCO2 vanwege de specifieke vloeistofkarakteristieken.

De auteurs stellen expliciet dat het huidige werk niet de transfereerbaarheid over het gehele fasediagram vaststelt, maar de methode valideert bij een representatief nabij-kritisch staatpunt. Zij concluderen dat dit raamwerk een fundament biedt voor het uitbreiden van MDFT naar complexere scenario's, zoals interfaces en opsluiting (confinement), en dat toekomstige ontwikkelingen machine learning-gebaseerde functionalen kunnen inzetten om transfereerbaarheid over verschillende thermodynamische condities te bereiken. De geïdentificeerde directe bruikbaarheid is het potentieel om de rationele design en optimalisatie van scCO2-gebaseerde industriële processen te faciliteren door middel van snelle, nauwkeurige screening van solvatie-eigenschappen.