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

Resumo Técnico: Estudo do Mecanismo de Dinâmica Molecular da Adsorção de Dodecil Benzeno Sulfonato de Sódio em Nanoporos de Sílica

Definição do Problema
O dodecil benzeno sulfonato de sódio (SDBS), um surfactante aniônico amplamente utilizado no desenvolvimento de campos de petróleo e gás, apresenta desafios significativos para o tratamento de águas residuais devido à sua persistência na água produzida, o que interfere nos processos de separação óleo-água. Embora os nanoporos de sílica sejam reconhecidos como adsorventes polares eficazes para a remoção de tais substâncias, o mecanismo de adsorção microscópica dentro de ambientes confinados em nanoescala permanece pouco compreendido. Especificamente, há uma carência de estruturas teóricas sistemáticas relativas às interações dinâmicas entre superfícies de sílica e moléculas de surfactantes aniônicos, bem como aos mecanismos de formação e estabilidade de camadas de adsorção secundárias. Os métodos experimentais tradicionais possuem limitações para resolver esses detalhes em nível atômico.

Metodologia
O estudo emprega simulações de dinâmica molecular (MD) utilizando o pacote de software Materials Studio 2020 para investigar o comportamento de adsorção do SDBS em nanoporos de sílica.

• Campo de Força: O campo de força COMPASS III foi utilizado para calcular a energia potencial, incluindo interações de ligação, ângulo, diedro e não ligantes (eletrostáticas e de van der Waals).
• Construção do Modelo: Um modelo microscópico 3D de uma solução aquosa de SDBS foi construído usando o módulo Amorphous Cell, dissociando o SDBS em ânions dodecil benzeno sulfonato (SDB-) e íons de sódio (Na+). Modelos de nanoporos de sílica foram gerados e otimizados.
• Modificação de Superfície: Para explorar os efeitos das propriedades superficiais, uma técnica de "hidroxilação completa" foi aplicada aos canais nanoestruturados de sílica. Isso envolveu a remoção de átomos de silício e oxigênio não ligados dentro de um raio de 10 Å e a adição de átomos de hidrogênio para criar uma superfície coberta por grupos silanol.
• Parâmetros de Simulação: As simulações foram conduzidas sob o ensemble NVT por 5 ns usando o módulo Forcite. O controle de temperatura foi gerenciado via termostato Nose-hoover. As interações eletrostáticas foram calculadas usando o método de soma de Ewald, enquanto as interações de van der Waals utilizaram o método baseado em átomos.
• Variáveis: O estudo analisou os efeitos da hidroxilação da superfície, temperatura do sistema (288 K, 303 K, 323 K, 343 K) e densidade da solução (variando de 1,016 a 1,025 g/cm³).
• Métricas de Análise: Os indicadores principais incluíram a Função de Distribuição Radial (RDF), o Deslocamento Médio Quadrático (MSD), o raio de rotação (raio de giração), distribuição de densidade e contornos de potencial eletrostático.

Principais Resultados

1. Impacto da Hidroxilação: A hidroxilação completa da superfície aumentou significativamente a capacidade de adsorção dos nanoporos de sílica para o SDBS. A área de superfície acessível aumentou aproximadamente 1190 Ų (superfície do solvente) e 2823,75 Ų (superfície de van der Waals).
   • Mecanismo: Em superfícies não hidroxiladas (hidrofóbicas), as moléculas de SDBS tendem a repousar de forma plana, com as cadeias alquílicas em contato com a superfície. Em superfies hidroxiladas (hidrofílicas), as moléculas adotam uma conformação mais vertical, com os grupos de ácido sulfônico voltados para a superfície e as caudas hidrofóbicas estendendo-se para a solução.
   • Dinâmica: A hidroxilação aumentou a flexibilidade e a capacidade de difusão das longas cadeias de SDB-. A análise de MSD revelou um modo de "ancoragem-oscilação", onde o grupo cabeça permanece firmemente fixo enquanto a cadeia da cauda mantém a mobilidade.
2. Efeitos da Temperatura: A camada de adsorção foi mais completa e ordenada entre 303 K e 323 K.
   • Abaixo ou acima desta faixa, a perturbação térmica ou a energia cinética insuficiente reduziu a ordem.
   • O aumento da temperatura aumentou a flexibilidade e a diversidade conformacional das longas cadeias de SDB-, mas reduziu sua estabilidade dinâmica, causando transições frequentes entre estados enrolados e esticados.
   • Crucialmente, a formação de uma segunda camada de adsorção é dependente da temperatura; ela é estável e distinta a 303 K e 323 K, mas desaparece a 288 K e 343 K.
3. Efeitos da Densidade da Solução: O aumento da densidade da solução (até 1,025 g/cm³) aumentou o número de moléculas adsorvidas tanto na primeira quanto na segunda camada, sem alterar as posições dos picos da RDF (indicando posições de camada fixas). A maior densidade promoveu a extensão das longas cadeias de SDB- e melhorou sua estabilidade cinética ao forçar uma conformação mais unificada e estendida.
4. Características da Segunda Camada de Adsorção: A segunda camada de adsorção é menos ordenada que a primeira, caracterizada por picos assimétricos na RDF e uma distribuição molecular difusa.
   • Em superfícies hidroxiladas, a segunda camada forma uma camada dupla ordenada, onde as moléculas se prendem via interações hidrofóbicas de extremidade de cauda à primeira camada vertical.
   • Em superfícies não hidroxiladas, a segunda camada é desordenada, formada por interações hidrofóbicas cadeia-cadeia de moléculas que repousam de forma plana.

Significância e Alegações
Os autores afirmam que este estudo fornece uma base teórica microscópica para compreender a remoção de SDBS e substâncias similares de soluções. Ao utilizar a dinâmica molecular, a pesquisa supera as limitações dos métodos experimentais tradicionais na visualização de interações dinâmicas em nanoescala.

O artigo postula que:

• A hidroxilação da superfície é uma estratégia viável e eficaz para aumentar a eficiência de adsorção dos nanoporos de sílica para surfactantes aniônicos.
• Existe uma janela termodinâmica específica (303–323 K) onde a adsorção é ótima e a formação da camada secundária é estável.
• O modo dinâmico de "ancoragem-oscilação" e as mudanças conformacionais específicas das cadeias de surfactante são fatores críticos no mecanismo de adsorção.
• A compreensão desses fatores permite a otimização das condições de processo e o desenvolvimento de novos materiais funcionais para o tratamento de águas residuais na indústria de petróleo e gás.

O estudo conclui que, embora a adsorção primária seja dominada pela monocamada, a formação e a estabilidade de uma camada secundária são altamente sensíveis à química da superfície, temperatura e densidade da solução, oferecendo alvos específicos para o design de materiais e controle de processos.