On mathematical characterization of a Bessel functions-based passive element in electronic circuits

This paper proposes a novel passive electronic element defined by modified Bessel functions that provides a physically interpretable, stable, and analytically tractable alternative to fractional-order models for capturing the broadband dispersive behavior of complex media like biological tissues.

The Gist

Imagine you are trying to model how a sponge absorbs water, or how a piece of muscle tissue reacts to an electrical pulse. In the world of science and engineering, these are called "relaxation phenomena."

When you poke a jelly dessert, it doesn't just snap back instantly like a metal spring; it wobbles, settles, and slowly returns to its original shape. Modeling that "wobble" is incredibly hard because it isn't just one simple movement—it’s a complex, messy mix of many different speeds and rhythms.

This paper introduces a new mathematical "tool" to describe this messiness more accurately and easily. Here is the breakdown of how it works.

1. The Problem: The "Clunky" Old Tools

For a long time, scientists have used two main ways to model these complex materials:

• The Simple Way (The Single Note): Like trying to describe a whole symphony using only one musical note. It’s easy to calculate, but it misses all the beauty and complexity of the actual sound.
• The Complex Way (The Infinite Sheet Music): Scientists use "fractional calculus," which is like trying to describe a symphony by writing down every single vibration of every single string. It’s incredibly accurate, but it’s a nightmare to actually use in a computer simulation or build into a real electronic circuit because it requires "infinite memory"—the math assumes the system remembers everything that ever happened to it since the beginning of time.

2. The Solution: The "Bessel" Instrument

The authors propose a new way to model this using something called Bessel functions.

Think of a Bessel function not as a single note, but as a perfectly tuned chord. A chord contains several notes that work together to create a rich, complex sound, but it is still much easier to play than a full orchestra.

By using these special mathematical functions, the researchers have created a "passive element"—a theoretical electronic component—that can mimic the complex "wobble" of biological tissues (like skin or muscle) without needing the impossible "infinite memory" of the old methods.

3. Why is this a big deal? (The Three Superpowers)

The paper highlights three reasons why this "Bessel tool" is better:

• Superpower 1: It’s "Real-World Ready" (Physical Realizability).
  Imagine trying to build a machine that requires a part that doesn't exist in nature. That’s what some old models do. The Bessel model, however, behaves exactly like real electrical components (resistors and capacitors). You could actually build a circuit that mimics this math.
• Superpower 2: It’s "Fast and Efficient" (Computational Tractability).
  Because the math is "closed-form" (meaning it has a neat, finished equation), a computer doesn't have to struggle with infinite calculations. It’s like the difference between trying to draw a circle by counting a billion tiny dots versus just using a compass. One is much faster and just as good.
• Superpower 3: It’s "Smart" (Spectral Control).
  The researchers found that by turning a single "dial" (a parameter called $\nu$), they can change how the model behaves. They can make it act like a smooth, gradual transition or a sharp, sudden one. This makes it perfect for "tuning" the model to match specific things, like the difference between dry skin and a piece of muscle.

4. The "Real World" Test: The Skin and Muscle Test

To prove it works, they tested their math against real data from biological tissues.

• They modeled dry skin (which is quite resistant).
• They modeled muscle tissue (which is more complex and has multiple layers of "wobble").

The math matched the real-world biological data almost perfectly. It captured the way electricity flows through these tissues across different frequencies, proving that this "mathematical chord" is a much better way to describe the "symphony" of life.

Summary

In short: Scientists found a way to describe the complex, messy way that materials (especially living ones) respond to energy. Instead of using math that is too simple to be true or too complex to be used, they found a "Goldilocks" solution: a mathematical tool that is complex enough to be accurate, but simple enough to be used in real computers and real electronic devices.

Technical Summary

Technical Summary: Mathematical Characterization of a Bessel Functions-Based Passive Element

1. Problem Statement

Modeling relaxation phenomena (the time-dependent return to equilibrium in complex media) is essential in fields like materials science, bioengineering, and condensed matter physics. Current modeling approaches face several critical limitations:

• Fractional-Order Models (e.g., Cole-Cole, Havriliak-Negami): While flexible and capable of capturing broadband dispersion, they often lack physical interpretability, require non-local (infinite-memory) operators, and lack closed-form time-domain expressions, making them difficult to implement in circuit-level simulations or real-time hardware.
• Classical Debye Models: These assume a single relaxation time (purely exponential decay), failing to capture the complex, asymmetric, and broad relaxation spectra observed in heterogeneous systems like biological tissues.

2. Methodology

The authors propose a novel passive electrical element whose impedance is defined analytically using modified Bessel functions of the first kind ($I_\nu$). The methodology follows these steps:

• Electro-Mechanical Analogy: The model is derived from linear viscoelasticity (specifically hemodynamic modeling). By exploiting the correspondence between mechanical stress/strain and electrical voltage/current, the authors translate viscoelastic relaxation functions into an electrical impedance framework.
• Mathematical Definition: The impedance $e_{ZB}(s)$ in the Laplace domain is defined as:
  $$e_{ZB}(s) = R_\infty \frac{I_\nu(\sqrt{s\tau})}{I_{\nu+2}(\sqrt{s\tau})}$$
  where $R_\infty$ is the high-frequency resistance, $\tau$ is the relaxation time, and $\nu > -1$ is a dimensionless parameter controlling spectral dispersion.
• Time-Domain Derivation: Unlike fractional models, the authors derive an exact analytical impulse response $z_B(t)$ using an infinite series of exponentially decaying modes based on the zeros of the Bessel function $J_{\nu+2}$.
• Hybrid Kernel Approximation: To improve computational efficiency, they introduce a "hybrid kernel" that uses a finite sum of exact modes for early-time response and a compact rational-exponential surrogate for the long-time "tail."

3. Key Contributions

• New Passive Element Class: Introduction of a mathematically rigorous element that captures distributed relaxation without the need for fractional calculus.
• Analytical Tractability: Provides closed-form expressions for both frequency-domain impedance and time-domain impulse responses.
• Hybrid Kernel Formulation: A novel method to balance numerical accuracy and computational speed, making the model suitable for real-time time-stepping simulations and finite-element methods.
• Mittag-Leffler Expansion of Admittance: A mathematical proof showing that the admittance can be represented as a Prony series (a sum of simple fractions), which is highly useful in system theory.

4. Results and Validation

• Mathematical Properties: The element is proven to be analytic, causal, monotonic, and BIBO stable. It is strictly passive, meaning it only stores or dissipates energy, satisfying the second law of thermodynamics.
• Spectral Control: The parameter $\nu$ allows for precise "tuning" of the frequency response. As $\nu \to -1$, the model recovers a classical RC circuit behavior; as $\nu$ increases, the transition from capacitive to resistive behavior becomes smoother and more gradual.
• Circuit Modeling: The authors demonstrated that a parallel configuration of a DC resistor ($R_0$) and multiple Bessel elements can model complex, multi-scale relaxation processes.
• Biological Validation: The model was successfully applied to experimental data for dry skin and muscle tissue. The Bessel-based model provided a high-quality fit to the magnitude and phase of the bioimpedance spectra, outperforming the descriptive capabilities of simpler models.

5. Significance

This work bridges the gap between phenomenological flexibility (the ability to fit complex data) and structural interpretability (the ability to link parameters to physical mechanisms).

By providing a model that is both mathematically elegant and computationally efficient, the authors offer a "drop-in replacement" for fractional-order models in electronic circuit simulators and bioimpedance spectroscopy hardware. This has significant implications for the development of more accurate medical diagnostic tools and the simulation of complex soft-matter environments.