Propriedades de Superfície Governam Transporte de Óleo em Nano-canais: Insights de Dinâmica Molecular

Contexto

Com demanda global de energia aumentando e recursos convencionais esgotando, óleo de xisto não-convencional atraiu intenso interesse de pesquisa por causa de suas vastas reservas [1,2]. Apesar do fato de que recursos de óleo de xisto excedem óleo convencional por mais de três vezes [2], a fração explorável permanece limitada devido ao seu confinamento dentro de poros nanoscópicos variando de 2 a 100 nm [4,5]. Este confinamento induz efeitos de superfície pronunciados que alteram propriedades do fluido — como densidade, molhabilidade e difusão — levando a comportamentos de transporte que diferem marcadamente de sistemas macroscópicos [10-12]. Por exemplo, água flui mais rápido em nanotubos de carbono que em tubos macroscópicos equivalentes [6,7], e o caráter hidrofílico/hidrofóbico de nanotubos de carbono muda com diâmetro [9]. Trabalhos prévios de dinâmica molecular mostraram que tensão de cisalhamento e tamanho de poro influenciam dinâmica de água e óleo em nanofendas de carbono e sílica [13-15]. No entanto, estudos sistemáticos sobre como química de superfície e topografia governam transporte de óleo de xisto são escassos. Entender estes efeitos é crucial para otimizar tecnologias de extração e tem implicações mais amplas para separações de membrana, transporte biomolecular e dispositivos de armazenamento de energia [16-20].

Métodos

Todas as simulações foram realizadas usando o módulo Discover do Materials Studio (Accelrys Inc.) com o campo de força COMPASS all-atom, que captura accuradamente interações ligadas, cross-term e não-ligadas em fases condensadas:

$$ {E}_{\mathrm{total}}={E}_{\mathrm{valence}}+{E}_{\mathrm{cross}\text{-}term}+{E}_{\mathrm{nonbond}} \tag{1} $$

Expressões completas para cada termo são dadas na publicação original. Forças de van der Waals foram modeladas com potencial Lennard-Jones (cut-off 15,5 Å), e eletrostática foi tratada via método Ewald. Simulações rodaram no ensemble NVT a 298 K usando termostato Andersen. Condições de contorno periódicas foram aplicadas em todas as direções. Dados foram amostrados a cada 5 ps.

O substrato de sílica, representando o mineral dominante em xisto, foi construído da superfície (001) e hidroxilado para mimetizar condições geológicas [26]. Duas superfícies opostas, cada uma 1,5 × 7 × 0,85 nm³, criaram um canal de fenda com larguras de 2, 4 e 6 nm. Octadecano, piridina e fenol — constituintes comuns de óleo de xisto — foram introduzidos para alcançar densidade de massa de 0,8 g cm⁻³ (40 moléculas no canal de 2 nm; 407, 344, 80 e 120 moléculas para piridina, fenol e octadecano nos canais de 4 e 6 nm, respectivamente). Após minimização de energia e equilíbrio de 500 ps, uma força de corpo constante equivalente a 3,1×10⁻¹⁴ N por átomo foi aplicada ao longo do eixo do canal para impulsionar as moléculas, uma abordagem padrão em estudos de transporte MD [27-29].

Resultados e Discussão

Efeito da Largura do Canal

Sob a força imposta, o número de moléculas cruzando o canal aumentou constantemente (Informação Suplementar, Fig. S1). A Figura 2 demonstra que canais mais largos produzem maiores deslocamentos COM após 2 ns. A posição COM é definida como

$$ {z}_{\mathrm{COM}}={\displaystyle \sum_i\frac{m_i}{M}{r}_i} \tag{5} $$

No canal de 2 nm, o COM moveu apenas 0,85 nm, enquanto no canal de 6 nm o deslocamento excedeu 30 nm. Esta diferença pronunciada surge de adsorção mais forte em poros mais estreitos, como quantificado pela energia média de interação

$$ {E}_{\mathrm{average\ interaction}}=\frac{E_{\mathrm{total}}-(E_{\mathrm{oil}}+E_{\mathrm{substrate}})}{N} \tag{6} $$

