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 plabody of resma motivation is older than the model. A 2020 letter readearch on the Hala attractor addresses the treatment of chaos in nonlinear dynamical systems. It departs from the conventional view of chaos from probe I–V characteristicsas an intrinsic and uncontrollable feature of certain systems, reframing it instead as a tunable property accessible through a discrete recurrence fitted to the electron-saturation regioninternal feedback. This research direction, now spanning more than six studies, indicates that chaoticity is not a binary state but a gradient ^[4],that cand a 2025 preprint developed the electric probe as an instrument for that purpose be systematically regulated. Through successive theoretical refinements and empirical comparison in plasma physics, the work shows that chaotic behaviour responds ^[5].to Aspecific recurrence obtained by discretising a one-dimensional I–V curve is a property of the chosen discrparameters, to external forcing, to spatial context, and to measurement perturbation. The framework also connects dissipative and Hamiltonian descriptions of chaotic systems, providing a route to interpret and control complex dynamical behaviour.

2. History and Theoretical Development

Theti progressation. It is not by itself evidence of chaosion of the Hala attractor framework is marked by a series of cumulative theoretical advances. The initial study ^[1] introduced the plasma, and the entry treats it as motiv Hala operator by modifying a classical Lorenz-type ordinary differential equation system. The key innovation was the incorporation only.

Sf a state-dependent nonlince January 2026 the ar feedback term has been studied as an opera, quantified by a quenching factor γ. The equation its own right, separas of the Hala attractor chaotic system are

ẋ = σ(y − x)

ẏ = x(ρ − z) − y

ż = αxy − βz + ℍ(x, y, z)

where the from the Lorenz host. The body of wHala operator acts along the z direction and is defined as

ℍ(x, y, z) = −γ(x² + y²)z

Thrork nughow has seven layers: precuut this entry article the classical parameter values σ = 10, ρ = 28, α = 1, β = 8/3 are asors, the attractor, tsumed unless stated otherwise. The operator, control protocols, a log-frequency is quadratic in the transverse coordinates and linear in z, extension, applications, and methods for auditing claimo it is negligible near the origin and dominant at large phase-space excursions. 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

Tis the property that allows it to stabilize the unstable saddle-foci of the host system without altering the local behaviour near the origin.

Numequrilibria C± are unstable scal simulation and analysis of the largest Lyapunov exponent revealed that as γ increaddlse-foci at γ = 0 and become stables the system undergoes a continuous transition, through a Hopseries of 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

Ss, from chaotic strange-attrable C± do not end tctor behaviour to stable fixed points. The chaos at once. Between γ_H and about 0.0016 thetic regime collapses under small changes in a single parameter, establishing that chaotic attraity is not an intractable attribute but a gradient property.

Figure 1. Symmetry breaking and manifold collapse. Phase-space projections of the Hala attractor at three operating points. γ = 0.0: the fully developed chaotic attractor; while superficially resembling classical dual-lobe flows, it lacks rotational symmetry and exhibits localized state-dependent phase compression. γ = 0.0001: micro-feedback initiation; the operator breaks the spatial topology, pinching trajectories along a directional vector and compressing the multi-layered manifold toward a dense nodal origin. γ = 0.01: operating within the negative Lyapunov regime (λ₁ ≈ −0.73), the operator has reduced the effective dimensionality; the chaotic attractor is suppressed and trajectories map onto a single stable focus.

Figure 2. Quantification of phase-space stabilization via the Hala operator. The largest Lyapunov exponent λ₁ as a function of the internal control feedback parameter γ. For γ ∈ [0.0, 0.0065] the system maintains a robust chaotic regime (λ₁ > 0), exhibiting a slow, controlled degradation of structural divergence. At the critical threshold γ_λ ≈ 0.0068 the exponent drops sharply and discontinuously across the zero boundary (red dashed line). For γ > γ_λ it is strictly negative, marking the loss of chaotic stretching and the convergence of trajectories onto a stable invariant set.

Three distinct threshor coexists with the two stable lds govern the controlled system and should be distinguished. The reverse-Hopf stabilization of the non-trivial equilibria,a C± occurs at γ_H ≈ 0.00176 (Tandble 1); the largest Lyapunov exponent λ₁ stays nearcrosses zero at γ_λ ≈ 0.00685 (Figure 2). The attractor is then dest; and the equilibrium branch itself terminates at the fold threshold γ_c = 1/54 ≈ 0.018519, deroyived in a boundary crisisSection 3.1. These satisfy ^[9]γ_H < γ_λ < ^[10]γ_c. BTheyond it, chaos survives only as a transient whose lifetime intermediate interval, in which the equilibria are already linearly stable while a chaotic attractor persists, is directly analogous to the coexistence window of the classical Lorenz system near ρ ≈ 24.06 – 24.74.

