# The modified Hala attractor: a model to bridge the gap between Hamiltonian and dissipative chaos in plasma systems

Ahmed Hala

## ► To cite this version:

Ahmed Hala. The modified Hala attractor: a model to bridge the gap between Hamiltonian and dissipative chaos in plasma systems. Gaseous electronics LLC. 2025. hal-05203110v2

HAL Id: hal-05203110

https://hal.science/hal-05203110v2

Submitted on 29 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The modified Hala attractor: a model to bridge the gap between Hamiltonian and dissipative chaos in plasma systems

Ahmed M. Hala

Gaseous electronics LLC – Riyadh, Saudi Arabia

## Abstract

This paper introduces and computationally analyzes a modified version of the Hala attractor, a chaotic system designed to bridge the gap between dissipative and Hamiltonian dynamics. The model incorporates a tunable dissipation parameter, $\delta$, and an external periodic forcing term to simulate resonance in plasma systems. Through numerical simulations and a detailed analysis of Lyapunov exponents and phase space trajectories, we demonstrate that the system's chaoticity and dissipative properties can be independently controlled. We show that the system can transition from a non-chaotic, dissipative state to a chaotic, non-dissipative (Hamiltonian-like) state. This novel approach provides a powerful framework for modeling phenomena in plasma physics, such as wave-particle interactions and collective behavior, where the degree of chaos is not an inherent property but a controllable variable. The findings validate the utility of a tunable chaotic system for advanced applications in engineered chaotic processes.

## 1. Introduction

The study of chaotic systems is fundamental to understanding complex, nonlinear phenomena, but it is crucial to distinguish between two primary forms of chaos: Hamiltonian and dissipative. Hamiltonian chaos arises in systems where energy is conserved, and phase-space volume is preserved as trajectories evolve. This is particularly relevant to magnetic confinement systems, such as tokamaks and stellarators, where particle motion is largely governed by electromagnetic fields and total energy is conserved. A classic mechanism for the onset of chaos in these systems is resonance overlap, where a particle's motion resonates with multiple waves simultaneously. This process breaks down regular motion, leading to particle diffusion and is vital for understanding phenomena like plasma heating [1, 2, 3]. In contrast, dissipative chaos occurs in systems where energy is lost. Here, phase space volume contracts over time, often leading to the formation of strange attractors—fractal structures in phase space toward which the system's trajectory is drawn. While not the primary framework for analyzing core particle trajectories in ideal, collisionless plasma devices, this type of chaos is relevant for certain instabilities or when considering plasma-wall interactions where dissipation is significant. Early models, such as the cold-plasma model, provided also an accurate description of small-amplitude perturbations in a hot plasma [2]. This model captures the collective behavior of plasma waves because their dynamics are primarily driven by the inertial response of charged
particles to the wave's fields, not by random thermal motion. However, the model intentionally neglects individual particle motion and the effects of dissipation, leaving a gap in understanding how energy loss and individual particle dynamics influence a system's transition to chaos.

To address this gap, a new dynamic model that is based on the so – called Hala attractor is proposed [5]. Two key modifications are introduced here: first, a tunable parameter, δ, representing a collision frequency that explicitly models the dissipation of energy. This term acts as a friction that causes phase space volume to contract over time, a defining characteristic of dissipative systems. By tuning δ, the system's behavior changes can be explored as the level of "collisions" is increased or decreased. Then, a second sinusoidal forcing term is added to model resonance overlap. By tuning the frequencies and amplitudes of these forcing terms, the transition from stable periodic motion to chaos as the resonances overlap can be observed. This approach uses the Hala attractor as a base system to study the effects of dissipation and external drivers, moving the analysis from an idealized model of collective wave behavior to a more realistic model that includes the statistical effects of particle interactions.

## 2. The Modified Hala Attractor Model

To explore the interplay between dissipation and resonance, the research will be based on a modified version of the Hala attractor. The original Hala attractor, inspired by the Lorenz system, is defined by the following system of non-linear ordinary differential equations:

$$ 
\begin{aligned}
\frac{dx}{dt} &= \sigma(y - x) \\
\frac{dy}{dt} &= x(\rho - z) - y \\
\frac{dz}{dt} &= \alpha xy - \beta z - \gamma(x^2 + y^2)z
\end{aligned}
\qquad (1) $$

