Derivation of nonlinear time-dependent macroscopic conductivity for an electropermeabilization model via homogenization

Resumo Técnico: Derivação da Condutividade Macroscópica Não Linear e Dependente do Tempo para um Modelo de Eletropereabilização via Homogeneização

Enunciado do Problema
O artigo aborda o desafio de modelar a eletroporação reversível (eletropereabilização) em escala de tecido. Embora os mecanismos ao nível celular de criação e fechamento de poros sejam relativamente bem compreendidos, modelos dinâmicos para campos elétricos em escala de tecido são escassos. Os modelos macroscópicos existentes frequentemente dependem de suposições ad hoc, como a prescrição de uma condutividade sigmoide dependente da intensidade do campo elétrico, sem derivar esse comportamento a partir de primeiros princípios. O problema específico é derivar um modelo macroscópico rigoroso para tecido biológico, representado como um meio periódico de células, partindo de um modelo fenomenológico ao nível celular (Kavian et al. [29]), onde a condutividade da membrana é acoplada a uma equação diferencial ordinária (EDO) para o grau de porosidade. Uma dificuldade fundamental reside no fato de que a condutividade microscópica é constante, embora dados experimentais sugiram que a condutividade efetiva do tecido é não linear, dependente do tempo e exibe efeitos de memória.

Metodologia
Os autores empregam a teoria de homogeneização periódica para unir as escalas micro e macro. A metodologia procede através dos seguintes passos:

1. Análise Dimensional e Escalonamento: Os autores realizam uma análise dimensional para identificar um parâmetro de escala relevante $\varepsilon$, definido como a razão entre o tamanho típico da célula e o tamanho da amostra de tecido. Esta análise revela um escalonamento de $\varepsilon^{-1}$ na condição de transmissão dinâmica através da membrana celular, o que é crucial para capturar o comportamento assintótico correto.
2. Formulação do Problema Microscópico: O modelo microscópico consiste em equações eletrostáticas nos domínios intra e extracelulares, acopladas a uma condição de transmissão dinâmica na membrana celular oscilante. Esta condição envolve um salto no potencial elétrico e um termo de condutividade superficial não linear acoplado a uma EDO para o grau de porosidade $X$.
3. Bem-postura do Problema Microscópico: Os autores provam a existência e unicidade das soluções do problema microscópico. Um obstáculo técnico significativo é que o termo de condutividade superficial não linear não é nem monótono nem Lipschitz contínuo no sentido padrão. Para superar isso, eles utilizam uma não linearidade modificada (truncamento) e provam uma condição de Lipschitz de um lado (monotonicidade generalizada) ao longo das soluções, permitindo a aplicação de argumentos de monotonicidade.
4. Expansões Assintóticas Formais: Expansões de duas escalas formais são derivadas para identificar a estrutura do problema limite. Isso revela um sistema acoplado onde o termo de ordem principal depende de um termo corretor envolvendo efeitos de memória.
5. Homogeneização Rigorosa: As expansões formais são rigorosamente justificadas usando convergência de duas escalas. Um grande desafio é passar para o limite nos termos não lineares definidos nas superfícies oscilantes (membranas celulares). Os autores estabelecem a convergência forte de duas escalas para o salto no potencial e para o grau de porosidade, combinando resultados de compacidade com os argumentos de monotonicidade mencionados acima.
6. Simulação Numérica: Um algoritmo numérico é desenvolvido para resolver a resultante equação integral-diferencial macroscópica. O algoritmo envolve a resolução de problemas celulares auxiliares (estacionários e dependentes do tempo) para computar coeficientes efetivos que dependem do histórico do gradiente de potencial.

Principais Contribuições

• Derivação Rigorosa de Condutividade Não Linear: A principal contribuição é a derivação matemática mostrando que uma condutividade efetiva não linear e dependente do tempo emerge naturalmente de condutividades microscópicas constantes e de um modelo de porosidade dinâmica. Isso fornece uma justificativa teórica para a dependência de forma sigmoide da condutividade do tecido em relação à intensidade do campo elétrico observada em experimentos.
• Efeitos de Memória: O modelo macroscópico derivado é uma equação integral-diferencial não linear que exibe efeitos de memória. A densidade de corrente efetiva depende não apenas do campo elétrico instantâneo, mas também de seu histórico, mediado pela evolução do grau de porosidade.
• Tratamento de Não Linearidades Não Monótonas: O artigo fornece um arcabouço rigoroso para homogeneizar problemas com condições de transmissão dinâmica envolvendo não linearidades que não são globalmente monótonas ou Lipschitz, provando a bem-postura via monotonicidade generalizada ao longo das trajetórias das soluções.
• Validação Quantitativa: Os autores computam numericamente a condutividade efetiva e demonstram que ela captura a evolução não monotônica característica observada em experimentos: uma queda inicial (devido ao carregamento capacitivo da bicamada lipídica) seguida por um aumento acentuado (devido à formação de poros e transporte de íons).

Resultos

• Modelo Macroscópico: O problema limite é formulado como um sistema envolvendo um potencial macroscópico $u_{-1}$ e um corretor $u_0$. A densidade de corrente efetiva é dada por um termo de divergência contendo uma parte instantânea e uma parte dependente do tempo envolvendo uma integral sobre o histórico do gradiente de potencial.
• Dinâmica da Condutividade: Simulações numéricas em um cenário bidimensional (tecido entre placas paralelas) confirmam que a condutividade efetiva $\sigma_{eff}$ evolui ao longo do tempo. Para campos elétricos altos, $\sigma_{eff}$ inicialmente diminui (refletindo a fase de carregamento capacitivo onde os íons são atraídos para a membrana) e depois aumenta bruscamente conforme os poros se formam.
• Dependência de Campo: As simulações mostram que a condutividade efetiva em estado estacionário depende não linearmente da intensidade do campo elétrico aplicado, exibindo uma tendência sigmoide. Isso se alinha com observações experimentais e modelos fenomenológicos anteriores, mas aqui derivado de uma fundação microscópica.
• Convergência: O algoritmo numérico demonstra convergência de segunda ordem em relação ao tamanho da malha e convergência de primeira ordem em relação aos passos de tempo, validando a abordagem computacional.

Significância e Alegações
O artigo alega fornecer uma explicação matemática rigorosa para a dinâmica observada experimentalmente da condutividade do tecido durante a eletroporação. Ao derivar o modelo macroscópico a partir de um modelo fenomenológico ao nível celular, os autores demonstram que o comportamento complexo, não linear e dependente do tempo da condutividade do tecido é uma consequência direta da física microscópica do carregamento da membrana e da evolução dos poros. O trabalho preenche a lacuna entre os mecanismos celulares e as observações em escala de tecido sem recorrer ao ajuste empírico da função de condutividade. Os autores enfatizam que seu modelo captura a evolução não trivial da condutividade efetiva, incluindo a queda característica e o subsequente aumento, o que concorda qualitativamente com dados experimentais. A significância reside na justificativa rigorosa da dependência não linear e de forma sigmoide da condutividade do tecido em relação ao campo elétrico, oferecendo uma base teórica para modelos usados na estimulação elétrica de células biológicas.