Quasi-two-dimensional dispersions of Brownian particles with competitive interactions: Dynamical clustering, non-Gaussianity and hydrodynamic correlations

Resumen Técnico: Dispersiones cuasi-bidimensionales de partículas brownianas con interacciones competitivas

Planteamiento del problema
El artículo aborda el comportamiento dinámico de dispersiones cuasi-bidimensionales (Q2D) de partículas brownianas sujetas a interacciones competitivas de corto alcance atractivas (SA) y de largo alcance repulsivas (LR), conocidas como sistemas SALR. Si bien las propiedades estructurales de estos sistemas —específicamente la formación de cúmulos de tamaño finito y microfases— están bien documentadas tanto en contextos tridimensionales (3D) como Q2D, la comprensión sistemática de su dinámica sigue incompleta. Las cuestiones clave abiertas incluyen:

1. ¿Cómo influyen las interacciones hidrodinámicas (HIs) en la estabilidad de los cúmulos y en el transporte colectivo en geometrías confinadas Q2D?
2. ¿Cómo se manifiesta la interacción entre la formación de cúmulos y la difusión colectiva anómala característica de los sistemas Q2D (donde $D(q)$ y $H(q)$ divergen cuando $q \to 0$)?
3. ¿Cuáles son los mecanismos que gobiernan la dinámica no gaussiana y la heterogeneidad dinámica en estas fases agrupadas?
4. ¿Cómo afectan los efectos inerciales de corto tiempo y el retraso hidrodinámico a la formación de correlaciones colectivas en las fases de cúmulos en comparación con las referencias de esferas duras?

Metodología
Los autores emplean un enfoque comparativo utilizando dos técnicas de simulación distintas para aislar los efectos de la hidrodinámica:

1. Dinámica de Langevin (LD): Utilizada para modelar sistemas donde las interacciones hidrodinámicas están ausentes. Esto sirve como referencia para la dinámica sobreamortiguada.
2. Dinámica de Colisión Multipartícula (MPC): Utilizada para resolver explícitamente la hidrodinámica fluctuante, capturando tanto los regímenes inercial como el sobreamortiguado. El solvente se modela como partículas puntuales que experimentan rotación estocástica (SRD) dentro de celdas de colisión, asegurando la conservación de la masa, el momento y la energía.

El sistema consiste en $N_p$ partículas esféricas idénticas confinadas a una monocapa planar dentro de un solvente Newtoniano 3D. El potencial de interacción combina una atracción de Lennard-Jones generalizada de pendiente pronunciada con una repulsión electrostática (Yukawa) de pantalla. El estudio varía sistemáticamente la fuerza de atracción ($\epsilon^*$) y la fracción de área ($\phi_{2D}$) para mapear el comportamiento de fase.

Herramientas analíticas clave:

• Funciones de correlación de enlaces ($B(t)$) y análisis de la vida media del cúmulo ($\tau_b$) para caracterizar el arresto dinámico.
• Parámetro de orden hexagonal local ($q_6$) y su función de correlación temporal ($g_6(t)$) para evaluar la relajación estructural.
• Parámetros no gaussianos ($\alpha_2$) y la función de van Hove de autocorrección para cuantificar la heterogeneidad dinámica.
• Funciones hidrodinámicas: Tanto la función estática $H(q)$ como la dependiente del tiempo $H(q,t)$, descompuestas en partes propia (self) y distinta (distinct), para investigar la difusión colectiva y las correlaciones de momento.
• Funciones de correlación de corriente longitudinal ($J(q,t)$) para investigar la propagación de los modos de sonido.

Contribuciones clave y resultados

• Comportamiento de fase y termalización: El estudio confirma cuatro fases generalizadas: Fluido Disperso (DF), Cúmulos en Equilibrio (EC), Percolación Aleatoria (RP) y Percolación de Cúmulos (CP). La termalización se verificó, mostrando que los estados estacionarios se alcanzan independientemente de las condiciones iniciales, excepto a temperaturas muy bajas donde ocurren derivas de energía lentas.

