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 A rspecurrence obtained by discretising a one-dimensional I–V curve is a property of the chosen discretisation. It is not by itself evidencecific parameters, to external forcing, to spatial context, and to measurement perturbation. The framework also connects dissipative and Hamiltonian descriptions of chaos in the plasma, and the entry treats it as motivation only.

Stic systems, providing a route to interpret and control complex dynamince Januaal behaviour.

2. History and Theoretical Development

The progryession 2026 the feedback term has been studof the Hala attractor framework is marked by a series of cumulative theoretical advances. The initial study ^[1] introduced as anthe Hala operator in its own right, separate from the Lorenz hostby modifying a classical Lorenz-type ordinary differential equation system. The body of work now has seven layers: precursors, the attractor,key innovation was the incorporation of a state-dependent nonlinear feedback term, quantified by a quenching factor γ. tThe operator, controequations of the Hala attractor chaotic system are

ẋ = σ(y − x)

ẏ = x(ρ − z) − y

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

where the Hala protocols, a log-frequency extension, applications, and methods for auditing claimoperator acts along the z direction and is defined as.

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

Throughout Tthis entry covers the first fourarticle the classical parameter values inσ = 10, ρ = 28, α = 1, β = 8/3 daretail and summarises the rest assumed unless stated otherwise.

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 opequilibria C± are unstable saddle-focirator is quadratic in the transverse coordinates and linear in atz, γso = 0 and become stable through a Hopf bifurcation at 

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

3.2. Boundary crisis γ_cr

Stit is negligible near the origin and dominant at lablrge C± do not end the chaos at once. Between γ_H and about 0.0016 the chaotic attractor coexists with the two stable eqphase-space excursions. This is the property that allows it to stabilize the unstable saddle-foci of the host system without altering the local behaviour near the origin.

Numerilibria, andcal simulation and analysis of the largest Lyapunov exponent λ₁ stays nearrevealed that as 0.85γ (Figure 2). The attractor is then destroyed in a boundary crisisncreases the system undergoes a continuous transition, through a series ^{[9] [10]}. Beyondf it, chaos survives only as a transient whose lifetime shortens as γ grows. In an ensemble of eight initial conditions, all remainebifurcations, from chaotic strange-attractor behaviour to stable fixed points. The chaotic regime collapses under small changes in a single parameter, establishing that chaoticity 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 d chaoistic after 6000 time units at γ = 0.00164, three at γ = 0.00184, and nonenct thresholds govern the controlled system and should be distinguished. The reverse-Hopf stabilization of the non-trivial equilibria C± occurs at γ_H ≥≈ 0.00196. A gated analysis places the canonical quench76 (Table 1); the largest Lyapunov exponent crosses zero at γ_cr_λ ≈ 0.001968 ^[9].

P(Figure 2); astnd the crisis, λ₁ equals the real part of the lequilibrium branch itself terminates at the fold threshold γ_c = 1/54 ≈ 0.018519, deriveading eigenvalue at C± in Section 3.1. These satisfy γ_H < γ_λ < γ_c. AThe 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 timeintermediate 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 ρ ^[11]≈ 24.06 – 24.74.

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

The sequilibriacond work C± r^[2] ecxpandede to infinity and cease to exist a this foundation by adding a dissipation paramet

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

Foer γδ > γ_c no attracting equilibrium remains and the trajectory is unbounded. The escape is in x and y, not in zand an external periodic forcing term of adjustable amplitude and frequency, in order to explore the interplay among dissipation, resonance, and chaos. 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 growparameter δ 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. The sum of the Lyapunov 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 Js equals the time-averaged divergence of the vector field, Σλᵢ = ⟨∇·F⟩, and governs the evolution of phase-space volume; the largest exponent λ₁ meansuary 2026 the feedback term wasres the rate of separated from the Lorenz system and applied to the Sprott A flow ^[12]ion along the unstable direction. For a three-dimensional autonomouse unmodified host with state (x₁, x₂, x₃), the operator on axis kthe divergence is the constant −σ − 1 − β, fixed by the linear damping coefficients alone, whereas λ₁ is

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

The a outcodyname 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.

A ical property that those coefficients do not directly set. The two can therefore be varied independegative largestntly, and the analysis of the Lyapunov exponent is therefore necessaryspectra, phase-space trajectories, and volume contraction established that they 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

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

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 fothir qud study ^[3] connenching 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.ted these advances to experiment. Current–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 terminating in a chaotic regime. A spatial dependence was observed, with chaos manifesting more Oustcome 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

