🔢 mathematics

Stability Analysis of Pantograph Delay Differential Equations

This paper establishes analytic stability criteria for pantograph delay differential equations to partition the parameter plane into distinct stability regions, validates these findings with numerical simulations, and explores the chaotic dynamics of a proposed proportional-delay analogue of the Mackey-Glass equation.

Original authors: Sachin Bhalekar

Published 2026-05-22

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sachin Bhalekar

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 you are trying to predict how a ball will roll down a hill. Usually, you look at the hill right now and the ball's current speed to guess where it will be next. But in the world of Pantograph Delay Differential Equations, the rules are a bit stranger.

In this paper, the author, Sachin Bhalekar, studies a specific type of mathematical equation where the "hill" the ball is rolling on depends not just on where the ball is now, but on where it was at a proportional time in the past.

Think of it like this: If it is currently 10:00 AM, the equation doesn't look at 9:00 AM (a fixed delay). Instead, it looks at a time that is a fraction of the current time. If the "proportion" is 0.5, at 10:00 AM, it looks at 5:00 AM. At 2:00 PM, it looks at 1:00 PM. The "look-back" window grows as time goes on.

Here is a breakdown of what the paper discovers, using simple analogies:

1. The Problem: Why is this hard?

Standard math tools for predicting stability (whether a system settles down or flies apart) usually rely on a "magic formula" called a characteristic equation.

The Analogy: Imagine trying to solve a puzzle where the pieces keep changing size as you try to fit them together. Because the time delay in these equations changes as time passes, the standard "magic formula" breaks down. You can't just plug in a number and get an answer; the math gets messy because the delay is constantly stretching.

2. The Goal: Will the system explode or settle?

The paper asks a simple question: If we start with a small push, will the system eventually calm down and stop moving (Stability), or will it grow uncontrollably until it breaks (Instability)?

The author divides the "control knobs" (parameters $a$ and $b$ ) into different zones:

The "Explosion" Zone (Instability):
- If the current force ( $a$ ) and the past force ( $b$ ) are both pushing in the same direction (positive), the system is like a snowball rolling down a hill that gets steeper every second. It will grow forever.
- Even if one is negative, if the total push ( $a+b$ ) is positive, the system is unstable. It's like trying to stop a car by pressing the gas pedal just a little bit harder than you press the brake.
The "Calm" Zone (Stability):
- If the current force is pulling back (negative) and the past force isn't too strong, the system acts like a shock absorber. It wobbles a bit but eventually comes to a rest.
- The paper proves mathematically that if the "pull back" is strong enough and the "proportion" of the delay is within a certain range, the system will definitely settle down.

3. The "Ultra-Slow" Mystery

One of the most interesting findings is what happens in the "gray area" between total chaos and total calm.

The Analogy: Imagine a pendulum that is supposed to stop. In a normal system, it stops quickly. In these equations, under certain conditions, the pendulum slows down so incredibly slowly that it looks like it's frozen, but it's actually still moving.
The paper calls this "ultra-slow convergence." The system takes a massive amount of time to settle, often oscillating (swinging back and forth) with a rhythm that gets longer and longer as time goes on. It's like a clock that ticks once a year, then once a decade, then once a century.

4. The Chaos Experiment (The Mackey-Glass Analogue)

The author takes a famous chaotic system (the Mackey-Glass equation, often used to model blood cell production) and swaps the fixed time delay for this "proportional" delay.

The Result: By tweaking the parameters, they found that the system can behave in three ways:
1. Stable: It settles to a steady value.
2. Delay-Dependent: It is stable only if the "proportion" factor is just right. If you change the proportion slightly, it becomes unstable.
3. Chaotic: It starts swinging wildly and unpredictably, never settling down. The paper shows that even with this changing delay, the system can produce "strange attractors"—complex, swirling patterns that look like chaos but follow hidden rules.

5. How They Figured It Out

Since the math was too hard to solve with a single formula, the author used a three-pronged approach:

Series Solution: Writing the answer as an infinite list of numbers to see if the numbers get huge (unstable) or small (stable).
Energy Functions (Lyapunov): Imagine a ball in a bowl. If the ball always rolls to the bottom, the system is stable. The author built a mathematical "energy bowl" to prove the ball would settle.
Computer Simulations: They wrote a computer program (a "predictor-corrector" method) to simulate the equations step-by-step. This allowed them to map out exactly where the "explosion" zones and "calm" zones are on a graph.

Summary

The paper doesn't just say "this is stable" or "this is unstable." It draws a detailed map of the parameter plane. It shows that while we can't use the old, simple formulas, we can still predict the behavior of these "growing delay" systems.

The key takeaway is that proportional delay creates unique behaviors: systems that are stable only under specific conditions, and systems that take an incredibly long time to settle down, often oscillating in a way that gets slower and slower as time marches on. The author also provides a roadmap for where chaos begins in these specific types of equations.