The shecond wortk ^[2] expans as γ grows. In an ensemble of eight initial conditions, all remaineded this foundation by adding a dissipation parameter δ and an external periodic forcing term of adjustable amplitude and frequency, in order to explore the interplay among dissipation, resonance, and chaotic after 6000 time units at γ = 0.00164, three at γ = 0.00184, and none at γ ≥ 0.00196. A gateds. The parameter δ scales the linear damping of the host, allowing the system to be tuned continuously between a strongly dissipative regime and a conservative one. Two quantities must be distinguished in this analysis places the canonical quench at γ_cr ≈ 0.00198. The sum of the Lyapunov exponents equals the time-averaged divergence of the vector field, Σλᵢ = ⟨∇·F⟩, ^[9].

Pand governst the crisis, λ₁ equal the evolution of phase-space volume; the largest exponent λ₁ measures 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 timeate of separation along the unstable direction. For the unmodified host the divergence is the constant −σ − 1 − β, fixed by the linear damping coefficients alone, whereas λ₁ is a dynamical property that those coefficients do not directly set. The two can therefore be varied independently, and the analysis of the Lyapunov spectra, phase-space trajectories, and volume contraction established that they ^[11]are.

As δ was reduced toward zero the phase-space contraction vanished, Σλᵢ → 0, while the largest Lyapunov exponent remained positive. In three dimensions, with the exponent along the flow direction vanishing, this forces λ₃ = −λ₁: stretching along the unstable direction is balanced exactly by contraction along the stable one, with no net loss of volume. Chaotic behaviour therefore persists in the conservative limit, in which Liouville phase-space volume is preserved. The term Hamiltonian-like is used here in that restricted sense — the flow becomes volume-preserving rather than derived from a Hamiltonian function, since the host admits no natural energy function or symplectic structure. The distinction matters because the two limits support chaos of qualitatively different kinds. In the dissipative regime volumes contract to zero and trajectories collapse onto a measure-zero strange attractor with a well-defined basin, so that initial conditions are forgotten. In the conservative limit no attracting set can exist, and chaos instead occupies a positive-measure region of phase space, typically interleaved with regular islands, in which initial conditions are never forgotten. The largest Lyapunov exponent remains positive across the transition while the geometric object supporting the chaos changes entirely, and it is this passage between the two descriptions that the modified formulation was constructed to bridge. The Hala operator occupies a specific place in this picture. Its contribution to the divergence is −γ(x² + y²), so that

∇·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.

= −σ − 1 − β − γ(x² + y²)

3.3. Fold γ_c

Unlike the linear damping, which contributes a constant, the operator supplies contraction that vanishes on the z-axis and grows quadratically with distance from it. The controlled flow is accordingly near-conservative close to the axis and increasingly dissipative at large excursions — neither uniformly dissipative nor uniformly volume-preserving, but a state-dependent combination of the two. A consequence follows for the framework as a whole: in the limit δ → 0 the operator becomes the only remaining source of volume contraction, so that γ alone carries the system from conservative chaos to a quenched state. The quenching mechanism described throughout this entry does not depend on background dissipation for its action. External forcing yielded distinct resonant phenomena. In a strongly damped system a periodic drive produces a bounded response set by the balance between energy input and dissipation; as δ → 0 that balance is removed, and a drive tuned near a natural frequency of the flow accumulates energy limited only by the nonlinearity. This regime is directly relevant to radio-frequency plasma heating, where a nearly collisionless plasma is close to conservative on the heating timescale and energy is deposited through wave–particle resonance rather than through collisional dissipation.

The equilthibriard study C±^[3] rconnecede to infinity and ceaseted these advances 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

Foperiment. Curr γ > γ_c eno 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 tht–voltage traces obtained from Langmuir probes in quiescent plasma were analysed in the electron saturation region. Quadratic fitting of that region yielded a discrete recurrence relation whose bifurcation diagram displayed a classic period-doubling cascade 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

Iinating in a chaotic regime. A spatial dependence was observed, with chaos manifesting more strongly at the magnetically confined plasma boundary and greater regularity in the quiescent bulk, and temporal transitions appeared when the diagnostic probe actively perturbed the system, a form of measurement January 2026 the feedback term was separated from the Lorenz system and applied to the Sprott A flowintrusion. A suitably parameterized Hala attractor model reproduced these behaviours, including the fixed point → limit cycle → period-doubling → chaos ^[12]. Frorute a three-dimensional autonomous host with state (x₁, x₂, x₃), the operator on axis k ind the hysteresis observed in the plasma I–V characteristics curves, providing evidence for the model’s ability to interpret and predict real physical plasma phenomena.

