REVIEW 3 major objections 4 minor 39 references
Lyapunov stability under $q$-dilatation and $q$-contraction of coordinates
T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Replacing the ordinary Jacobian by a q-deformed one, this paper finds that coordinate contraction lowers positive Lyapunov exponents for chaotic orbits, dilatation raises them, and both tend to destabilize regular KAM-tori orbits.
desk verdict The analytical bifurcation curves are a real new result, but the headline q-trend in the Lyapunov exponents is mostly built into the definition, and the physical interpretation needs a lot more support. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The carrier of the argument is the $q$-derivative, defined by $\partial^{(q)}_x f(x) = (f(qx)-f(x))/((q-1)x)$, which measures change under a $q$-dilatation of the variable rather than a translation. For a function that is homogeneous of degree $\psi$ in a coordinate, this derivative reduces to $[\psi]_q f(x)/x$ with $[\psi]_q = (q^{\psi}-1)/(q-1)$, a $q$-number that equals $\psi$ at $q=1$. Applying this derivative componentwise to the equations of motion gives the $q$-Jacobian $J^{(q)}$, whose products or eigenvalues replace the ordinary Jacobian in the Lyapunov exponent formula, yielding $\lambda^{(q)}$. Since $[1]_q=1$, the $q$-deformation acts only through nonlinear terms, and the homogeneity degrees of those terms determine the direction and size of the effect; this is what makes the Hénon map's quadratic term, the Hénon–Heiles cubic terms, and the diamagnetic Kepler quartic terms respond differently under contraction versus dilatation.
What would settle it
Take a chaotic orbit of the Hénon map, rescale the coordinates literally by $q$ in the equations of motion ($x \to qx$, $y \to qy$), and compute the standard maximal Lyapunov exponent of the rescaled map; smooth coordinate changes leave Lyapunov exponents invariant, so if the standard exponent stays equal to the $q=1$ value while the $q$-Jacobian $\lambda^{(q)}$ changes, the interpretation of $\lambda^{(q)}$ as the Lyapunov stability under coordinate $q$-dilatation is refuted.
Extended reading notes
Core claim
The central claim is that the $q$-deformed Jacobian, obtained by replacing ordinary derivatives with the $q$-derivative in the linearized equations, yields a $q$-dependent Lyapunov exponent $\lambda^{(q)}$ that tracks whether coordinates are contracted or dilated. For $q>1$ the positive $\lambda^{(q)}$ of chaotic orbits increases relative to $q=1$; for $q<1$ it decreases; and for regular orbits both transformations tend to push the $q$-Lyapunov exponent upward, signalling destabilization. In the Hénon map the analytical stability curves $P_{12}$, $P_{1'2'}$, and $P_{24}$, which mark the period-1 to 2, period-1' to 2', and period-2 to 4 bifurcations, vary with $q$, and the two curves that are degenerate at $q=1$ become distinct at $q\neq 1$, so the $q$-deformation resolves a degeneracy invisible in the standard limit. In the conservative systems the same qualitative rule holds for chaotic initial conditions, with exceptions clustered around specific KAM tori, and the maximal $q$-Lyapunov exponent follows a power law in energy in the Hénon–Heiles system and a linear law in the diamagnetic Kepler problem.
Load-bearing premise
The load-bearing premise is that the $q$-deformed Jacobian measures stability under coordinate contraction and dilatation: unlike the ordinary Jacobian, it is not the derivative of any physical perturbation of the trajectory, so the claim that $\lambda^{(q)}$ quantifies coordinate-deformation stability rests on interpreting the algebraic rescaling $[\psi]_q$ as a genuine coordinate effect.
Editorial extensions
If this is right
- In the Hénon map, the analytical curves $P_{12}$ and $P_{1'2'}$ coincide at $q=1$ and split for $q\neq 1$, resolving a degeneracy in the period-doubling structure that the standard analysis cannot see.
- Bifurcation points that are stable at $q=1$ acquire positive $q$-Lyapunov exponents under $q>1$ dilatation, meaning coordinate dilatation destabilizes them while contraction stabilizes them.
- For chaotic initial conditions in both conservative systems, contraction lowers the positive $q$-Lyapunov exponent and dilatation raises it, making $\lambda^{(q)}$ a monotone stability indicator along the $q$ direction for chaos.
- Regular orbits near KAM tori do not follow the simple rule: both contraction and dilatation can increase the $q$-Lyapunov exponent, so the $q$-deformation acts as a generic perturbation for tori, consistent with KAM stability theory.
- The energy dependence of the maximal $q$-Lyapunov exponent remains a power law in the Hénon–Heiles system and a linear law in the diamagnetic Kepler problem, with $q$ changing the fitted exponents or the slope and intercept.
Reading between the lines
- A direct testable extension is to literally rescale coordinates in the equations of motion, $x \to qx$, compute the standard Lyapunov exponent of the rescaled system, and compare it with $\lambda^{(q)}$; the paper's interpretation predicts a match only if the $q$-derivative genuinely captures coordinate dilatation, and a mismatch would isolate the algebraic part of the effect.
- The $q=1$ degeneracy splitting suggests that $q$ can be used as a symmetry-breaking parameter to unfold other hidden coincidences in bifurcation diagrams, potentially serving as a numerical continuation tool for stability curves in maps with higher-degree terms.
- Since $[\psi]_q$ is increasing in $\psi$ for $q>1$ and decreasing for $q<1$, a heuristic prediction is that chaotic orbits whose instability is dominated by higher-degree homogeneous terms show the strongest response to dilatation; classifying orbits by the dominant degree in the $q$-Jacobian could turn the observed trend into a per-term prediction.
- The exceptions near KAM tori may reflect that tori with different rotation numbers respond differently to the algebraic weighting of the $q$-Jacobian; connecting the size of the $q$-effect to the torus's winding number would be a natural next step.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a q-deformed Lyapunov exponent obtained by replacing ordinary derivatives in Jacobians with q-derivatives, and applies it to the Hénon map, the Hénon-Heiles system, and the diamagnetic Kepler problem. The central reported findings are that q-contraction (q<1) decreases positive Lyapunov exponents for chaotic trajectories, q-dilatation (q>1) increases them, and both transformations tend to increase the exponents of regular orbits. For the Hénon map the authors also derive analytical bifurcation curves P1, P12, P1'2', and P24, which split at q≠1 and are compared with numerically computed maximal q-Lyapunov exponents. The paper claims that these results demonstrate Lyapunov stability under coordinate contraction and dilatation and reveal a hidden degeneracy in the usual q=1 bifurcation structure.
Significance. If the proposed interpretation were established, the work would provide a general formalism for assessing stability under nonlinear coordinate deformations and would extend q-deformation methods to Lyapunov analysis in both dissipative and conservative systems. The analytical Hénon bifurcation curves and their numerical verification are self-contained and potentially useful, as is the observation that the q=1 degeneracy of P12 and P1'2' is lifted for q≠1. However, the central physical claim is not currently supported: the q-Jacobian is an algebraic rescaling of the ordinary Jacobian evaluated on the same q=1 trajectory, not the linearization of any defined q-deformed dynamics, and the main monotonic trend follows directly from the definition of [d]_q. The manuscript also contains a concrete formal error in the Hénon sum rule. The numerical results may be reproducible, but the interpretation needs either a dynamical derivation or an explicit reframing.
major comments (3)
- [Section II, Eq. (9)] The definition of the q-Lyapunov exponent is not derived from a variational equation. The text states that 'it is possible to construct the corresponding linearized equations dy(q)/dt = J(q) y(q)' using the q-Jacobian, but no q-deformed map or vector field is defined whose tangent dynamics would produce this equation. The trajectories entering J(q) are still generated by the original q=1 equations, and J(q) is obtained by pointwise rescaling of the ordinary Jacobian entries. Consequently, the claim that λ(q) measures 'Lyapunov stability under dilatation and contraction of coordinates' is an interpretation without derivation. The manuscript should either define a concrete q-deformed dynamics whose variational equation is dy(q)/dt = J(q)y(q), or clearly state that λ(q) is a diagnostic of the rescaled Jacobian along unchanged trajectories.
- [Section III A, text after Eq. (11)] The statement 'we have two q-LEs, λ1(q) > λ2(q), which satisfy λ1(q) + λ2(q) = -b' is incorrect. The determinant of J(q) in Eq. (11) is det J(q) = -b, independent of q, so the sum of the Lyapunov exponents is ln|b|, not -b. This is a concrete formal error in a central example and should be corrected; it also indicates that the q-LE is being treated as if it were an ordinary Lyapunov exponent of the matrix product, for which the standard sum rule would apply.
- [Sections III B 1, III B 2, and IV] The main qualitative trend—contraction lowers positive q-LEs and dilatation raises them—is essentially a consequence of the definition. Each nonlinear monomial derivative is multiplied by [d]_q, and [d]_q is strictly increasing in q for every d≥2, e.g. [2]_q = q+1. Since the trajectory is unchanged and the q-LE is a time average of log norms of these rescaled matrices, the stated monotonicity is expected from the rescaling alone and does not by itself establish a dynamical effect. The paper acknowledges exceptions for KAM tori in the abstract and conclusion, further showing that the trend is not universal. To make the central claim load-bearing, the authors should provide a quantitative argument separating the rescaling effect from any genuinely new stability property, or reframe the contribution as an analysis of rescaled Jacobian spectra.
minor comments (4)
- [Figure 3 caption] The caption labels the bottom row as 'q = 0.5', but the text and the subsequent discussion indicate that the bottom row should be q = 1.5; this typo should be fixed.
- [Eq. (8) and Eq. (9)] The notation |J(t)| and |Ji| is ambiguous for matrices; the authors should specify the matrix norm used in the numerical computations, for example the induced 2-norm or singular value norm.
- [Reference [8]] The author name 'Ayse' appears incomplete or incorrect; the reference should be verified and completed, since it is the source of the q-derivative definition in Eq. (1).
- [Section III B 1, Fig. 5] The power-law fits report many significant digits without uncertainties or goodness-of-fit measures, which would be helpful for assessing the claim that the exponent changes with q.
Circularity Check
The main q-LE stability trend is a rescaling artifact of the q-Jacobian definition; the analytical bifurcation curves remain self-contained.
-
other
[Abstract; Section II Eq. (9); Section III A Eq. (11)]
"Simulations show that q-contraction (q-dilatation) generally decreases (increases) positive Lyapunov exponents relative to the q = 1 case... J_n^(q) = ... -[ψ_x^(2)]_q x_n ... b; 1; 0, [ψ_x^(2)]_q = (q^2-1)/(q-1), and λ^(q) = lim_{n→∞} (1/n) ln(∏_i |J_i^(q)|)."
For the Hénon map, [2]_q = q+1, so J^(q) is the ordinary Jacobian with the (1,1) entry multiplied by (q+1)/2, evaluated on the same q=1 trajectory {x_n}. More generally, the q-Jacobian is the standard Jacobian with each monomial derivative multiplied by the q-number [d]_q, which is strictly increasing in q for every degree d≥2. Since λ^(q) in Eq. (9) is a time-average of logarithms of norms of these pointwise rescaled matrices, the reported trend that contraction decreases and dilatation increases positive Lyapunov exponents follows by construction from the monotonicity of [d]_q. No q-deformed evolution law is defined: the orbit feeding J^(q) is still generated by the ordinary q=1 equations. Thus the central numerical claim is equivalent to the definitional rescaling in Eqs.
full rationale
Most of the paper is self-contained numerical computation: the analytical bifurcation curves P1, P12, P1'2', P24 are derived from the q-Jacobian determinant and match the zero-LE contours, and these mathematical consequences are not circular. The central interpretive claim, however, is the q-dependence of positive Lyapunov exponents. That claim reduces by construction to the definition of the q-derivative: for the Hénon map, [2]_q = q+1, so J^(q) is the ordinary Jacobian with the (1,1) entry multiplied by (q+1)/2, evaluated on the same q=1 trajectory; more generally, every monomial derivative in the q-Jacobian is multiplied by [d]_q, a strictly increasing function of q. Because λ^(q) is the time-average of log norms of these pointwise rescaled matrices, the stated 'contraction decreases, dilatation increases' pattern for chaotic orbits follows immediately from that monotonicity rather than from an emergent property of the systems. No q-deformed tangent-space dynamics is defined, so the paper's assertion that this measures Lyapunov stability under coordinate transformations is an interpretation imported from Ref. [12] (same author group) and from the q-derivative formalism, not a derived result. The analytical content and the PSS simulations retain independent value, but the main stability conclusion is a rescaling artifact. The inconsistent sum rule λ1^(q)+λ2^(q) = -b (whereas det J^(q) = -b would give ln|b| for two-dimensional maps) is a correctness issue rather than a circularity, but it reinforces that the q-LE is treated formally. Overall score 6: the central 'prediction' for positive Lyapunov exponents reduces by construction, while the analytical curves and the regular-orbit exceptions are not circular.
Assumptions & free parameters
free parameters (2)
- power-law exponent β for Hénon-Heiles =
β ≈ 2.99 (q=0.5), 3.45 (q=1.0), 2.38 (q=1.5)
- linear fit coefficients (a,b) for diamagnetic Kepler =
q=0.5: (1.12334, -0.290927); q=1.0: (1.62469, -0.480846); q=1.5: (2.41462, -0.612911)
assumptions (3)
- ad hoc to paper The q-derivative defined in Eq. (1) and its homogeneous extension in Eq. (4) correctly represent the effect of coordinate dilatation/contraction on stability.
- domain assumption The q-Lyapunov exponent computed from the q-Jacobian via Eq. (9) captures the stability of trajectories under the q-deformed perturbation.
- domain assumption The standard numerical algorithm [14] for Lyapunov exponents remains valid when the Jacobian is replaced by the q-Jacobian.
Cite this review
Pith. "Pith review of Lyapunov stability under $q$-dilatation and $q$-contraction of coordinates." pith.science (2026). https://pith.science/paper/P2ZGAZBV
@misc{pith2026241118691,
author = {Pith},
title = {Pith review of: Lyapunov stability under $q$-dilatation and $q$-contraction of coordinates},
year = {2026},
howpublished = {\url{https://pith.science/paper/P2ZGAZBV}},
note = {Machine review of arXiv:2411.18691}
}
abstract
This study examines the Lyapunov stability under coordinate $q$-contraction and $q$-dilatation in three dynamical systems: the discrete-time dissipative H\'enon map, and the conservative, non-integrable, continuous-time H\'enon-Heiles and diamagnetic Kepler problems. The stability analysis uses the $q$-deformed Jacobian and $q$-derivative, with trajectory stability assessed for $q > 1$ (dilatation) and $q < 1$ (contraction). Analytical curves in the parameter space mark boundaries of distinct low-periodic motions in the H\'enon map. Numerical simulations compute the maximal Lyapunov exponent across the parameter space, in Poincar\'e surfaces of section, and as a function of total energy in the conservative systems. Simulations show that $q$-contraction ($q$-dilatation) generally decreases (increases) positive Lyapunov exponents relative to the $q = 1$ case, while both transformations tend to increase Lyapunov exponents for regular orbits. Some exceptions to this trend remain unexplained regarding Kolmogorov-Arnold-Moser (KAM) tori stability.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[2]
The angular momentum Lz around z is a conservative quantity, and we consider Lz = 0
The diamagnetic problem The Hamiltonian of the diamagnetic Kepler problem with total energy ED is given by [26] HD = 1 2 p2 x + p2 y + p2 z + γ x2 + y2 − e r , (16) with γ = e2B2 8mc2 , r = p x2 + y2 + z2, e and m are electron’s charge and mas m, and the intensity of the static magnetic field in the z-direction is B = 4.72 × 105T . The angular momentum Lz...
work page 2021
-
[1]
The H´ enon-Heiles problem The H´ enon-Heiles (H-H) problem is a classic example of a dynamical system that exhibits chaotic behavior. It was introduced by M. H´ enon and C. Heiles in 1964 [20] as a simplified model to study the motion of stars in a galaxy. The Hamiltonian is given by HH = 1 2 p2 x + p2 y + 1 2 Ax2 + By 2 + Dx2y − C 3 y3, (12) where A, B,...
work page 1964
-
[3]
C. F. Gauss, reprinted in Werke 3, 123 (1876)
- [4]
- [5]
-
[6]
Heine, Handbuch der Kugelfunktionen, Theorie und Anwendung , Vol
E. Heine, Handbuch der Kugelfunktionen, Theorie und Anwendung , Vol. 1 (Reimer, Berlin, 1878)
-
[7]
G. Gasper and M. Rahman, Hypergeomeytic Series (University Press: Cambridge, 1997)
work page 1997
-
[8]
C. Tsallis, A. R. Plastino, and W. M. Zheng, Chaos Solitons Fractals 8, 885 (1997)
work page 1997
Show all 39 references
-
[9]
U. M. S. Costa, M. Lyra, A. R. Plastino, and C. Tsallis, Phys. Rev. E 56, 245 (1997)
1997
-
[10]
Ayse, Phys
E. Ayse, Phys. Rev. Lett. 78, 3245 (1997)
1997
-
[11]
Lyra and C
M. Lyra and C. Tsallis, Phys. Rev. Lett. 80, 53 (1998)
1998
-
[12]
Jaganathan and S
R. Jaganathan and S. Sinha, Phys. Lett. A 338, 277 (2005)
2005
-
[13]
Patidar and K
V. Patidar and K. K. Sud, Nonlinear Sci. Numer. Simulat. 14, 827 (2009)
2009
-
[14]
P. C. Rech and M. W. Beims, J. Phys. : Conf. Ser. 246, 012006 (2010)
2010
-
[15]
Barreira and Y
L. Barreira and Y. B. Pesin, Lyapunov exponents and smooth ergodic theory (Providence, Rhode Island, 2022)
2022
-
[16]
A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, Physica D 16, 285 (1985)
1985
-
[17]
H´ enon, Commun
M. H´ enon, Commun. Math. Phys. 50, 69 (1976)
1976
-
[18]
J. A. C. Gallas, Physica A 222, 125 (1995)
1995
-
[19]
D. L. Hitzl and F. Zele, Physica D 14, 305 (1985)
1985
-
[20]
M. W. Beims and J. A. C. Gallas, Physica A 238, 225 (1997)
1997
-
[21]
J. A. C. Gallas, Phys. Rev. Lett. 70, 2714 (1993)
1993
-
[22]
H´ enon and C
M. H´ enon and C. Heiles, The Astr.Journal69, 73 (1964)
1964
-
[23]
Y. F. Chang, M. Tabor, J. Weiss, and C. Corliss, Phys. Lett. A 85, 211 (1981)
1981
-
[24]
Y. F. Chang, M. Tabor, and J. Weiss, J. Math. Phys. 23, 531 (1982)
1982
-
[25]
Grammativos, B
B. Grammativos, B. Dorizzi, and R. Padjem, Phys. Lett. 89A, 11 (1982)
1982
-
[26]
Sahadevan
R. Sahadevan. and M. Lakshmanan, J. Phys. A: Math. Gen. 19, L949 (1986)
1986
-
[27]
I. I. Shevchenko and A. V. Mel’nikov, JETP Letters 77, 642 (2003)
2003
-
[28]
Friedrich and D
H. Friedrich and D. Wintgen, Phys. Rep. 183, 37 (1989)
1989
-
[29]
Holle, G
A. Holle, G. Wiebusch, J. Main, B. Hager, H. Rottke, and K. Welge, Phys. Rev. Lett. 56, 2594 (1986)
1986
-
[30]
Welch, M
G. Welch, M. Kash, C. Iu, L. Hsu, and D. Kleppner, Phys. Rev. Lett. 62, 893 (1989)
1989
-
[31]
M. Kash, L. Hsu, and D. Kleppner, Phys. Rev. Lett. 63, 1133 (1989)
1989
-
[32]
Delande, A
D. Delande, A. Bommier, and J. Gay, Phys. Rev. Lett. 66, 141 (1991)
1991
-
[33]
J. M. Mao and J. B. Delos, Phys. Rev. A 45, 1746 (1992)
1992
-
[34]
M. W. Beims and G. Alber, Phys. Rev. A 48, 3123 (1993)
1993
-
[35]
M. W. Beims and G. Alber, J. Phys. B 29, 4139 (1996)
1996
-
[36]
M. W. Beims, Phys. Rev. A 56, R2503 (1997)
1997
-
[37]
M. W. Beims and J. A. C. Gallas, Phys. Rev. A 62, 043410 (2000)
2000
-
[38]
A. J. Lichtenberg and M. A. Lieberman, Regular and Chaotic Dynamics (Springer-Verlag, New York, 1992)
1992
-
[39]
Manchein and M
C. Manchein and M. W. Beims, Caos, Solitons and Fractals 39, 2041 (2009)
2009
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.