1. Introduction

Introduction

The Hala attractor was introduced in August 2025 as a Lorenz-like system in which the degree of chaos is a tunable quantity, not a fixed property ^[1]. Two companion studies followed in the same month. One adds a dissipation parameter and periodic forcing to connect dissipative and volume-preserving regimes ^[2]. The other relates the model to Langmuir-probe current–voltage traces from a quiescent multidipole plasma ^[3].

The body of work on the Hala attractor, developed by Dr. Ahmed M. Halasma motiva, represents a contribution is older than the model. A 2020 letter reafield of nonlinear dynamics and chaos from probe I–V characteristics through a discrete recurrence fitted to the electron-saturation region ^[4], and a 2025 prtheory. It highlights a different treatment than the conventional vieprintw developed the electric probeof chaos as an instrument for that purpostrinsic, uncontrollable feature ^[5].of A rcecurrence obtained by discretising a one-dimensional I–V curve is a rtain systems, instead reframing it as a tunable, controllable property of the chosen discretisation. It is not by itself evidence of. This cohesive research direction, spanning three studies, demonstrates that chaos in the plasma, and the entry treats it as motivation only.

Sincs not a binary, all-or-nothing state Janbuary 2026 the feedback term has been studied as an operator in its own right, separate from the Lorenz host. The body of work now has seven layers: precursors, the attractor, the operator, control protocols, a log-frequency extension, applications, and methods for auditing claims. This entry covers the first four in detail and summarises the rest.

2. Definition of the system

The host is the Lorenz system ^[6] with the Hala operator added to the z-equation:

Where the Hala operator is defined as:

The reference parameters are σ = 10, ρ = 28, α = 1, β = 8/3, and γ ≥ 0 is the quench factor. Three properties follow directly ^[7].

• The origin is untouched. The operator is cubic, so the linearisation at the origin is that of the Lorenz system for every γ.
• Contraction depends on the state. The divergence of the flow is −(σ + 1 + β) − γ(x² + y²), which is negative everywhere and largest in magnitude far from the z-axis.
• The nontrivial equilibria move outward. The equilibria C± keep y = x and z* = ρ − 1 = 27, with

Figure 1. Phase-space projections (x, z) of the Hala attractor at three quench factors, from the initial condition (0.1, 0.1, 0.1). At γ = 0 the Lorenz attractor is recovered. At γ = 10⁻⁴ the attractor is visually unchanged. At γ = 0.01 the trajectory spirals into the stable equilibrium C+ (marked), where the largest Lyapunov exponent is −1.23.

A modified formulation adds a dissipation parameter δ and external periodic forcing ^[2]. As δ tends to zero the phase-space volume is preserved while the largest Lyapunov exponent stays positive. “Hamiltonian-like” is used for that system only in this restricted sense: it means volume preservation, and no Hamiltonian structure is claimed.

3. Threshold structure

Three thresholds organise the reference system, in the order γ_H < γ_cr < γ_c. Two are known in closed form.

3.1. Hopf threshold γ_H

Thet rather a gradient that can be systematically regulated. Through a blend of successive theoretical refinements and empirical validation, particularly in the realm of plasma physics, the work equilibria C± are unstable saddle-foci at γ = 0 and become stable through a Hopf bifurcation at 

This follows from the Routh–Hurwitz condition on the Jacobian at C± ^[8]. At γ_H the leading eigenvalue pair is ±10.18i. Above γ_H the real part is negative: −0.027 at γ = -0.00176, −0.57 at γ = 0.0068 and −1.23 at γ = 0.01.

3.2. Boundary crisis γ_cr

Stastablishes that chaotic behavior can ble C± do not end the chaos at once. Between γ_H and about 0.0016 the chaotic attractor coexists with the two stable equilibriacontrolled via specific parameters, external forcing, spatial context, and the largest Lyapunov exponent λ₁ stays near 0.85 (Figure 2). The attractor is then destroyed in a boundary crisieven measurement perturbations. In addition, this innovative perspective reconciles ^{[9] [10]}. Beyond it, chaos survives only as a transient whose lifetime shortens as γ grows. In an ensemble of eight initial conditions, all remained chaotic after 6000 time units at γ = 0.00164, three at γ = 0.00184, and none at γ ≥ 0.00196. A gated analysis places the canonical quench at γ_cr ≈ 0ssipative and Hamiltonian views of chaotic systems, offering a new perspective through which to interpret and control complex dynamical behaviors.00198 ^[9].

