1. 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 plasma motivation is older thanof work on the Hala attractor and the model. A 2020 letter read chaos from probe I–V characteristics through a discrete recurrence fitted to the electron-saturation region ^[4], ified Hala attractor and its hybrid version present a coherent and progressive research trajectory in which chaos is refranmed a 2025 preprint developed the electric probe as an instrument for that purposfrom being an uncontrollable, intrinsic feature ^[5].of A rcecurrence obtained by discretising a one-dimensional I–V curve is rtain systems to instead a property of the chosen discretisation. It is not by itself evidence of chaos in the plasma that can be regulated, tuned, and modulated via specific parameters, external forcing, spatial context, and the entry treats it as motivation only.

Smeasurement perturbations; this is accomplincshe January 2026 the feedback term has been studied as an operator in its own right, separate from the Lorenz hostd through successive theoretical refinements and empirical validation, especially in plasma physics. The body of work now has seven layers: precursors, the afirst study introduces the Hala 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

The by modifying a Lorenz-type ordinary differential equation system, incorporating a nonlinear feedback term equilibria C± are unstable saddle-foci at γ = 0 and become stableantified by a parameter γ, and demonstrates 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

Stable C± do numerical simulations anot end the chaos at once. Between γ_H and about 0.0016 the chaotic attractor coexists with the two stable equilibria, and the laalculation of the largest Lyapunov exponent λ₁ stays near 0.85 (Figure 2). The attractor is then destroyed in a boundary crisisthat as γ increases, the system undergoes ^[9]a ^[10]. Beycond tit, chaos survives only as a nuous 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 ≈ 0.00198 ^[9].

Ption—through bifurcations—from chaotic strange attractor behavior to stable fixed points; thus the chaotic regime can be collapsed by only modest changes in γ, indicating that chasot the crisis, λ₁ equals the real part 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 γ_cr. The estimate must therefore state its discard timicity is not an all-or-nothing attribute but a gradient property. The second work expands the Hala attractor framework by adding a dissipation parameter δ and an external periodic forcing term with adjustable amplitude ^[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

Theand frequilibria C± recede to infinity and ceaseency so as 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

Fplore γ > γ_c no attracting equilibrium remains and the trajectory is unbounded. 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 rate of about 11 per time unit (Figure 4). The calibration rule for any pulse is therefore γ_on < γ_c.

Figure 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

In Jthe interplay among dissipation, resonance, and chaos, thereby enabling the system to interpolate between strongly dissipative behavior and Hamiltonian-like (nearly volume preserving) behavior; by analyzing the Lyapunov spectra (at least the lanuary 2026 the feedback term was separated from the Lorenz system and applied to the Sprott A flow ^[12].gest exponent), phase space trajectories, and measures of volume contraction For a three-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.

Thdivergence, the authors show that as dissipation is re odutcome is specific to the host (Figures 5 and 6, Table 1). With the operator 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.

Aced toward zero (or to a critical threshold) the system’s phase space contraction vanishes while negativhe largest Lyapunov exponent is therefore necessary for quenching but does not guarantee a 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.remains positive, preserving chaotic behavior even in the Hamiltonian limit; external forcing yields resonant phenomena and distinct resonance sweeps, which helps in modeling Outphysicome 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

Threeal processes such as plasma heating where energy exchange via resonance protocolis use the operator. They differ in what happens to the information carried by the trajectory: collapse discards it, deferral postpones the loss, and release returns the system to its attractor without restoring the earlier trajectory.

5.1. Irreversible: Successive Controlled Collapse (SCC)

SCC switchimportant. The third study ties these theoretical advances to experimental data: the authors examine empirical current-voltage (I-V) traces obthe operator on in pulses ^{[20] [21] [22]}. Eaained by Langmuir probes in quieschent pulse raises γ smoothly to an amplitude γ_on < γ_c, holds it, and returns it to zero. Pulses are separated by a free-evolution interval T_free. During a pulse the trajectory contracts towards C±(γ_on); the state at the end of the pulse is the landing point.

The lasma, notably in the electron saturation region, and via quadratic fitting of that region extract a discrete recurrence relation whose bifurcation diagram displays a period-doubling collapse is irreversiblescade culminating 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.005chaotic regime; they observe spatial dependence, with chaos manifesting more strongly at the magnetically confined plasma boundary and relative order or stability in the quiescent bulk, and backward integration from the collapsed state diverges.

Twotemporal transitions induced when the diagnostic probe actively fpeatures 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).

