Propiedades de transporte de co-solutos en hidrogeles termosensibles

Abstract

En este trabajo presentamos un estudio acerca de la transferencia de masa de diferentes cosolutos a través de partículas de hidrogeles de poli N-Isopropilacrilamida (pNIPAM) en soluciones acuosas. Estos cosolutos son: fenol, 5- Fluorouracilo (5FU), metano y metanol. El estudio se aborda desde distintas perspectivas: combinando cálculos teóricos, simulaciones computacionales de dinámica molecular átomo-por-átomo y grano-grueso, y comparando con reportes experimentales. Las partículas de hidrogel son diseñadas con estructuras rellenas o huecas, y en particular estudiamos los procesos de encapsulación/liberación de fenol hacia/desde ellas. En la primera etapa del proyecto diseñamos un modelo original de tipo grano-grueso para geles de pNIPAM, el cual captura de forma adecuada la transición bobina-glóbulo y glóbulo-bobina del polímero. Con este modelo construimos arreglos poliméricos complejos aproximados a un hidrogel, como lo es una membrana, con los cuales se observa una mejora en el consumo de recursos de las simulaciones de hasta 54 veces. Por otra parte, los cálculos teóricos son abordados mediante la teoría del funcional de la densidad dinámico (DDFT, por sus siglas en inglés), a partir de la cual obtenemos los perfiles de densidad del cosoluto y su evolución en el tiempo. La DDFT es una técnica que permite el estudio a nivel molecular de la dinámica de sistemas Fickianos en los que los modos hidrodinámicos no son relevantes, mediante la única hipótesis de que podemos aproximar las correlaciones dinámicas por las que existirían en un sistema equivalente (con igual perfil de densidad) en equilibrio. A diferencia de otros modelos que consideran una difusión ideal, la DDFT toma en cuenta la interacción efectiva entre el polímero y el cosoluto (U ef) y entre moléculas de cosoluto, para calcular el potencial químico de la carga dentro del hidrogel, y a partir de esto, su flujo. Además, la DDFT requiere de una definición explícita del coeficiente de difusión efectivo del cosoluto dentro de la matriz de hidrogel (D ef ). Estos parámetros se extraen de simulaciones computacionales tipo átomo-por-átomo. La interacción U ef se obtiene a partir de la energía libre de trasporte, la cual se calcula con el método de relación de aceptación de Bennett (BAR en inglés). Por su parte, el D ef se calcula a partir del desplazamiento cuadrático medio de las partículas de prueba embebidas en el hidrogel y la relación de Stokes-Einstein. La validez y la implementación del método de resolución de la DDFT es probada comparando sus resultados contra simulaciones de dinámica molecular átomo-por-átomo en condiciones análogas, esto es, para la absorción de fenol en un hidrogel hueco de pNIPAM. Las predicciones teóricas y las simulaciones muestran resultados casi idénticos para los perfiles de distribución de carga y velocidades de absorción. Finalmente, realizamos una comparación entre los resultados de la DDFT con reportes experimentales obtenidos de la literatura para laliberación de fenol desde hidrogeles completamente rellenos. La longitud típica de escalado con que normalizamos las unidades en la DDFT es ajustada a las dimensiones del microgel experimental, al igual que las concentraciones iniciales. Las curvas de liberación de cosolutos obtenidas con DDFT coinciden satisfactoriamente (dentro de la misma escala temporal) con los experimentos.