Phase-Field Models for Particle-Stabilised Emulsions

Resumo Técnico: Modelos de Campo de Fase para Emulsões Estabilizadas por Partículas

1. O Problema

As emulsões estabilizadas por partículas (como emulsões de Pickering e géis de emulsão interfacialmente travados bicontínuos, conhecidos como bijels) são fundamentais na ciência de materiais macios. Embora simulações computacionais ofereçam insights cruciais sobre sua formação, os métodos atuais (como Dinâmica Molecular Coarse-Grained, Dinâmica de Partículas Dissipativas e Lattice-Boltzmann) enfrentam limitações severas de escala.

• Limitação Principal: Essas técnicas resolvem partículas individualmente. Para sistemas macroscópicos relevantes, como a formação de bijels via Separação de Fases Induzida por Transferência de Solvente (STrIPS), que envolvem bilhões de nanopartículas e ocorrem em escalas de tempo de segundos e centenas de micrômetros, a resolução individual torna-se computacionalmente inviável.
• Desafio Específico: Métodos de campo de fase existentes que tentam incluir partículas frequentemente resolvem-nas como objetos discretos ou acoplam suas densidades diretamente ao potencial químico do líquido de forma que gera deformações não físicas e instabilidades numéricas na interface líquido-líquido.

2. Metodologia

Os autores desenvolveram uma nova estratégia baseada na teoria de campo de fase, tratando as nanopartículas como um campo contínuo em vez de objetos discretos, permitindo a simulação de grandes escalas espaciais e temporais de forma eficiente.

• Funcional de Energia Livre: O modelo baseia-se em um funcional de energia livre que descreve uma mistura heterogênea de líquidos imiscíveis ($\phi$) e nanopartículas ($\psi$). Inclui termos de Flory-Huggins para os líquidos, uma aproximação de mistura ideal para as partículas e um termo de acoplamento interfacial que torna as partículas ativas na superfície.
• Equações Dinâmicas Revisadas:
  • As equações iniciais de Cahn-Hilliard acopladas mostraram-se inadequadas, pois a acumulação de partículas distorcia a interface líquida e gerava tensões interfaciais negativas (não físicas).
  • Solução: Os autores propuseram um desacoplamento parcial das equações dinâmicas. O potencial químico do líquido é calculado independentemente da densidade local de partículas (assumindo que partículas muito menores não alteram o potencial químico do líquido macroscopicamente), enquanto as partículas ainda respondem ao gradiente da interface. Isso elimina as distorções da interface.
• Mecanismo de "Jamming" (Travamento): Para modelar a estabilização da emulsão, foi introduzida uma mobilidade dependente da densidade. Quando a densidade de partículas na interface atinge um limiar crítico ($\psi_c$), a mobilidade das espécies cai abruptamente. Isso simula o "travamento" (jamming) das partículas na interface, impedindo o transporte de massa e, consequentemente, a coalescência e o amadurecimento de Ostwald.
• Aplicação ao STrIPS: O modelo foi integrado a um framework existente de STrIPS, onde a difusão de um solvente altera o parâmetro de interação entre os líquidos, induzindo a separação de fases.

3. Contribuições Principais

• Novo Formalismo de Campo de Fase: Desenvolvimento de um conjunto de equações diferenciais parciais que descreve a dinâmica acoplada da separação de fases líquida e adsorção de nanopartículas sem a necessidade de resolver partículas individualmente.
• Estabilidade Numérica: Resolução do problema de instabilidade e deformação de interface encontrado em abordagens anteriores, permitindo simulações estáveis com altos parâmetros de acoplamento.
• Modelagem de Escala Macroscópica: Capacidade de simular a formação de bijels em escalas de tempo e espaço relevantes para experimentos reais (milímetros/micrômetros e segundos), algo impossível com métodos de partículas discretas.
• Mecanismo de Estabilização: Demonstração clara de como o "travamento" interfacial inibe o amadurecimento de Ostwald através da redução da mobilidade interfacial.

4. Resultados

• Formação de Emulsões: O modelo reproduz com sucesso a formação de dois tipos de morfologias:
  • Emulsão de Pickering: Gotas discretas estabilizadas por partículas (composição inicial não crítica).
  • Bijel: Rede bicontínua de canais líquidos interconectados (composição inicial crítica).
• Supressão do Amadurecimento de Ostwald: Simulações mostram que, na ausência de partículas, as emulsões sofrem coalescência e amadurecimento de Ostwald (crescimento de gotas grandes às custas das pequenas). Com a presença de partículas e o mecanismo de jamming, a estrutura é travada em um estado metaestável, impedindo o crescimento posterior.
• Influência da Concentração de Partículas:
  • Maior concentração inicial de nanopartículas ($\psi_0$) leva a um travamento mais rápido da separação de fases.
  • Isso resulta em domínios líquidos menores e mais monodispersos.
• Simulação de STrIPS:
  • O modelo capturou o gradiente característico de tamanho de domínio observado experimentalmente em bijels fabricados via STrIPS (domínios menores perto da interface com o solvente e maiores no substrato).
  • O gradiente surge devido à diferença nas taxas locais de separação de fases mediada pelo solvente e adsorção de partículas.
  • A relação inversa entre concentração de partículas e tamanho médio do domínio foi confirmada, alinhando-se com evidências experimentais.

5. Significado e Implicações

Este trabalho representa um avanço significativo na modelagem computacional de sistemas de matéria macia complexa.

• Viabilidade Computacional: Permite o estudo de fenômenos de formação de emulsões em escalas que antes eram inacessíveis, preenchendo a lacuna entre simulações microscópicas e observações macroscópicas.
• Projeto de Materiais: A ferramenta oferece um meio eficiente para investigar como parâmetros de processo (como concentração de partículas, taxa de remoção de solvente e propriedades de molhabilidade) influenciam a morfologia final de materiais como membranas, catalisadores e sistemas de entrega de fármacos.
• Limitações e Futuro: O modelo atual não inclui hidrodinâmica (fluxo convectivo) nem flutuações térmicas explícitas, o que pode ser relevante em estágios tardios de coalescência ou em sistemas STrIPS específicos. No entanto, a abordagem fornece uma base sólida para futuras extensões que incorporem esses efeitos, mantendo a eficiência computacional.

Em suma, o modelo de campo de fase proposto oferece uma abordagem promissora e robusta para a investigação morfológica de emulsões estabilizadas por partículas, validando mecanismos teóricos e fornecendo insights para o controle de materiais complexos.