History and Theoretical Development

Past tThe crisis, λ₁ equals the real parthistorical progression of the leading eigenvalue at C±. A Lyapunov estimate taken without discarding the transient, or over a short window, reports a positive value well past γ_crHala attractor framework is marked by a series of logical and cumulative theoretical advances. The estimate mustinitial study ^[1] intherefore state its discard time ^[11].

Figure 2. Largest Lyapunov exponent λ₁ against the quench factor γ for the reference system, after discarding a transient of 6000 time units and averaging over 600. Dashed line: real part of the leading eigenvalue at C±. Dotted line: Hopf threshold γ_H ≈ 0.001393. Shaded band: boundary-crisis window, 0.00165 to 0.00198, in which chaos persists only as a long transient. Right panel: detail of the collapse.

3.3. Fold γ_c

Toduced the equiHalibria C± recede to infinity and cease to exist at

The value is independent of σ and β (Figure 3) ^[7]. The same construction gives γ_c = 1/(2z*) for other Lorenz-family hosts: 1/42 for the Chen system and 1/40 for the Lü system at their standard parameters.

Figure 3. Continuation of the equilibria C± in the quench factor. The coordinate z* = 27 does not change. The coordinates x* = y* diverge as γ approaches the fold γ_c = 1/54. Dotted lines mark γ_H and γ_c.

3.4. Over-driving

Fa attractor by modifying a classic Lorenz-type or γ > γ_c dino attracting equilibrium remains and the trajectory is unboundedary differential equation system. The escape is in x and y, not in z. The operator pins z to a small constant, z∞ ≈ 2.77 at γ = 0.14, so the term x(ρ − z) no longer saturates and x and y grow exponentially at a ratekey innovation was the incorporation of a nonlinear feedback term, quantified by a parameter of$\gamma$. aboNut 11 per time unit (Figure 4). The calibration rule for any pulse is thereforemerical simulations and analysis of the γ_on < γ_c.

Filargest Lyapure 4. Over-driven case, γ = 0.14, well above γ_c. Left: |x| and |y| grow exponentially, reaching 10⁶ by t ≈ 1.5. Right: z settles on the constant z∞ ≈ 2.77.

4. The Hala operator on other hosts

Innov exponent January 2026 the feedback term wevealed that as s$\gamma$ increparated from the Lorenz ses, the system and applied to the Sprott A flow ^[12]. For a undergoes a conthree-dimensional autonomous host with state (x₁, x₂, x₃), the operator on axis k is

Here “host” means the mathematical system that carries the operator; it does not mean a physical device. The operator has been surveyed on eleven benchmark hosts ^{[13] [14]} : Lorenz, Chen  ^[15], Lü ^[16], Dadras, Thomas, Lorenz-84, Sprott A ^[17] (the Nosé–Hoover oscillator ^{[18] [19]}), Aizawa, Halvorsen, Rössler and Rabinovich–Fabrikant.

The ounuous transitcome is specific to the host (Figures 5 and 6, Table 1). With the opeon—through a series of bifurcations—frator on the z-equation, four patterns appear.

• Fixed point, then loss of boundedness. Lorenz, Chen, Lü and Dadras reach a stable equilibrium in a finite window of γ and become unbounded above it. For the Lorenz family the window closes at γ_c = 1/(2z*).
• Fixed point at large γ only. Thomas and Lorenz-84 pass through periodic states and reach a fixed point only for γ of order 0.2 to 0.6. They stayed bounded up to γ = 1000.
• Limit cycle. Sprott A and Aizawa settle on periodic orbits and never reach a fixed point. Sprott A has no equilibria, so none can be stabilised. Halvorsen becomes periodic and then unbounded.
• No bounded quench. Rössler shows narrow periodic windows before it becomes unbounded. Rabinovich–Fabrikant stays chaotic until it becomes unbounded.