Throngly 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 intrusion. A suitably parameterized Hala attractor model ree protocols 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 toproduced these behaviours, including the fixed point → limit cycle → period-doubling → chaos route and 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

The itsHala attractor without restoring the earlier trajectory.

5.1. Irreversible: Successive Controlled Collapse (SCC)

SCC admits a deterministic control protocol termed Successwitches the operator on in pulsesve Controlled Collapse (SCC). Brief activation of the feedback factor ^[20]γ ^[21]> ^[22].0 Eacth pulse raises γ rough a smoothly to an sigmoidal ramp of amplitude γ__on fonrces < γ_c, holds it, and returns it to zero. Pulses are chaotic 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._free Dduring a pulse the trajectory contracts towards C±(γ_on); thwhich γ returns to zero, generates families of discrete landing points whose distribution traces the geometry of the underlying strange attractor. Because the collapse phase is globally contracting, distinct pre-collapse state at the end ofs are carried onto common landing points, and the pulse is there-image of a landing point.

Th cannot be recollapse is vered from the point itself. SCC is therefore an irreversible in a practical sense. A smooth flow is invertprotocol, in contrast with the reversible 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 unitsnd partially reversible mechanisms of Section 5. A uniqueness argument 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

The per axis converge to within 2.5 × 10⁻⁹ of C± underquilibria onto which the collapse occurs follow in closed form. Setting γẋ = 0 =gives 0.005y = x, and subackward intstituting into ẏ = 0 yiegrlds x(ρ − z − 1) = 0, so that either x = 0, corresponding from the collapsto the origin, or

z* = ρ − 1 = 27

Thed state diverges.

Twvertical coordinate of the non-trivial features set the landing-point diequilibria is therefore independent of γ. Substributiotuting y = x and z = z* (Finto ż = 0 gurives 7).x*²(α − 2γz*) = βz*, hence

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

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

At γ = 0 this γ_oreturns approacx* = ±6√2 ≈ ±8.485, the clas γ_c thesical Lorenz equilibria recede and the plane z = 27 becomes a slow regC±. 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

ion, because ẏdependent of σ =and 0β. whBeyond γ_c no real non-triver y = x on it. Tial equilibrium exists within the symmetric ansatz, and trajectories then linger without landing (a ghost bottleneck), which is a second reason to keep γ_on well below γ_cno longer relax to static equilibria but are governed by transient dynamics. This value bounds the admissible pulse amplitude for the protocol, as developed in Section 3.3.

Figure 73. 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.

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.

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

3.2 Free-evolutigure 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.

on interval and the landing-point distribution

5.2. Regime reversal: collapse and release

WiThe conthdrawing the operator after a collapse lets the trajectory leave the equilibrium and return to the chuous trajectory evolution during the uncontrolled phase is evaluated by varying the free-evolution interval, 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 cotic attractorndensation at larger T_free represents physical relaxation ^[23]. Oalong the Lorenz host the recovered attrace stable manifold, confirming that longer evolution times allow the controlled 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 amplitorude matches the original statistically: mean z of 23.4 against 23.6, and standard deviust 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 x of 7.85 against 7.91 (Figure 9).

TheSection 3.1 gives x*² = 72/(1 − 7.56) < 0, so no reversal applies to the regimeequilibrium exists for the trajectory to collapse onto. The trajectory after release is not the trajectory before collapse, because all initial conditions that landed on the same equilibrium share one futureresulting 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 operator also removes the time-reversal symmetry of the Nosé–Hoover host: under (x, y, z, t) → (x, −y, −z, −t) the host is invariant and the operator leaves a residual 2γz(x² + y²)ional requirement for maintaining bounded discretization across successive switching cycles is accordingly a calibration constraint on the pulse amplitude, γ_on < γ_c, and not a solver setting.

O

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 Sprott A, plain release does not restore chaos. The operator drives the trajectory onto a limit cycle that lies in a regularpoints recorded at short free-evolution intervals form a tight cluster near z = 27, well inside the region vacated by the equilibrium branch. This clustering is a ghost bottleneck: a localized region of the conservative host. After releasephase space in which the flow retains the slow character of a nearby equilibrium structure that the trajectory stays on a regular orbit there, with λ₁ ≈ 0.0004 against 0.018 for the original chaotic orbit (Figure 10). Recovery on this host needs an additional expansion stepcan 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, the plasne z = ρ − 1 = 27 carrin the Adaptive Reversible Expansion and Harmonic Anchoring scheme,es the equilibrium branch for every value of γ, and on that plane the second component of the vector field vanishes identically, ẏ = x(ρ − z) − y = 0 whenever y = x. The plane ics th reports recovery with moderate phase sliperefore a slow surface of the flow irrespective of where the equilibria themselves lie. ^[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

DSeferralcond, adds a periodic modulation that keeps restoring local instability as the operator's dissipation grows, so that collapse occurs at a larger γ ^[24] s γ approaches γ_c from below the equilibria recede toward large |x| along 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. Thraje figures here use one explicit realisation, a resoctories adiabatically tracking a continuous sigmoidal ramp enter this region and traverse it slowly.

Under ant modulation of short free-evolution interval (T_free = 5.0) inthe 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).

