🔢 mathematics

Energy, Passivity, and Diffusive Representations in Anisotropic Fractional Maxwell--Cole--Cole Media

This paper establishes a rigorous energy formulation and diffusive representation for linear anisotropic fractional Maxwell–Cole–Cole media, demonstrating that fractional polarization power decomposes into nonnegative storage and dissipation components to prove electromagnetic passivity, global energy decay, and uniqueness while characterizing the broadened dielectric-loss bands and integrable initial layers inherent to fractional-order relaxation.

Original authors: Nikko John Leo S. Lobos

Published 2026-10-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Nikko John Leo S. Lobos

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

Materials around us are rarely perfect insulators. When an electric field pushes against a substance like glass, water, or biological tissue, the atoms inside do not simply snap into place instantly. Instead, they take a moment to rearrange, a process known as dielectric relaxation. This delay means the material's response depends not just on the electric force applied right now, but also on how that force has changed in the recent past. For decades, scientists have used a standard model to describe this behavior, but it often fails to capture the messy, complex reality found in nature, where relaxation happens over a wide range of speeds rather than a single, uniform pace. To bridge this gap, researchers have turned to fractional calculus, a mathematical tool that allows for "memory" in physical laws, describing how a system remembers its history with a fading, non-uniform intensity. Understanding exactly how energy flows, gets stored, and is lost in these complex, memory-laden materials is crucial for designing better sensors, communication devices, and medical imaging tools, yet the precise accounting of energy in such systems has remained a difficult puzzle.

A new study by Nikko John Leo S. Lobos at the University of Santo tackles this puzzle by developing a rigorous energy framework for a specific class of materials that exhibit this fractional memory. The research focuses on linear, anisotropic media, meaning materials where the electrical properties differ depending on the direction of the electric field, and which contain multiple channels of relaxation. The author constructs a detailed energy balance sheet that tracks exactly how much energy is stored in the material's memory and how much is dissipated as heat. The key achievement is proving that even though the mathematical description of these materials involves a sharp, singular spike in power right at the moment an electric field is turned on, the total energy remains finite and well-behaved. The study demonstrates that the energy supplied to the material is always positive, confirming that these complex media are passive and do not spontaneously generate energy, a fundamental requirement for any physical system.

The researchers began by defining the rules for how these materials behave when an electric field is applied starting from a state of rest. They found that when the field is switched on, the internal polarization of the material does not rise smoothly. Instead, it forms a thin, intense layer of activity right at the beginning, where the rate of change is extremely high, theoretically infinite for an instant, but mathematically manageable. This initial burst is followed by a slower, long-lasting tail of activity. The study separates the total power into two distinct parts: a storage component that represents the energy held in the material's history, and a dissipation component that represents the energy lost to friction-like processes within the material. Both of these components are shown to be non-negative, meaning energy is either stored or lost, never created. This decomposition holds true even when the material has multiple different relaxation channels, each with its own speed and strength, and even when the material is not uniform throughout.

To ensure this energy accounting was robust, the author compared two different ways of describing the material's memory. The first method, called the history representation, calculates energy based on the difference between the material's current state and its past states. The second method, known as the diffusive representation, imagines the material as a collection of countless simple, exponential decay processes happening simultaneously. The study proves that while both methods correctly account for the total energy, the history method always predicts a slightly higher amount of stored energy than the diffusive method. The difference between them is not a new form of energy, but rather a measure of how much information is lost when simplifying the complex history into a set of simpler decay modes. This finding clarifies that there is no single, unique way to define the "free energy" of such a material; rather, multiple valid descriptions exist, and they are consistent with one another.

The paper also explores how changing the fractional order, a parameter that controls the breadth of the material's memory, affects its behavior. When this parameter is lowered, the material's response becomes less like a sharp, quick reaction and more like a broad, slow spread. In terms of frequency, this means the material absorbs energy over a much wider range of frequencies rather than at a single peak. For example, in a single-channel model, reducing the fractional order from the standard limit down to a lower value broadens the bandwidth of energy absorption from roughly one decade to nearly four decades. In the time domain, this change manifests as a stronger initial burst of power followed by a much slower, algebraic decay, meaning the material takes significantly longer to settle back to equilibrium after being excited. This behavior explains why certain materials, like biological tissues, can absorb electromagnetic waves across a wide spectrum, a property that is critical for applications ranging from radar to medical diagnostics.

Finally, the study verifies these theoretical findings against exact solutions for specific scenarios, such as a sudden step in voltage or a continuous wave. These checks confirm that the energy balance holds true, that the initial singularities are integrable, and that the long-term behavior matches the expected stationary response. The work establishes a solid foundation for understanding how energy moves through complex, memory-dependent materials without relying on approximations that might fail at the critical moment of initialization. By providing a clear, mathematically sound picture of energy storage and loss, this research offers a reliable tool for engineers and physicists working with advanced dielectric materials, ensuring that their models respect the fundamental laws of thermodynamics and passivity, even in the most complex, fractional regimes.

Technical Summary: Energy, Passivity, and Diffusive Representations in Anisotropic Fractional Maxwell–Cole–Cole Media

Problem Statement
This study addresses the formulation of an energy balance for linear, anisotropic Maxwell–Cole–Cole media characterized by multiple fractional polarization channels, passive conductivity, and prescribed external currents. A central challenge in hereditary dielectric models is defining admissible storage and dissipation functionals when the material possesses several fractional relaxation channels, memory begins at a finite time ( $t=0$ ), and the system retains Maxwell boundary supplies. Previous literature has established constitutive models, time-domain solvers, and regularity theories for fractional electrodynamics; however, a separate finite-time question remains regarding how to decompose polarization power into storage and dissipation for initialized anisotropic systems without assuming infinite prehistory. The study focuses on the continuum energy structure rather than proposing new fractional Maxwell equations or time-stepping algorithms.

