📄 pharmacology and toxicology

Tunneling Effect with Time-Dependent Effective Potential Barrier: A Semiclassical (WKB) Reinterpretation of Drug Release Kinetics in Polymeric Nanocapsules

This paper proposes a semiclassical reinterpretation of drug release kinetics in polymeric nanocapsules by modeling the potential barrier as a time-decaying function rather than a static one, demonstrating through global optimization of ex-vivo data that an exact transmission formula combined with an exponential decay hypothesis provides a superior and more parsimonious fit than previous multifractal or rational decay models.

Original authors: de Albuquerque, D. F., de Albuquerque, M. A. S.

Published 2026-08-05

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: de Albuquerque, D. F., de Albuquerque, M. A. S.

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

Imagine you are trying to sneak a secret note out of a locked room. In the real world, if the door is locked and you don't have a key, you're stuck. But in the weird, tiny world of quantum physics, particles like electrons have a superpower called "tunneling." It's as if the particle can magically phase through the solid wall of the door and pop out the other side without ever actually opening it. Scientists have been using this idea to understand how medicine escapes from tiny capsules inside our bodies. They imagine the capsule wall as a "potential barrier"—a digital wall that the drug molecules try to tunnel through to get out.

For a long time, researchers tried to explain how fast this happens by treating the wall as a static, unchanging object, and then using some very fancy, complex math involving "fractals" (patterns that repeat themselves infinitely) to describe how time passes. But this paper asks a simple question: Is that really how it works? The authors argue that instead of using magic math to turn space into time, we should just look at the wall itself. In the real world, these drug capsules are made of polymers that slowly break down, swell up, or dissolve in water. This means the "wall" isn't solid forever; it gets weaker and thinner as time goes on. The paper proposes a new, simpler way to model this: imagine the barrier is a door that is slowly rusting away. As it rusts, it becomes easier for the drug to slip through. By tracking how fast the door rusts, we can predict exactly when the medicine will be released, without needing to invent complicated new rules for how time works.

The Story of the Rusting Door

In this study, the authors, Douglas and Maria Antonia de Albuquerque, decided to take a fresh look at how drugs escape from polymeric nanocapsules. They started by pointing out a problem with the current "fractal" models. Those models treat the barrier as a fixed wall and then use a mathematical trick (called a Wick rotation) to pretend that space is actually time. The authors argue this is like trying to explain why a cake is baking by saying the oven is actually a clock. It might look like it works on paper, but it doesn't make physical sense.

Instead, they propose a much more grounded idea: the barrier height changes over time. Think of the drug capsule as a fortress. The drug molecules are prisoners trying to tunnel out. In the old models, the fortress walls were made of indestructible steel. In this new model, the walls are made of ice. As time passes, the ice melts. The "barrier" gets lower and lower until, eventually, the drug doesn't need to tunnel at all; it can just walk right over the top.

The team tested two different ideas for how this "melting" happens:

Exponential Decay: The wall melts quickly at first and then slows down, like a snowball rolling down a hill.
Rational (Hill-type) Decay: The wall stays solid for a moment, then starts melting slowly, like a slow-cooking stew.

They used a very precise mathematical formula (the "exact transmission formula") to calculate how much medicine gets out at every single moment. This is a big upgrade from older methods that used an approximation (the WKB method) which worked fine for thick walls but created weird, jagged "steps" in the data when the wall got thin. The new formula is smooth, like a real curve, and handles the moment the wall collapses perfectly.

The Great Fit-Off

To see if their idea actually worked, the authors went hunting for real data. They couldn't use the data from the original fractal studies because those only reported a single number at 24 hours (too little info to test a time-based model). Instead, they found a different study that measured how a drug called 5-Fluorouracil (5-FU) moved through chicken skin over time. They carefully digitized the graphs from that study, turning the lines into thousands of data points.