A Figura 3 mostra uma correlação inversa clara: maior energia de interação leva a deslocamento reduzido, com aumento triplo na energia causando queda >30 vezes no deslocamento COM. Efeitos de interação diminuem à medida que o canal alarga, ressaltando a natureza dependente de tamanho do transporte nanoscópico.

Perfis de densidade em camadas (Fig. 2) revelam uma camada interfacial de 5-Å de espessura adjacente às paredes. Orientação do octadecano dentro destas camadas é predominantemente paralela à superfície (θ ≈ 80°-90°), enquanto moléculas no centro do canal exibem orientações aleatórias. Coeficientes de difusão, derivados de deslocamentos quadráticos médios (Eq. 7), variam com posição radial: são mais altos no centro e diminuem em direção à parede (Fig. 5a). O canal de 6 nm exibe perfil parabólico íngreme, enquanto o canal de 2 nm mostra difusão quase uniforme, baixa, assemelhando-se a um pistão. Estas tendências podem ser atribuídas à competição entre adsorção na parede (que suprime movimento na região de "sombra" próxima à superfície) e transferência de momento inter-molecular.

Efeito da Polaridade

Simulações em canal de sílica de 4 nm com números comparáveis de fenol, piridina e octadecano ilustram o papel da polaridade molecular. Apesar de forças idênticas, o deslocamento COM do octadecano é ~16× maior que o de fenol ou piridina (Fig. 8). Energias totais de interação (vdW + eletrostática) são marcadamente maiores para espécies polares, com eletrostática contribuindo aproximadamente metade do total (Fig. 8). Momentos dipolares calculados via métodos de primeiros princípios confirmam que octadecano (0,032 Debye) é virtualmente não-polar, enquanto fenol (1,306 Debye) e piridina (2,245 Debye) possuem dipolos significativos. Consequentemente, moléculas polares experimentam ligação mais forte à superfície de sílica, reduzindo mobilidade.

Efeito da Composição do Material

Comparações entre canais de sílica, ouro e calcita (Fig. 9) mostram que ouro imobiliza fortemente octadecano, enquanto calcita permite deslocamento moderado. A Figura 10 indica que o deslocamento COM correlaciona com energia média de interação: sílica exibe menor energia e maior deslocamento, ouro o oposto, e calcita intermediário. Energias similares entre sílica e calcita ainda com deslocamentos diferentes sugerem que tipo de átomo de superfície e textura também influenciam transporte.

Efeito da Rugosidade de Superfície

Introduzindo cavidades na nanoescala (3 Å ou 6 Å profundas) na superfície de sílica altera fluxo local. Cavidades abrigam moléculas de óleo aprisionadas, que reduzem velocidade de camadas de fluido adjacentes (Fig. 11). Deslocamentos COM para cavidades de 3-Å e 6-Å são 3,95 Å e 3,07 Å, respectivamente — ligeiramente maiores que o valor de superfície plana (3,17 Å) para cavidade rasa, mas menores para cavidade mais profunda. O realce com cavidades rasas provavelmente resulta de alargamento efetivo do canal, enquanto cavidades mais profundas impõem adsorção mais forte e impedem fluxo.

Conclusões

Este estudo de dinâmica molecular demonstra que largura de nano-canal, polaridade do óleo, material do substrato e rugosidade de superfície ditam coletivamente transporte de óleo. O canal de 6 nm proporciona deslocamento COM >30 vezes maior que o canal de 2 nm, e o coeficiente de difusão no centro do canal quase dobra aquele próximo à parede. Moléculas polares e substratos de alta interação suprimem marcadamente mobilidade, enquanto cavidades superficiais rasas podem modestamente melhorar transporte. Estes insights são essenciais para projetar sistemas de extração nano-estruturados e para aplicações mais amplas onde transporte de fluido em geometrias confinadas é crítico.

Abreviações

COM:

Centro de massa

COMPASS:

Potencial molecular otimizado para fase condensada para estudos de simulação atomística

MD:

Dinâmica molecular

vdW:

van der Waals