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

Resumen Técnico: Propiedades de Solvatación Precisas en CO2 Supercrítico con la Teoría de la Densidad Funcional Molecular

Planteamiento del Problema
El dióxido de carbono supercrítico (scCO2) es un disolvente altamente eficiente y ambientalmente benigno para aplicaciones industriales que van desde la extracción hasta la síntesis de nanopartículas. La optimización de estos procesos requiere la predicción precisa de las propiedades de solvatación, específicamente las energías libres de solvatación y las estructuras de solvatación, a través de diversas condiciones termodinámicas. Si bien las simulaciones de Dinámica Molecular (MD) proporcionan una alta precisión predictiva y una visión microscópica, sufren de costos computacionales prohibitivos, particularmente al calcular propiedades termodinámicas mediante estrategias de muestreo avanzado (por ejemplo, inserción de Widom, integración termodinámica). Por el contrario, los métodos paramétricos existentes (ecuación de estado, basados en densidad o modelos de solvatación) ofrecen costos computacionales más bajos pero sufren de una transferibilidad limitada, requiriendo a menudo la re-parametrización para nuevos solutos o modificaciones del sistema. Existe la necesidad de un método que combine la precisión de MD con la eficiencia computacional requerida para estudios paramétricos a gran escala.

Metodología
Este estudio emplea la Teoría de la Densidad Funcional Molecular (MDFT), un formalismo de la teoría de la densidad clásica (cDFT) adaptado para fluidos moleculares. El enfoque minimiza un funcional de solvatación, $\Delta\Omega_{solv}[\rho]$, para determinar la densidad molecular de equilibrio $\rho(\mathbf{r}, \omega)$ y la energía libre de solvatación. El funcional se descompone en tres componentes:

1. Contribución Ideal ($\Delta\Omega_{id}$): Contabiliza la entropía de la densidad del disolvente.
2. Contribución Externa ($F_{ext}$): Representa las interacciones soluto-disolvente derivadas de campos de fuerza estándar (EPM2 para CO2, OPLS-AA para solutos moleculares, SPC/E para agua).
3. Contribución de Exceso ($\Delta\Omega_{exc}$): Representa las interacciones disolvente-disolvente. Este término se aproxima utilizando la aproximación de Fluido de Referencia Homogéneo (HRF), que asume que las fluctuaciones de densidad cerca del soluto se asemejan a las del bulbo homogéneo. El término HRF es una expansión cuadrática dependiente de las funciones de correlación de pares directas (DCF) del disolvente puro en el bulbo, obtenidas de simulaciones de MD previas y de la teoría de ecuaciones integrales.

Para abordar las limitaciones de la aproximación HRF (que descuida fluctuaciones de orden superior), los autores investigaron tres niveles de aproximación para el funcional de puente ($\Delta\Omega_{bridge}$), que representa la diferencia entre el funcional de exceso exacto y el término HRF:

• Sin Puente: $\Delta\Omega_{bridge} = 0$.
• Puente de Esferas Duras (HSB): Aproxima el puente utilizando un funcional de exceso de esferas duras (FMT) calibrado con la energía libre del bulbo.
• Puente de Gran Escala (CGB): Modela el puente como un funcional de un campo de densidad de gran escala, derivado de una expansión polinómica de la ecuación de estado del bulbo.

Los cálculos de MDFT se realizaron en una caja cúbica ($L=30$ Å) con una rejilla espacial de $0.3$ Å y un índice angular máximo de $n_{max}=4$. Estos resultados se compararon con simulaciones de MD utilizando 4000 moléculas de CO2 y el método de la relación de aceptación de Bennett para cálculos de energía libre.

Contribuciones Clave

• Aplicación del Formalismo: El estudio aplica con éxito un formalismo de cDFT molecular, dependiente tanto de las posiciones espaciales como de las orientaciones moleculares, al scCO2.
• Parametrización del Funcional de Puente: Los autores construyeron y parametrizaron dos funcionales de puente distintos (HSB y CGB) para el scCO2. Demostraron que el funcional CGB, con una longitud de gran escala de $\sigma_g = 1.45$ Å, reproduce con precisión la contribución de puente homogénea derivada de la ecuación de estado.
• Eficiencia Computacional: El estudio cuantifica la ventaja computacional de MDFT, señalando que un solo cálculo requiere menos de un minuto en una computadora portátil (aprox. 0.01 CPU·horas), en comparación con las aproximadamente 640 CPU·horas de la simulación de MD equivalente.

Resultados
El marco de MDFT fue probado contra referencias de MD para dos clases de solutos: solutos esféricos (tipo argón) y nueve solutos moleculares (incluyendo agua, etanol, tolueno y el propio CO2) en un punto de estado cercano al crítico ($n_0 = 0.8 n_c$, $T = 1.05 T_c$).

• Propiedades Estructurales: Para todos los solutos, las funciones de distribución radial (RDF) predichas por MDFT mostraron un traslape esencialmente perfecto con los resultados de MD. Este acuerdo se mantuvo independientemente de la aproximación del funcional de puente utilizada (HRF, HSB o CGB).
• Energías Libres de Solvatación:
  • Para solutos esféricos, la aproximación HRF (sin puente) sobreestimó sistemáticamente las energías libres de solvatación para solutos grandes ($\sigma > 7$ Å). Sin embargo, tanto las correcciones HSB como las CGB produjeron resultados casi indistinguibles de las referencias de MD en un amplio rango de tamaños y energías de interacción.
  • Para solutos moleculares, las tres variantes de MDFT (HRF, HSB, CGB) reprodujeron las energías libres de solvatación de MD con un error relativo promedio de aproximadamente el 1%.
  • Notablemente, para solutos moleculares en scCO2, se encontró que la contribución del funcional de puente era insignificante. La aproximación HRF por sí sola fue suficiente para capturar la termodinámica, contrastando con estudios previos de MDFT en agua líquida donde los términos de puente eran esenciales. Esto se atribuye a la ausencia de una red de enlaces de hidrógeno fuertemente direccional en el scCO2.
• Precisión vs. Costo: El enfoque de MDFT logró una precisión casi cuantitativa tanto en las propiedades estructurales como en las termodinámicas, reduciendo el esfuerzo computacional en un factor de aproximadamente 60,000 en comparación con MD.

Significancia y Reivindicaciones
El artículo afirma que la MDFT ofrece una alternativa robusta y computacionalmente eficiente a MD para modelar la solvatación en scCO2. La principal significancia radica en demostrar que la aproximación de Fluido de Referencia Homogéneo (HRF) es suficiente para recuperar los resultados de referencia de MD para la termodinámica de solvatación en scCO2, siempre que las correlaciones del bulbo se conozcan con precisión. Este hallazgo sugiere que las complejas correcciones de puente requeridas para otros fluidos (como el agua) pueden ser innecesarias para el scCO2 debido a sus características de fluido específicas.

Los autores declaran explícitamente que el presente trabajo no establece la transferibilidad a través de todo el diagrama de fases, sino que valida el método en un punto de estado representativo cercano al crítico. Concluyen que este marco proporciona una base para extender la MDFT a escenarios más complejos, como interfaces y confinamiento, y que desarrollos futuros podrían aprovechar funcionales basados en aprendizaje automático para lograr la transferibilidad a través de diferentes condiciones termodinámicas. La utilidad inmediata identificada es el potencial para facilitar el diseño racional y la optimización de procesos industriales basados en scCO2 mediante el cribado rápido y preciso de las propiedades de solvatación.