They then tried to fit their "melting wall" model to this real data. The results were impressive. For two of the three systems they tested (a free drug in a gel and a drug in nanocapsules), their model fit the data almost perfectly. The math showed that the "exponential" melting model was the clear winner. It described the release of the drug with high precision, and the numbers they found for the "melting speed" were very reliable.

However, the third system (nanocapsules inside a gel) was a bit trickier. The data wasn't quite enough to pin down all the numbers perfectly. The authors had to make a choice: should they assume this system has the same "energy ratio" (a measure of how hard it is for the drug to get out) as the free drug? When they forced the numbers to match, the fit was still good. But when they looked at the nanocapsules alone, they found a different "energy ratio" that was significantly different. This suggests that the environment inside the gel changes how the drug behaves, and there isn't just one universal rule for every single capsule type.

What This Means for the Future

The paper concludes that we don't need to use complex fractal geometry or weird space-time tricks to understand drug release. We just need to acknowledge that the barrier changes over time. By treating the capsule wall as something that degrades or swells, the authors created a model that is not only mathematically cleaner but also physically realistic.

They found that the "exponential" model is the best fit for the data they had, offering a clear, smooth prediction of how the drug releases. While they couldn't prove that every single drug system behaves exactly the same way (since the data for some systems was a bit fuzzy), their method successfully ruled out the need for the complicated fractal approach. They showed that a simple, time-dependent barrier is enough to explain the complex dance of drug release, providing a tool that is both accurate and easier to understand. The authors suggest that future experiments should gather even more detailed data to confirm these findings, but for now, the "melting wall" theory stands as a strong, realistic alternative to the old "magic math" models.

Technical Summary: Tunneling Effect with Time-Dependent Effective Potential Barrier

Problem Statement
Recent models describing drug release from polymeric nanoparticles have utilized an analogy with the quantum tunneling effect, treating the delivery system as a static rectangular potential barrier. These models often employ a multifractal formalism where the time evolution of drug release is obtained via a formal Wick rotation ( $x \to t$ ), effectively substituting a spatial coordinate with a temporal one without direct physical justification. The authors argue that this approach lacks rigorous physical grounding, as the time dependence should arise from the physical degradation or swelling of the polymeric matrix rather than an artificial coordinate transformation. Furthermore, previous applications of the WKB (Wentzel-Kramers-Brillouin) approximation in this context introduce artificial step-like transitions at the point of barrier collapse, failing to provide a smooth description of the release kinetics as the barrier height drops below the particle energy.

Methodology
The paper proposes a semiclassical reinterpretation of drug release kinetics that preserves the standard structure of the Schrödinger equation for a rectangular barrier but introduces time dependence through a physically motivated effective barrier height, $U_{\text{eff}}(t) = U_0 f(t)$ .

Time-Dependent Barrier: The barrier height is modeled as a function of time reflecting the progressive degradation or swelling of the polymeric matrix. Two hypotheses for the decay function $f(t)$ are tested:
- Exponential Decay: $f(t) = f_\infty + (1 - f_\infty)e^{-t/\tau}$ .
- Rational (Hill-type) Decay: $f(t) = f_\infty + \frac{1 - f_\infty}{(1 + t/\tau)^n}$ .
  These hypotheses predict distinct dependencies for the finite barrier collapse time $t^*$ on the energy ratio $m = U_0/E$ .
