Topological Signatures of Diffusive Release in Porous Media

Resumo Técnico: Assinaturas Topológicas de Liberação Difusiva em Meios Porosos

Definição do Problema
A liberação difusiva em meios porosos é governada pela geometria em escala de poros, incluindo conectividade, tortuosidade e o arranjo espacial de obstáculos sólidos. Embora simulações numéricas mecanísticas (ex: métodos de elementos finitos) sejam o padrão para estudar essas relações, elas são computacionalmente dispendiosas, particularmente ao realizar a triagem de um grande número de estruturas candidatas para aplicações que variam desde a entrega de fármacos até o transporte de contaminantes. Há uma necessidade de descritores estruturais interpretáveis e baseados em geometria que possam prever perfis de liberação qualitativos sem o custo de simulações completas de difusão.

Metodologia
O estudo utiliza a homologia persistente, uma ferramenta da análise de dados topológicos (TDA), para quantificar a topologia multiescala da fase sólida em meios porosos.

1. Configuração da Simulação:

   • Meios porosos sintéticos foram gerados dentro de um cubo unitário utilizando obstáculos esféricos com variados tamanhos e arranjos (seis classes estruturais: aleatório monodisperso, aleatório polidisperso, tipo canal, tipo camada, tipo gargalo e tipo agrupado).
   • As amostras foram geradas para corresponder a porosidades alvo de 40%, 60% e 80%.
   • A liberação difusiva foi modelada como um problema misto de Dirichlet–Neumann em um espaço de poros voxelizado. A curva de liberação cumulativa $Q(t)$ foi computada usando um método de elementos finitos (FEM) conforme, derivado da grade de voxels.

2. Extração de Características Topológicas:

   • Em vez de analisar o espaço de poros diretamente, os autores analisaram a fase sólida ($S$) através de um processo de filtração. A fase sólida foi gradualmente expandida (offset) por um raio $r$, criando uma sequência de conjuntos $F_r = \{x \in \Omega : \text{dist}(x, S) \leq r\}$.
   • A homologia persistente foi computada nesta filtração para rastrear o nascimento e a morte de características topológicas: componentes conectados ($H_0$), loops ($H_1$) e vazios/cavidades ($H_2$).
   • Seis estatísticas de resumo foram extraídas para cada dimensão de homologia (contagem, contagem finita, soma de persistência, persistência média, persistência máxima e escala de nascimento média), resultando em 18 características topológicas por amostra.

3. Análise e Classificação:

   • Definição de Regime: As curvas de liberação foram categorizadas em três regimes baseados em métricas de pontos temporais específicos: "rápida inicial" (alta liberação inicial), "cauda longa" (alta fração não liberada ao final) e "liberação tardia" (aumento significativo na liberação durante a metade final da simulação).
   • Associação Estatística: Testes de Kruskal–Wallis foram realizados para determinar se os resumos topológicos diferiam significativamente entre os regimes de liberação, estratificados por porosidade alvo para isolar efeitos topológicos de efeitos de porosidade.
   • Classificação: Um modelo de regressão logística multinomial foi treinado utilizando apenas os 18 atributos topológicos para classificar as amostras nos três regimes de liberação.
   • Tempo Computacional: Os tempos de execução (wall-clock times) foram comparados entre o solver de difusão FEM e o pipeline de extração de características de homologia persistente.

Principais Resultados

• Assinaturas Topológicas Independentes de Porosidade: Mesmo controlando a porosidade total, amostras com diferentes regimes de liberação exibiram assinaturas topológicas estatisticamente distintas.

  • A liberação rápida inicial foi associada a contagens de persistência menores e somas de persistência menores (topologia da fase sólida mais simples).
  • A liberação de cauda longa foi associada a contagens de persistência maiores e somas de persistência maiores (organização multiescala mais complexa de componentes sólidos, loops e cavidades).
  • A "soma de persistência" emergiu como o sinal mais estável e discriminativo em todos os níveis de porosidade.
  • O poder discriminativo das características topológicas diminuiu à medida que a porosidade aumentou (ex: a 80% de porosidade, a fase sólida é mais esparsa, tornando as características de maior dimensão menos distintas).

• Desempenho de Classificação: Utilizando apenas os 18 atributos topológicos, o modelo de regressão logística alcançou acurácias de teste entre 0,64 e 0,76 (significativamente acima do nível de chance de 0,33 para uma tarefa de três classes balanceadas). O desempenho foi maior a 40% de porosidade e declinou em porosidades mais altas, consistente com os resultados de associação estatística.

• Eficiência Computacional: A extração de características por homologia persistente foi substancialmente mais rápida que as simulações de difusão FEM.

  • Para resolução de grade $N=16$, o FEM foi de 6,6 a 64,7 vezes mais lento que a extração de PH.
  • Para resolução de grade $N=32$, o FEM foi de 81,7 a 631,0 vezes mais lento.
  • Os tempos de extração de PH permaneceram abaixo de 0,018 segundos por amostra, enquanto os tempos de FEM variaram até ~3 segundos.

Significância e Alegações
O artigo afirma que a homologia persistente fornece um descritor leve, interpretável e baseado em geometria para triagem de comportamento de liberação difusiva. As principais contribuições são:

1. Demonstrar que o comportamento de liberação depende não apenas do volume do espaço de poros (porosidade), mas criticamente da organização multiescala da fase sólida.
2. Mostrar que resumos topológicos podem distinguir efetivamente diferentes regimes de liberação (rápida inicial, liberação tardia, cauda longa) sem resolver a equação de difusão.
3. Estabelecer que essas características topológicas podem ser extraídas a uma fração do custo computacional das simulações tradicionais de FEM, tornando-as adequadas para a priorização inicial de estruturas candidatas antes de simulações detalhadas ou avaliação experimental.

Os autores posicionam este método como uma ferramenta complementar para triagem, e não como um substituto para simulações baseadas em física, destacando potenciais aplicações em sistemas farmacêuticos (ex: comprimidos de liberação controlada) onde a arquitetura de poros dita os perfis de liberação de fármacos.