Molecular Dynamics Mechanism Study of Sodium Dodecyl Benzene Sulfonate Adsorption on Silica Nanopores

Resumen Técnico: Estudio del Mecanismo de Dinámica Molecular de la Adsorción de Dodecil Benceno Sulfonato de Sodio en Nanoporos de Sílice

Planteamiento del Problema
El dodecil benceno sulfonato de sodio (SDBS), un surfactante aniónico ampliamente utilizado en el desarrollo de campos de petróleo y gas, plantea desafíos significativos en el tratamiento de aguas residuales debido a su persistencia en el agua producida, lo que interfiere con los procesos de separación de aceite y agua. Aunque se reconoce que los nanoporos de sílice son adsorbentes polares eficaces para la eliminación de tales sustancias, el mecanismo de adsorción microscópica dentro de entornos confinados a escala nanométrica sigue siendo poco comprendido. Específicamente, existe una falta de marcos teóricos sistemáticos con respecto a las interacciones dinámicas entre las superficies de sílice y las moléculas de surfactante aniónico, así como a los mecanismos de formación y estabilidad de las capas de adsorción secundarias. Los métodos experimentales tradicionales presentan limitaciones para resolver estos detalles a nivel atómico.

Metodología
El estudio emplea simulaciones de dinámica molecular (MD) utilizando el paquete de software Materials Studio 2020 para investigar el comportamiento de adsorción del SDBS en nanoporos de sílice.

• Campo de Fuerza: Se utilizó el campo de fuerza COMPASS III para calcular la energía potencial, incluyendo interacciones de enlace, ángulo, diedro y no ligadas (electrostáticas y de van der Waals).
• Construcción del Modelo: Se construyó un modelo microscópico 3D de una solución acuosa de SDBS utilizando el módulo Amorphous Cell, disociando el SDBS en aniones de dodecilbencenosulfonato (SDB-) e iones de sodio (Na+). Los modelos de nanoporos de sílice fueron generados y optimizados.
• Modificación de la Superficie: Para explorar los efectos de las propiedades superficiales, se aplicó una técnica de "hidroxilación completa" a los nanocanales de sílice. Esto implicó la eliminación de átomos de silicio y oxígeno no ligados dentro de un radio de 10 Å y la adición de átomos de hidrógeno para crear una superficie cubierta con grupos silanol.
• Parámetros de Simulación: Las simulaciones se realizaron bajo el ensamble NVT durante 5 ns utilizando el módulo Forcite. El control de temperatura se gestionó mediante un termostato de Nose-hoover. Las interacciones electrostáticas se calcularon mediante el método de sumación de Ewald, mientras que las interacciones de van der Waals utilizaron el método basado en átomos.
• Variables: El estudio analizó los efectos de la hidroxilación de la superficie, la temperatura del sistema (288 K, 303 K, 323 K, 343 K) y la densidad de la solución (en un rango de 1.016 a 1.025 g/cm³).
• Métricas de Análisis: Los indicadores clave incluyeron la Función de Distribución Radial (RDF), el Desplazamiento Cuadrático Medio (MSD), el radio de rotación (radio de giro), la distribución de densidad y los contornos de potencial electrostático.

Resultados Clave

1. Impacto de la Hidroxilación: La hidroxilación completa de la superficie aumentó significativamente la capacidad de adsorción de los nanoporos de sílice para el SDBS. El área superficial accesible aumentó aproximadamente 1190 Ų (superficie del solvente) y 2823.75 Ų (superficie de van der Waals).
   • Mecanismo: En superficies no hidroxiladas (hidrofóbicas), las moléculas de SDBS tendían a permanecer planas con las cadenas alquílicas en contacto con la superficie. En superficies hidroxiladas (hidrofílicas), las moléculas adoptaron una conformación más erguida con los grupos de ácido sulfónico orientados hacia la superficie y las colas hidrofóbicas extendiéndose hacia la solución.
   • Dinámica: La hidroxilación aumentó la flexibilidad y la capacidad de difusión de las largas cadenas de SDB-. El análisis de MSD reveló un modo de "anclaje-oscilación" donde el grupo de cabeza permanecía firmemente fijo mientras la cadena de la cola mantenía la movilidad.
2. Efectos de la Temperatura: La capa de adsorción fue más completa y ordenada entre 303 K y 323 K.
   • Por debajo o por encima de este rango, la perturbación térmica o la insuficiencia de energía cinética redujeron el orden.
   • El aumento de la temperatura aumentó la flexibilidad y la diversidad conformacional de las largas cadenas de SDB-, pero redujo su estabilidad dinámica, causando transiciones frecuentes entre estados enrollados y estirados.
   • Crucialmente, la formación de una segunda capa de adsorción depende de la temperatura; era estable y distinta a 303 K y 323 K, pero desapareció a 288 K y 343 K.
3. Efectos de la Densidad de la Solución: El aumento de la densidad de la solución (hasta 1.025 g/cm³) incrementó el número de moléculas adsorbidas tanto en la primera como en la segunda capa sin alterar las posiciones de los picos del RDF (lo que indica posiciones de capa fijas). Una mayor densidad promovió la extensión de las largas cadenas de SDB- y mejoró su estabilidad cinética al forzar una conformación más unificada y extendida.
4. Características de la Segunda Capa de Adsorción: La segunda capa de adsorción es menos ordenada que la primera, caracterizada por picos asimétricos en el RDF y una distribución molecular difusa.
   • En superficies hidroxiladas, la segunda capa forma una doble capa ordenada donde las moléculas se adhieren mediante interacciones hidrofóbicas de extremo de cola a la primera capa erguida.
   • En superficies no hidroxiladas, la segunda capa es desordenada, formada por interacciones hidrofóbicas cadena-cadena de moléculas que yacen de forma plana.

Significancia y Reivindicaciones
Los autores afirman que este estudio proporciona una base teórica microscópica para comprender la eliminación de SDBS y sustancias similares de la solución. Al utilizar la dinámica molecular, la investigación supera las limitaciones de los métodos experimentales tradicionales para visualizar las interacciones dinámicas a escala nanométrica.

El artículo postula que:

• La hidroxilación de la superficie es una estrategia viable y eficaz para mejorar la eficiencia de adsorción de los nanoporos de sílice para surfactantes aniónicos.
• Existe una ventana termodinámica específica (303–323 K) donde la adsorción es óptima y la formación de la capa secundaria es estable.
• El modo dinámico de "anclaje-oscilación" y los cambios conformacionales específicos de las cadenas de surfactante son factores críticos en el mecanismo de adsorción.
• La comprensión de estos factores permite la optimización de las condiciones del proceso y el desarrollo de nuevos materiales funcionales para el tratamiento de aguas residuales en la industria del petróleo y el gas.

El estudio concluye que, si bien la adsorción primaria está dominada por la monocapa, la formación y estabilidad de una capa secundaria son altamente sensibles a la química de la superficie, la temperatura y la densidad de la solución, ofreciendo objetivos específicos para el diseño de materiales y el control de procesos.