Exact Transmission Formula: Instead of relying on the thick-barrier WKB approximation, which breaks down near the collapse point ( $U_{\text{eff}} \approx E$ ), the authors utilize the exact transmission formula derived from the Schrödinger equation. This formula is real-analytic in $f(t)$ and smoothly transitions into the resonance (over-barrier) regime via analytic continuation ( $\sinh(iz) = i \sin z$ ), avoiding artificial discontinuities.
Traversal Time Definition: The model adopts the kinetic definition of traversal time ( $\tau = \int |\psi|^2 dx / J$ ) to quantify the time scale of the release process and validate the adiabatic approximation, ensuring the barrier changes slowly relative to the particle's transit time.
Numerical Validation: The model was fitted against digitized ex-vivo permeation kinetics of 5-Fluorouracil (5-FU) across chicken skin, sourced from Rata et al. [1]. Three systems were analyzed: NCA-1-5-FU (nanocapsules), G-NCA-1-5-FU (nanocapsules in gel), and G-5-FU (free drug in gel).
- Optimization: Global optimization (differential evolution) followed by local refinement (scipy.optimize.curve_fit) was employed to avoid local minima, a necessity identified due to the oscillatory nature of the exact formula past the collapse point.
- Comparison: The results were compared against the original multifractal model parameters reported by Băcăită et al. [2] and tested for parameter identifiability and sensitivity to the energy ratio $m$ .

Key Results

Fit Quality: The exact transmission formula substantially improved the fit quality ( $R^2 > 0.999$ ) for all three systems compared to the WKB approximation, producing genuinely smooth theoretical curves that correctly model the barrier collapse region.
Model Selection: The exponential decay model was found to be well-identified for the G-5-FU and NCA-1-5-FU systems, with all parameter uncertainties below 11% and 6%, respectively. In contrast, the rational (Hill) model remained poorly identified for all systems, with parameter uncertainties often exceeding the estimates themselves. By the principle of parsimony, the exponential model is preferred.
Identifiability Issues: The G-NCA-1-5-FU system remained non-identified when the energy ratio $m$ was free, a finding attributed to the limited number of digitized data points (13) without replicate error bars.
Energy Ratio ( $m$ ) Sensitivity: A sensitivity check fixing $m$ to the value derived from the well-identified G-5-FU system ( $m=1.377$ ) resolved the non-identifiability for G-NCA-1-5-FU with negligible loss in fit quality. However, applying this shared value to NCA-1-5-FU (which has its own well-identified $m=1.188 \pm 0.007$ ) degraded the fit significantly (maximum residual increased by a factor of ~4). The difference between the two independently fitted values is formally significant ( $\sim 4\sigma$ ), leading to the rejection of a single universal energy ratio $m$ across all studied systems.
Comparison with Multifractal Model: While a direct numerical comparison with the multifractal model [2] was impossible due to the lack of time-resolved data in that study (only 24-hour endpoints were available), the proposed model successfully fitted independent kinetic data with physically interpretable parameters, addressing the methodological gaps of the fractal approach.

Significance and Claims
The paper claims that the potential barrier framework for drug release does not require a multifractal formalism or artificial space-time transformations. The authors assert that:

Physical Motivation: Introducing time dependence via matrix degradation ( $U_{\text{eff}}(t)$ ) provides a direct physical mechanism absent in the Wick rotation approach.
Mathematical Consistency: Using the exact transmission formula eliminates the artificial step-like transitions of the WKB approximation, providing a smooth, real-analytic description of the release kinetics through the resonance regime.
Predictive Power: The distinct functional forms of the exponential and rational decay hypotheses offer a testable criterion (the dependence of barrier collapse time $t^*$ on the energy ratio $m$ ) to distinguish between degradation mechanisms.
Interpretability: The model parameters ( $\tau$ , $f_\infty$ , $T_0$ , $X_n$ ) have direct physical interpretations related to matrix degradation timescales and barrier geometry, unlike the phenomenological constants of fractal models.
Limitations: The authors acknowledge that the adiabatic hypothesis relies on order-of-magnitude estimates of particle energy and that the current validation is limited by the sparsity of the digitized data, particularly regarding the non-identifiability of parameters in the G-NCA-1-5-FU system. They emphasize that the rejection of a universal $m$ is a specific finding for the tested systems and that further experimental validation (e.g., replicate measurements, additional time points) is required to generalize the findings.

Technical Summary: Tunneling Effect with Time-Dependent Effective Potential Barrier

More like this