3. Deterministic Control and Discretization via Successive Controlled Collapse

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 ouHala attcome is specific to the host (Figures 5 and 6, Table 1). With the operator on the z-equation, four patterns appearactor admits a deterministic control protocol termed Successive Controlled Collapse (SCC). Brief activation of the feedback factor γ > 0 through a smooth sigmoidal ramp of amplitude γ_on for.

• 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.

Aces negchative largest Lyapunov exponent is therefore necessary for qotic trajectories to collapse onto the stable equilibria of the controlled flow in a globally contracting manner. Repeating this pulse, separated by a free-evolution interval T_free duerinching but does not guarantee a fixedg which γ returns to zero, generates families of discrete landing point. The operator's axis matters as well: the z-axis is a poor match for hosts whose instability iss whose distribution traces the geometry of the underlying strange attractor. Because the collapse phase is globally contracting, distinct pre-collapse states are carried by another variable, and axis choice for Rössler, Rabinovich–Fabrikant and Halvorsen remains openonto common landing points, and the pre-image of a landing point cannot be recovered from the point itself. SCC is therefore an irreversible protocol, in contrast with the reversible and partially reversible mechanisms of Section 5.

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. A uniqueness Oargutcome 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 decadement establishes that no two landing points coincide under non-identical pulse timing sequences, so that the protocol replaces the continuous chaotic flow with a programmable discrete set that remains faithful to the geometry of the original attractor.

3.1 Equilibrium structure and the calibration threshold

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

5. Control protocols

Three protocols use the operator. Th equilibria onto which the collapse occurs follow in closed form. Setting ẋ = 0 giveys y = x, andiffer in substituting into ẏ = 0 yields wx(ρ − z − 1) = 0, so that eithaer x = 0, corresppeondinsg to the inforigin, orma

z* = ρ − 1 = 27

The vertioncal carried by the trajectory: collapse doordinate of the non-trivial equilibria is therefore independent of γ. Substitutiscng y = x arnds z = z* into ż = 0 gives x*²(α − 2γz*) = βz*, dhefnce

x*² = βz* / (α − 2γz*) = 72 / (1 − 54γ)

At γ = 0 this returns x* = ±6√2 ≈ ±8.485, the classical postpones the loss, and release returns the system to Lorenz equilibria C±. As γ increases the branch grows monotonically and diverges when the denominator vanishes, giving the closed-form fold threshold

γ_c = α / [2(ρ − 1)] = 1/54 ≈ 0.018519

independents of σ attnd β. Beyond γ_c no reactor without restoring the earlierl non-trivial equilibrium exists within the symmetric ansatz, and trajectory.

5.1. Irreversible: Successive Controlled Collapse (SCC)

SCC switches the operator on in pulsesno longer relax to static equilibria but are ^[20]governed ^[21]by ^[22].transient Edynach pulse raises γ smoothly to an amics. This value bounds the admissible pulse amplitude γ_on < γ_c, holds it,for the protocol, as developed in Section 3.3.

Figure 3. Equilibrium continuation of the Hala–Lorenz system as a function of the control parameter γ. The fixed-point branches (x*, y*, z*) preserve the structural symmetry x* = y*, with z* = 27 remaining invariant across the sweep. As the control parameter approaches the fold threshold γ_c = 1/54 ≈ 0.0185 (red dotted line), x* and y* asymptote toward infinity, marking the critical limit past which real non-trivial equilibrium branches cease to exist within the symmetric ansatz.

3.2 Free-evolution interval and the landing-point distribution

The acontind returns it to zero. Pulses are separuous trajectory evolution during the uncontrolled phase is evaluated by a varying the free-evolution interval T_free. During a pulse the trajectory, which reveals the intrinsic dissipative properties of the flow under periodic switching:

• Short evolution interval (T_free = 5.0). Trajectories retain residual displacement from the preceding parameter pulse, resulting in a broader spatial dispersion of landing points that reflects transient manifold structure.
• Extended evolution interval (T_free = 50.0). Trajectories undergo prolonged phase-space volume contraction, drawing landing points into significantly tighter clusters along the local manifold.

This spatial condensatracts towards C±(γ_on); the state at the end of the pulse is the landing point.

Tion at larger T_free represents physical relaxation along the stable manifold, confirming that longer evolution times allow the controllapse is irreversible in a practical sened system to settle closer to its asymptotic invariant sets.

Figure 4. Phase-space volume contraction under the Successive Controlled Collapse protocol. Side-by-side comparison of landing-point constellations evaluated at short (T_free = 5.0, left) and extended (T_free = 50.0, right) free-evolution intervals, alongside the target region (cyan marker). Shorter evolution intervals capture transient dispersion caused by residual displacement after the pulse, whereas prolonged free evolution permits phase-space volume contraction, drawing trajectories into tighter clusters along the invariant manifold.

