🔢 mathematics

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

This paper rigorously derives a nonlinear, time-dependent macroscopic conductivity model for electropermeabilization in periodic biological tissues via homogenization, mathematically explaining the experimentally observed sigmoidal conductivity dynamics and memory effects arising from cell membrane porosity evolution.

Original authors: Tobias Gebäck, Ioanna Motschan-Armen, Irina Pettersson

Published 2026-02-03

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Tobias Gebäck, Ioanna Motschan-Armen, Irina Pettersson

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: From Tiny Cells to Tissue Behavior

Imagine a piece of biological tissue (like a slice of liver) not as a solid block, but as a massive crowd of tiny, individual cells packed together. Each cell is like a small house with a wall (the membrane) that usually keeps things out.

The scientists in this paper wanted to understand what happens when you zap this tissue with a strong electric pulse. This process is called electroporation. It's like using electricity to temporarily poke tiny holes in the cell walls so that drugs or DNA can get inside.

The big mystery they solved is this: Why does the whole tissue act like a different kind of material than the individual cells?

The Micro View: Inside the lab, the fluid inside and outside the cells has a constant ability to conduct electricity (like a steady stream of water).
The Macro View: When you measure the whole tissue, its ability to conduct electricity changes wildly. It doesn't just stay steady; it drops, then shoots up, and it changes depending on how strong the electric "zap" is.

The paper provides a rigorous mathematical "recipe" to explain how the steady, simple behavior of individual cells turns into the complex, changing behavior of the whole tissue.

The Analogy: The Crowd of Sponges

To understand their math, imagine a room filled with thousands of sponges (the cells) floating in water.

The Walls: Each sponge has a rubber skin (the membrane).
The Zapping: You turn on a giant electric fan (the electric field) blowing through the room.
The Reaction:
- Phase 1 (The Drop): At first, the electric wind pushes all the water molecules (ions) toward the sponges. They get stuck trying to get through the rubber skin. Because the water is stuck at the edges, the flow through the room slows down. The "conductivity" (how easily electricity flows) drops.
- Phase 2 (The Popping): If the wind is strong enough, the rubber skin stretches and pops tiny holes (pores). Now, the water can rush through the sponges. The flow speeds up dramatically. The conductivity shoots up.
- Phase 3 (The Shape): If you turn the wind up even stronger, the conductivity doesn't just keep going up forever; it levels off. It forms an "S" shape (a sigmoid curve).

The paper proves mathematically that this "S" shape and the initial drop happen naturally just because of how the cells are arranged and how their skins react, even if the water inside and outside never changes its properties.

The Mathematical Journey

The authors did three main things to prove this:

1. The Micro-Macro Zoom Out (Homogenization)

They started with a very detailed model of a single cell. It was like looking at a single brick in a wall. The math for one brick is simple, but the wall is made of billions of them.

The Challenge: You can't simulate billions of cells on a computer; it would take forever.
The Solution: They used a mathematical technique called homogenization. Think of it like taking a high-resolution photo of a brick wall and zooming out until you can't see the individual bricks anymore. You just see the "average" texture of the wall.
The Twist: Usually, when you zoom out, the math stays simple. But here, because the cell walls are "smart" (they react to the voltage by opening holes), the zoomed-out math becomes nonlinear and time-dependent. The "average" wall behaves differently than the sum of its parts.

2. The Memory Effect

The resulting formula for the tissue isn't a simple snapshot; it has memory.

Analogy: Imagine a door that doesn't just open or close based on who is standing in front of it right now. Instead, the door remembers how hard people have been pushing on it in the last few seconds.
In the paper's model, the tissue's ability to conduct electricity depends on the history of the electric field. The "holes" in the membranes don't open instantly; they take time to form and close. This creates a "lag" or memory in the system.

3. The Proof of Stability