A m chaotic stranegative largest Lyapunov exponent is therefore necessary for quenching but does not guarantee a e attractor behavior to stable fixed point. The operator's axis matters as well: the z-axis is a poor match for hosts whose instability is carried by another variable, and axis choice for Rössler, Rabinovich–Fabrikant and Halvorsen remains open.

Figure 5. The Hala operator on eleven benchmark hosts. Grey: chaotic attractor at γ = 0. Blue: periodic orbit under the operator. Orange dot: fixed point under the operator. The quench factor used is given above each panel. The Aizawa cycle lies at constant z and is seen edge-on. Rabinovich–Fabrikant has no bounded quenched state with the operator on z.

Figure 6. Long-time outcome on each host against the quench factor, operator on the z-equation, from 31 values of γ between 10⁻⁴ and 10³. Grey: chaos. Blue: periodic or quasi-periodic. Green: fixed point. Orange: unbounded.

Table 1. Outcos. This finding demonstrated that the chaotic regime could be collapsed with only small changes in a single parame of the Hala operator on the z-equation of each host. Ranges are read from the grid of Figure 6 and are accurate to about a quarter of a decade.

Periodic windows inside the chaotic range (Dadras, Thomas, Rössler) are listed where the grid resolved them.

5. Control protocols

Threeer, proving that chaoticity is not an intractable attribute but a gradient protocols use the operator. ty. 

They differ in whatsecond work ^[2] hexpappens to the informnded this foundation carried by the trajectory: collapse discards it, deferral postpones the lossby adding a crucial dissipation parameter, δ, and release returns the system to its attractor without restoring the earlier trajectory.

5.1. Irreversible: Successive Controlled Collapse (SCC)

SCC san external periodic forcing term witches the operator on in pulses ^{[20] [21] [22]}. E adjustach publse raises γ smoothly to an ame amplitude γ_on < γ_c, holds it, and returns it to zero. Pulses are separated by a free-evolution interval T_free. During a pulsand frequency. This theoretical expansion was designed to explore the trajectory contracts towards C±(γ_on); the state at the end of the pulse is the landing point.

Tintricate interplay among dissipation, resonance, and chaos. By manipulating the collapse is irreversible in a practical sense. A smooth flow is invertible at every finite time, so the map is non-injective only in the limit of long pulses or below numerical resolution. On the reference system, 50 initial conditions spread over ten units per axis converge to within 2.5 × 10⁻⁹ of C± under γ = 0.005, and backward integration fromparameters, the system could be made to interpolate between strongly dissipative behavior, where phase space volumes contract significantly, and Hamiltonian-like (nearly volume-preserving) behavior. The analysis of the collLyapunov spectra, phapsed state diverges.

Two space trajectories, and fmeatures set the landing-point distribution (Figure 7).

• Pulse length. A pulse long enough for full collapse erases the history, so every cycle then starts from the same point and the protocol becomes exactly periodic. A spread of landing points requires incomplete collapse or non-identical timing.
• Free interval. With T_free = 5 the released trajectory has not left the neighbourhood of the equilibrium when the next pulse arrives, and every landing falls on the same point. With T_free = 50 the trajectory returns to the full attractor, and landings divide between C+ and C− (102 and 98 of 200).

Assures of volume contraction showed a γ_on remapproaches γ_c the equilibria recede and the plane z = 27 becomes a slow region, because ẏ = 0 wherever y = x on it. Trajectories then linger without landing (a ghost bottleneck), which is a second reason to keep γ_on well below γ_c.

Figure 7. SCC landing points on the reference system for two free intervals, with γ_on = 0.01, pulse length 3 and sigmoid ramp width 0.25, over 200 pulses. Grey dots: state at pulse onset. Orange circles: landing point at pulse end. Crosses: C±(γ_on). Light grey: the Lorenz attractor.

The pulse ramp can be generated autonomously by making γ a fourth state variable that obeys a logistic equation, γ̇ = kγ(1 − γ/γ_target), decoupled from (x, y, z). The local growth rate lags the threshold: with k = 0.3 and γ_target = 0.01, γ crosses γ_H at t ≈ 25.6 but the growth rate stays negative only from t ≈ 32, when γ ≈ 0.006 (Figure 8).

