Dynamic behavior of a two-dimensional discrete-time HIV infection model: stability and bifurcation analysis

Human immunodeficiency virus, known as HIV, is a relentless pathogen that attacks the body's immune system, specifically targeting the cells responsible for coordinating the defense against illness. Without treatment, this infection can progress to a stage where the body can no longer fight off common diseases. For decades, scientists have used mathematical models to understand how the virus behaves inside the human body. Traditionally, these models have been built on the idea of continuous time, treating the infection as a smooth, unbroken flow where changes happen at every single instant. However, in the real world, medical data is rarely collected in a smooth stream. Doctors and researchers take measurements at specific intervals—days, weeks, or months apart. This gap between the smooth theory of continuous time and the step-by-step nature of real-world data has led a team of researchers at Alexandria University to explore a different approach: a model that moves in discrete steps, mirroring the way we actually observe the disease.

The researchers, Sarah Ramadan, Sanaa Salman, and Ahmed El-Sayed, set out to build a simplified but powerful computer simulation of HIV infection. They started with a known continuous model that describes the interaction between uninfected immune cells and the virus, then translated it into a step-by-step system. In this new framework, the state of the infection is updated in distinct jumps rather than a smooth glide. The team focused on two main populations: the healthy immune cells that the virus targets, and the infected cells that are actively producing new virus particles. By adjusting the rate at which the virus infects healthy cells, they watched to see how the system responded over time. Their goal was not just to see if the virus would die out or persist, but to uncover the hidden patterns that emerge when the infection rate changes, looking for signs of complex, unpredictable behavior that simpler models might miss.

What they found was a world of surprising complexity hidden within the simple rules of infection. When the infection rate is low, the system behaves predictably, settling into a steady state where the virus is either cleared or held at a constant, manageable level. However, as the researchers increased the infection rate in their simulations, the system began to change in dramatic ways. Instead of settling down, the numbers representing the virus and the immune cells began to oscillate. They did not just fluctuate randomly; they entered into rhythmic cycles, rising and falling in a regular pattern. In some scenarios, these cycles doubled in length, creating a pattern where the virus peaks twice as often before the cycle repeats. In other cases, the behavior became even more intricate, with the system tracing out a smooth, closed loop that never quite repeats the exact same path, a state known as quasiperiodic motion.

Perhaps most strikingly, the team discovered that under certain conditions, the system could tip into chaos. In this state, the viral load and immune cell counts fluctuate in a way that appears random and is impossible to predict over the long term, even though the rules governing the system are fixed and deterministic. This means that the erratic swings in viral levels sometimes seen in patients might not be caused by external noise or random chance, but could be an inherent feature of the biological interaction itself. The researchers used a tool called the Lyapunov exponent to measure this behavior, a method that acts like a thermometer for chaos. When the reading was negative, the system was stable or rhythmic; when it turned positive, the system had entered a chaotic regime. Their simulations showed that the transition from order to chaos happened through a specific route, where the system first began to oscillate, then the oscillations doubled, and finally, the pattern broke down into chaos.

The study also identified specific thresholds where these dramatic shifts occur. By varying the infection rate parameter, the team pinpointed the exact moments when the system loses its stability. They found that the model undergoes two distinct types of transitions. The first is a flip bifurcation, where the system jumps from a steady state to a cycle that repeats every two steps, and then to cycles of four, eight, and so on, until it becomes chaotic. The second is a Neimark-Sacker bifurcation, where the system moves from a steady state directly into a smooth, looping oscillation that never repeats exactly. These findings suggest that the complex, sometimes erratic behavior of HIV in the body is not necessarily a sign of a failing treatment or an unpredictable environment, but rather a natural consequence of the nonlinear feedback between the virus and the immune system.

The implications of this work extend beyond the mathematics. The researchers suggest that the "viral blips"—sudden, temporary spikes in viral load that occur even in patients who are otherwise stable on medication—might be a deterministic result of these internal dynamics rather than random errors. If the system is naturally prone to these oscillations and chaotic states, then the timing of viral surges could be an intrinsic property of the infection itself. This perspective shifts the focus from trying to eliminate all fluctuations to understanding the underlying rhythms of the disease. The study does not claim to have solved the problem of HIV or to offer a new cure, but it provides a robust theoretical framework for understanding why the virus behaves the way it does. By showing that discrete-time models can capture rich, complex behaviors like periodicity and chaos, the work offers a new lens through which to view the persistence of the virus and the potential strategies for controlling it. The results, confirmed through extensive computer simulations, highlight that the battle between the virus and the immune system is a dynamic, ever-shifting dance of numbers that can tip into unpredictability with the slightest change in conditions.