Gradient dynamics model for chemically driven running drops

1. Problema e Contexto

O artigo aborda o fenômeno de gotas líquidas que se auto-propelam sobre um substrato sólido devido a reações químicas. Especificamente, o foco está em gotas que contêm uma espécie química capaz de alterar a molhabilidade da superfície através de adsorção reversível.

• Limitações de Modelos Anteriores: Modelos existentes para gotas correntes frequentemente assumem um acoplamento "unilateral" (a concentração afeta a molhabilidade, mas a dinâmica da gota não afeta a cinética de reação ou difusão) ou tratam a cobertura do substrato como irreversível. Essas abordagens não são termodinamicamente consistentes, o que significa que não respeitam rigorosamente as leis da termodinâmica de não-equilíbrio (como a produção de entropia e o balanço de energia).
• O Desafio: Desenvolver um modelo que descreva a dinâmica de gotas em um regime de não-equilíbrio sustentado, onde a gota se move indefinidamente sem parar, mantendo um gradiente de molhabilidade local através de reações químicas reversíveis e troca com o meio ambiente.

2. Metodologia

Os autores desenvolveram um modelo baseado na estrutura de Dinâmica de Gradiente (Gradient Dynamics), que garante a consistência termodinâmica.

• Sistema Físico: Consideram gotas sessais em um meio ambiente confinado em uma fenda vertical entre duas placas paralelas. O sistema inclui três espécies químicas:
  1. A: Espécie dissolvida na gota.
  2. B: Espécie adsorvida no substrato.
  3. C: Espécie dissolvida no meio ambiente (ar/líquido circundante).
• Reações Químicas: O modelo considera reações reversíveis:
  • $A \rightleftharpoons B$ (adsorção sob a gota).
  • $B \rightleftharpoons C$ (dessorção para o meio ambiente fora da gota).
• Formulação Termodinâmica:
  • Define-se um funcional de energia livre ($F$) que inclui energia de molhamento, tensão superficial, energia de mistura e contribuições entrópicas das concentrações.
  • As equações de evolução são derivadas da forma padrão de dinâmica de gradiente: $\partial_t \psi = \nabla \cdot (Q \nabla \frac{\delta F}{\delta \psi}) - M \frac{\delta F}{\delta \psi} + R$, onde $Q$ e $M$ são matrizes de mobilidade e $R$ representa as correntes de reação.
  • A consistência termodinâmica é garantida através do balanço de detalhes locais (local detailed balance) nas correntes de reação, conectando as taxas de reação direta e inversa aos potenciais químicos.
• Acionamento de Não-Equilíbrio (Chemostats): Para obter movimento sustentado, o sistema é aberto. As espécies $A$ (na gota) e $C$ (no meio) são tratadas como chemostats (reservatórios com potenciais químicos fixos e distintos, $\mu_a$ e $\mu_c$).
  • Quando $\mu_a \neq \mu_c$, o sistema é forçado para fora do equilíbrio termodinâmico.
  • Isso cria um fluxo contínuo de partículas: adsorção sob a gota e dessorção fora dela, gerando um gradiente de molhabilidade que impulsiona a gota.

3. Contribuições Principais

• Modelo Termodinamicamente Consistente: É a primeira descrição de dinâmica de gradiente para o caso de cobertura reversível do substrato que inclui explicitamente a troca com o meio ambiente de forma termodinamicamente consistente.
• Derivação Reduzida: Os autores partem de um sistema fechado (com conservação de massa total) e, através da aproximação de reservatórios (chemostats), derivam um modelo reduzido para a altura do filme ($h$) e a densidade de adsorvato ($b$).
• Novo Termo de Reação: O modelo reduzido difere de abordagens anteriores (como a do caso irreversível) pela inclusão de um termo de reação adicional ($-R_{bc}$) que captura a troca de adsorvato com o fluido ambiente, essencial para o movimento sustentado.
• Mecanismo de Auto-Propulsão Sustentada: Demonstra que, diferentemente do caso irreversível (onde a gota eventualmente para ao encontrar sua própria "esteira" menos molhável), o caso reversível com chemostats distintos permite uma auto-propulsão indefinida.

4. Resultados Numéricos e Análise

Os autores realizaram simulações numéricas usando o método de elementos finitos (biblioteca oomph-lib) em um domínio periódico unidimensional.

• Transição de Estados:
  • Caso Passivo ($\mu_a = \mu_c$): O sistema relaxa para um estado de equilíbrio termodinâmico, resultando em uma gota parada. Gradientes iniciais de concentração são dissipados.
  • Caso de Não-Equilíbrio ($\mu_a > \mu_c$): A gota desenvolve um gradiente de molhabilidade assimétrico e inicia o movimento. A adsorção ocorre sob a gota (aumentando a densidade de $B$) e a dessorção ocorre atrás dela, mantendo o gradiente.
• Bifurcação Drift-Pitchfork: A emergência do movimento é caracterizada por uma bifurcação drift-pitchfork supercrítica.
  • Existe um limiar crítico na diferença de potenciais químicos ($\mu_a - \mu_c$) ou nas taxas de reação.
  • Abaixo do limiar, a gota está parada. Acima do limiar, a velocidade da gota ($v$) escala com a raiz quadrada da distância ao limiar ($v \propto |\Delta \mu|^{1/2}$).
  • O movimento só ocorre se houver uma diferença finita de molhabilidade entre a região sob a gota e a região fora dela. Isso contrasta com gradientes externos impostos, onde qualquer gradiente minúsculo já causa movimento.
• Dependência das Taxas de Reação:
  • O movimento é maximizado para taxas de reação intermediárias.
  • Se as reações forem muito lentas, o gradiente não se forma suficientemente.
  • Se forem muito rápidas, o perfil de concentração torna-se simétrico (saturação), eliminando o gradiente de força motriz.
  • A difusividade do adsorvato ($D_b$) influencia o tipo de bifurcação (supercrítica vs. subcrítica), podendo levar a histerese e multistabilidade.

5. Significado e Impacto

• Fundamental: O trabalho estabelece uma base teórica rigorosa para a hidrodinâmica de filmes finos reativos, corrigindo inconsistências termodinâmicas de modelos anteriores. Ele demonstra como a quebra de simetria espacial e a manutenção de estados de não-equilíbrio podem ser descritas de forma unificada.
• Aplicação: O modelo é relevante para o entendimento de matéria ativa (active matter), incluindo gotas catalíticas, condensados biomoleculares e sistemas de entrega de fármacos autônomos.
• Experimental: O uso de chemostats permite traduzir parâmetros de controle experimental (como concentrações de reagentes no banho) diretamente para o modelo, facilitando a comparação com dados experimentais de gotas correntes.
• Futuro: O modelo é extensível para incluir substratos líquidos, efeitos de Marangoni e reações de ordem superior, abrindo caminho para o estudo de padrões dinâmicos complexos em sistemas reativos.

Em resumo, o artigo fornece uma descrição matemática robusta e termodinamicamente correta de como gotas podem se mover autonomamente em superfícies sólidas através de ciclos de reação reversíveis, identificando as condições exatas (bifurcações) para o início e a sustentação desse movimento.