3.3 Pulse calibration

The pulse. 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 integramplitude must satisfy γ_on < γ_c. When this condition is violated the protocol has no target: at γ_on = 0.14, far beyond γ_c = 1/54, the equilibrium relation of Section 3.1 gives x*² = 72/(1 − 7.56) < 0, so no real equilibrium exists for the trajectory to collapse onto. The resulting unbounded vertical excursion, in which the coordinate z exceeds 10⁵, is the correct behaviour of the vector field under an over-driven pulse rather than an artifact of the integrator; imposing state clipping or domain bounding would conceal a genuine dynamical escape rather than correct a numerical fault. The operational requirement for maintaining bounded discretization from the collapsed state diverges.

Twacross successive switching cycles is accordingly a calibration constraint on the pulse amplitude, γ_on features set the l< γ_c, and not a solver setting.

Figure 5. Trajectory escape under an over-driven control pulse, γ_on = 0.14 ≫ γc, across T_free = 5.0 (left) and T_free = 50.0 (right). With the pulse amplitude driven far beyond the fold boundary, no real equilibrium exists within the symmetric ansatz, and the nonlinear feedback term −γ(x² + y²)z drives trajectories out of the bounded phase space along vertical manifold towers (z > 10⁵), demonstrating the necessity of pulse calibration.

3.4 Transient trapping and the ghost bottleneck

Landing- point distribs recorded at short free-evolution (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).

Aintervals form a tight cluster near z = 27, well inside γ_on approaches γ_cthe region vacated by the equilibria recede and the plane z = 27 becomes a slowum branch. This clustering is a ghost bottleneck: a localized region, because ẏ = 0 wherever y = x on it. T of phase space in which the flow retains the slow character of a nearby equilibrium structure that the trajectories then linger without landing (a ghost bottleneck), which is a second reason toy can no longer reach, so that passing trajectories stall as though an attractor were still present. The mechanism has two components in the Hala–Lorenz system.

First, ktheep plane z = ρ − 1 = 27 carries γ_onthe 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

Wiequilibrium branch for every value of γ, and on thdrat plawing the operator after a collne the second component of the vector field vanishes identically, ẏ = x(ρ − z) − y = 0 whenever y = x. The plapne ise lets the trajectory leav therefore a slow surface of the flow irrespective of where the equilibrium and return to the chaotic attractora themselves lie. Second, as γ approaches γ_c from below the equilibria recede toward large ^[23].|x| Oalon the Lorenz host the recovered attractor matches the original statisg that plane, and past γ_c they cease to exist altogether; what remains at moderate |x| is a region of near-vanishing phase-space velocity with no equilibrium to terminate it. Trajectories adiabatically: mean z of 23.4 against 23.6, and sta tracking a continuous sigmoidal ramp enter this region and traverse it slowly.

Under ard deviation of x short free-evolution interval (T_free of= 7.85 against 7.91 (Figure 9).

T5.0) integration ends while reversal applies to the regime. The trajectory after release is not the trajectory before collapse,trajectories are still crossing the bottleneck, and the recorded landing points accumulate there. The structure is termed a ghost because all initial conditions that landed on the same equilibrium share one future. The operator also removes the time-rit haunts, figuratively speaking, the coordinates of a state the system can no longer occupy: it produces the visual and structural illusion of a static attractor over short observation windows, and it disperses when trajectories are given sufficient time (T_free = 50.0) to eversal symmetry of the Nosé–Hoover host: under (x, y, z, t) → (x, −y, −z, −t) the host is cape the bottleneck and relax onto true asymptotic invariant sets. The contrast between the two panels of Figure 4 is the direct signature of this effect.

4. Universality of the Hala Operator

The uninvariant and the ersality of the Hala operator leaves a residual 2γz(x² + y²).

Ois established by direct computation across a benchmark set Sprott A, plain release does not restore chaos. Trather than by appeal to a single host system. Applied in its canonical z-directed form to eleven canonical chaotic systems, the operator drives the trajectory onto a limit cyclelargest Lyapunov exponent negative in ten. The long-horizon behaviour behind that lies in a regular region of the conservative host. After release the trajectory stays on a regular orbit there, with λ₁ ≈ 0.0004count is more differentiated, and is summarised in Table 1. Clean, sustained convergence to a fixed point is obtained in four hosts — Lorenz, Chen, Lü, and Dadras. Two further hosts reach zero-dimensional states only at large multiples of their thresholds: Thomas at approximately 2.8 × 10⁴ γ_c and Lorenz-83 agt ainst 0pproximately 8 γ_c.018 The remainder terminate for the original chaotic orbit (Figure 10). Recovery on this host needs an additional eon sustained one-dimensional limit cycles (Sprott A, Aizawa, Halvorsen) or diverge (Rössler, Rabinovich–Fabrikant). A negative largest Lyapunov exponent is therefore necessary but not sufficient for quenching, and the two criteria are reported separately here.

