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

Resumo Técnico: Propriedades de Solvatação Precisas em CO2 Supercrítico com Teoria do Funcional de Densidade Molecular

Definição do Problema
O dióxido de carbono supercrítico (scCO2) é um solvente altamente eficiente e ambientalmente benigno para aplicações industriais, que variam desde a extração até a síntese de nanopartículas. Otimizar esses processos requer a previsão precisa das propriedades de solvatação, especificamente energias livres de solvatação e estruturas de solvatação, através de diversas condições termodinâmicas. Embora as simulações de Dinâmica Molecular (MD) forneçam alta precisão preditiva e insights microscópicos, elas sofrem com custos computacionais proibitivos, particularmente no cálculo de propriedades termodinâmicas via estratégias de amostragem avançadas (ex: inserção de Widom, integração termodinâmica). Por outro lado, os métodos paramétricos existentes (equação de estado, baseados em densidade ou modelos de solvatação) oferecem menores custos computacionais, mas sofrem de transferibilidade limitada, muitas vezes exigindo reparametrização para novos solutos ou modificações no sistema. Há uma necessidade de um método que combine a precisão da MD com a eficiência computacional necessária para estudos paramétricos de larga escala.

Metodologia
Este estudo emprega a Teoria do Funcional de Densidade Molecular (MDFT), um formalismo de teoria do funcional de densidade clássica (cDFT) adaptado para fluidos moleculares. A abordagem minimiza um funcional de solvatação, $\Delta\Omega_{solv}[\rho]$, para determinar a densidade molecular de equilíbrio $\rho(\mathbf{r}, \omega)$ e a energia livre de solvatação. O funcional é decomposto em três componentes:

1. Contribuição Ideal ($\Delta\Omega_{id}$): Contabiliza a entropia da densidade do solvente.
2. Contribuição Externa ($F_{ext}$): Representa as interações soluto-solvente derivadas de campos de força padrão (EPM2 para CO2, OPLS-AA para solutos moleculares, SPC/E para água).
3. Contribuição de Excesso ($\Delta\Omega_{exc}$): Representa as interações solvente-solvente. Este termo é aproximado usando a aproximação do Fluido de Referência Homogêneo (HRF), que assume que as flutuações de densidade próximas ao soluto se assemelham às do bulk homogêneo. O termo HRF é uma expansão quadrática dependente das funções de correlação direta de pares (DCF) do solvente puro em bulk, obtidas de simulações de MD prévias e da teoria de equações integrais.

Para abordar as limitações da aproximação HRF (que negligencia flutuações de ordem superior), os autores investigaram três níveis de aproximação para o funcional de ponte ($\Delta\Omega_{bridge}$), que representa a diferença entre o funcional de excesso exato e o termo HRF:

• Sem Ponte: $\Delta\Omega_{bridge} = 0$.
• Ponte de Esferas Rígidas (HSB): Aproxima a ponte usando um funcional de excesso de esferas rígidas (FMT) calibrado pela energia livre do bulk.
• Ponte de Granularidade Grossa (CGB): Modela a ponte como um funcional de um campo de densidade de granularidade grossa, derivado de uma expansão polinomial da equação de estado do bulk.

Os cálculos de MDFT foram realizados em uma caixa cúbica ($L=30$ Å) com uma grade espacial de $0,3$ Å e um índice angular máximo de $n_{max}=4$. Estes resultados foram comparados com simulações de MD usando 4000 moléculas de CO2 e o método da razão de aceitação de Bennett para cálculos de energia livre.

Principais Contribuições

• Aplicação do Formalismo: O estudo aplica com sucesso um formalismo de cDFT molecular, dependente tanto de posições espaciais quanto de orientações moleculares, ao scCO2.
• Parametrização do Funcional de Ponte: Os autores construíram e parametrizaram dois funcionais de ponte distintos (HSB e CGB) para o scCO2. Eles demonstraram que o funcional CGB, com um comprimento de granularidade grossa de $\sigma_g = 1,45$ Å, reproduz com precisão a contribuição de ponte homogênea derivada da equação de estado.
• Eficiência Computacional: O estudo quantifica a vantagem computacional da MDFT, observando que um único cálculo requer menos de um minuto em um laptop (aprox. 0,01 CPU·horas), comparado a aproximadamente 640 CPU·horas para a simulação de MD equivalente.

Resultos
O framework de MDFT foi testado contra benchmarks de MD para duas classes de solutos: solutos esféricos (tipo argônio) e nove solutos moleculares (incluindo água, etanol, tolueno e o próprio CO2) em um ponto de estado próximo ao crítico ($n_0 = 0,8 n_c$, $T = 1,05 T_c$).

• Propriedades Estruturais: Para todos os solutos, as funções de distribuição radial (RDF) previstas pela MDFT mostraram uma sobreposição essencialmente perfeita com os resultados de MD. Este acordo manteve-se independentemente da aproximação do funcional de ponte utilizada (HRF, HSB ou CGB).
• Energias Livres de Solvatação:
  • Para solutos esféricos, a aproximação HRF (sem ponte) superestimou sistematicamente as energias livres de solvatação para solutos grandes ($\sigma > 7$ Å). No entanto, tanto as correções HSB quanto as CGB produziram resultados quase indistinguíveis dos benchmarks de MD em uma ampla gama de tamanhos e energias de interação.
  • Para solutos moleculares, as três variantes de MDFT (HRF, HSB, CGB) reproduziram as energias livres de solvatação de MD com um erro relativo médio de aproximadamente 1%.
  • Notavelmente, para solutos moleculares em scCO2, a contribuição do funcional de ponte mostrou-se negligenciável. A aproximação HRF sozinha foi suficiente para capturar a termodinâmica, contrastando com estudos anteriores de MDFT em água líquida onde os termos de ponte eram essenciais. Isso é atribuído à ausência de uma rede de ligações de hidrogênio fortemente direcional no scCO2.
• Precisão vs. Custo: A abordagem MDFT alcançou precisão quase quantitativa tanto em propriedades estruturais quanto termodinâmicas, reduzindo o esforço computacional em um fator de aproximadamente 60.000 em relação à MD.

Significância e Alegações
O artigo afirma que a MDFT oferece uma alternativa robusta e computacionalmente eficiente à MD para modelar a solvatação em scCO2. A principal significância reside em demonstrar que a aproximação do Fluido de Referência Homogêneo (HRF) é suficiente para recuperar os resultados de benchmark de MD para a termodinâmica de solvatação do scCO2, desde que as correlações do solvente bulk sejam conhecidas com precisão. Este achado sugere que as complexas correções de ponte frequentemente exigidas para outros fluidos (como a água) podem ser desnecessárias para o scCO2 devido às suas características fluídicas específicas.

Os autores declaram explicitamente que o presente trabalho não estabelece a transferibilidade para todo o diagrama de fases, mas sim valida o método em um ponto de estado próximo ao crítico representativo. Eles concluem que este framework fornece uma base para estender a MDFT para cenários mais complexos, como interfaces e confinamento, e que desenvolvimentos futuros poderiam utilizar funcionais baseados em aprendizado de máquina para alcançar a transferibilidade através de diferentes condições termodinâmicas. A utilidade imediata identificada é o potencial para facilitar o design racional e a otimização de processos industriais baseados em scCO2 através da triagem rápida e precisa de propriedades de solvatação.