Bifurcations and Hyperchaos in Mathematical Models of Sleep
Adriano Scibilia1, Luigi Fortuna2,3
1Institute of Intelligent Industrial Systems and Technologies for Advanced Manufacturing, National Research Council, 20133 Milan, Italy.
Bioengineering (Basel, Switzerland)
|July 28, 2026
Summary
This study uses nonlinear dynamics to model sleep regulation, revealing complex behaviors like hyperchaos. The findings offer a new framework for understanding irregular sleep patterns and developing digital sleep technologies.
Area of Science:
- Neuroscience
- Complex Systems Science
- Computational Biology
Background:
- Sleep regulation involves intricate nonlinear feedback loops, including homeostatic pressure, circadian rhythms, and neuromodulation.
- Existing sleep models often simplify these complex dynamics, necessitating advanced analytical approaches.
Purpose of the Study:
- To present a nonlinear-dynamical interpretation of sleep regulation, emphasizing bifurcation structure and Lyapunov stability.
- To analyze sleep-stage dynamics using advanced mathematical frameworks and computational models.
- To explore the potential of bifurcation theory in understanding irregular sleep transitions and informing future sleep technologies.
Main Methods:
- Summarized representative sleep models and analyzed sleep-stage dynamics literature.
- Combined bifurcation analysis and Lyapunov exponent calculations on sleep-population models.
- Developed and analyzed a four-dimensional slow-fast feedback model incorporating cortical activation, sleep-promoting activity, homeostatic drive, and circadian pacemaker.
Main Results:
- Identified routes from regular oscillations to period-doubling, chaotic bands, and positive Lyapunov regions using baseline diagrams.
- Lyapunov maps of the four-dimensional model revealed parameter regions supporting hyperchaotic dynamics (positive largest and second-largest Lyapunov exponents).
- Demonstrated how chronic sleep debt and insomnia can be modeled as shifts in homeostatic drive, inhibitory gain, and circadian coupling.
Conclusions:
- Bifurcation theory provides a robust framework for interpreting complex and irregular sleep transitions.
- The nonlinear-dynamical approach offers insights into the fragmentation and intermittency observed in sleep.
- This research paves the way for novel sleep-technology and digital twin applications based on a deeper understanding of sleep dynamics.
Related Concept Videos
Sleep-Wake Cycles
Sleep is an essential physiological process vital to maintaining overall well-being. The reticular activating system (RAS), a network of neurons in the brainstem, regulates wakefulness and sleep. While it may seem passive, sleep consists of distinct cycles, each with its unique characteristics and functions. Two key sleep phases are non-rapid eye movement (NREM) and rapid eye movement (REM).
NREM Sleep
NREM sleep comprises four progressive stages that seamlessly merge:
NREM Sleep
NREM sleep comprises four progressive stages that seamlessly merge:
Limits with Oscillating Discontinuities
An oscillating discontinuity is a type of discontinuity in which a function’s values fluctuate infinitely often as the input approaches a particular point. Unlike jump discontinuities, where the function suddenly shifts between two values, or infinite discontinuities, where the function diverges without bound, an oscillating discontinuity arises from rapid back-and-forth variation. Because the function never stabilizes toward a single value, no finite limit exists at that point.One of the most...
Stages of Sleep
Sleep progresses through distinct stages, each characterized by specific brain wave patterns and physiological responses ranging from wakefulness to stages of non-rapid eye movement, known as non-REM, to rapid eye movement, referred to as REM. Understanding these stages helps in recognizing how sleep supports various bodily and cognitive functions.
Before sleep begins, in wakefulness, the brain exhibits primarily beta waves, which are high in frequency and low in amplitude, indicating alertness...
Before sleep begins, in wakefulness, the brain exhibits primarily beta waves, which are high in frequency and low in amplitude, indicating alertness...
The Squeeze Theorem
Certain mathematical functions exhibit unpredictable or highly variable behavior near specific input values, making direct evaluation of their limits challenging. This complexity may arise from rapid oscillations or irregular patterns that obscure the function’s trend. In such cases, the Squeeze Theorem offers a reliable method for determining limits.According to the Squeeze Theorem, if a function is confined between two other functions near a particular point, and both outer functions approach...
BIBO stability of continuous and discrete -time systems
System stability is a fundamental concept in signal processing, often assessed using convolution. For a system to be considered bounded-input bounded-output (BIBO) stable, any bounded input signal must produce a bounded output signal. A bounded input signal is one where the modulus does not exceed a certain constant at any point in time.
To determine the BIBO stability, the convolution integral is utilized when a bounded continuous-time input is applied to a Linear Time-Invariant (LTI) system.
To determine the BIBO stability, the convolution integral is utilized when a bounded continuous-time input is applied to a Linear Time-Invariant (LTI) system.
Geometry of Hyperbolas
A hyperbola consists of all points where the absolute difference of distances to two fixed points, called foci, remains constant. The standard equation isEach branch extends infinitely and approaches two asymptotes, which guide the curve’s behavior. The parameters a and b define key features: a measures the distance from the center to each vertex along the transverse axis, while b influences the slopes of the asymptotes. The asymptotes have equationsA rectangle centered at the origin with...