Three hosts exphibit ansion step, as in the Adaptive Reversible Expansion and Harmonic Anchoring scheme, which reports recoveryxis mismatch. For Rössler the appropriate damping direction is x rather than z; for Rabinovich–Fabrikant it is y, since z-damping decouples from the term bounding y and produces catastrophic growth; and for Halvorsen the correct axis remains unresolved, with moderate phase slipz evidently unsuitable. This dependence ^[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

Dgefnerral adds a periodic modulation that keeps restoring local instability as tralises into a concrete per-host axis-selection procedure applicable to any future host, and is the precise sense in which the operator is universal.

Table 1. Divergence structure, equilibrium persistence, and verified long-horizon behaviour under the operator across all eleven canonical hosts: not merely whether λ₁ crosses zero, but what the trajectory actually does once it has.

The operator'’s dissipation grows, so that collapse occurs at a larger γ ^[24] sigmoidal control schedule, previously imposed as an external timetable, is the exact solution of a logistic equation. The fconstructigures here use one explicit realisation, a resonant modulation of the drive on can therefore be recast as a single autonomous four-dimensional system in which γ is a state variable rather than a 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).

N ^[4]. Framed this war the deferred collapse the chaos is intermittent. At ε = 0.04 and γ = 0.0029 the y, the Hala attractor is no longer an operator applied to a host system with an externally chosen γ; it is a self-contained four-dimensional dynamical system whose own trajectory stays on one wcarries it from strange-attractor dynamics to a quenched register.

Figure 6. Quenching fbehavior up to 84 time units between bursts, ur of the Hala operator across benchmark systems.

Becagainust 20 for unmodulated chaos at γ = 0.0015. Past e γ evolves autonomously and is decoupled from (x, y, z), a single realisation of the collapse it locks onto a forced oscillation about one equilibrium (Figures 12 and 13).

On Sproupled flow sweeps through the chaotic-to-quenched transition along its own trajectory. Whether that transition is detectable from within the dynamics was tested using a sliding-window, multi-direction local Lyapunov estimat A the unmodulated operator gives a limit cycle of period 5.58 (Figure 14). Modulating the constant drive, 1 → 1 + ε sin(ω_b t) withe computed continuously as the trajectory evolves: three independent perturbation directions, renormalised every 1.6 time units, on the Lorenz host with γ _base ε= 0, γ _target = 0.301, and ω_b = 0.41, turns the cycle into a two-frequency torus with λ₁ ≈ 0 (Figure 15). Its Poincaré sek = 0.3. The estimate does detect the transition, settling reliably negative once γ(t) is past threshold. As Figure 7 shows, however, it does so only well after the crossing, and this detection is a closed curvelag is a quantified limitation of the self-regulation claim.

Figure 7. Top: the autonomous logistic state γ(t) compared with the externally determined static threshold γ_H ≈ 0.00176. Bottom: the online local Lyapunov estimate (raw, three-direction average, and smoothed), which settles reliably negative only well after γ(t) has crossed the static threshold.

This is travhersed once in about 5.2 crossings (Figure 16). The mod sense in which universality and self-regulation, used loosely in the original formulation here restores quasi-periodicity; it does not restore chaos.

“B, can be made precise: universality is a per-host axis-diagnosis procedure rather than a fixed rule, and self-regulaiding” is used ition is an autonomous property of a single augmented flow rather than a descriptive sense foron of externally swept behaviour.

5. Reversible and Partially Reversible Chaos Control

Whereas SCC is thirrevese orbits. The braid types of periodic orbitrsible, two further mechanisms provide reversible and partially reversible control. They act at different stages of the Lorenz equations were classified by Birman and Williamscollapse and should be distinguished: the first recovers an attractor after collapse ^[25]; brhaids words for the modulated orbits have not yet been computedoccurred, the second defers collapse before it occurs.

Figur5.1 Recove 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.

Figurey: 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.

Figurve 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).

Fsigurble 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

Fracontional ordrol by collapse and releaser.

The first Rmeplacing the time derivatives by Caputo derivatives of order q gives the fractional-order Hala chanism contracts the chaotic trajectory onto a lower-dimensional invariant set — a stable fixed point in the Lorenz system. Its equilibria are those of the integer-order system,, a stable limit cycle in the Sprott A system — and then withdraws the control, reintroducing phase-space expansion so the Matignon stability boundaryat the trajectory evolves back toward its original attractor q_c(γ)^[5]. Fis known in closed form from the Hopf thresholds: (293 − √65761)/26244 for thgures 8 and 9 show that the reconstructed Lorenz attractor closely matches the original butterfly geometry, whereas the recovered Sprott A attractor preserves its characteristic multi-loop topology but exhibits phase slip and incomplete temporal reconstruction.

