🔢 mathematics

Exploring Multi-Rhythmicity in the light of Hilbert's 16th Problem : A Liénard--Levenson-Smith Perspective and the Potency of Perturbative Methods

This paper synthesizes mathematical theory and perturbative methods to classify nonlinear oscillators within the generalized Liénard--Levenson--Smith framework, offering insights into Hilbert's 16th problem and the dynamics of multi-rhythmic biological systems while critically addressing the framework's limitations regarding the uniqueness of limit cycles.

Original authors: Sandip Saha, Gautam Gangopadhyay, Deb Shankar Ray

Published 2026-07-23

📖 8 min read🧠 Deep dive

CC BY 4.0

Original authors: Sandip Saha, Gautam Gangopadhyay, Deb Shankar Ray

Original paper licensed under CC BY 4.0 (https://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 universe as a giant, cosmic dance floor. Some dancers spin in perfect, predictable circles, while others wobble, speed up, and slow down in chaotic patterns. In the world of science, this dance is called "dynamics," and the most fascinating dancers are the ones that settle into a steady, repeating rhythm on their own, without anyone pushing them. Scientists call these self-sustaining loops "limit cycles." You can find them everywhere in real life: the steady thump of a human heart, the daily rise and fall of sleepiness (circadian rhythms), or the rhythmic flashing of fireflies.

But here is the tricky part: what happens when a system tries to dance to two or three different rhythms at once? Can a single system have multiple stable loops, like a dancer who can instantly switch between a slow waltz and a fast salsa, both of which are equally stable? This question touches on a famous, century-old puzzle in mathematics known as Hilbert's 16th Problem. It asks a simple but stubborn question: for a specific type of mathematical equation, what is the maximum number of these rhythmic loops possible? For over 100 years, mathematicians have been trying to count these loops, but the answer has remained elusive, especially when the equations get complicated. The paper explicitly notes that this problem "remains an ongoing and formidable challenge."

This paper, written by Sandip Saha, Gautam Gangopadhyay, and Deb Shankar Ray, dives into this puzzle by looking at a specific family of equations called the Liénard–Levenson–Smith (LLS) system. Think of the LLS system as a master blueprint for building oscillators. The authors show how many different biological and physical systems—from the way fish get their spots to how cells divide—can be translated into this blueprint. They then use a set of mathematical "magnifying glasses" (called perturbative methods) to zoom in on these equations to estimate the upper bounds on how many rhythmic loops they can support and to design specific examples.

The researchers found that by tweaking the "knobs" (parameters) in these equations, they could systematically design systems that dance to two, three, or even more distinct rhythms simultaneously. They created new versions of famous oscillators (like the Van der Pol and Rayleigh models) that can hold multiple stable limit cycles, providing theoretical insights and concrete examples for these specific cases. However, they also hit a snag: they discovered that the standard rules used to predict these rhythms sometimes fail. In some cases, the math says a rhythm shouldn't exist, yet simulations show it clearly does. The paper concludes that while we can now build complex multi-rhythmic models and estimate their limits, the old mathematical rules need an update to explain these surprising new dances.

The Big Picture: Counting the Rhythms

At its heart, this paper is about organizing chaos. The authors start by taking messy, real-world equations that describe things like sugar breaking down in yeast (glycolysis) or how fish get their patterns, and they simplify them into a standard format. They call this the "LLS form." It's like taking a thousand different recipes for cakes and realizing they all use the same basic ingredients: flour, sugar, and eggs. Once you know the ingredients, you can predict how the cake will rise.

In this case, the "ingredients" are the forces that push the system (like a restoring force) and the forces that slow it down (like friction or damping). The paper classifies these systems into three main families:

The Liénard family: Where the slowing down depends on how far you are from the center.
The Rayleigh family: Where the slowing down depends on how fast you are moving.
The Extended family: A mix of both, which covers even more complex scenarios like the body's internal clock.

By converting real-world problems into these families, the authors can use powerful mathematical tools to predict the behavior of the system without having to simulate every single second of its life.

The Magic Tools: How They Count the Loops

To figure out how many rhythms a system can hold, the authors use two clever techniques called the Renormalization Group (RG) and the Krylov–Bogoliubov (KB) methods. Imagine you are trying to predict the path of a leaf floating in a river. The water is turbulent, and the leaf spins wildly. Instead of tracking the leaf every millisecond, these methods let you look at the "average" flow. They strip away the tiny, messy wiggles and focus on the big picture: is the leaf spiraling inward, outward, or settling into a perfect circle?

Using these tools, the authors derived "amplitude equations." Think of these as a scorecard that tells you the size of the rhythm. If the scorecard says the rhythm can be size 1, size 2, or size 3, then the system has three possible stable loops. The paper shows that the number of loops depends on the "powers" of the variables in the equation. For example, if the equation involves terms like $x^2$ or $x^4$ , the complexity increases, allowing for more potential rhythms.

The Discovery: Designing Multi-Rhythmic Systems

The most exciting part of the paper is the "Systematic Designing" section. The authors didn't just count existing rhythms; they built new ones. They took the famous Van der Pol oscillator (a classic model for heartbeats) and the Rayleigh oscillator (often used for sound waves) and added extra terms to them.

By carefully choosing the numbers (parameters) in these equations, they demonstrated that they could create systems with:

Bi-rhythmicity: Two stable rhythms co-existing.
Tri-rhythmicity: Three stable rhythms co-existing.

They showed that for a specific setup, a system could have five concentric loops (like rings on a tree trunk). Three of these rings would be stable (the system would happily stay there), and two would be unstable (the system would quickly move away from them). This is a big deal because it proves that nature could theoretically support complex, multi-layered rhythms, which might explain how biological systems manage to be so robust.

For instance, they looked at a model for glycolysis (sugar processing in cells) and showed that by adding a little bit of extra nonlinearity, you could turn a single-rhythm system into a double-rhythm one. This suggests that cells might have a hidden "backup" rhythm or a way to switch modes, which could be crucial for survival.

The Conflict: When the Rules Break

However, the paper also highlights a significant problem. There is a long-standing rule in this field that says: "If the damping force at the center is zero or positive, you can't have a limit cycle." It's like saying, "If there's no friction at the start, the dancer can't spin."

The authors found a conflict. They created a system where the math says the damping at the center is zero (which should mean no rhythm), but when they ran computer simulations, the system did develop a stable rhythm. It was as if the dancer started spinning even though the floor was perfectly slippery.

The paper argues that the old rule ( $F(0,0) < 0$ ) is too strict. It suggests that we need a new way to look at these systems, perhaps by checking how energy changes over a full cycle rather than just looking at the center point. They propose that even if the math at the very center looks "boring," the system can still find a rhythm if the energy balance over time works out. This is a crucial correction because it means we might have been missing valid rhythms in our models all along.

What This Means for the Future

The paper doesn't claim to have solved Hilbert's 16th Problem entirely—that's a mountain too high for one climb. Instead, it provides a better map for a specific part of the mountain. It shows that for these generalized systems, we can now estimate the upper bounds on the number of rhythms based on the complexity of the equation and construct specific examples.

The authors suggest that this framework could be a game-changer for biology and engineering. If we can design systems with multiple stable rhythms, we might be able to create better pacemakers that adapt to different heart rates, or understand how the brain switches between different states of consciousness. They even hint at using machine learning to find these rhythms in real data, moving from theory to practice.

In short, this paper takes a dry, abstract math problem and turns it into a toolkit for understanding the complex, multi-layered rhythms of life. It shows us that while the universe might be chaotic, it follows rules we can learn, and sometimes, those rules allow for more than one perfect dance.

Technical Summary: Exploring Multi-Rhythmicity in the Light of Hilbert's 16th Problem

Problem Statement
This systematic review addresses the second part of Hilbert's 16th problem, which concerns the existence and maximum number of limit cycles in planar polynomial vector fields. Specifically, the paper investigates the dynamics of nonlinear oscillatory kinetic equations in two variables, aiming to establish a systematic classification of these systems and determine the upper bounds on the number of limit cycles they can support. The authors focus on the transition from mono-rhythmic (single limit cycle) to multi-rhythmic (bi- and tri-rhythmic) behaviors, a phenomenon observed in various biological systems such as glycolysis, circadian rhythms, and cell cycles. A central challenge addressed is the limitation of the standard Liénard–Levenson–Smith (LLS) framework, which traditionally assumes the existence of a unique limit cycle, and the need to reconcile this with systems exhibiting multiple concentric limit cycles or cases where standard stability conditions (e.g., $F(0,0) < 0$ ) appear to fail yet limit cycles persist.

Methodology
The paper employs a multi-faceted approach combining mathematical reduction, perturbative analysis, and numerical simulation:

Reduction Scheme: The authors demonstrate how general autonomous kinetic equations (two-variable ODEs) can be transformed into a second-order differential equation in the generalized LLS form: $\ddot{\xi} + F(\xi, \dot{\xi})\dot{\xi} + G(\xi) = 0$ . This is achieved through linear transformations shifting the steady state to the origin.
Classification: Based on the polynomial structures of the damping function $F(\xi, \dot{\xi})$ $F (ξ, \dot{ξ})$ and the restoring force $G(\xi)$ $G (ξ)$ , systems are classified into three families:
- Liénard Family: Where the restoring force is a function of displacement only ( $G(\xi)$ ).
- Rayleigh Family: Where the damping and restoring forces depend primarily on velocity ( $F(\dot{\xi})$ and $G(\dot{\xi})$ ).
- Extended Liénard Family: Where the restoring force is velocity-dependent ( $G(\xi, \dot{\xi})$ ).
Perturbative Methods: To analyze amplitude-phase dynamics and estimate the number of limit cycles, the paper utilizes two primary techniques:
- Renormalization Group (RG) Method: Used to remove secular terms from perturbation series, deriving flow equations for amplitude and phase.
- Krylov–Bogoliubov (KB) Averaging Method: Applied to weakly nonlinear systems to derive amplitude equations. The authors analyze the roots of these amplitude equations to determine the number of limit cycles.
Systematic Design: The authors construct generalized models for Van der Pol (VdP) and Rayleigh oscillators by introducing higher-order polynomial terms in the damping and restoring forces. They systematically tune parameters to achieve specific numbers of limit cycles (bi-rhythmicity and tri-rhythmicity).
Energy Analysis: To address conflicts where standard LLS conditions (like $F(0,0) < 0$ ) are not met, the paper employs an energy consumption analysis per cycle to verify the existence of limit cycles in generalized systems.

Key Contributions and Results

Unified Framework: The paper establishes that diverse biological oscillators (glycolysis, fish pigmentation, bacterial respiration, cell cycles, circadian rhythms, and Ca $^{2+}$ oscillations) can be universally represented within the LLS framework, providing a common mathematical language for these phenomena.
Upper Bound Estimation: Using the KB method, the authors derive a relationship between the polynomial degrees of the damping ( $M$ ) and restoring ( $N$ ) functions and the maximum number of limit cycles. They provide tables (Table 1 and Table 2) estimating the maximum number of distinct conjugate roots (and thus limit cycles) for various combinations of $N$ and $M$ . For instance, for a Rayleigh system with $N=1$ , the maximum number of limit cycles is determined by the degree of the velocity polynomial.
Design of Multi-Rhythmic Models:
- VdP Family: The authors successfully designed generalized VdP oscillators (Eq. 97 and 99) capable of exhibiting three stable limit cycles (tri-rhythmicity) by incorporating higher-order terms ( $x^8, x^{10}$ ) in the damping function.
- Rayleigh Family: Similarly, generalized Rayleigh oscillators (Eq. 101 and 103) were constructed to support three stable limit cycles by extending the velocity-dependent damping terms.
- Parameter Spaces: The study delineates specific parameter regions (e.g., $\alpha, \beta, \gamma, \delta$ ) where mono-, bi-, and tri-rhythmic behaviors occur, supported by bifurcation diagrams and phase portraits showing concentric stable and unstable limit cycles.
Resolution of LLS Limitations: The paper identifies a conflict where certain generalized systems (e.g., Eq. 105 and 106) possess stable limit cycles despite violating the standard condition $F(0,0) < 0$ (where $F(0,0)=0$ or is undefined). Through energy analysis, the authors demonstrate that these systems still exhibit limit cycle behavior, suggesting the need for a generalized condition beyond the traditional LLS theorem.

Significance and Claims
The authors claim that this work provides a comprehensive framework for understanding nonlinear oscillations by bridging mathematical theory (Hilbert's 16th problem), perturbative analysis, and real-world biological applications. The significance lies in:

Systematic Construction: Offering a recipe for constructing minimal models with a desired number of limit cycles, moving beyond the traditional assumption of a unique limit cycle.
Theoretical Insight: Contributing to the ongoing discourse on Hilbert's 16th problem by providing specific cases and upper bound estimates for polynomial vector fields in the context of LLS systems.
Interdisciplinary Applicability: Demonstrating the universality of the LLS condition in modeling complex rhythmic processes in physics, biology, and engineering.
Future Directions: The paper modestly suggests that these generalized LLS systems serve as testbeds for exploring complex dynamics, including stochastic effects, time delays, and potential integration with machine learning for equation discovery, though it does not present these as completed results.

The review concludes that while the LLS framework has limitations regarding the uniqueness of limit cycles, its extension through perturbative methods and higher-order generalizations offers a powerful tool for analyzing and designing multi-rhythmic systems.

The Big Picture: Counting the Rhythms

The Magic Tools: How They Count the Loops

The Discovery: Designing Multi-Rhythmic Systems

The Conflict: When the Rules Break

What This Means for the Future

More like this