🔢 mathematics

Fractal-Fractional HIV Dynamics with Mittag-Leffler Kernel: Analysis, Stability, and Numerical Simulations

This paper proposes and analyzes a fractal-fractional HIV model using the Atangana-Baleanu-Caputo operator with a Mittag-Leffler kernel to capture memory effects, establishing solution existence, Hyers-Ulam stability, and numerical reliability through advanced approximation methods and sensitivity analyses.

Original authors: Niaz Ali Shah, Samad Noeiaghdam

Published 2026-07-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Niaz Ali Shah, Samad Noeiaghdam

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

Imagine the human body as a bustling city under siege by an invisible enemy: the HIV virus. For decades, scientists have tried to map out how this battle plays out using math. Usually, they use "standard" math, which assumes that what happens today depends only on what is happening right now. It's like driving a car where you only look at the road directly in front of you, ignoring the turns you just made.

But the authors of this paper, Niaz Ali Shah and Samad Noeiaghdam, argue that the body is more like a smart, memory-keeping driver. The immune system doesn't just react to the present; it remembers past infections, past treatments, and how it has responded over time.

Here is a simple breakdown of what they did, using everyday analogies:

1. The New "Memory" Map (The Model)

The researchers built a new, more sophisticated map of the HIV battle. Instead of a flat, standard map, they used a Fractal-Fractional Model.

The Analogy: Think of a standard map as a flat piece of paper. A "fractal" map is like a crumpled piece of paper or a coastline—it has twists, turns, and details at every level of zoom. This allows the model to capture the messy, complex reality of how the virus spreads and how the immune system fights back.
The "Memory" Ingredient: They used a special mathematical tool called the Mittag-Leffler kernel (part of the Atangana-Baleanu operator). Imagine this as a "memory bank." In their model, the virus and the immune cells don't just react to the current moment; they carry the "weight" of their history. If the immune system was strong yesterday, that strength lingers and influences today, just like how a muscle remembers a workout from last week.

2. The Five Key Players

The model tracks five specific groups in this biological city:

Healthy Soldiers (T): The healthy CD4+ T-cells (the city's main defenders).
Infected Soldiers (I): The CD4+ cells that the virus has taken over.
The Enemy Army (V): The free-floating HIV viruses.
Special Forces (Z): The CD8+ T-cells, which are the elite troops sent to hunt down infected cells.
Activated Special Forces (Za): The CD8+ cells that have been "woken up" and are actively fighting.

3. Proving the Map Works (Stability and Existence)

Before showing the results, the authors had to prove their map wasn't broken.

Existence: They proved that if you start with a specific number of soldiers and viruses, the math guarantees a solution exists. It's like proving that if you drop a ball, physics guarantees it will fall somewhere, rather than vanishing into thin air.
Uniqueness: They showed there is only one correct path for the battle to take given the starting conditions. You won't get two different outcomes from the same starting point.
Stability (The "Wobble" Test): They performed a "Hyers-Ulam stability" test. Imagine you are balancing a stack of plates. If you nudge the stack slightly (a small error in data or a tiny change in the virus), does the whole tower collapse? Their math says no. The system is stable; small mistakes won't ruin the entire prediction.

4. Running the Simulation (The Computer Experiment)

Since these equations are too complex to solve with a pencil and paper, they used a special computer method called the Newton polynomial approximation combined with the Atangana-Toufik scheme.

The Analogy: This is like using a high-speed drone to simulate a storm. Instead of guessing the wind, the computer calculates the path step-by-step, looking back at previous steps (the memory) to predict the next move.

5. What the Results Showed

When they ran the simulations, they tweaked the "knobs" (parameters) to see what happened.

The "Memory" Effect: When they turned up the "memory" (using the fractal dimension), the system behaved differently. The immune response became smoother and more persistent. It was as if the body remembered how to fight better over time, leading to a slower, more controlled decline in health rather than a sudden crash.
The Virus: The virus population dropped faster when the "memory" was stronger, suggesting that accounting for history helps predict how quickly the virus can be suppressed.
The Soldiers: The healthy soldiers (T-cells) recovered better in the memory-inclusive model compared to the standard model.

6. The Sensitivity Check (The Heatmap and Tornado)

To understand which factors mattered most, they used two visual tools:

The Heatmap: Imagine a weather map where red means "very sensitive" and blue means "not sensitive." They colored their model to show that the viral load (how much virus is in the body) is extremely sensitive to how fast the virus infects cells and how fast it is produced.
The Tornado Diagram: This is a bar chart that looks like a tornado, ranking the factors from most dangerous to least dangerous. It showed that if you want to control the disease, you need to focus on the parameters that cause the biggest "winds" in the chart—specifically, the rates of infection and viral production.

The Bottom Line

The paper claims that by adding memory and complexity (fractals) to the math, we get a much more realistic picture of how HIV behaves. The standard models are like watching a movie in black and white; this new model is like watching it in 4K with surround sound. It captures the lingering effects of the past, showing that the immune system's history is crucial for understanding the future of the infection.

The authors conclude that this new mathematical framework is a reliable tool for understanding the short-term and long-term behavior of HIV, and they suggest it could be used for other diseases like Malaria or COVID-19 in the future.

Technical Summary: Fractal-Fractional HIV Dynamics with Mittag-Leffler Kernel

Problem Statement
The paper addresses the limitations of classical integer-order mathematical models in capturing the complex, memory-dependent, and hereditary properties of HIV transmission dynamics. While traditional models describe biological processes using standard differential equations, they often fail to account for the historical dependence and fractal characteristics inherent in immune responses and viral progression. The authors propose that a fractal-fractional framework is necessary to more realistically model the transient and long-term behaviors of HIV, specifically focusing on the interactions between healthy CD4+ T-cells, infected CD4+ T-cells, free virions, and the immune response mediated by CD8+ T-cells (both resting and activated).

Methodology
The study employs a rigorous mathematical and numerical approach centered on the Atangana-Baleanu-Caputo (ABC) operator with a Mittag-Leffler kernel.

Model Formulation: The authors extend a classical five-compartment HIV model (Healthy T-cells, Infected T-cells, Virions, CD8+ T-cells, and Activated CD8+ T-cells) into a fractal-fractional differential equation system. This formulation incorporates two parameters: a fractional order ( $\alpha$ ) to capture memory effects and a fractal dimension ( $\gamma$ ) to describe the heterogeneity of the system.
Theoretical Analysis:
- Existence and Uniqueness: Using fixed-point theorems and functional analysis within a Banach space, the authors reformulate the system into an integral form. They prove the existence of solutions by demonstrating that the system's kernels satisfy Lipschitz conditions and are contractions.
- Stability Analysis: The paper conducts a Hyers-Ulam (HU) stability analysis to verify the robustness of the solutions against small perturbations, ensuring that approximate solutions remain close to exact solutions.
Numerical Scheme: To solve the non-linear fractal-fractional system, the authors utilize a numerical technique combining the Newton polynomial approximation method with the Atangana-Toufik numerical scheme. This approach approximates the continuous integrand using historical data points to handle the non-local nature of the fractional operator.
Visualization and Sensitivity: The study employs sensitivity heatmaps and tornado diagrams to analyze the influence of key epidemiological parameters on the system's output, specifically the viral load and cell populations.

Key Contributions

Novel Framework: The proposal of a specific fractal-fractional HIV model using the ABC operator with a Mittag-Leffler kernel, which allows for the simultaneous incorporation of memory effects and fractal characteristics.
Mathematical Rigor: A comprehensive proof of the existence and uniqueness of the model's solutions, alongside a formal Hyers-Ulam stability analysis, which establishes the theoretical reliability of the proposed system.
Advanced Numerical Implementation: The application of the Newton polynomial approximation within the Atangana-Toufik scheme to generate approximate solutions for complex parameter settings, providing a computational tool for analyzing fractal-fractional biological systems.
Parameter Sensitivity: A detailed visualization of how variations in transmission rates, clearance rates, and immune activation parameters affect the dynamics of the HIV system, identifying dominant factors through tornado diagram analysis.

Results
The numerical simulations and theoretical analysis yielded the following findings:

Memory and Fractal Effects: The fractional order ( $\alpha$ ) and fractal dimension ( $\gamma$ ) significantly influence the system's dynamics. Higher values of $\alpha$ generally lead to faster dynamics and quicker clearance of infected cells and viral loads, while reducing $\gamma$ (e.g., from 1 to 0.8) enhances memory effects, resulting in smoother transitions and prolonged immune responses.
Dynamics of Compartments:
- Healthy T-cells ( $T$ ): Show recovery over time, with the rate of recovery dependent on the fractional order.
- Infected Cells ( $I$ ) and Virions ( $V$ ): Exhibit rapid decay, with the fractional model showing faster suppression of viral load compared to integer-order models under certain conditions.
- Immune Response ( $Z, Z_a$ ): The activated CD8+ T-cells demonstrate coupled oscillatory behavior with viral loads, indicating a strong feedback mechanism. The fractal dimension influences the persistence of these immune cells.
Stability: The system was confirmed to be Hyers-Ulam stable, indicating that the model is robust against small perturbations in initial conditions or parameters.
Sensitivity: Sensitivity analysis revealed that viral load and infected cell populations are highly sensitive to infection transmission and production rates, while immune compartments are more responsive to activation and clearance parameters. Phase portraits (2D and 3D) revealed complex, non-closed, and oscillatory trajectories, suggesting rich transient dynamics driven by nonlinear coupling and memory effects.

Significance and Claims
The paper claims that the proposed fractal-fractional framework provides an "effective and reliable approach" for analyzing HIV transmission dynamics. By integrating memory and hereditary properties, the model offers a more realistic description of the disease's progression compared to classical integer-order models. The authors assert that their work highlights the critical role of parameter values in shaping disease control strategies.

The significance of the study is framed within the context of improving mathematical modeling for biological systems. The authors suggest that this framework can be extended to other complex biological systems, such as Malaria and COVID-19 infections, in future investigations. The work emphasizes the utility of combining advanced fractional calculus operators with robust numerical schemes and visualization tools to gain deeper insights into epidemiological dynamics.