🌀 nonlinear sciences

Extreme events in MLC circuit

This study elucidates the mechanisms behind extreme events in the Murali-Lakshmanan-Chua (MLC) circuit, revealing that an external periodic force drives large chaotic deviations via PM intermittency and manifold dynamics, while statistical analysis confirms these rare occurrences follow generalized Pareto and extreme value distributions.

Original authors: Tapas Kumar Pa, Dibakar Ghosh

Published 2026-04-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Tapas Kumar Pa, Dibakar Ghosh

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 watching a tiny, chaotic electronic circuit called the MLC circuit. Think of this circuit not as a boring box of wires, but as a wild, unpredictable dancer spinning on a stage. Usually, this dancer moves in a complex but contained pattern, staying within a specific area of the stage (this is called a "chaotic attractor").

However, sometimes, without warning, this dancer suddenly leaps wildly off the stage, travels a huge distance, and then snaps back. In the world of science, these massive, rare jumps are called Extreme Events (EEs). They are like a sudden, massive tsunami in a calm ocean or a stock market crash on an otherwise stable day.

This paper by Tapas Kumar Pal and Dibakar Ghosh investigates why this electronic dancer makes these terrifying leaps. They found that the circuit doesn't just randomly jump; there is a specific choreography behind the chaos.

Here is the breakdown of their discovery using simple analogies:

1. The Setup: The Dancer and the Wind

The circuit is being pushed by an external force, like a periodic wind blowing on our dancer. The researchers turned the "wind speed" (the amplitude of the force) up and down slightly.

The Discovery: They found that when the wind speed hits a very specific, narrow range, the dancer stops doing small, contained spins and starts making massive, dangerous leaps.
The Trigger: This happens through a process called PM Intermittency. Imagine the dancer is walking a tightrope. Suddenly, the tightrope breaks, and the dancer falls into a huge, wild loop before catching the rope again. This "breaking of the tightrope" is the precursor to the extreme event.

2. Why Does the Dancer Jump? (Three Explanations)

The authors didn't just say "it happens." They explained how it happens using three different perspectives, like looking at the same storm through three different windows.

A. The Force Field (The Invisible Hand)

Imagine the stage isn't flat; it has invisible "wind zones."

The Analogy: Most of the time, the dancer is in a gentle breeze (low-intensity force field). But suddenly, the dancer drifts into a zone where a massive hurricane (high-intensity force field) is waiting.
The Result: The hurricane grabs the dancer and violently pulls them far away from their usual spot. This "pull" is the extreme event. The external force acts like a magnet that suddenly snaps the trajectory away from its home.

B. The Stable and Unstable Roads (The Manifold Map)

Imagine the stage has two types of roads painted on the floor:

The Stable Road (Red): A smooth, sticky path that keeps the dancer moving in a circle.
The Unstable Road (Green): A slippery, repelling path that pushes the dancer away.
The Mechanism: Usually, the dancer stays on the Stable Road. But because the two roads are right next to each other, the dancer sometimes accidentally steps onto the edge of the Unstable Road. The moment they touch it, the road acts like a springboard, launching them far away. The researchers used math (Floquet multipliers) to map exactly where these roads are and how they interact to cause the launch.

C. The Slow-Fast Dance (The Rollercoaster)

The circuit has parts that move fast and parts that move slow.

The Analogy: Think of a rollercoaster. The "slow" part is the climb up the hill (the stable path). The "fast" part is the drop.
The Mechanism: The dancer slowly climbs up the "stable" hill. But at the very top, there is a cliff (the unstable path). Once the dancer reaches the edge of the cliff, they can't stay there anymore. They are ejected violently down the other side. This "ejection" is the extreme event.

3. The Statistics: Predicting the Rare

The researchers also looked at the data like a weather forecaster. They asked: "How often do these jumps happen, and how big are they?"

The Findings: They found that these extreme jumps follow specific mathematical rules (called Generalized Pareto and Generalized Extreme Value distributions).
What this means: Even though the jumps look random, they aren't totally random. They follow a pattern. If you know the rules of the game, you can calculate the probability of a "tsunami" happening, even if you can't predict the exact second it will occur.

The Big Picture

This paper is important because it shows us that chaos isn't just messy; it has a structure.

By understanding the "wind zones" (force fields) and the "roads" (manifolds) in this simple electronic circuit, scientists hope to understand how extreme events happen in much bigger, more complex systems—like how a power grid might crash, how a stock market might collapse, or how a weather system might create a super-storm.

In short: The MLC circuit is a laboratory for studying how small, regular systems can suddenly explode into chaos, and the authors have mapped the exact mechanics of that explosion.

1. Problem Statement

The paper investigates the emergence of Extreme Events (EEs) within the Murali-Lakshmanan-Chua (MLC) circuit, a well-known nonautonomous, dissipative, and nonlinear electronic system. While EEs are recognized in various natural and anthropogenic systems (e.g., tsunamis, stock crashes), their specific dynamical genesis in the MLC circuit had not been fully elucidated. The study aims to:

Determine if EEs occur in the MLC circuit under varying external forcing amplitudes.
Identify the specific dynamical mechanisms responsible for the sudden, large deviations of chaotic trajectories that characterize EEs.
Statistically validate the rarity and distribution of these events using Extreme Value Theory (EVT).

2. Methodology

The researchers employed a combination of numerical simulation, dynamical systems analysis, and statistical theory:

System Model: The study utilizes the MLC circuit equations:
$\dot{x} = f \sin(\omega t) - x - y - g(x)$
$\dot{y} = \beta x$
Where $g(x)$ is a piecewise linear function representing the Chua diode.
Parameter Variation: The amplitude ( $f$ ) of the externally applied periodic force was varied in the range $[0.1312, 0.1315]$ while keeping other parameters fixed ( $a = -1.121, b = -0.6047, \beta = 0.05, \omega = 0.265$ ).
Event Definition: Two types of events were analyzed based on the observable $x$ $x$ :
- $x_{max}$ : Local maxima exceeding a threshold $H_{th} = m + 6\mu$ (where $m$ is the mean and $\mu$ is the standard deviation).
- $x_{min}$ : Local minima falling below a threshold $H_{th} = m - 6\mu$ .
Dynamical Analysis Techniques:
1. Bifurcation Analysis: To map the transition from periodic to chaotic states and identify the onset of EEs.
2. Force Field Visualization: Analyzing the synergistic effect of the external periodic force on the phase space trajectory.
3. Manifold Decomposition (Floquet Multipliers): Transforming the nonautonomous system into an autonomous one to calculate stable and unstable manifolds using Floquet theory.
4. Slow-Fast Dynamics: Decomposing the critical manifold to understand the interaction between stable (attracting) and unstable (repelling) segments.
Statistical Analysis:
- Threshold Excess Values: Analyzed using the Generalized Pareto Distribution (GPD).
- Inter-Extreme-Spike-Intervals (IESI): Analyzed using the Generalized Extreme Value (GEV) distribution.
- Goodness of Fit: Validated using Kolmogorov-Smirnov (K-S) tests, P-P plots, and Q-Q plots.

3. Key Contributions

The paper makes three primary theoretical contributions regarding the mechanism of EE genesis in nonautonomous systems:

Identification of the Precursor: The study confirms that EEs in the MLC circuit emerge via the Pomeau-Manneville (PM) intermittency route. A sudden expansion of the chaotic attractor occurs near $f \approx 0.131451$ , transitioning from a one-periodic orbit to large-amplitude chaos.
Dynamical Mechanism - Force Field Synergy: The authors propose that the external periodic force creates a "force field" in the phase space. When a chaotic trajectory moves from a low-intensity force region to a high-intensity region, it experiences a sudden "pull" or deflection, causing the large excursion characteristic of an EE.
Geometric Mechanism - Manifold Interaction:
- Floquet Multipliers: The coexistence of stable (attracting) and unstable (repelling) manifolds creates a dynamic environment where trajectories are ejected from the stable region.
- Slow-Fast Dynamics: The chaotic trajectory wanders on the stable manifold but is abruptly repelled when it approaches the unstable manifold (or fold curves), leading to the extreme spike. This mechanism was successfully applied to both $x_{max}$ and $x_{min}$ events.

4. Results

Dynamical Behavior:
- For $f$ values in the non-extreme chaotic region (e.g., $f=0.131228$ ), no EEs were observed.
- For $f$ values in the extreme region (e.g., $f=0.1314, 0.131417$ ), the time series exhibited rare, high-amplitude spikes crossing the $6\sigma$ threshold.
- Phase space analysis revealed that these spikes correspond to trajectories being ejected from the main chaotic attractor due to the interaction with the external force field and unstable manifolds.
Statistical Validation:
- Threshold Excess Values: The data for both $x_{max}$ and $x_{min}$ events followed the Generalized Pareto Distribution (GPD) with high confidence (K-S test p-values indicated failure to reject the null hypothesis).
- Inter-Extreme-Spike-Intervals (IESI): The time intervals between consecutive extreme events followed the Generalized Extreme Value (GEV) distribution.
- The statistical analysis confirmed that these events are indeed "extreme" (rare and significantly deviated from the mean) and not merely noise.

5. Significance

Mechanistic Insight: The paper provides a comprehensive geometric and dynamical explanation for how EEs are generated in nonautonomous systems, moving beyond mere observation to mechanistic understanding involving force fields and manifold interactions.
Methodological Framework: It demonstrates a robust framework for analyzing EEs in nonautonomous systems by transforming them into autonomous forms to apply Floquet theory and slow-fast analysis.
Predictive Potential: By identifying the PM intermittency route as a precursor and characterizing the statistical distributions (GPD/GEV), the study offers potential pathways for early warning systems and risk mitigation in nonlinear electronic circuits and analogous physical systems.
Future Directions: The authors suggest that understanding the magnitude and direction of the "force field" and the specific nature of non-Gaussian distributions near unstable manifolds are promising avenues for future research.

In conclusion, the study successfully links the macroscopic observation of extreme events in the MLC circuit to microscopic dynamical mechanisms (intermittency, force fields, and manifold geometry), providing a holistic view of EE genesis in driven nonlinear systems.