The math involved some tricky, non-smooth functions (like a switch that flips on and off). The authors had to prove that their equations actually have a solution and that the solution is unique (there's only one correct answer). They did this by showing that even though the math is messy, the physical quantities (like voltage) stay within realistic bounds, preventing the model from "exploding" into nonsense.

The Results: What the Computer Said

After doing the heavy math, they ran computer simulations to see what the "zoomed-out" tissue would do.

The Drop: They saw the conductivity dip at the very beginning of the pulse. They explained this as the "charging phase," where ions are gathering at the cell walls like traffic jamming at a toll booth before the gates open.
The Surge: Once the gates (pores) opened, the conductivity spiked.
The S-Curve: When they tested different strengths of electric pulses, the final conductivity followed a perfect "S" curve.
- Low voltage? No change.
- Medium voltage? A sharp increase.
- High voltage? It levels off.

This matched real-world experiments perfectly.

The Bottom Line

The paper's main achievement is a mathematical explanation for a phenomenon scientists have seen in the lab for years.

They proved that you don't need to assume the tissue is weird or that the fluids change properties. You just need to look at how a crowd of tiny, reactive cells interacts with electricity. The complex, changing behavior of the tissue is a natural consequence of the simple rules governing the individual cells.

In short: They built a bridge between the microscopic world of cell membranes and the macroscopic world of tissue conductivity, showing that the tissue's "personality" (its changing conductivity) is just the collective behavior of its tiny, hole-making cells.

Technical Summary: Derivation of Nonlinear Time-Dependent Macroscopic Conductivity for an Electropermeabilization Model via Homogenization

Problem Statement
The paper addresses the challenge of modeling reversible electroporation (electropermeabilization) at the tissue scale. While cellular-level mechanisms of pore creation and resealing are relatively well understood, dynamic models for tissue-scale electric fields remain scarce. Existing macroscopic models often rely on ad-hoc assumptions, such as prescribing a sigmoid-like conductivity dependent on the electric field strength, without deriving this behavior from first principles. The specific problem is to derive a rigorous macroscopic model for biological tissue, represented as a periodic medium of cells, starting from a phenomenological cell-level model (Kavian et al. [29]) where membrane conductivity is coupled with an ordinary differential equation (ODE) for the degree of porosity. A key difficulty lies in the fact that the microscopic conductivity is constant, yet experimental data suggests the effective tissue conductivity is nonlinear, time-dependent, and exhibits memory effects.

Methodology
The authors employ periodic homogenization theory to bridge the micro- and macro-scales. The methodology proceeds through the following steps:

Dimensional Analysis and Scaling: The authors perform a dimensional analysis to identify a relevant scaling parameter $\varepsilon$ , defined as the ratio between the typical cell size and the tissue sample size. This analysis reveals an $\varepsilon^{-1}$ scaling in the dynamic transmission condition across the cell membrane, which is crucial for capturing the correct asymptotic behavior.
Microscopic Problem Formulation: The microscopic model consists of electrostatic equations in the intra- and extracellular domains, coupled with a dynamic transmission condition on the oscillating cell membranes. This condition involves a jump in electric potential and a nonlinear surface conductivity term coupled with an ODE for the porosity degree $X$ .
Well-posedness of the Microscopic Problem: The authors prove the existence and uniqueness of solutions to the microscopic problem. A significant technical hurdle is that the nonlinear surface conductivity term is neither monotone nor Lipschitz continuous in the standard sense. To overcome this, they utilize a modified nonlinearity (truncation) and prove a one-sided Lipschitz condition (generalized monotonicity) along the solutions, allowing for the application of monotonicity arguments.
Formal Asymptotic Expansions: Formal two-scale expansions are derived to identify the structure of the limit problem. This reveals a coupled system where the leading order term depends on a corrector term involving memory effects.
Rigorous Homogenization: The formal expansions are rigorously justified using two-scale convergence. A major challenge is passing to the limit in the nonlinear terms defined on the oscillating surfaces (cell membranes). The authors establish strong two-scale convergence for the jump in potential and the porosity degree by combining compactness results with the monotonicity arguments mentioned above.
Numerical Simulation: A numerical algorithm is developed to solve the resulting macroscopic integro-differential equation. The algorithm involves solving auxiliary cell problems (stationary and time-dependent) to compute effective coefficients that depend on the history of the potential gradient.

Key Contributions

Rigorous Derivation of Nonlinear Conductivity: The primary contribution is the mathematical derivation showing that a nonlinear, time-dependent effective conductivity emerges naturally from constant microscopic conductivities and a dynamic porosity model. This provides a theoretical justification for the sigmoid-shaped dependence of tissue conductivity on electric field strength observed in experiments.
Memory Effects: The derived macroscopic model is a nonlinear integro-differential equation exhibiting memory effects. The effective current density depends not only on the instantaneous electric field but also on its history, mediated by the evolution of the porosity degree.
Handling Non-Monotone Nonlinearities: The paper provides a rigorous framework for homogenizing problems with dynamic transmission conditions involving nonlinearities that are not globally monotone or Lipschitz, by proving well-posedness via generalized monotonicity along the solution trajectories.
Quantitative Validation: The authors numerically compute the effective conductivity and demonstrate that it captures the characteristic non-monotonic evolution observed in experiments: an initial drop (due to capacitive charging of the lipid bilayer) followed by a sharp increase (due to pore formation and ion transport).

Results

Macroscopic Model: The limit problem is formulated as a system involving a macroscopic potential $u_{-1}$ and a corrector $u_0$ . The effective current density is given by a divergence term containing an instantaneous part and a time-dependent part involving an integral over the history of the potential gradient.
Conductivity Dynamics: Numerical simulations in a two-dimensional setting (tissue between parallel plates) confirm that the effective conductivity $\sigma_{eff}$ evolves over time. For high electric fields, $\sigma_{eff}$ initially decreases (reflecting the capacitive charging phase where ions are pulled toward the membrane) and then increases sharply as pores form.
Field Dependence: The simulations show that the steady-state effective conductivity depends nonlinearly on the applied electric field strength, exhibiting a sigmoid trend. This aligns with experimental observations and previous phenomenological models, but is here derived from a microscopic foundation.
Convergence: The numerical algorithm demonstrates second-order convergence with respect to mesh size and first-order convergence with respect to time steps, validating the computational approach.

Significance and Claims
The paper claims to provide a rigorous mathematical explanation for the experimentally observed dynamics of tissue conductivity during electroporation. By deriving the macroscopic model from a cell-level phenomenological model, the authors demonstrate that the complex, nonlinear, and time-dependent behavior of tissue conductivity is a direct consequence of the microscopic physics of membrane charging and pore evolution. The work bridges the gap between cellular mechanisms and tissue-scale observations without resorting to empirical fitting of the conductivity function. The authors emphasize that their model captures the nontrivial evolution of effective conductivity, including the characteristic drop and subsequent rise, which agrees qualitatively with experimental data. The significance lies in the rigorous justification of the nonlinear, sigmoid-shaped dependence of tissue conductivity on electric field strength, offering a theoretical basis for models used in electrical stimulation of biological cells.