Figure 8. Autonomous logistic ramp of the quench factor on the Lorenz host: k = 0.3, γ_target = 0.01, seed γ(0) = 10⁻⁶. Top: γ(t), with γ_H marked. Middle: z(t). Bottom: local growth rate of a tangent vector over windows of 1.6 time units.

5.2. Regime reversal: collapse and release

Withdrrkable result: as dissipation was reduced toward zero, the system's phase space contraction vanished, yet the largest Lyapunov exponent remawing the operator after a collapse lets the trajectory leave the equilibrium and return to theed positive. This demonstrated the persistence of chaotic attractbehavior ^[23].even Oin the LoHamiltonian limit, where enz host the recovered attractor matches the original statistically:ergy is conserved. External forcing, meanwhile, yielded distinct resonant phenomean z of 23.4 against 23.6, and standard deviation of x of 7.85 against 7.91 (Figure 9).

Ta, which are vital for modeling pheysical reversal applies to the regime. The trajectory after release is not tprocesses such as plasma heating, where energy exchange via resonance is paramount.

The thirajectory before collapse, because all initial conditionsd and arguably most significant study ^[3] that landied on the same equilibrium share one future. The operator also removes the time-reversal symmetry of the Nosé–Hoover host: under (x, y, z, t) → (x, −y, −z, −t) the host is invariant and the operator leaves a residual 2γz(x² + y²).

On Sp these theoretical advances to real-world phenomena through empirical validation. Dr. Hala examined experimental current-voltage (I-V) traces obtained frottm A,Langmuir probes plain release does not restore chaos. The operator drives the trajectory onto a limit cycle that lies in a regular region of the conservative host. After release the trajectory stays on a regular orbit there, withquiescent plasma, specifically in the electron saturation region. Through quadratic fitting of this region, he extracted a discrete recurrence relation whose λ₁bifurcation diagram ≈displayed 0.0004 againsta classic 0.018period-doubling cascade for culminathe originaling in a chaotic orbit (Figure 10). Recovery on this host needs an additional expansion step, as in the Adaptive Reversible Expansion and Harmonic Anchoring scheme, which reports recovery with moderate phase slipregime. He also observed a spatial dependence, with chaos manifesting more strongly at the magnetically ^[23].

Figure 9. Collapse and release on the Lorenz host. Top: (x, z) projections during chaos (γ = 0, t < 40), collapse (γ = 0.01, 40 ≤ t < 80) and after release (γ = 0, t ≥ 80). Bottom: x(t), with the collapse interval shaded.

Figure 10. Collapse and release on the Sprott A host from the initial condition (0, 5, 0). Top: (x, y) projections during chaos (γ = 0, t < 300), on the limit cycle (γ = 0.05, 300 ≤ t < 700) and after release (γ = 0, t ≥ 700). Bottom: z(t), with the operator interval shaded. The released trajectory stays on a regular orbit.

5.3. Deferral: topological braiding

Deconfinerral adds a periodic modulation that keeps restoring local ind plasma boundary and relative order or stability as the operator's dissipation grows, so that collapse occurs at a larger γ ^[24] in the quiescent bulk. The figFures here use one explicit realisation, a resonant modulation of the drive parameter:

The frequency ω_b is the rotation frequency about C± at the Hopf threshold. Modulation at that frequency pumps the weakly damped spiral and ejects the trajectory from the neighbourhood of the equilibrium. The last γ with sustained chaos rises from 0.00165 at ε = 0 to 0.00235, 0.00295 and 0.0034 at ε = 0.02, 0.04 and 0.06, an extension of 42%, 79% and 106% (Figure 11).

Nethermore, temporal transitions were observed when the diagnostic probe ar cthe deferred collapse the chaos is intermittent. At ε = 0.04 and γ = 0.0029 the trajectory stays on one wing for up to 84 time units between bursts, against 20 for unmodulated chaos at γ = 0.0015. Past the collapse it locks onto a forced oscillation about one equilibrium (Figures 12 and 13).