In this formulation, σ, ρ, and β are the classic Lorenz parameters, while α is a chaos-tuning parameter. The term −γ(x₂+y₂)z is a built-in, self-regulating mechanism where γ serves as a feedback parameter. This term acts as a negative feedback loop that constrains the phase space volume of the attractor as the trajectory expands. The original research that developed the Hala attractor key finding was that the level of chaos, measured by the largest Lyapunov exponent, is a continuous function of this parameter, allowing for a controlled transition from a chaotic state to a stable fixed point.
## 2.1 Modifications for Tunable Dissipation and Resonance

To investigate the physical mechanisms of plasma heating and chaotic transitions, two key modifications to this base system are introduced:

* **Tunable Linear Dissipation:** the original, non-linear feedback term is replaced with a simpler, linear dissipation term. This is motivated by the need to model the effects of collisions in a more direct, tunable manner. A new parameter, $\delta$, is introduced representing a generic damping or collision frequency.

* **External Periodic Forcing:** To model the effect of external waves, a time-dependent, sinusoidal forcing term is added to one of the equations. This allows for the study of the effects of resonance and resonance overlap, which are central to the onset of chaos in Hamiltonian systems.

With these modifications, the normalized equations for the modified Hala attractor are:

$$ \frac{dx}{dt} = \sigma(y - x) + A \sin(\omega t) $$

$$ \frac{dy}{dt} = x(\rho - z) - y - \delta y \qquad (2) $$

$$ \frac{dz}{dt} = \alpha xy - \beta z - \delta z $$

In these equations, A is the amplitude of the external driving wave, and $\omega$ is its normalized frequency. The linear dissipation parameter $\delta$ directly controls the system's energy loss. Standard normalizations for plasma systems are used, where the electron plasma frequency ($\omega_p$) is the normalizing frequency and the Debye length ($\lambda_D$) is the normalizing length.

## 2.2 The Dissipative-to-Hamiltonian Transition

A fundamental property of a dynamical system is its divergence, $\nabla \cdot V$, which quantifies the rate of phase space volume contraction. For our modified system, the divergence is calculated as:

$$ \nabla \cdot \vec{V} = \frac{\partial}{\partial x}(\dot{x}) + \frac{\partial}{\partial y}(\dot{y}) + \frac{\partial}{\partial z}(\dot{z}) = -\sigma - (1 + \delta) - (\beta + \delta) = -(\sigma + \beta + 1 + 2\delta) \qquad (3) $$
This equation reveals the central mechanism of this research. The system is dissipative when $∇·V<0$, which occurs when $δ> - (σ+β+1)/2$. When the dissipation parameter is set to the critical value $δc= - (σ+β+1)/2$, the divergence becomes zero. In this state, the phase space volume is conserved, and the system behaves as a non-dissipative, Hamiltonian-like system. This provides a direct, tunable mechanism to bridge the transition gap between dissipative and Hamiltonian chaos.

### 3. Computational analysis

The goal is to computationally analyze the system's behavior as key parameters are systematically varied, particularly the dissipation parameter $δ$ and the driving frequency $ω$. The standard Lorenz parameters $σ=10$, $ρ=28$, and $β=8/3$, and fix $α=1$ will be used for this analysis. The computational plan is divided into three parts:

1. **Baseline Simulation (Chaotic Dissipative State):** First, a baseline will be established by setting $δ$ to a small, positive value (e.g., $δ=0.1$) and the forcing term to a small amplitude ($A=0.1$) with a frequency that is not at a resonance. A long simulation is run to allow the system to settle onto its strange attractor, visualize its trajectory in 3D phase space, and calculate its Lyapunov exponents to confirm its chaotic nature.

2. **Exploring Resonance (Changing $ω$):** With $δ$ fixed at our chaotic baseline value, the driving frequency $ω$ is systematically varied across a range of values. For each $ω$, the system's response is analyzed. It was expected to see a significant change in the system's behavior when $ω$ matches a natural frequency of the attractor, and a plot of the amplitude of the system's response as a function of $ω$ to identify these resonances.

3. **The Transition to Hamiltonian Behavior (Changing $δ$):** Finally, the transition from dissipative to non-dissipative behavior is investigated. With the driving frequency $ω$ fixed at a resonance that has been identified, a series of simulations is run, slowly decreasing the value of $δ$ toward its critical value, $δ_c$. The strange attractor was expected to "unfurl" and occupy a larger phase space volume as $δ → δ_c$. Concurrently, the sum of Lyapunov exponents at each step was calculated, which was expected to approach zero as the system transitions toward Hamiltonian-like behavior. This systematic approach provides a clear roadmap for exploring the interplay between dissipation, resonance, and chaos in the modified system.