Negration ends while trajectories area still cr 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 unmodulated chaos at γossing the bottleneck, and the recorded landing points accumulate there. The structure is termed a ghost because it 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.0015. Past the collapse it locks) to escape the bottleneck and relax onto a forced oscillation about one equilibrium (Figures 12 and 13)true asymptotic invariant sets. The contrast between the two panels of Figure 4 is the direct signature of this effect.

O

4. Universality of the Hala Operator

The un Spiverott A the unmodulatedsality of the Hala 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). The modulation here restores quasi-periodicity; it does not restore chaos.

“Bis established by direct computation across a benchmark set rather than by appeal to a single host system. Applied in its canonical z-directed form to eleven canonical chaotic systems, the operator drives the largest Lyapunov exponent negative in ten. The long-horizon behaviour behind that count 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 reaiding” is used in a descriptive sense for these orbits. The braid types of periodic orbits ofch zero-dimensional states only at large multiples of their thresholds: Thomas at approximately 2.8 × 10⁴ γ_c theand Lorenz equations-83 at approximately 8 γ_c. wThere classified by Birman and William remainder terminate on sustained one-dimensional limit cycles ^[25];(Sprott A, brAizawaid words for the modulated orbits have not yet been computed.

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

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

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

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

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

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

6. Extensions

Fractional order., Halvorsen) or diverge (Rössler, Rabinovich–Fabrikant). A negative largest Lyapunov exponent Replacing the time derivatives by Caputo derivatives of order q gives the fractional-order Hala systems therefore necessary but not sufficient for quenching, and the two criteria are reported separately here.

Three Ihosts equilibria are those of the integer-order system, so the Matignon stability boundary q_c(γ) is known in closed form frxhibit axis 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 Hopf 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 ^[8] term bounding y and produces catastrophic growth; and for Halvorsen the correct axis remains unresolved, with z evidently unsuitable. This dependence generalises into a concrete per-host axis-selection procedure applicable to any future host, and is the precise sense in which the operator is universal.

-

Timable-dependent quench 1. factor.

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 quopench factor may follow a waveform γ(t): linear, exponential, sigmoidal or oscillatingrator’s sigmoidal control schedule, previously imposed as an external timetable, is the exact solution of a logistic equation. The landing points of SCC are reporteconstruction can therefore be recast as a single autonomous four-dimensional system in which γ is a state variable rather than a parameter ^[4]. Framed tohis way, be insensitive to the waveform as long as the plateau stays below γ_cthe 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 carries it from strange-attractor dynamics to a quenched ^[21]register.

- M

Figulti-axisre operator6. Quenching Tbehaviour of the Hala operator can act on all thracross benchmark systems.

Because equations with separate factors (γx, γy, γz). A response-surface study on the Lorenz and Sprott A hosts found the zγ evolves autonomously and is decoupled from (x, y, z), a single realisation of the coupled 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 dominantlocal ^[24].

- Other lines. Lyater work carries the operator into a punov estimate computed continuum analogue (a state-ously as the trajectory evolves: three independent bulk viscosity), a Nosé–Hoover–Hala host for a Lawson stabilityperturbation directions, renormalised every 1.6 time units, on the Lorenz host with γ _base = 0, γ _target = manifold0.01, and a log-frequency reformulation of wave dispersion. These lie outside the scopek = 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 lag is a quantified limitation of this entrye self-regulation claim.

7. Context and applications

7.1. Relation to established chaos control

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 the OGY method ^[26]sense in which universality and dself-relayed feedbackgulation, used loosely in the original formulation, can ^[27]be stmabilise unstable periodic orbits that already exist in the uncontrollde precise: universality is a per-host axis-diagnosis procedure rather than a fixed rule, and self-regulation is an autonomous property of a single augmented flow, using small pertu rather than a description of externally swept behaviour.

5. Reversible and Partially Reversible Chaos Control