Oively perturbed the system, acting as a form of measurement intrusion. The author then showed that a suitably parameterized Hala attractor model could reproduce all these behaviors, including Sprott A the unmodulated operator gives athe fixed point → limit cycle of→ period 5.58 (Figure 14). Modulating the constant drive, 1 → 1 + ε sin(ω_b t) with ε = 0.3 and ω_b = 0.41, turn-doubling → chaos route, as well as thhysteresis e cycle into a two-ffrequency torus with λ₁ ≈ 0 (Figure 15). Its Poincaré section is a closed curve traversed once in about 5.2 crossings (Figure 16). Tects in the I-V curves. This empirical match provided compelling evidence for the modulation here restores quasi-periodicity; it does not restore chaos.

“Braiel's utility and the theoretical framework's ability to interpret anding” is used in a descriptive sense for these orbits. The braid types of periodic orbits of the Lorenz equations wpredict real-world plasma phenomena.

Applications and Broader Implications

There classified by Birman and Williams ^[25]; braid ollective wordsk for the modulated orbits have not yet been computed.

Figure 11. Largest Lyapunov exponent against the quench factor for four modulation depths ε on the Lorenz host (ω_b = 10.18; transient of 4000 discarded). The collapse moves to larger γ as ε increases.

Figure 12. Time series x(t) on the Lorenz host. Top: unmodulated chaos, ε = 0, γ = 0.0015. Middle: intermittent chaos near the deferred collapse, ε = 0.04, γ = 0.0029. Bottom: past the collapse, ε = 0.04, γ = 0.0031, a forced oscillation on one wing.

Figure 13. Modulated Lorenz host at ε = 0.06 under increasing dissipation: (x, z) projections at γ = 0.0015, 0.0025 and 0.0033, with λ₁ = 0.86, 0.83 and 0.52.

Figure 14. Sprott A under the operator alone, γ = 0.05. Left: transient in the (x, y) plane. Centre: the limit cycle. Right: z(t).

  Figure 15. Sprott A under the operator and drive modulation, γ = 0.05, ε = 0.3, ω_b = 0.41: (x, y) and (y, z) projections of the torus and a segment of z(t).

Figure 16. Poincaré sections at x = 0 with ẋ > 0 for Sprott A at γ = 0.05. Left: without modulation the limit cycle gives one point. Right: with modulation the 1878 crossings fill a closed curve.

6. Extensions

Fractional order. Replon the Hala attractor opens up an importacing the time derivatives by Caputo derivatives of order q gives the fractional-order Hala system. Its equilibria are those of the integer-ordert implications and applications. The central tenet—that one can engineer system, so the Matignon stability boundary q_c(γ) is known in closed form from the Hopf thresholds: (293 − √65761)/26244 for the Lorenz host and (181 − √631)/7938 for the Chen host. Chaos has been maps whose chaoticity is under control—has relevance in numerous domains beyond plasma physics.

• Plasma Physics: The ability to control chaotic behavior is fundamental to understanding and manipulating wave-particle interactions, plasma heating mechanisms, and instabilities. It also holds promise for mitigating boundary phenomena in fusion devices, a critical challenge in the quest for clean energy.

• Secure Communications: Chaos can be harnessed for encryption or masking signals due to its inherent unpredictability and sensitivity to initial conditions. Chaotic signals can be used as carriers for secure data transmission, a field known as chaos-based communication. Related literature supports this potential, with optomechanical systems showing that modulated coupling can enable a switch between chaotic and regular behavior, suggesting applications in low-power optical secure communication.

• Random Number Generation: Micromechanical and nanomechanical resonators under nonlinear driving can be pushed into chaotic oscillation regimes through amplitude or frequency modulation. This chaotic output can be converted into high-quality True Random Number Generators (TRNG), which are essential for cryptography and computational simulations.

• Mixing and Fluid Dynamics: In industrial and chemical engineering, controlled chaos can be used to optimize mixing processes, ensuring uniform distribution of substances and enhancing reaction efficiency.

• Biological and Biomedical Systems: The principles of tunable chaos could be applied to models of biological systems, from heart rhythms to neural networks, to better understand and potentially control pathological states that involve erratic or chaotic dynamics.

Limitations and Future Directions

Despited over the (γ, q) plane with a full-memory predictor–corrector scheme these important contributions, ^[8].

