Phase-Field Models for Particle-Stabilised Emulsions

Aquí presento un resumen técnico detallado del artículo "Phase-Field Models for Particle-Stabilised Emulsions" (Modelos de campo de fase para emulsiones estabilizadas por partículas), basado en el documento proporcionado.

1. El Problema

Las emulsiones estabilizadas por partículas (como las emulsiones Pickering y los geles de emulsión interfacialmente atascados o bijels) son fundamentales en la ciencia de la materia blanda. Sin embargo, su simulación computacional enfrenta limitaciones severas debido a las escalas espaciotemporales involucradas:

• Limitación de escalas: Los métodos tradicionales (como la dinámica molecular o la dinámica de partículas disipativas) resuelven individualmente cada partícula. Para sistemas con nanopartículas, esto restringe las simulaciones a escalas sub-micrométricas y sub-microsegundos.
• Complejidad del proceso STrIPS: En la formación de bijels mediante separación de fases inducida por transferencia de solvente (STrIPS), los fenómenos ocurren en escalas macroscópicas (cientos de micrómetros y segundos), involucrando miles de millones de nanopartículas. Resolver interacciones individuales en estas escalas es computacionalmente inviable.
• Fallas en modelos previos: Los intentos anteriores de aplicar modelos de campo de fase a estas emulsiones a menudo seguían resolviendo partículas individuales o, al eliminarlas, no lograban capturar adecuadamente la influencia de las partículas en la morfología y la tensión interfacial, generando inestabilidades numéricas (tensiones interfaciales negativas).

2. Metodología

Los autores proponen un nuevo marco teórico basado en la teoría de campo de fase que trata a las nanopartículas como un campo continuo en lugar de objetos discretos, permitiendo simulaciones eficientes a gran escala.

• Funcional de Energía Libre: Se define un funcional de energía libre que incluye contribuciones de los líquidos inmiscibles ($\phi$), las nanopartículas ($\psi$) y la interfaz. Se utiliza una expresión de Flory-Huggins para los líquidos y una aproximación de mezcla ideal para las partículas.
• Desacoplamiento de Potenciales Químicos:
  • Se identificó que las ecuaciones de Cahn-Hilliard estándar, donde el potencial químico del líquido depende directamente de la densidad de partículas, causan deformaciones no físicas de la interfaz y tensiones interfaciales negativas.
  • Solución: Se introdujeron ecuaciones dinámicas desacopladas parcialmente. El potencial químico del líquido se mantiene independiente de la densidad de partículas (asumiendo que las partículas no alteran la termodinámica del líquido a gran escala), mientras que las partículas se adsorben en la interfaz basada en el gradiente de concentración del líquido.
• Movilidad Dependiente de la Densidad (Atascamiento/Jamming):
  • Para modelar la estabilización, se introdujo una función de movilidad dependiente de la densidad de partículas ($\tilde{M}(\psi)$).
  • Cuando la densidad de partículas en la interfaz supera un umbral crítico ($\psi_c$), la movilidad cae drásticamente. Esto simula el "atascamiento" (jamming) de la capa de partículas, fijando la interfaz y deteniendo la separación de fases.
• Simulación de STrIPS: El modelo se acopló a un campo de solvente que difunde fuera del sistema, aumentando el parámetro de interacción $\chi$ y desencadenando la descomposición espinodal, replicando el proceso experimental STrIPS.

3. Contribuciones Clave

• Formalismo de Campo de Fase para Emulsiones: Desarrollo de un conjunto de ecuaciones diferenciales parciales acopladas que describen la dinámica de la separación de fases de líquidos y la adsorción/atascamiento de nanopartículas sin resolver partículas individuales.
• Resolución de Inestabilidades Numéricas: La propuesta de desacoplar los potenciales químicos evita la distorsión de la interfaz y las tensiones interfaciales negativas, permitiendo simulaciones estables con altas concentraciones de partículas.
• Modelado de la Cinética de Atascamiento: Implementación de una movilidad no lineal que captura la transición de un estado fluido a un estado "atascado" (jammed), lo cual es crucial para la formación de estructuras estables como los bijels.
• Integración con STrIPS: Extensión exitosa del modelo para incluir la dinámica de difusión de solvente, permitiendo estudiar la formación de bijels en condiciones experimentales realistas.

4. Resultados Principales

• Formación de Emulsiones: El modelo reproduce exitosamente tanto la formación de emulsiones Pickering (gotas dispersas) como de bijels (redes bicontinuas), dependiendo de la composición inicial del líquido.
• Mecanismo de Estabilidad: Se demostró que la adsorción de nanopartículas y el subsiguiente atascamiento de la interfaz suprimen eficazmente el maduración de Ostwald (crecimiento de gotas grandes a expensas de las pequeñas). La movilidad reducida en la interfaz impide el transporte de masa a través de ella.
• Influencia de la Concentración de Partículas:
  • Tamaño de Dominio: Existe una relación inversa entre la concentración inicial de nanopartículas y el tamaño promedio de los dominios. Mayores concentraciones de partículas provocan un atascamiento más rápido, deteniendo la separación de fases en etapas tempranas y produciendo estructuras con dominios más pequeños y monodispersos.
  • Gradiente en Bijels STrIPS: Las simulaciones 3D de STrIPS replicaron el gradiente característico de tamaño de dominio observado experimentalmente: los dominios son más pequeños cerca de la interfaz con el solvente (donde la difusión es rápida y el atascamiento ocurre temprano) y más grandes cerca del sustrato sólido (donde la difusión es lenta y el coarsening avanza más tiempo).
• Validación: Los resultados cualitativos coinciden con la evidencia experimental previa y simulaciones de partículas resueltas, validando el enfoque de campo de fase.

5. Significado e Impacto

Este trabajo representa un avance significativo en la simulación de sistemas de materia blanda complejos:

• Escalabilidad: Permite estudiar fenómenos en escalas macroscópicas (micrómetros a milímetros) que eran inaccesibles para métodos de resolución de partículas individuales.
• Diseño de Materiales: Proporciona una herramienta computacional eficiente para predecir y optimizar la morfología de emulsiones estabilizadas por partículas, crucial para aplicaciones en catálisis, almacenamiento de energía, separación de membranas y farmacéutica.
• Comprensión Mecanística: Aclara el papel cinético de las nanopartículas en la estabilización de interfaces, diferenciando entre la disponibilidad de partículas y la cinética de su adsorción/atascamiento.
• Futuro: Aunque el modelo actual omite la hidrodinámica (asumiendo dominancia difusiva), establece una base sólida que puede ser extendida para incluir efectos de convección y propiedades específicas de las partículas (forma, ángulo de contacto), abriendo nuevas vías de investigación en el diseño de materiales complejos.

En resumen, el modelo de campo de fase presentado ofrece un equilibrio óptimo entre precisión física y eficiencia computacional, permitiendo la exploración de la formación y morfología de emulsiones estabilizadas por partículas a escalas relevantes para la aplicación industrial y científica.