Arturbs the system; the authors γ_othen approaches γ_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

Withdshow that a suitably parameterized Hala attractor model reproduces these behaviors, including the fixed point → limit cycle → period-doubling → chaos route, as well as hysteresis effects in I-V curawing thve operator after a collapse lets the trajectory leave the equilibrium and return to thes. Together, these three contributions establish that chaotic attractbehavior ^[23]. Oin systhe Lorenz host the recovered attractor matches the original statistically: mean z of 23.4 against 23.6, and standard deviation of x of 7.85 against 7.91 (Figure 9).

Thems of the sort exemplified by plasmas is not an unavoidable given, but depends in systematic ways on fee reversdbal applies to the regime. The trajectory after release is not the trajectory before collapse, because all initial conditions that landed 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 invariantck strength, dissipation magnitude, external forcing, spatial positioning, and measurement intrusion; that transitions between stability, periodicity, and chaos follow well-known dynamical systems routes (bifurcations, period-doubling cascades); and the operator leaves a residual 2γz(x² + y²).

Oat even systems approaching Hamiltonian (non S-dissiprott A, plain release does not restoreative, volume preserving) regimes can retain 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, with λ₁ ≈ 0.0004 against 0.018 for the original chaotic orbit (Figure 10). Recovery on (positive Lyapunov exponent) despite absence or near absence of dissipation. The empirical match between model predictions and Langmuir probe data underscores the theoretical framework’s utility for interpreting real plasma phenomena. Moreover, this host needs an additional expansion step, as in the Adaptive Reversible Expansion and Harmonic Anchoring schemebody of work opens up important implications and applications: one may engineer systems whose chaoticity is under control, which reports recovery with moderate phase slip ^[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

Deferral addhas relevance in plasma physics (wa periodic modulation that keeps restoring local ve-particle interactions, plasma heating, instability as the operator's dissipation grows, so that collapse occurs at a larger γ ^[24] . The fies, boundary phenomena in fusion devices), secure communications (using chaos for encryptigureson 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).

Nor masking), mixing applications, and possibly in any engineereard the 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 unsystems where unpredictability and sensitivity are both liabilities and assets. Additional examples from related literature support this potential: for instance, optomechanical systems with modulated chaos at γ = 0.0015. Past the collapse it locks onto a forced oscillation about one equilibrium (Figures 12 and 13).

Ooupling and synthetic magnetic fields have been shown to alternate between chaotic and regular behavior depending on Sprott A the ununable boundary phase modulated operator gives a limit cycle of period 5.58 (Figure 14). Modulating the constant drive, 1 → 1 + ε sin(ω_b t) with ε = 0.3 and ω_b = 0.41, turns the cycle into a two-frequency torus with λ₁ ≈ 0 (Figure 15). Its Poincaré section is a closed curve traversed once in about 5.2 crossings (Figure 16). Theion, suggesting that chaos can be used in low-power optical secure communication devices or integrated chaotic light sources; likewise, micromechanical and nano-mechanical resonators under nonlinear driving have been converted into true random number generators when pushed into chaotic oscillation regimes through amplitude or frequency modulation . There restores quasi-periodicity; it does not restore chaos.

“Brmain, however, several limitations and unresolved questions: most aiding” is used in a descriptive sense for these orbits. The braid types of periodic orbits of the Lorenz equations were classified by Birman and Williamsalyses so far focus on the largest Lyapunov exponent rather than full Lyapunov spectra (all exponents), and detailed fractal dimension (e.g., Kaplan-Yorke dimension, ^[25]; brHaiusd words for the modulated orbits havorff dimension) of the strange attractors under varying parameters are 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.fully mapped; Replacing the time derivatives by Caputo derivatives of order q gives the fractional-order Haldditionally, the models are low-dimensional and idealized, whereas real plasma system. Its equilibria are those of the integer-order system, so the Matignon sts are spatially extended, noisy, heterogeneous, and subject to measurement back-action and perturbations, meaning that robustness of tunability boundary q_c(γ) is known in closed form from the Hopfunder realistic noise and boundary conditions still needs deeper study; further, the precise thresholds: (293 − √65761)/26244 for the Lorenz host and (181 − √631)/7938 for the Chen host. Chaos has been mapped over the (γ, q) plane with a full-memory predictor–corrector scheme at which bifurcations occur (in terms of γ, δ, forcing amplitude or frequency, etc.) are primarily determined ^[8].

- Time-dependent quench factor. The qnumench factor 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. Trically, and rigorous theoretical/analytical proofs are lacking; from the empirical side, more extensive measurement of spatial and temporal variation, higher operator can act on all three equations with separate factors (γx, γy, γz). A response-surface study oresolution diagnostics, and controlled experiments varying the same parameters as in the Lorenz and Sprott A hosts found the z-direction dominantmodel would strengthen the validation. Looking ahead, ^[24].

- Other lines. Lafuturer 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 lie outside the scope of this entry.

7. Context and applications

7.1. Relation to established chaos control

Theshould aim to compute the full Lyapunov spectrum across parameter space and relate that to attractor dimensionality and complexity; analyze the fractal dimensions and metric entropy; OGY method ^[26] axtend dthelayed feedback ^[27] stabi modelise unstable periodic orbits that already exist in the uncontrolled flow, using small perturb to spatially extended, high-dimensional systems (partial differential equations 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.

Nonl, spatial coupling) to better mimic realistic plasma behavior; study the effects of stochastic perturbations, measurement noise, and probe-inducear state feedback that drives the Lorenzd perturbations more system to fixed points or limit cycles is an established class ^{[28] [29] [30]}. Thatically; investigate specific contribution here is the cubic form −γ(x² + y²)z, which leaves the origin's linearisation unchanged, makes the contraction state-dependent, and yields closed-form thresholds.

7.2. Plasma physics and gaseous electronics

Chontrol strategies for chaos (how to suppress it or trigger it reliably) in engineered systems; explore thermodynaosmic control in discharges has been demonstrated experimentally without a model, by small modulation of the discharge voltage in a thermionic argon discharge ^[31] a/ energy / entropy considerations especially near Hamiltonian limits; refine understanding of measuremendt by feedback and synchronisation in plasma diodes, ionisation waves and drift waves ^[32]. A hack-action and the limits of observability (inspired by quantum-chaost model carrying the operator adds predicted thresholds, end states and margins before hardware is built.

To literature showing that measurement strealise the operator in a device, z must be a slow, tegy can itself be a controllable variable, such as discharge current or anode bias, and x² + y² a measured fluctuation power. As an external loop this is a feedback multiplier. An internal realisl parameter affecting emergence or suppression of chaos); and finally pursue practical application requires a physical loss channel with the same algebraic form, and that channel has to be identified for each prototypes in controlled plasma heating, boundary instability mitigation in fusion device.

Thes, Hala_Opserator_Skill ^[33] pacure communickages thistions, signal procedure as a method with no built-inssing, mixing, and random numbers. It selects hosts by structural match to the experiment, 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 generation. Collectively, the work on the Hala attractor offers a paradigm shift: chaos can be conceived not as a binary defect but as a continuum, something that can be engineered, tuned, and harnessed, in theory, simulation, and no result on a surrogate host is claimed for a device ^[34]. In that experiment, thereby reconciling disstudy ipower-matched linear feedback outperformed the operator at 7 of 9 registered gains on the decisive host.

7.3. Other proposed uses

SCC hative and Hamiltonian views of chaotic systems and providing as been proposed for the parametric control of fission reactors ^[35], fpowerful new lens throrugh logic primitives built from Sprott A and Rössler hosts (a clock and a switch), and for secure communication and random-number generation. These proposals rest on numerical modelswhich to interpret and control complex dynamical behavior—foundational, versatile, illuminating, and deeply promising.

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.