- Time-dependent quench factor. Tthe quench factbor may follow a waveform γ(t): linear, exponential, sigmoidal or oscillating. The landing points of SCC are reported to be insensitive to the waveform as long as the plateau stays below γ_c ^[21].

- Multi-axis operator. Thdy of work acknowledges several limitations and outlines a clear path for future research. Most analyses to date have focused primarily on the largest Lyapunov exponent, operator can actrather than the full Lyapunov spectra. A mon all three equations with separate factors (γx, γy, γz). A response-surface study on the Lorenz and Sprott A hosts found the z-direction dominantre complete understanding of the system's dynamics requires mapping all exponents and relating them to the attractor's dimensionality. ^[24].

- OThe detailed fractheral lines.dimension Lat(er work carries the operator into a continuum analogue (a state-dependent bulk viscosity), a Nosé–Hoover–Hala host for a Lawson stability manifold, and a log-frequency reformulation of wave dispersion. These li.g., Kaplan-Yorke dimension, Hausdorff dimension) of the strange attractors under varying parameters has not yet been fully mapped, which is essential for characterizing their complexity.

More outside the scope of this entry.

7. Context and applications

7.1. Relation to established chaos control

Thver, the current models are OGY low-dimethod ^[26]nsional and idelayed feedbackalized. Real ^[27] plastmabilise unstable periodic orbits that already exist in the uncontrolled flow, using small perturbations that vanish on the target orbit. The Hala operator changes the vector field: it creates stability at C± and moves them. Its control effort does not vanish at the target.

No systems are spatially extended, noisy, heterogeneous, and subject to measurement back-action and perturbations. A deeper study is nlinear state feedback that drives the Lorenz system to fixed points or limit cycles is an established class ^{[28] [29] [30]}. Theded to assess the robustness of chaos tunability under these sprecifalistic contribution here is the cubic form −γ(x² + y²)z, which leaves the origin's linearisationditions. The precise thresholds for bifurcations in terms of $\gamma$, δ, unchanged, makes the contraction state-dependent, and yields closed-form thresholds.

7.2. Plasma physics and gaseous electronics

Chaos cod external forcing parameters are currentrol in discharges has been demonstrated experimentally without a model, by small modulation of the discharge voltage in a thermionic argon discharge ^[31] ay determined numerically, and rigorous theoretical or and by feedback and synchronisation in plasma diodes, ionisation waves and drift waves ^[32]. A hlytical proofs are lacking. Frost model carrying the operator adds predicted thresholds, end states and margins before hardware is built.

T the empirical side, more extensive measurements of realise the operator in a device, z must be a slow, controllablespatial and temporal variable, such as discharge current or anode biation, higher-resolution diagnostics, and x² + y² a measured fluctuation power. As an external loop this is a feedback multiplier. An internal realisation requires a physical loss channel with the same algebraic form, and that channel has to be identified for each devicecontrolled experiments that systematically vary the parameters corresponding to the model would greatly strengthen validation.

The Hala_OperatLoor_Skillking ^[33] packheages this procedure as a method with no built-in numbers. It selects hosts by structural match to the experimed, future work should aim to:

• Advanced Theoretical Analysis: Computing the full Lyapunov spectrum across the parameter space and relating it to attractor dimensionality and complexity.

• Higher-Dimensional Models: Extending the models to spatially extended, high-dimensional systems (e.g., partial differential equations) to better mimic realistic plasma behavior.

• Noise and Perturbation Studies: Systematically investigating the effects of stochastic perturbations, measurement noise, and probe-induced perturbations. This could be inspired by quantum-chaos literature, which shows that the measurement strategy itself can act as a control parameter.

• Control Strategies: Investigating robust control strategies for chaos—how to reliably suppress it or trigger it—in engineered systems.

• Thermodynamic Considerations: Exploring thermodynamic and entropy considerations, especially near the Hamiltonian limits, to understand the energy implications of controlling chaos.

• Practical Prototype Development: Pursuing the development of practical prototypes in fields like controlled plasma heating, boundary instability mitigation in fusion devices, secure communications, and random number generation to translate the theoretical findings into tangible technologies.

