{"id":"60a63c43-c88a-4231-bbac-7f21122500f7","arxiv_id":"1908.01467","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"The regime classification is built on a time-evolution formula that is algebraically inconsistent with the model's energy eigenvalues.","lead":"This paper claims that the expectation values of position and momentum in a q-deformed harmonic oscillator show periodic, quasi-periodic, or chaotic behavior depending on the deformation parameter q and the coherent amplitude alpha_q. The time series used for this classification contain a derivation error, so the reported dynamics do not correspond to the stated model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model cannot support genuine chaos: point-spectrum unitary evolution makes every bounded-observable expectation value almost periodic, so the reported positive Lyapunov exponents are numerical artifacts; Eq. (10) also has a separate phase error.","rationale":"The reader's rejection is correct, and our independent analysis makes it more robust. The Eq. (10) phase typo is real: substituting Eq. (9) into E_n gives E_n=(1+q^2[n])/2, so the time-evolution phase should be exp(-it/2)exp(-it q^2[n]/2); the printed exp(-it(1+q^2)[n]/2) and the exponents in Eqs. (12)-(13) describe a different dynamics. But even a corrected calculation could not support the headline claim, since the model's unitary evolution with discrete, accumulating spectrum yields almost periodic expectation values; positive Lyapunov exponents cannot be genuine. The paper's time-series diagnostics are heuristic and, in the absence of code, data, or convergence checks, cannot overturn this. We therefore keep the REJECT verdict, while identifying the almost-periodicity objection rather than the phase typo as the deepest load-bearing problem. Our agreement with the reader is partial because their weakest assumption is valid but not the fundamental one.","tokens_in":9500,"tokens_out":15759,"duration_ms":157495,"concrete_test":"Recompute Fig. 8 using the corrected phase factor exp(-it/2)exp(-it q^2[n]/2) and with data lengths T=10^4, 10^5, 10^6 while varying the embedding dimension; if the Rosenstein/Wolf lambda_max decreases systematically toward zero as T grows (the expected behavior for an almost periodic trigonometric series), the chaotic regime is a numerical artifact. If lambda_max stabilizes at a positive value, check the point-spectrum/almost-periodicity argument directly, because a positive value would contradict it.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing condition for the paper's central claim is that the reported positive Lyapunov exponents indicate genuine chaos. That condition cannot hold for this model. For 0<q<1 the Hamiltonian (1) has a pure point spectrum: E_n=(1+q^2[n])/2, with [n]=(1-q^{2n})/(1-q^2), and the eigenvalues accumulate at the finite value 1/(2(1-q^2)). Under unitary evolution, the state is a superposition of eigenstates with discrete phases, so any bounded-observable expectation value (and <X(t)>_q, <P(t)>_q are of this type because X,P are bounded here) is an absolutely convergent trigonometric series in the frequencies E_m-E_n. Such a function is almost periodic. Almost periodic functions have zero largest Lyapunov exponent and cannot exhibit exponential separation of nearby trajectories. The positive lambda_max in Sec. III E and Fig. 8 therefore must be a finite-time or reconstruction artifact of the Rosenstein/Wolf algorithms, not a property of the q-deformed oscillator. This objection is independent of the phase-factor error in Eq. (10), which is another real defect: using Eq. (9) the phase should be exp(-it/2)exp(-it q^2[n]/2), not exp(-it/2)exp(-it(1+q^2)[n]/2), and this propagates into Eqs. (12)-(13). Either way, the claimed chaotic regimes are unsupported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the q-deformed harmonic oscillator with Hamiltonian H_q = (AA† + A†A)/2 and deformed commutation relation AA† − q^2 A†A = I. It constructs a q-deformed coherent state, derives the time-evolved expectation values ⟨X(t)⟩_q and ⟨P(t)⟩_q, and classifies their time series as periodic, quasi-periodic, or chaotic on the basis of recurrence plots, power spectra, first-return-time distributions, and Lyapunov exponents. The central claim is an explicit regime classification in the (q, α_q) plane (Fig. 10), with positive largest Lyapunov exponents in the chaotic regions.","tokens_in":9896,"tokens_out":17543,"duration_ms":164014,"significance":"If the central claim were correct, the paper would be a useful extension of earlier expectation-value-based studies of quantum dynamics to a q-deformed oscillator, and the regime classification could interest the quantum-optics community. The manuscript is explicit about the time-evolution formulas and provides a concrete set of falsifiable predictions, which are strengths. However, the model's energy spectrum is not consistently defined, and, independently of that ambiguity, the point-spectrum unitary evolution of this bounded-observable system cannot produce positive Lyapunov exponents. The reported numerical diagnostics therefore do not support the claimed chaotic regimes, and the classification in Fig. 10 is not reliable.","major_comments":[{"comment":"The energy spectrum is defined inconsistently. From Eq. (7), A†A|n⟩ = [n]|n⟩ and AA†|n⟩ = [n+1]|n⟩, so Eq. (1) gives E_n = ([n]+[n+1])/2 = [n] + q^{2n}/2 = (1+(1+q^2)[n])/2. The printed Eq. (3), ([n]+q^{2n})/2, does not follow from the Hamiltonian and its q→1 limit is (n+1)/2, not the stated E_n = n+1/2. The phase factor in Eq. (10), exp(−it/2)exp(−it[n](1+q^2)/2), corresponds to the former spectrum, not to the printed Eq. (3). Thus the manuscript uses two different spectra. If Eq. (3) is intended, Eq. (10) should contain exp(−it q^2[n]/2) instead of exp(−it(1+q^2)[n]/2); if Eq. (10) is intended, Eq. (3) must be corrected. All subsequent expressions, Eqs. (11)–(13), and all numerical results inherit this ambiguity, so this is a load-bearing inconsistency.","section":"II A, Eqs. (1), (3), (7), (9), (10)"},{"comment":"Positive Lyapunov exponents are impossible for this model. For 0<q<1 the Hamiltonian has a pure point spectrum with eigenvalues accumulating at a finite value, and X and P defined in Eqs. (5)–(6) are bounded operators because A and A† are bounded (sup_n [n] = 1/(1−q^2)). Consequently ⟨X(t)⟩_q and ⟨P(t)⟩_q are absolutely convergent trigonometric series in the bounded frequency differences E_m−E_n, i.e., they are almost periodic functions. Almost periodic functions do not exhibit exponential separation of nearby trajectories, so the positive λ_max values reported in Sec. III E and Fig. 8 must be finite-time or reconstruction artifacts of the Rosenstein/Wolf algorithms, not genuine properties of the oscillator. This objection is independent of the spectrum ambiguity in the previous comment and invalidates the central claim of chaotic regimes.","section":"III E, Fig. 8; Eqs. (1), (5), (6)"},{"comment":"The first-return-time analysis is circular as evidence for chaos. The paper defines ergodic behavior through the exponential first-return-time distribution F_1(τ) = (1/τ)e^{−t/τ} and then uses the fact that certain time series fit this distribution as a verification of ergodicity and as support for chaos. This applies the same criterion both as definition and as confirmation. Moreover, exponential first-return statistics also occur for non-chaotic stochastic and multi-frequency quasi-periodic signals, so this diagnostic cannot compensate for the absence of genuine Lyapunov exponents.","section":"II C 3 and III D"}],"minor_comments":[{"comment":"In the second terms of Eqs. (12) and (13), the power of |α_q| appears to be off by one factor: the A† contribution should carry α_q^* |α_q|^{2n} (equivalently α_q^{−1}|α_q|^{2n+2}) rather than α_q^{−1}|α_q|^{2n+1}.","section":"Eqs. (12)–(13)"},{"comment":"The caption of Fig. 2 specifies q = 0.95 and |α_q|^2 = 1, while Sec. III A discusses the autocorrelation function peaking at random intervals for α_q = 1 and q = 0.9; the caption and the text should be made consistent.","section":"Fig. 2 and Sec. III A"},{"comment":"The abstract's claim that detailed analysis of chaotic quantum systems based on expectation values 'has not been reported in the literature' is overstated, because Refs. [6–8] and [31] already perform such analyses; the novelty should be stated as the extension to q-deformed systems.","section":"Abstract and Introduction"},{"comment":"The regime diagram in Fig. 10 is presented as an 'approximate demarcation' but the text does not explain how the boundaries were determined, what the axes' scaling is, or how the diagram was constructed from the numerical diagnostics; this makes the central classification difficult to reproduce.","section":"Fig. 10"},{"comment":"The normalization of the q-deformed coherent state in Eq. (8) and the bound in Eq. (15) are cited from Ref. [24] but not derived; since the q-exponential does not satisfy e_q(x)e_q(−x)=1, the normalization convention should be stated explicitly.","section":"II A, Eq. (8) and Eq. (15)"},{"comment":"The parameter choices for the Lyapunov computations (embedding dimension, delay, and tolerance) are not fully specified for all figures; Fig. 9 shows m = 4, 6, 8 but the delay and the length of the linear fitting region are not given, which is important because the positive exponents are likely numerical artifacts.","section":"III E, Fig. 9"}],"recommendation":"reject","confidential_remarks":"The almost-periodic argument in the second major comment is, in my view, decisive: a bounded observable under point-spectrum unitary evolution cannot have a positive Lyapunov exponent, so the paper's numerical classification is not a property of the q-deformed oscillator. The additional spectrum inconsistency between Eq. (3) and Eq. (10) makes the manuscript unreliable even as a numerical study. I do not see a fix that stays within the paper's current scope, so rejection is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, this one is a clean reject, and the quickest path to that verdict isn't the numerical evidence—it's the model. H has discrete eigenvalues E_n = (1+q^2[n])/2 accumulating at a finite value, and X and P are bounded operators for q<1. Under unitary evolution, <X(t)> and <P(t)> are absolutely convergent trigonometric series in the differences E_m - E_n. That makes them almost periodic. Almost periodic functions do not have exponential separation of nearby trajectories, so a positive largest Lyapunov exponent from Rosenstein or Wolf is a finite-time or reconstruction artifact, not chaos. This objection does not depend on any detail of their time series.\n\nThere is also a separate algebraic error. Eq. (10) writes the phase factor as exp(-it/2) exp(-it(q^2+1)[n]/2), but Eqs. (3)-(4) and identity (9) give E_n = (1+q^2[n])/2. The phase should be exp(-it/2) exp(-it q^2[n]/2). The extra exp(-it[n]/2) propagates into Eqs. (12)-(13) and into every numerical result. So even before the almost-periodicity objection, the time series they analyze is not the one generated by their model.\n\nCredit where due: the paper is clearly written; applying the expectation-value toolbox (recurrence plots, power spectra, return times, Lyapunov exponents) to a q-deformed oscillator is a natural extension; the regime diagram in Fig. 10 is a genuine attempt at classification; and they cite relevant prior work, including Buzek, whose phase-space plots overlap with theirs. But no code or data are given, and the classification rests on the two defects above. The first-return-time analysis is also self-referential—the exponential fit is used both to define and to verify ergodicity—though that is minor compared with the main issue.\n\nFor whom is this paper? Someone curious about q-deformed oscillator dynamics might glance at the figures, but as a research claim it does not hold. I would not send it to peer review; a referee would return the same two objections. If the authors fix the phase factor and re-examine the numerics, they will likely find that the positive Lyapunov exponents disappear.","headline":"The q-deformed oscillator has a pure point spectrum and bounded X, P, so the reported positive Lyapunov exponents are numerical artifacts; Eq. (10) also has a phase-factor error. The paper is readable but the central claim doesn't survive.","tokens_in":10380,"tokens_out":3176,"would_cite":false,"duration_ms":29778,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The q-deformed harmonic oscillator can be periodic, quasi-periodic, or chaotic depending on q and the coherent amplitude.","keywords":["q-deformed harmonic oscillator","quantum chaos","expectation values","coherent states","Lyapunov exponent","recurrence plot","power spectrum","first-return-time distribution"],"falsifier":"Numerically integrate the time-dependent Schrödinger equation for $H_q$ using the matrix elements in Eq. (7) for $q = 0.9$ and $\\alpha_q = 1$, extract $\\langle X(t)\\rangle_q$, and compare its recurrence plot and largest Lyapunov exponent with those reported here; any significant difference shows the analytic time evolution in Eq. (10) needs correction.","tokens_in":9342,"feed_emoji":"🌀","tokens_out":11687,"duration_ms":98628,"temperature":0.7,"pith_summary":"The paper claims that the q-deformed harmonic oscillator, despite obeying a linear Schrödinger equation, can show periodic, quasi-periodic, or chaotic dynamics in the expectation values of its deformed position and momentum operators, depending on the deformation parameter $q$ and the coherent amplitude $\\alpha_q$. This matters because quantum systems are often expected to be quasi-periodic; here, the non-linearly spaced deformed energy levels produce genuinely irregular time series. The authors classify the regimes with recurrence plots, power spectra, first-return-time distributions, and Lyapunov exponents, and find positive Lyapunov exponents in the chaotic region. The central claim is that the dynamics of observables in this simple deformed system can be chaotic without invoking a classical limit.","feed_headline":"q-deformed oscillator can be chaotic, periodic, or quasi-periodic","feed_subtitle":"Expectation values of position and momentum reveal a full phase diagram in the q–α_q plane.","key_machinery":"The engine of the paper is the q-number $[n] = (1 - q^{2n})/(1 - q^2)$, which replaces $n$ in the ladder-operator relations $A|n\\rangle_q = \\sqrt{[n]}\\,|n-1\\rangle_q$ and $A^\\dagger|n\\rangle_q = \\sqrt{[n+1]}\\,|n+1\\rangle_q$. This q-number makes the energy spectrum non-linear and produces multi-frequency time dependence in the coherent-state expectation values. The other load-bearing piece is the time-evolved deformed coherent state in Eq. (10), whose phase factors $e^{-it/2} e^{-it[n](1+q^2)/2}$ generate the sums in Eqs. (12)-(13) for $\\langle X(t)\\rangle_q$ and $\\langle P(t)\\rangle_q$; the resulting time series are then fed into four diagnostics: recurrence plots, power spectra, first-return-time distributions, and the standard largest-Lyapunov-exponent algorithms.","core_discovery":"Starting from the q-deformed Hamiltonian $H_q = \\frac{1}{2}(AA^\\dagger + A^\\dagger A)$ with $AA^\\dagger - q^2 A^\\dagger A = I$, the energy levels are $E_{q,n} = ([n] + q^{2n})/2$, with $[n] = (1-q^{2n})/(1-q^2)$. Because $[n]$ grows non-linearly with $n$, the spectrum is not equally spaced, unlike the ordinary oscillator. The paper shows that for coherent states $|\\alpha\\rangle_q$, the expectation values $\\langle X(t)\\rangle_q$ and $\\langle P(t)\\rangle_q$, computed from the time-evolved state in Eq. (10), are periodic for $q \\lesssim 0.1$ (at $\\alpha_q = 1$), quasi-periodic for $0.1 < q \\lesssim 0.2$ and again near $q \\to 1$, and chaotic for roughly $0.2 < q < 0.99$, with the chaotic window widening as $\\alpha_q$ grows. Positive largest Lyapunov exponents and exponential first-return-time distributions confirm the chaotic classification.","pith_inferences":["A direct consequence of the paper's own Eqs. (3) and (9) is that $E_{q,n} = 1/2 + q^2[n]/2$, so the phase factor in Eq. (10) would read $e^{-it/2} e^{-it q^2[n]/2}$. Recomputing the time series with this factor would test whether the regime boundaries in Fig. 10 move.","The autocorrelation decay in Fig. 2 suggests a quantum-information probe: in the chaotic band, the fidelity of the deformed coherent state should decay irregularly, and quantifying that decay could give a state-based signature of the transition.","Because the chaotic range depends on $\\alpha_q$, the system is a tunable source of complex time series; engineering the deformation could serve as a controlled testbed for randomness generation or for time-series classification methods."],"forward_implications":["Chaos in a quantum system can be diagnosed from expectation-value time series of just two observables, without needing a classical or semiclassical limit.","The regime map in Fig. 10 gives a parameter-by-parameter recipe: for fixed $\\alpha_q$, sweeping $q$ moves the system through periodic, quasi-periodic, chaotic, and back toward quasi-periodic and periodic phases.","Positive Lyapunov exponents in the chaotic band imply exponential divergence of nearby initial coherent states, so the deformation parameter controls sensitivity to initial conditions.","Because the chaotic range widens with $\\alpha_q$, larger coherent amplitudes at fixed deformation are more chaotic, while strong deformation restores order.","The same diagnostic toolbox can be applied to other deformed algebras to test whether the periodic-to-chaotic regime sequence is universal."],"supporting_citations":[{"why":"Defines the q-deformed Hamiltonian, its energy spectrum, and the normalization bound on $|\\alpha_q|^2$.","marker":"[24]"},{"why":"Shows that expectation values in the non-deformed oscillator are periodic, providing the baseline comparison.","marker":"[25]"},{"why":"Introduces the expectation-value time-series approach and first-return-time analysis used here.","marker":"[6]"},{"why":"Supplies the exponential first-return-time distribution criterion for ergodic chaotic behaviour.","marker":"[7]"},{"why":"Provides the standard algorithm used to estimate the largest Lyapunov exponents.","marker":"[29]"},{"why":"Provides an independent algorithm used to cross-check the Lyapunov-exponent results.","marker":"[30]"},{"why":"Gives the recurrence-plot formalism and the diagonal-line criteria for classifying periodic, quasi-periodic, and chaotic series.","marker":"[26]"},{"why":"Supplies the forms of the deformed position and momentum operators $X$ and $P$ used in Eqs. (5)-(6).","marker":"[20]"},{"why":"Provides the deformed coherent-state and ladder-operator relations used in the time-evolution calculation.","marker":"[23]"}],"fun_headline_variants":["q-deformed oscillator: periodic, quasi-periodic, or chaotic","Chaos emerges in q-deformed oscillator expectation values","Mapping chaos in a q-deformed harmonic oscillator","q-oscillator dynamics: from periodic to chaotic","Observables reveal chaos in q-deformed quantum system"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation of every time series and every Lyapunov exponent rests on the time-evolved coherent state in Eq. (10); if the phase factor in that expression does not follow from the Hamiltonian (1), the computed regime boundaries are not reliable.","fun_headline_variants_meta":{"raw":{"variants":["q-deformed oscillator: periodic, quasi-periodic, or chaotic","Chaos emerges in q-deformed oscillator expectation values","Mapping chaos in a q-deformed harmonic oscillator","q-oscillator dynamics: from periodic to chaotic","Observables reveal chaos in q-deformed quantum system"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000779,"raw_usage":{"total_tokens":3441,"prompt_tokens":943,"completion_tokens":2498,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":559,"completion_tokens_details":{"reasoning_tokens":2430}},"tokens_in":559,"tokens_out":2498,"duration_ms":19255,"temperature":1.0,"reasoning_tokens":2430,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:13:10.987608+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the time-dependent Schrödinger equation for $H_q$ using the matrix elements in Eq. (7) for $q = 0.9$ and $\\alpha_q = 1$, extract $\\langle X(t)\\rangle_q$, and compare its recurrence plot and largest Lyapunov exponent with those reported here; any significant difference shows the analytic time evolution in Eq. (10) needs correction.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the q-deformed Hamiltonian, its energy spectrum, and the normalization bound on $|\\alpha_q|^2$."},{"cited_title":"Sivakumar, J","cited_arxiv_id":null,"evidence_quote":"Shows that expectation values in the non-deformed oscillator are periodic, providing the baseline comparison."},{"cited_title":"Casati, B","cited_arxiv_id":null,"evidence_quote":"Introduces the expectation-value time-series approach and first-return-time analysis used here."},{"cited_title":"Casati and L","cited_arxiv_id":null,"evidence_quote":"Supplies the exponential first-return-time distribution criterion for ergodic chaotic behaviour."},{"cited_title":"Sudheesh, S","cited_arxiv_id":null,"evidence_quote":"Provides the standard algorithm used to estimate the largest Lyapunov exponents."},{"cited_title":"Marwan, M","cited_arxiv_id":null,"evidence_quote":"Provides an independent algorithm used to cross-check the Lyapunov-exponent results."},{"cited_title":"Rcamier, M","cited_arxiv_id":null,"evidence_quote":"Gives the recurrence-plot formalism and the diagonal-line criteria for classifying periodic, quasi-periodic, and chaotic series."},{"cited_title":"Guha and P","cited_arxiv_id":null,"evidence_quote":"Supplies the forms of the deformed position and momentum operators $X$ and $P$ used in Eqs. (5)-(6)."},{"cited_title":"Batouli, M","cited_arxiv_id":null,"evidence_quote":"Provides the deformed coherent-state and ladder-operator relations used in the time-evolution calculation."}],"review_version":1}