Geometric 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 ^[8]recovery alone is a weak criterion. The attractor is an invariant set determined by the vector field, so once the control is withdrawn any initial condition in the basin regenerates the same geometry; the discriminating quantity is the trajectory’s phase on the attractor, which is not recovered in the Sprott A case.

- Time-dependent quench factor. The qudistinction betweench factor may follow a waveform γ(t): linear, exponential, sigmoidal or oscillating the two hosts is structural. The Lorenz system possesses equilibria that serve both as collapse targets and as positional references on release, whereas Sprott A is a hidden attractor with no equilibria at all, so its controlled state carries no comparable positional memory. The landing points of SCC are reported to be insensitive to the waveform as long as the plateau stays below γ_c ^[21]attainable degree of reversibility is governed accordingly by the equilibrium structure of the host, with equilibrium-based systems recovering more faithfully than hidden attractors.

-

Figure 8. MReversible control of the Lorenz attractor showing the transition from chaos to a stable fixed point and the subsequent recovery of the butterfly attractor with high geometric fidelity.

Figure 9. Partial reversible control of the hidden Sprott A attractor showing the transition from chaos to a stable limit cycle and the subsequent recovery of the attractor with preserved topology and observable phase slip.

5.2 Deferral: ti-axis operator.

opological braiding

The opserator can act on all three equations with separate factors (γx, γy, γz). A response-sucond mechanism has a different objective. Rather than reconstructing a previously collapsed attractor, it counteracts the collapse mechanism itself by periodically restoring local instability as nonlinear dissipation increases ^[6]. The geometrfaice study on the Lorenz and Sprott A hosts foun modulation, applied with braiding strength ε, repeatedly reorients trajectories toward the z-expanding direction dominants of the ^[24].

- Other lines.flow, Ldelater work carries the operator into a continuum analoguying the loss of chaotic stretching and broadening the critical transition between chaos and order.

The (aeffect state-dependent bulk viscosity), a Nosé–Hoover–Hala host foris quantified by the displacement of the zero-crossing of the largest Lyapunov exponent with braiding strength: from γ ≈ 0.0146 at Lawsε = 0 ton γ ≈ 0.0188 sat ε = 0.02, γ ≈ 0.0231 ability ε = 0.04, manifoldd γ ≈ 0.0277 at ε = 0.06, and a log-frequency reformulation of wave dispersion. These lie outside the scopeextension of the chaotic range by roughly ninety per cent at the strongest setting. Figures 10 and 11 show the corresponding upward shift of the exponent and the widening of this entrye critical transition region.

7. Context and applications

Figure 10. Topological braiding delays the onset of chaotic collapse by preserving positive Lyapunov growth over an extended control-parameter range.

7.1. Relation to established chaos control

Figure 11. Braiding-induced restoration of local chaotic stretching broadens the critical transition region and produces intermittent recovery of chaotic dynamics near collapse.

The OGY mcorrethodsponding ^[26]phase-space trand delayed feedbackjectories in Figures 12–15 demonstrate that the Lorenz system retains a braided chaotic ^[27]attractor, whereastabilise unstable Sprott A evolves into a braided quasi-periodic orbits that already exist in the uncontrolled flow, using small perturbations that vanish on the target orbit. The Hala ope on a torus rather than collapsing directly to a lower-dimensional state. Poincaré sections (Figure 15) show a thickened closed curve with a dominant return period of approximately fifteen crossings, indicating quasi-periodicity rather than periodicity and confirming that the rich topology is not fully collapsed.

Braiding tor changeherefore defers the vector field: it creates stability at C± and moves them. Its control effort does not vanish at the target.

Nonset of stabilization rather than reversing a collapse that has already occurred. To keep the two mechanisms distinct, this entry refers to that of Section 5.1 as recovery and to that of Section 5.2 as deferralinear state feedback that driv, reserving the term reversible control for the former.

Figure 12. Deferral of collapse in the Lorenz attractor through topological braiding, preserving a braided chaotic state under increasing nonlinear dissipation.

Figure 13. Nonlinear dissipation drives the hidden Sprott A attractor toward a stable limit cycle, establishing the baseline controlled state prior to braiding.

Figure 14. Topological braiding defers collapse in the Sprott A system by transforming the controlled limit cycle into a braided quasi-periodic attractor.