Methodology
The analysis employs a finite-terminal Caputo derivative to model the temporal polarization law for each anisotropic channel, while retaining standard macroscopic Maxwell equations for spatial evolution. The model assumes zero prehistory (quiescent initial polarization) and allows for weak initial singularities in the polarization rate, specifically $O(t^{\alpha-1})$ , which are characteristic of Caputo evolution equations.

The methodology proceeds through the following steps:

Mathematical Formulation: The problem is defined within a weighted temporal function class ( $X_\alpha$ ) that admits the integrable fractional initial layer generated by quiescent data. The constitutive law separates the vacuum magnetic response, high-frequency electric displacement, relaxing polarization, Ohmic conduction, and external currents.
History-Based Decomposition: Using the Caputo convolution, the fractional polarization power is decomposed into a non-negative history-storage functional ( $U_{\alpha}$ ) and a non-negative constitutive dissipation rate ( $Q_{\alpha}$ ). This derivation utilizes squared polarization increments and does not rely on a generalized fractional product rule.
Global Energy Theorem: The local power decomposition is integrated with the Poynting theorem to derive a global field-plus-memory energy identity. This identity accounts for instantaneous field energy, hereditary storage, conductive loss, and fractional dissipation.
Diffusive Representation: A second admissible energy decomposition is constructed using a positive diffusive representation of the Caputo kernel. This approach resolves the power-law kernel into a continuum of exponential relaxation modes, generating auxiliary memory states.
Comparative Analysis: The history and diffusive storage functionals are compared analytically. The difference between them is identified as an exact, non-negative weighted variance of polarization increments.
Limiting Cases and Parameter Studies: The analysis examines convergence to Debye, instantaneous-dielectric, vanishing-susceptibility, and vacuum limits. Additionally, a reproducible single-channel parameter study quantifies the effects of the fractional order $\alpha$ on dielectric loss, wave attenuation, and transient step responses.

Key Contributions and Results

Energy Decomposition: The paper establishes that for linear anisotropic Maxwell–Cole–Cole media with zero prehistory, the fractional polarization power can be strictly separated into a non-negative history-storage functional and a non-negative dissipation rate. The instantaneous dissipation may diverge at initialization (scaling as $t^{\alpha-1}$ ) but remains integrable, while the stored history energy vanishes initially.
Passivity and Uniqueness: Combining the energy decomposition with the Poynting theorem yields a global passivity inequality. For closed systems with zero boundary supply, the total field-plus-memory energy is non-increasing. This estimate provides a proof of conditional uniqueness for weak solutions under homogeneous passive boundary conditions (e.g., PEC or periodic).
Diffusive vs. History Storage: A second energy representation based on diffusive variables is derived. The study proves that the diffusive storage never exceeds the history storage. Their difference is an exact non-negative weighted variance that vanishes if and only if the polarization history is zero. This clarifies that while multiple storage functionals are compatible with the same constitutive law, they are not thermodynamically equivalent in terms of magnitude, though neither is claimed to be a unique thermodynamic free energy.
Limiting Behavior: Under stated convergence assumptions, the model recovers the standard Debye form as $\alpha \to 1^-$ , the instantaneous dielectric limit as $\tau \to 0^+$ , and the vacuum Maxwell energy when all polarization channels are removed.
Fractional Order Effects: The parameter study demonstrates that the fractional order $\alpha$ $α$ controls the distribution of loss across time and frequency rather than acting as a uniform loss factor.
- Frequency Domain: Decreasing $\alpha$ broadens the dielectric-loss band (increasing the half-maximum width from ~1.14 decades at $\alpha=1$ to ~3.68 decades at $\alpha=0.4$ ) and lowers the peak loss height.
- Time Domain: A smaller $\alpha$ produces a stronger integrable initial power layer and a slower algebraic relaxation tail (replacing Debye exponential settling).
- Attenuation: The exact attenuation curves cross; a smaller $\alpha$ increases loss at frequencies far from the characteristic frequency but decreases it near the peak, indicating that attenuation is not globally monotone with respect to $\alpha$ .

Significance and Scope
The paper claims to provide a rigorous energy formulation that connects fractional initial layers, hereditary storage, electromagnetic passivity, diffusive relaxation, dispersion, and absorption within a spatially local constitutive model. It establishes that the finite-terminal Caputo derivative, when applied to initialized systems with zero prehistory, yields a consistent energy balance where the "memory" is strictly hereditary and the dissipation is non-negative.

The study explicitly limits its scope to linear, nonmagnetic, spatially local media with time-independent passive tensors. It excludes nonlinear polarization, magnetic dispersion, spatial nonlocality, aging coefficients, and thermal transport. The author notes that the derived non-negative loss is an electromagnetic energy-dissipation rate, not a thermodynamic entropy-production rate; a full thermodynamic interpretation would require a coupled thermal balance and Clausius–Duhem analysis. Furthermore, while the paper validates the exact constitutive responses and energy identities, it does not claim to validate the convergence of fully discrete Maxwell solvers or time-stepping algorithms, noting that such numerical convergence requires separate analysis of mesh, quadrature, and spectral nodes.

The results serve to clarify the physical meaning of stored and dissipated energy in initialized fractional media, offering exact consistency checks (via step responses) for future history and spectral discretization schemes.

More like this