Thermodynamics of hard-sphere fluids in polydisperse random porous media: Extended scaled particle theory

Resumo Técnico: Termodinâmica de Fluidos de Esferas Duras em Meios Porosos Aleatórios Polidispersos

1. O Problema
O estudo foca na descrição termodinâmica de fluidos de esferas duras (HS - Hard-Sphere) confinados em meios porosos aleatórios que apresentam polidispersidade de tamanho (distribuição de tamanhos variados das partículas da matriz).

• Contexto: Sistemas de matéria mole (coloides, proteínas, polímeros) frequentemente exibem distribuição de tamanhos. A polidispersidade afeta significativamente o empacotamento, as propriedades termodinâmicas e os equilíbrios de fase.
• Limitação Atual: Embora existam teorias analíticas robustas para fluidos em matrizes monodispersas (tamanho único) e métodos numéricos (como simulações de Monte Carlo) para matrizes polidispersas, falta uma descrição totalmente analítica que lide simultaneamente com a polidispersidade da matriz e altas frações de empacotamento da matriz ($\eta_0$). Teorias existentes frequentemente falham ou divergem em regimes de alta densidade da matriz.

2. Metodologia
Os autores desenvolveram uma extensão da Teoria da Partícula Escalonada (SPT - Scaled Particle Theory), especificamente generalizando a abordagem SPT2 para sistemas polidispersos.

• Modelo: O fluido confinado é tratado como um sistema monodisperso (partículas de tamanho fixo), enquanto a matriz porosa é modelada como um conjunto "congelado" (quenched) de obstáculos de esferas duras com distribuição de tamanhos contínua.
• Generalização da SPT2:
  • A teoria calcula o trabalho reversível necessário para criar uma cavidade no meio poroso para inserir uma partícula de teste escalonada.
  • Introduziram conceitos de porosidade geométrica ($\phi_0$) e porosidade termodinâmica ($\phi_1$).
  • Para lidar com a polidispersidade, substituíram somatórios discretos por integrais sobre a função de distribuição de tamanho da matriz $f(\sigma_0)$, utilizando momentos da distribuição ($m_l$).
• Desenvolvimento de Novas Aproximações:
  • Revisaram aproximações anteriores (SPT2a, SPT2b, SPT2b1, SPT2b3*) e identificaram que elas apresentavam divergências ou perda de precisão em altas frações de empacotamento da matriz ($\eta_0 \geq 0.25$).
  • Inovação Principal: Propuseram uma nova aproximação, denominada SPT2b3****. Esta versão é derivada da expansão da SPT2b3* (equação 2.17), negligenciando o último termo da expansão logarítmica (equação 2.16). Esta simplificação elimina divergências indesejadas e melhora a estabilidade numérica em altas densidades.
• Correções: As expressões finais para o potencial químico e pressão incluem correções de Carnahan-Starling para garantir precisão em fluidos densos.

3. Contribuições Principais

• Descrição Analítica Completa: Derivação de expressões analíticas fechadas para o potencial químico excessivo e pressão de fluidos HS em matrizes polidispersas.
• Nova Aproximação SPT2b3****: Desenvolvimento de uma aproximação que supera as limitações das teorias anteriores, mantendo a precisão em uma faixa mais ampla de frações de empacotamento da matriz ($0.1 \leq \eta_0 \leq 0.3$).
• Versatilidade de Distribuições: A metodologia é aplicável a qualquer distribuição contínua de tamanhos. Os autores demonstraram isso utilizando distribuições Retangulares e Schultz-Zimm.
• Generalização: O framework pode ser estendido para fluidos polidispersos confinados em matrizes polidispersas e para sistemas com interações atrativas ou associativas.

4. Resultados
Os resultados teóricos foram validados comparando-os com dados de simulações de Monte Carlo (MC) da literatura (Leon et al., Ref. [45]) e teorias de equações integrais.

• Comparação com Simulações:
  • Em baixas e médias frações de empacotamento da matriz ($\eta_0 \leq 0.2$), as aproximações SPT2b3* e a nova SPT2b3** mostram excelente acordo com os dados de MC.
  • Em altas frações de empacotamento ($\eta_0 = 0.25$ e $0.3$), as aproximações anteriores (SPT2a, SPT2b1, SPT2b3*) começam a desviar-se significativamente dos dados de simulação.
  • A SPT2b3**** mantém um acordo notável com os dados de Monte Carlo mesmo em $\eta_0 = 0.3$, capturando corretamente a dependência do potencial químico com a densidade do fluido e a largura da distribuição de tamanhos da matriz.
• Efeito da Polidispersidade:
  • O aumento da polidispersidade da matriz reduz o potencial químico do fluido para uma dada densidade, indicando uma maior eficiência de empacotamento dos obstáculos de tamanhos variados.
  • A influência da polidispersidade torna-se mais pronunciada à medida que a fração de empacotamento da matriz aumenta.
• Distribuição de Tamanho de Poro: A teoria permite calcular a distribuição de tamanho de poro ($v_1$), mostrando que, em altas densidades da matriz, a posição do pico dessa distribuição é sensível à largura da distribuição de tamanhos das partículas da matriz.

5. Significado e Impacto

• Ferramenta de Referência: O trabalho fornece um modelo de referência analítico preciso para sistemas complexos de matéria mole confinada, substituindo a necessidade de simulações computacionais intensivas para estimativas termodinâmicas rápidas.
• Aplicações Práticas: O modelo é diretamente aplicável ao estudo de:
  • Soluções de proteínas e polímeros em ambientes biologicamente "aglomerados" (crowding).
  • Separação de fases e agregação em géis porosos.
  • Difusão de partículas em meios desordenados.
• Avanço Teórico: Ao resolver o problema da divergência em altas densidades de matriz, esta extensão da SPT preenche uma lacuna importante na teoria de fluidos confinados, permitindo a exploração de regimes de densidade anteriormente inacessíveis a métodos puramente analíticos.

Em suma, o artigo estabelece um novo padrão analítico para a termodinâmica de fluidos em meios porosos polidispersos, oferecendo uma ferramenta robusta e versátil para a física estatística de sistemas desordenados.