Figure 15. Poincaré analysis indicating the persistence of organized braided dynamics and quasi-periodic structure following deferral of collapse.

6. Applications and Broader Implications

Thes the Lorenz system to fixed points or limit cycles is an established clascollective work on the Hala attractor so far carries implications across several domains. The central proposition — that systems can be engineered whose chaoticity is under control — has ^[28]relevance ^[29]well ^[30]. Thbeyond spplasma physics.

• Plasma physics. Control of chaotic behaviour is fundamental to understanding and manipulating wave–particle interactions, plasma heating mechanisms, and instabilities. It also offers a route to mitigating boundary phenomena in fusion devices, a critical challenge in the pursuit of clean energy.
• Secure communications. Chaos can be harnessed for encryption or signal masking through its unpredictability and sensitivity to initial conditions, with chaotic signals serving as carriers for secure data transmission. Reported results on optomechanical systems, in which modulated coupling switches the response between chaotic and regular behaviour, suggest 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, and the resulting output converted into high-quality true random number generators for cryptography and computational simulation.
• Mixing and fluid dynamics. In industrial and chemical engineering, controlled chaos can optimize mixing processes, improving uniformity of distribution and enhancing reaction efficiency.
• Biological and biomedical systems. The principles of tunable chaos may apply to models of biological systems, from cardiac rhythms to neural networks, where pathological states involve erratic or chaotic dynamics.
• Equilibrium-based hosts of the Lorenz class. High-fidelity recovery enables switchable operation between chaotic exploration and stable equilibrium, of interest wherever the original chaotic dynamics must be reliably restored after temporary stabilization.
• Hidden-attractor hosts. Deferral provides a mechanism for delaying or modulating collapse while preserving essential attractor topology, offering strategies for critical-transition mitigation in complex systems where exact trajectory reconstruction is inherently unattainable.

7. Limitations and Future Directions

Sevecral lific contribution here is the cubic form −γ(x² +mitations qualify the results above and define the path for further work.

• Lyapunov characterization remains incomplete. Full spectra have been computed in specific cases but have not been mapped systematically across parameter space, and the fractal dimensions of the strange attractors under varying parameters have not been established. The Kaplan–Yorke dimension is determined directly by the full spectrum, so the two tasks are one.
• Threshold determination is only partially analytic. The equilibrium fold at γ_c = 1/54 follows in closed form from the symmetric ansatz, as shown in Section 3.1, and the reverse-Hopf threshold γ_H is accessible through the Routh–Hurwitz conditions evaluated at C±. The Lyapunov crossing γ_λ is presently numerical, and the three thresholds have not been assembled into a single analytic map over γ, δ, and the external forcing amplitude and frequency.
• The self-regulation claim carries a detection lag. The online local Lyapunov estimate of Section 4 settles negative only well after γ(t) has crossed the static threshold. The magnitude of that lag, its dependence on window length and renormalisation interval, and whether it can be reduced without loss of robustness remain open.
• Discretization robustness is unestablished. The SCC uniqueness argument holds in exact arithmetic. Under finite precision and timing jitter in T_free, distinct collapses may land within machine epsilon of one another, and the stability of the discrete landing set against round-off and jitter has not been quantified.
• The models are low-dimensional and idealized. Real physical plasma systems are spatially extended, noisy, heterogeneous, and subject to measurement back-action and perturbation. The robustness of chaos tunability under these conditions requires assessment. From the empirical side, more extensive measurements of spatial and temporal variation, higher-resolution diagnostics, and controlled experiments that systematically vary the parameters corresponding to the model would strengthen validation.
• The control mechanisms of Section 5 are demonstrated for two representative low-dimensional hosts, with reversibility evaluated primarily through qualitative phase-space reconstruction and Lyapunov-based indicators. Quantitative trajectory-recovery metrics, robustness to noise and parameter uncertainty, and validation on higher-dimensional or experimentally realized chaotic systems remain to be established.

Future y²)z,work which leaves the orshould accordingly aim to:

• Complete the Lyapunov characterization. Compute the full spectrum across parameter space and relate it to attractor dimensionality and complexity through the Kaplan–Yorke relation.
• Close the analytic threshold map. Derive γ_H in closed form and cross-check the numerically located γ_λ against the analytic fold condition.
• Extend to higher-dimensional models. Move to spatially extended systems described by partial differential equations, better representing realistic plasma behaviour.
• Characterize noise and perturbation response. Investigate stochastic perturbation, measurement noise, and probe-induced perturbation, including the possibility that the measurement strategy itself acts as a control parameter.
• Develop quantitative control metrics. Establish measures of recovery fidelity and deferral effectiveness for hyperchaotic, networked, and experimentally accessible systems, and investigate adaptive, data-driven implementations capable of real-time stabilization, recovery, and resilience enhancement.
• Examine thermodynamic implications. Explore entropy and energy considerations, especially near the Hamiltonian limit, to establish the cost of controlling chaos.
• Pursue prototype development. Translate the theoretical findings into practical devices in controlled plasma heating, boundary instability mitigation in fusion devices, secure communications, and random number generation.