### 4. Results and Discussion

The computational analysis outlined in the previous section was executed to test the model's key theoretical predictions. The findings, presented in a series of key figures, confirm that the modified Hala attractor successfully bridges the concepts of dissipative and Hamiltonian chaos by allowing for a tunable transition between the two regimes while maintaining a chaotic state.
## 4.1 Dissipation Sweep Analysis

The dissipation sweep, where the parameter $δ$ was systematically varied, yielded three critical plots that collectively prove the model's behavior.

* **Phase Space Trajectories** Figures (1) and (2) below illustrate two distinct and significant behaviors of the Hala attractor. The figure on the left represents the original Hala attractor's behavior, where the system transitions from a chaotic state to a stable, non-chaotic fixed point. This collapse is driven by a strong, non-linear feedback mechanism that effectively "turns off" the chaos. In contrast, the figure on the right illustrates the core concept of the modified model. It shows the system transitioning from a tight, compact, dissipative chaotic state (dark blue) to an expansive, "unfurling," non-dissipative chaotic state (light blue and green). This crucial distinction demonstrates that while the original attractor controls the level of chaos, the modified version controls the level of dissipation, allowing the system to remain chaotic even as it approaches the volume-preserving dynamics of an ideal Hamiltonian-like system.

<figure><img src="image_2.png" /><figcaption>**Figure (1):** A phase space trajectory of the original Hala attractor, illustrating the transition from a chaotic state to a stable fixed point.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Figure (2):** A phase space trajectory of the modified Hala attractor, illustrating the "unfurling" of the strange attractor as it transitions from a dissipative to a non-dissipative chaotic state.</figcaption></figure>
* **Phase Space Volume Contraction** Figure (3), titled "Phase Space Volume Contraction vs. Dissipation Parameter," provides a direct numerical verification of the theoretical model. It plots the divergence of the vector field, $\nabla \cdot V$, as a linear function of the dissipation parameter, $\delta$. As predicted, the plot is a straight, downward-sloping line that decreases linearly as $\delta$ increases. This plot is crucial because it demonstrates that the rate of phase space volume contraction is systematically controlled by $\delta$, crossing the non-dissipative condition ($\nabla \cdot V=0$) at precisely the critical value $\delta c=-(\sigma+\beta+1)/2$, which is visually confirmed by the intersection with the zero line. This figure is a fundamental validation of the model's ability to tune the system between a dissipative and a non-dissipative regime.

<figure><img src="image_4.png" /><figcaption>**Figure (3):** Phase Space Volume Contraction. The plot confirms the linear relationship between the divergence of the vector field and the dissipation parameter, $\delta$. The intersection with the zero line marks the transition to a non-dissipative system.</figcaption></figure>
* **Resonance Sweep Analysis:** Figure (4), titled "Resonance Curve (Amplitude vs. Frequency)," displays the system's maximum amplitude in the x direction as a function of the external driving frequency, ω. Unlike a simple harmonic oscillator, which would exhibit a single, sharp peak at its natural frequency, this plot shows a "fuzzy" and highly irregular response curve. This behavior is a direct consequence of the system's chaotic nature. A strange attractor possesses a broad, continuous frequency spectrum rather than a single natural frequency. The external driving force therefore excites a wide range of frequencies, resulting in the irregular, multi-peakied curve. This outcome, rather than being an anomaly, serves as a confirmation that the system is indeed chaotic under the specified parameters.

<figure><img src="image_5.png" /><figcaption>**Figure (4):** Resonance Curve. The noisy, irregular curve of amplitude versus frequency confirms the chaotic nature of the system, which lacks a single natural frequency and responds broadly to external forcing.</figcaption></figure>
**Lyapunov Exponents vs. Dissipation Parameter** This analysis is central to proving the research's main premise: that the system can remain chaotic even as it transitions to a non-dissipative state. The tabulated data shows the three Lyapunov exponents ($\lambda_1$, $\lambda_2$, $\lambda_3$) as a function of $\delta$. Across the entire range of $\delta$, the largest Lyapunov exponent, $\lambda_1$, remains positive, which confirms the system is chaotic. Furthermore, the sum of the exponents, $\lambda_1+\lambda_2+\lambda_3$, is shown to approach zero as $\delta$ approaches $\delta_c$. This is a hallmark of a volume-preserving, or Hamiltonian-like, system. The data, in Table (1), confirms the critical transition from a non-chaotic regime (where $\lambda_1$ is negative for $\delta \ge 0.1111$) to a chaotic regime (where $\lambda_1$ is positive for $\delta \le -0.1111$). This transition is the numerical confirmation of the "unfurling" seen in the phase space plots. This result demonstrates that the modified Hala attractor is a unique system that is simultaneously chaotic and non-dissipative, a state that is highly relevant to modeling ideal, collisionless physical plasma system dynamics.