Whereas SCC is irreversibations that vanish on the target orbit. The Hala operator changes the vector field: it creates stability at C± and moves them. Its control effort dole, two further mechanisms provide reversible and partially reversible control. They act at different stages of the collapse and should be distinguished: the first recovers an attractor after collapse has occurred, the second defers collapse before it occurs.

5.1 Recovery: reversible control by collapse and release

The firs not vanish at the target.

Nt mechanism contracts the chaotic trajectory onto a lower-dimensionalinear state feedback that drives invariant set — a stable fixed point in the Lorenz system to fixed points or , a stable limit cycles is an established class in the Sprott A system ^[28]— ^[29]and ^[30]. Tthen specificwithdraws the contribution here is the cubic form −γ(x² + y²)z, which leaves theol, reintroducing phase-space expansion so that the trajectory evolves back toward its origin's linearisational attractor ^[5]. Figunchanged, makes the contraction state-dependent, and yields closed-form thresholds.

7.2. Plasma physics and gaseous electronics

Chaosres 8 and 9 show that the reconstructed Lorenz attractor closely matches the original conbuttrol in discharges has been demonstrated experimentally without a model, by small modulation erfly geometry, whereas the recovered Sprott A attractor preserves its characteristic multi-loop topology but exhibits phase slip and incomplete temporal reconstruction.

Geof methe discharge voltage in a thermionic argon dischargeric recovery alone is a weak criterion. The attractor is an invariant set ^[31] adetermined by feedback and synchronisathe vector field, so once the control is withdrawn any initial condition in plasma diodes, ionisation waves and drift wave the basin regenerates the same geometry; the discriminating quantity is the trajectory’s ^[32].phase Aon host model carrying the operator adds predicted thresholds, end states and margins before hardware is built.

Tthe attractor, which is not recovered in the Sprott A case. The distinction between the two hosts is structural. The Lorenz system possesses equilibria that serve both as collapse targets and as positional references on realise the operator in a device, z must be a slow,lease, whereas Sprott A is a hidden attractor with no equilibria at all, so its controllable variable, such as discharge current or anode bias, and x² + y² a measured state carries no comparable positional memory. The 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. Reversible 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: topological braiding

The second fluctuation power. As an external loop this is a feedback multiplier. An internal realismechanism 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 requirincreases a^[6]. pThysical loss channele geometric modulation, applied with the same algebraic form, and that channel has to be idebraiding strength ε, repeatedly reorients trajectories toward the expanding directions of the flow, delaying the loss of chaotic stretching and broadening the critical transition between chaos and order.

The effect is quantified forby each device.

Tthe displacement of the Hala_Opzerator_Skillo-crossing of ^[33]the plackages this procergest Lyapunov exponent with braiding strength: from γ ≈ 0.0146 at ε = 0 to γ ≈ 0.0188 at ε = 0.02, γ ≈ 0.0231 at ε = 0.04, andure γ ≈ 0.0277 ast ε = 0.06, an method with no built-in numbers. It selects hosts by structural match toextension of the chaotic range by roughly ninety per cent at the strongest setting. Figures 10 and 11 show the corresponding upward shift of the experiment, rechecks the operator on all threonent and the widening of the critical transition region.

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

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

The axcorres, separates conclusionsponding phase-space trajectories in Figures 12–15 demonstrate that hold on every host from host-specific ones, offers the three protocols of Sthe Lorenz system retains a braided chaotic attractor, whereas Sprott A evolves into a braided quasi-periodic orbit on a torus rather than collapsing directly to a lower-dimensional state. Poincaré section 5, ans (Figure 15) show a thickened closesd curve with a reliability audit. A pre-registered null-model study sets the acceptance rule: a result counts only if it is host-dominant return period of approximately fifteen crossings, indicating quasi-periodicity rather than periodicity and confirming that the rich topology is not fully collapsed.

Brainvardiant, and no result on a surrogate host is claimed for a devicng therefore defers the onset of stabilization rather than reversing a collapse ^[34]. In that hastudy power-matched linear feedback outperformed the operator at 7 of 9 registered gains on the decisive host 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.

7.3. Other proposed uses

SCC2 has been proposed for the parametricdeferral, reserving the term reversible control offor 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

The fcollectission reactorsve work on the Hala attractor so ^[35], foar logic primitives built from Sprott A and Rössler hosts (a clock and a switch), and for scarries implications across several domains. The central proposition — that systems can be engineered whose chaoticity is under control — has relevance well beyond plasma 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

Secuvere communication and random-number generation. These proposals resal limitations 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 on numerical models.work should accordingly aim to:

8. Limitations and open questions

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

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

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