8. Conclusion

Claigms in's nonlinearisation unchanged, makes the contraction state-dependent, and yields closed-form thresholds.

7.2. Plasma physics and gaseous electronics

C dynamics carry an unusual verification burden, and the caveat belongs alongside the results rather than after them. A positive largest Lyapunov exponent computed from a finite trajectory in finite-precision arithmetic is an estimate, not a proof: transient chaos, control in discharges has been demonstrated experimentally without a model, by small modulation of the discharge voltage in a thermionic argon dischargequasi-periodicity of long period, and integrator artifact can each reproduce the signature over the observation window. Newly proposed chaotic systems frequently prove to be smooth coordinate transformations of systems already in the literature, so novelty is properly established through invariants rather than through the appearance of a phase portrait. Control claims carry ^[31]a pand by feedback and synchronisation in plasma diodes, ionisation waves and drift wavesrallel burden. Suppressing chaos is trivial in itself, since sufficient damping stabilizes anything; the substantive question is whether an intervention achieves what indiscriminate dissipation would not — whether it is minimal, ^[32].structured, Aand host model carrying the operator adds predicted thresholds, end statephysically realizable, and whether it stabilizes an orbit already embedded in the flow, defers a transition, or destroys the attractor outright. These are distinct claims and margins before hardware is built.

Toare best reported separately. Models translate imperealise the operator in a device, z must be a slow, controllable variable, fectly to the laboratory, where noise, delay, finite actuator bandwidth, and measurement back-action reshape the problem rather than merely perturb it.

Measuchred as discharge current or anode bias, and x² + y² a measured fluctuation power. As an external loop this is against these standards, the Hala attractor chaotic system occupies a defined position. It is a layered construct rather than a single modified equation. The Hala operator supplies state-dependent feedback multiplier. An internal realisation requires a physical loss channel with the same algebraic form,along a diagnosed axis; beneath it the quenching factor γ is not an externally swept parameter but an embedded logistic state, so that the complete system is autonomous and four-dimensional; and that channel has to be identified for each device.

Thee actions the system performs — one-way discretization under SCC, reversible Hala_Operator_Skillnd partially reversible control, ^[33]and topackages this procedure as a method with no built-in numbers. It selects hosts by structural match to the experiment, rechecks the operator on all three axes, separates conclusionsological braiding — are its purpose rather than a by-product of damping. The canonical chaotic systems on which those actions are exercised are hosts: parametrized substrates established and explored at length by their original authors, serving here as a benchmark set rather than as the object of study. This is why behaviour is reported per host, and why no claim is made that holdthey respond alike.

The appron every host from host-specific oach also sits deliberately outside the dominant lineage in chaos control. Methods in the tradition of Ott, Grebogi and Yorke ^[7], and the times-delayed feedback of Pyragas ^[8], offapply small pers the three protocols of Section 5, and closes with a reliability audit. A pre-registered null-model study setrturbations to stabilize unstable periodic orbits already embedded in an attractor, leaving the attractor itself intact; minimality of intervention is their governing principle and their strength. The Hala operator instead restructures the acceptance rule: a result counts only if it is host-invariantvector field. It alters the equilibrium continuation, introduces a fold available in closed form at γ_c = 1/54, and cano result on a surrogate host is claimed for a device carry the flow to a fixed point, onto a programmable discrete set, or back to its original geometry. The two families answer different ^[34].questions In— that study power-matched linear feedback outperformed the operawhat can be stabilized without disturbing an attractor, and what an attractor can be made to become — and their results are not interchangeable.

The criterion or at 7 of 9 registered gains on the decisive hostf a negative largest Lyapunov exponent is likewise separated throughout from verified long-horizon behaviour, the two coinciding in four of the eleven benchmark hosts.

7.3. Other proposed uses

SCC Whas been proposed for the parametric control of fission reactorst remains outstanding is stated plainly in Section 7: the models are low-dimensional, empirical contact is confined to Langmuir probe characteristics, ^[35],and fror logic primitives built from Sprott A abustness under noise, delay, and finite actuator bandwidth has not been established.

Nonde of Rössler hosts (a clock and a switch)this counsels against the work. It is the discipline that makes results durable, and for secure communication and random-number generation. These proposals rest on numerical modelsthe questions it leaves open — controlling chaos with minimal intervention, containing it where it is destructive, and harnessing it where it is useful — remain among the most consequential in the field.

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.