Technical Summary: Stability Analysis of Pantograph Delay Differential Equations

Problem Statement
This paper investigates the stability properties of the scalar pantograph delay differential equation (DDE):
$\dot{x}(t) = ax(t) + bx(qt), \quad 0 < q < 1$
where the delay is proportional to the current time, represented as $x(qt) = x(t - (1-q)t)$. Unlike constant-delay DDEs, pantograph equations do not yield a simple characteristic quasi-polynomial via the standard exponential ansatz ( $x(t) = e^{\lambda t}$ ), rendering classical spectral stability theory inapplicable. Furthermore, the explicit dependence on both $t$ and $qt$ prevents direct reduction to a characteristic equation, and standard numerical methods for ordinary differential equations require significant modification to handle the proportional delay.

Methodology
The study employs a hybrid approach combining analytical techniques with numerical simulations:

Series Representation: The authors utilize the known power-series solution of the equation to derive elementary sufficient conditions for instability.
Lyapunov Methods:
- Lyapunov–Krasovskii Theory: A functional is constructed to derive sufficient conditions for asymptotic stability.
- Lyapunov–Chetaev Theory: A function is defined to establish sufficient conditions for instability.
Numerical Simulation: A predictor-corrector method (detailed in the Appendix) is implemented to validate theoretical findings and explore regions of the parameter plane ( $a, b, q$ ) where analytical results are inconclusive. This method utilizes linear interpolation for the delayed term $x(qt)$ when the argument falls between discrete time steps.
Extension to Nonlinear Systems: The stability framework is applied to a $q$ -analogue of the Mackey–Glass equation to investigate bifurcations and chaotic behavior.

Key Contributions and Results

Analytical Stability Criteria:
- Instability: The system is proven unstable if $a > 0$ and $b > 0$ , or if $a > -b > 0$ . Additionally, if $a + b > 0$ , the system is unstable for all $q \in (0, 1)$ .
- Asymptotic Stability: A rigorous sufficient condition for asymptotic stability is derived: $a < 0$ , $a < b < -a$ , and $b^2/a^2 < q < 1$ . The authors note this condition is sufficient but not necessary.
- Parameter Plane Partitioning: The study maps the parameter plane into three distinct regions: delay-independent unstable, delay-independent stable, and delay-dependent stability regions.
Numerical Observations and Refinements:
- Refined Stability Region: Numerical computations suggest that the sufficient condition derived from the Lyapunov–Krasovskii functional is conservative. The system appears to be asymptotically stable for all $q \in (0, 1)$ within the broader region defined by $a < 0$ and $a < b < -a$ .
- Threshold Curve ( $b^*$ ): In the region where $a < 0$ and $b < a$ , the simulations indicate the existence of a threshold value $b^*(a)$ that separates parameters leading to instability for all $q$ from those exhibiting delay-dependent stability. This threshold is empirically approximated by the quadratic fit:
  $b^*(a) \approx 0.014455a^2 + 1.391065a - 0.03994898$
- Dynamical Behavior: The paper identifies "ultra-slow" convergence or divergence in boundary cases, often accompanied by oscillations with periods that increase over time due to the unbounded growth of the delay $\tau(t) = (1-q)t$ .
Mackey–Glass Analogue:
- The authors formulate a proportional-delay version of the Mackey–Glass equation: $\dot{x}(t) = -\beta x(t) + \frac{\alpha x(qt)}{1 + x(qt)^c}$ .
- By linearizing around non-zero equilibria, the problem reduces to the pantograph equation with specific coefficients.
- Two bifurcation thresholds for the parameter $\alpha$ are identified:
  $\alpha^*_1 = \frac{\beta c}{c - 2}, \quad \alpha^*_2 = \frac{\beta^2 c}{b^* + \beta(c - 1)}$
- These thresholds delineate regimes of delay-independent stability, delay-dependent stability, and instability. For $\alpha > \alpha^*_2$ , the system exhibits chaotic oscillations and coexisting chaotic attractors.

Significance and Limitations
The paper's primary significance lies in addressing the analytical challenges posed by proportional delays, which preclude standard characteristic equation methods. By combining series solutions, Lyapunov theory, and numerical validation, the work provides a comprehensive stability diagram for the linear pantograph equation and extends these insights to a nonlinear chaotic system.

The authors explicitly acknowledge the modesty of their claims regarding the stability boundaries. They state that the conjecture that the entire region $a < 0, a < b < -a$ is asymptotically stable for all $q$ requires a rigorous proof. Similarly, the threshold curve $b^*(a)$ is currently empirical; deriving it analytically or proving sharp bounds remains an open question. The work serves as a foundational step toward a more complete theory for pantograph-type delay systems, highlighting the unique dynamical features (such as ultra-slow divergence and time-varying oscillation periods) inherent to proportional delays.