Int, rechecks the operator on all three axes, separates conclusions that hold on every host from host-specific ones, offers the three protocols of Section 5, and closes with a reliability audit. A pre-registered null-model study sets the acceptance rule: a result counts only if it is host-invariant, and no result on a surrogate host is claimed for a deviconclusion, the work on the Hala attractor offers a foundational, versatile, and promising framework. It provides a new perspective on chaos, not as an uncontrollable force ^[34].of In that study power-matched linear feedback outperformed the operator at 7 of 9 registered gains on the decisive host.

7.3. Other proposed uses

SCC hature, but as a controlling continuum that can be engineered and harnes been proposed for the parametric control of fission reactors ^[35], for a wide range of applications, reconciloingic primitives built from Sprott A and Rössler hosts (a clock and a switch), and for secure communication and random-number different views of chaotic systems and paving the way for a new generation. These proposals rest on numerical model of controlled dynamical systems.

8. Limitations and open questions

• Evidence is numerical. All results are from low-dimensional ordinary differential equations. No experiment has yet realised the operator in a plasma device.
• The crisis value is not closed-form. γ_H and γ_c are exact; γ_cr is numerical and depends on how long a transient is tolerated.
• Reversal is regime reversal. Collapse-and-release recovers the attractor on the Lorenz host and fails on Sprott A without an added expansion step. No protocol recovers the pre-collapse trajectory.
• Deferral is realisation-dependent. The size of the shift depends on the modulation form and frequency; at ω_b = 8 the same depths gave little or no shift.
• Braid invariants are missing. The topological description of the modulated orbits is qualitative.
• Axis selection is unresolved for Rössler, Rabinovich–Fabrikant and Halvorsen.
• Full Lyapunov spectra have not been computed for all hosts.
• A linear baseline can win. Any claim of advantage needs a power-matched linear-feedback comparison.

9. Conclusion

• Definition. The Hala attractor is the Lorenz system with one term added to its third equation: the Hala operator H = −γ(x² + y²)z. The term is a damping on z whose strength grows with the squared distance of the state from the z-axis. A single non-negative number, the quench factor γ, sets how strong it is. At γ = 0 the system is the ordinary Lorenz system; every regime described in this entry is reached by raising γ alone.

• What the operator does to the flow. Three structural effects follow from the form of the term. Because it is cubic in the state, it vanishes to first order at the origin, so the origin keeps exactly the stability it has in the Lorenz system. Because it multiplies z by x² + y², the rate at which phase-space volume shrinks is no longer constant: it is −(σ + 1 + β) − γ(x² + y²), fastest far from the z-axis. And the two nontrivial equilibria C± keep their height z* = ρ − 1 = 27 but move outward along the line y = x as γ grows, with x*² = 72/(1 − 54γ).

• Threshold sequence. Raising γ takes the reference system (σ = 10, ρ = 28, β = 8/3) through four regimes in a fixed order. Below the Hopf threshold γ_H ≈ 0.001393 the equilibria C± are unstable and the motion is chaotic. Between γ_H and roughly 0.0016 the equilibria are stable but the chaotic attractor still exists alongside them. Near γ_cr ≈ 0.002 the attractor is destroyed in a boundary crisis and every trajectory eventually settles on C+ or C−. At the fold γ_c = 1/54 ≈ 0.0185 the equilibria recede to infinity, and above it trajectories are unbounded.

• Exact and numerical results. Two of the three thresholds are exact. The Hopf threshold is γ_H = (293 − √65761)/26244, obtained from the Routh–Hurwitz condition at C±, and the fold is γ_c = α/[2(ρ − 1)] = 1/54, obtained from the condition that C± exist. The crisis value γ_cr has no closed form. Just past the crisis, chaos persists as a transient that can last thousands of time units, so the measured value depends on how long a transient the observer discards. A Lyapunov exponent computed without discarding that transient stays positive well past the true end of chaos, which is why any reported value must state its discard time.

• Calibration rule. The quench factor must always stay below γ_c. Above the fold there is no equilibrium for the trajectory to settle on. The operator then holds z at a small constant value (about 2.77 at γ = 0.14), which removes the saturation that normally limits x and y in the Lorenz system, and x and y grow exponentially at a rate of about 11 per time unit. Applying the operator too strongly therefore produces divergence, not stronger quenching.