• Autodifusión y subdifusión: A medida que aumenta la fuerza de atracción, la autodifusión se ve significada y considerablemente suprimida.

  • En las fases fluidas (DF, EC, RP), las partículas exhiben una difusión browniana distinta de corto y largo tiempo, ralentizada por el agrupamiento transitorio en lugar del encajonamiento (caging).
  • En la fase de Percolación de Cúmulos (CP), emerge una relajación de dos pasos, que corresponde a dos regímenes subdifusivos distintos, coincidiendo con la formación de un gel.

• Vida media de los cúmulos y criterios de transición:

  • Se encontró que la vida media de los cúmulos ($\tau_b$) está influenciada más fuertemente por la fuerza de atracción que por la concentración de partículas.
  • Se identificaron dos criterios dinámicos para la transición de fases no agrupadas a agrupadas:
    1. Un umbral en la vida media media del cúmulo: $\tau_b / \tau_D \approx 10$ separa las fases DF y EC.
    2. El tiempo de relajación del orden hexagonal local ($\tau_6$): El cruce $\tau_6 \approx \tau_D$ marca la transición DF–EC.
  • En la fase EC, la vida media de los cúmulos se alinea con la teoría de cruce de barreras de Kramers, mientras que en el régimen DF, siguen la teoría de disociación de Smoluchowski.

• Dinámica no gaussiana: Los sistemas Q2D-SALR con cúmulos exhiben una pronunciada dinámica no gaussiana. La función de van Hove de autocorrección en la fase de cúmulos en equilibrio muestra una forma aproximadamente exponencial, consistente con un mecanismo de "difusión de la difusividad". Esto indica una heterogeneidad dinámica significativa donde las partículas en los cúmulos tienen una movilidad restringida, mientras que aquellas en la fase diluida permanecen móviles.

• Interacciones Hidrodinámicas (HIs) y Difusión Colectiva:

  • Las simulaciones MPC revelan que la difusión colectiva a gran escala anomalamente aumentada, característica de los sistemas Q2D con interacciones hidrodinámicas, se preserva cualitativamente en las dispersiones Q2D-SALR.
  • Sin embargo, la presencia de agrupamiento suprime el pico de orden de rango intermedio en la función hidrodinámica en comparación con sus contrapartes 3D.
  • Crucialmente, el análisis de las correlaciones hidrodinámicas dependientes del tiempo ($H(q,t)$) y los modos de sonido indica que el agrupamiento conduce a un inicio más temprano de las interacciones hidrodinámicas que en los sistemas de referencia de esferas duras Q2D. Esto implica que las HIs se vuelven relevantes en las escalas de tiempo inerciales en los sistemas SALR agrupados.

Significancia y Reivindicaciones
El artículo afirma proporcionar un análisis dinámico exhaustivo que une la separación de microfases estructurales con las propiedades de transporte en sistemas SALR confinados. Al combinar LD y MPC, los autores demuestran que:

1. La fuerza de atracción es el parámetro de control dominante para el agrupamiento dinámico, superando los efectos de la concentración de partículas.
2. Las interacciones hidrodinámicas juegan un papel crítico en la dinámica colectiva de las dispersiones Q2D-SALR, preservando la difusión anómala de largo alcance mientras modifican las características espectrales de la función hidrodinámica.
3. El agrupamiento acelera la relevancia de la hidrodinámica, causando que las HIs se manifiesten en las escalas de tiempo inerciales, un hallazgo que desafía la suposición de HIs instantáneas en las descripciones teóricas de tales sistemas.
4. El estudio ofrece criterios dinámicos (basados en la vida media del cúmulo y la relajación del orden hexagonal) para distinguir entre fases dispersas y agrupadas, complementando los criterios estructurales existentes.

El trabajo subraya la necesidad de considerar tanto las interacciones competitivas como el acoplamiento hidrodinámico para comprender plenamente los mecanismos microscópicos subyacentes al agrupamiento dinámico y al transporte colectivo en dispersiones coloidales confinadas.