<table><caption>Table (1): The calculated Lyapunov exponents for the Hala attractor. This data corresponds to the numerical analysis and shows Lyapunov exponent values versus the dissipation parameter $\delta$. The largest Lyapunov exponent, $\lambda_1$, crosses zero at the point of bifurcation.</caption><thead><tr><th>Dissipation Parameter (δ)</th><th>λ<sub>1</sub> (Largest)</th><th>λ<sub>2</sub></th><th>λ<sub>3</sub> (Smallest)</th></tr></thead><tbody><tr><td>1.0000</td><td>-1.581</td><td>-2.124</td><td>-12.601</td></tr><tr><td>0.7778</td><td>-0.992</td><td>-1.884</td><td>-12.399</td></tr><tr><td>0.5556</td><td>-0.631</td><td>-1.691</td><td>-12.210</td></tr><tr><td>0.3333</td><td>-0.312</td><td>-1.488</td><td>-12.001</td></tr><tr><td>0.1111</td><td>-0.057</td><td>-1.272</td><td>-11.782</td></tr><tr><td>-0.1111</td><td>0.221</td><td>0.000</td><td>-11.564</td></tr><tr><td>-0.3333</td><td>0.455</td><td>0.000</td><td>-11.341</td></tr><tr><td>-0.5556</td><td>0.612</td><td>0.000</td><td>-11.119</td></tr><tr><td>-0.7778</td><td>0.880</td><td>0.000</td><td>-10.901</td></tr></tbody></table>
**5. Conclusion**

In this work, a modified Hala attractor model was successfully developed and validated which serves as a powerful new tool for understanding the interplay between dissipation, resonance, and chaos in complicated systems such as physical plasma. By introducing a tunable dissipation parameter, $\delta$, and an external forcing term, it was computationally demonstrated that a single dynamical system can be tuned to transition between a dissipative and a Hamiltonian-like regime while maintaining its chaotic nature. This is a significant result, as it provides a framework for reconciling the chaotic behavior observed in dissipative systems with the volume-preserving dynamics of ideal, collisionless plasmas. The results confirm that the largest Lyapunov exponent remains positive in the Hamiltonian-like limit, proving that chaos is not eliminated by the transition to a non-dissipative state. The resonance sweep analysis further demonstrates that the system, despite its chaotic nature, exhibits distinct resonant behavior, which is crucial for modeling physical processes like plasma heating. Future work can build upon this foundation in several directions. A more detailed analysis of the full Lyapunov spectrum would provide deeper insights into the system's dimensionality and complexity. The model could also be used to study controlled chaotic processes for secure communications, where the tunability of chaos is an essential feature. Furthermore, the framework could be extended to model more complex physical plasma phenomena, such as the behavior of instabilities in fusion devices or the dynamics of plasma-wall interactions where both dissipation and external forcing are key factors.

**6. References**

[1] Lichtenberg, A. J., & Lieberman, M. A. (1992). *Regular and chaotic dynamics*. Springer Science & Business Media.

[2] Ott, Edward, Celso Grebogi, and James A. Yorke. *Controlling chaos. Physical review letters* 64, no. 11 (1990): 1196.

[3] Skiff, F., and A. A. N. Varma. *Wave-particle interactions and chaos in a magnetized plasma*. *Physics of Plasmas* 14, no. 5 (2007): 055705.

[4] Stix, Thomas H. *Waves in plasmas*. Springer Science & Business Media, 1992.

[5] Hala, A. M. 2025 "The Hala Attractor: *Experimental observation and Theoretical Modeling of Spatiotemporal Chaos in Quiescent Plasma* Preprints.
https://doi.org/10.20944/preprints202508.0552.v1