• Host dependence. The term can be added to any three-dimensional chaotic system, called a host, and the result depends on the host. With the operator on the z-equation of eleven benchmark systems, four (Lorenz, Chen, Lü, Dadras) reach a stable fixed point over a bounded range of γ and diverge above it. Two (Thomas, Lorenz-84) reach a fixed point only at γ of order 0.2 to 0.6. Three (Sprott A, Aizawa, Halvorsen) settle on a periodic orbit and never reach a fixed point; Sprott A cannot, since it has no equilibria. Two (Rössler, Rabinovich–Fabrikant) have no bounded quenched state. For the Lorenz family the upper end of the fixed-point range is γ_c = 1/(2z*), where z* is the height of the host's equilibria. The operator is therefore a general construction but not a universal quench: a suppressed Lyapunov exponent does not guarantee a fixed point.

• Irreversible collapse. Successive Controlled Collapse applies the operator in pulses separated by intervals of free chaotic motion. Each pulse contracts the trajectory onto one of the two equilibria, and the state at the end of the pulse is a landing point. The collapse discards the trajectory's history, but only in a practical sense: a smooth flow can in principle be run backward at any finite time, so the information is lost only for long pulses or below numerical precision. One consequence is that a pulse long enough for complete collapse makes every cycle start from the same point, so the protocol repeats exactly. A varied set of landing points requires incomplete collapse or irregular timing, together with a free interval long enough for the trajectory to return to the whole attractor.

• Deferral. Deferral, called topological braiding in the source papers, adds a small periodic modulation that repeatedly pushes the trajectory away from the equilibria, so that collapse happens at a larger γ than it otherwise would. In the realisation used here the drive parameter ρ is modulated at the frequency with which trajectories spiral around C±. With modulation depths of 2%, 4% and 6% the last γ with sustained chaos rose by 42%, 79% and 106%. Just before the deferred collapse the chaos becomes intermittent, with long stretches on one wing between bursts. The effect depends on the form and frequency of the modulation; at a frequency away from the spiral frequency the same depths gave little or no shift.

• Regime reversal. If the operator is switched off after a collapse, the Lorenz host leaves the equilibrium and returns to its chaotic attractor, with the same statistics as before. What comes back is the regime, not the earlier motion: every trajectory that collapsed onto the same equilibrium follows the same path afterward, so the state before collapse cannot be recovered. The result is also host-specific. On Sprott A the operator leaves the trajectory on a regular orbit of the underlying conservative system, and switching the operator off does not return it to chaos; an additional expansion step is needed.

• Relation to prior work. Nonlinear state feedback that drives the Lorenz system to a fixed point or a limit cycle has been studied since the late 1990s, and the established OGY and delayed-feedback methods control chaos by stabilising orbits that already exist in the uncontrolled system. The Hala operator differs from the latter in that it changes the system itself and its control effort does not vanish at the target. Its specific contribution within the former class is the cubic form, which preserves the origin, gives closed-form thresholds, and supports the three protocols above.

• Status in plasma physics and gaseous electronics. The system was motivated by chaos observed in discharge plasmas, where control has been demonstrated experimentally without a model. A host model carrying the operator offers predicted thresholds, end states and safety margins before hardware is built. No experiment has yet realised the operator in a device. Doing so requires three things: a slow variable the experimenter can act on, such as discharge current or anode bias, to play the role of z; a measured fluctuation power to play the role of x² + y²; and either an external feedback multiplier or an identified physical loss channel with the same algebraic form.

• What a valid claim requires. Because outcomes differ from host to host, a result is claimable only if it holds on every host that plausibly represents the device, and a result on a surrogate model is not a result for a device. Any claim that the operator outperforms conventional control also needs a comparison against linear feedback using the same actuator power. In the one pre-registered test so far, linear feedback outperformed the operator at 7 of 9 gains on the decisive host.

• Open work. Four items remain. A device-level test against a linear baseline would establish whether the operator has practical value in a plasma. Computing braid invariants for the modulated orbits would turn the topological description from qualitative to exact. Choosing the operator's axis for the Rössler, Rabinovich–Fabrikant and Halvorsen hosts may yield bounded quenched states that the z-axis does not. Full Lyapunov spectra across all hosts would complete the dimensional picture.