{"id":"01409cab-ac74-411d-9089-918b02c60885","arxiv_id":"2501.14507","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"In a PT-symmetric kicked harmonic oscillator, non-resonant driving yields directed momentum current with ballistic energy growth, while resonant driving yields damped cosine oscillations of momentum and energy with a common frequency.","lead":"This paper simulates a PT-symmetric version of the quantum kicked harmonic oscillator and finds two behaviors: non-resonant kicking produces a steady momentum current and ballistic energy growth, while resonant kicking produces damped oscillations of momentum and energy. A generalist might care because non-Hermitian Floquet systems are a candidate platform for controlling cold atoms and solid-state electron transport.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model is internally inconsistent: the quadratic harmonic potential cannot be periodic with the kick period, so the discrete momentum lattice used in Sec. II and Eq. (9) is not justified.","rationale":"The paper's central claim is a numerical phenomenology for a PT-symmetric kicked harmonic oscillator. The strongest_claim is restated accurately in the reader's verdict. The key condition for the model to be a kicked harmonic oscillator with a discrete momentum lattice is that [θ,p]=iħ_eff with θ on a circle; however, the quadratic potential term η²θ²/2 cannot be periodic on that circle. The text explicitly asserts Hp(X+b)=Hp(X) and uses integer momentum eigenvalues, and those two statements are incompatible for ω≠0. If θ is real, the momentum spectrum is continuous and the basis expansion is invalid; if θ is periodic, the quadratic potential is discontinuous at the boundary. Either way, the model analyzed is not the one described by Eq. (1). This is an internal inconsistency, not a disagreement with an external consensus. The non-resonant mechanism in Eq. (9) relies on nearest-neighbor hopping between discrete momentum states, which is only justified by the compact geometry that makes the quadratic potential ill-defined. Therefore the central claims are not adequately supported. The secondary issues, such as the resonant frequencies being imposed by fitting the same ω_c and the absence of numerical convergence details or code, reinforce the rejection but are not the primary driver. The proposed concrete test would determine whether the reported phenomenology survives a well-defined periodic-potential replacement; until then, the verdict of rejection stands.","tokens_in":11899,"tokens_out":3548,"duration_ms":36353,"concrete_test":"Re-run the split-step simulation with a genuinely periodic, smooth confining potential, e.g., V_per(θ)=η²(1−cosθ), over θ∈[−π,π), for the two parameter sets used in the paper (η=2π/e², λ=1 and η=2π, λ=1). If the linear ⟨p⟩ growth and damped cos(ω_c t) behavior survive with the same scalings and frequencies, the central phenomenology is robust to the model correction; if they change or disappear, the reported dynamics are artifacts of the inconsistent periodic quadratic potential.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is the periodicity of the harmonic potential. In Sec. II the paper states Hp(X+b)=Hp(X) for Hp=Mω²X²/2, which is algebraically false for ω≠0. This periodicity is not a harmless simplification: it is what makes the dimensionless coordinate θ=2πX/b compact and the momentum eigenbasis p|φ_m⟩=mħ_eff|φ_m⟩ discrete. If θ is instead an unbounded real coordinate, the momentum spectrum is continuous and the expansion used in Sec. II and in the effective Hamiltonian Eq. (9) is invalid. If θ is forced to be periodic, the term η²θ²/2 is discontinuous at θ=±π, so the split-step Fourier evolution in Eq. (6) is not the evolution of the stated Hamiltonian. Consequently the numerical results for both the non-resonant directed current (Eqs. 7-8) and the resonant damped oscillations (Eqs. 11-14) are obtained from a model that is not well-defined as written. The mechanism claimed for the current, nearest-neighbor hopping on a momentum lattice, presupposes the same compact geometry that makes the quadratic potential inconsistent.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript numerically studies a PT-symmetric kicked harmonic oscillator defined by Eq. (1). For the non-resonant frequency ratio η=2π/e² it reports linear growth of mean momentum, ⟨p⟩=Gt, and ballistic energy growth, ⟨E⟩≈(1/2)G²t²+C, together with Gaussian wave packets; for the resonant ratio η=2π it reports damped cosine oscillations of momentum and energy with a common frequency ω_c, described by Eqs. (11)-(14). The paper attributes the non-resonant behavior to nearest-neighbor hopping on a discrete momentum lattice, as expressed in Eq. (9), and the resonant behavior to resonant coupling between the non-Hermitian driving and the harmonic oscillator.","tokens_in":12195,"tokens_out":5071,"duration_ms":48884,"significance":"If the model were well defined and the fitted laws were robust, the distinction between irrational and rational frequency ratios in a non-Hermitian kicked oscillator would be a useful addition to the Floquet literature. The paper does provide explicit fitting functions and a qualitative wave-packet characterization, which are in principle falsifiable. However, the central modeling assertion is algebraically incorrect, and the main quantitative results are fitted functions rather than derived predictions. As written, the claimed phenomena are therefore not established, and the paper's central contribution is not sound.","major_comments":[{"comment":"The statement that the harmonic potential is periodic, Hp(X+b)=Hp(X) for Hp=Mω²X²/2, is false for ω≠0. This periodicity is load-bearing: it justifies the compact coordinate θ=2πX/b, the discrete momentum eigenbasis p|φ_m⟩=mħ_eff|φ_m⟩, and the nearest-neighbor hopping picture in Eq. (9). If θ is an unbounded real coordinate, the momentum spectrum is continuous and the expansion used in the paper is invalid; if θ is forced to be periodic, the term η²θ²/2 is discontinuous at θ=±π and the split-step Fourier operator in Eq. (6) does not evolve the stated Hamiltonian. The numerical results in Figs. 1-5 are therefore obtained from a model that is not well defined as written.","section":"Sec. II, after Eq. (1)"},{"comment":"The damped cosine laws are not derived from the Hamiltonian; they are the functions used to fit the numerics, with adjustable saturation values, amplitudes, decay times τ, frequency ω_c, and the time-dependent phase shift tc=t0+D exp(γt). The claim that ⟨p⟩, ⟨E_k⟩, and ⟨E_p⟩ share one frequency is enforced by choosing the same ω_c=4π/15 in the fits rather than demonstrated by the dynamics, and no goodness-of-fit measures or parameter uncertainties are reported. As written, these equations are a phenomenological description and cannot support the conclusion that the system behaves as a damped harmonic oscillator with a fixed frequency.","section":"Sec. III B, Eqs. (11)-(14)"},{"comment":"The linear growth law ⟨p⟩=Gt and the ballistic law ⟨E_k⟩≈(1/2)G²t² rest on a fitted slope G whose values are not tabulated and for which no uncertainty is given. The additional claim that the growth rate G increases to 2π with increasing λ is not supported by any quantitative analysis in the text or by a dedicated figure. Without such support, the non-resonant result is a curve fit rather than a quantitative prediction.","section":"Sec. III A, Eq. (7) and Fig. 1"}],"minor_comments":[{"comment":"The heading 'RESULTS AND DISSCUSION' contains a typo, and 'superlattices' is misspelled as 'supper-lattices' in the Introduction.","section":"Section heading and Introduction"},{"comment":"The density matrix defined as ρ=(1/N)Tr(|ψ⟩⟨ψ|) is not a standard notation; if the intent is normalization after each kick, write ρ=|ψ⟩⟨ψ|/⟨ψ|ψ⟩ and define N accordingly.","section":"Sec. II, after Eq. (6)"},{"comment":"The split-step parameter N is not specified, and no convergence study is reported; since Uω is obtained by a Trotter-type product, the error from finite N should be quantified.","section":"Sec. II, Eq. (6)"},{"comment":"The figure captions give the fitted functions but omit the numerical values of many constants (e.g., t0, D, ps, Eks, Eps, A1, A2) and the fitting ranges, which makes it impossible to reproduce the fits from the text alone.","section":"Sec. III, figure captions"}],"recommendation":"reject","confidential_remarks":"I agree with the stress-test assessment: the periodicity assumption for the harmonic potential is not a harmless simplification but the foundation of the discrete momentum lattice used throughout. Because the model as written is inconsistent and the resonant laws are fitted rather than derived, the manuscript's central claims are not supported. A reformulation of the model and a derivation of the reported laws would be needed, which would constitute a substantially new paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has a clear narrative and asks a natural question, but the central model is not well-defined. In Sec. II the author states Hp(X+b)=Hp(X) for Hp=Mω^2X^2/2, which is simply false for any nonzero ω. That false periodicity is what makes θ compact and the momentum eigenstates discrete; without it, the split-step evolution with η^2θ^2/2 on a periodic grid is not the evolution of the stated Hamiltonian. The stress-test note is right: this is load-bearing, not a cosmetic simplification. The nearest-neighbor hopping picture in Eq. (9) and the directed-current mechanism both presuppose that lattice.\n\nWhat is genuinely new is the damped oscillation of momentum and energy under resonant conditions, but it is characterized only by fitted functions with many adjustable parameters—the same ω_c is imposed for p, Ek, and Ep, so the 'identical frequencies' claim is enforced, not discovered. No error bars, convergence tests, or simulation details are given. The non-resonant directed current and ballistic diffusion were already reported for the PT-symmetric kicked rotor in the cited refs [34,35], and the author acknowledges the similarity.\n\nThe paper is clearly written and the split-step method is standard; I give credit for an honest literature comparison. But the soft spots are not minor. The model needs to be redefined—for example, a genuinely periodic quadratic potential on a circle, with the discontinuity addressed, or a different regularization—and then the oscillation frequency and damping should be derived or at least predicted, not just fitted. As written, the claims are not supported.\n\nThis is a paper for someone studying non-Hermitian Floquet dynamics, but I would not send it to peer review in its current state. A referee would immediately hit the same inconsistency. If the author fixes the model and provides reproducible data, a resubmission could be worthwhile. As is, I would not cite it and would not bring it to the reading group.","headline":"The model is inconsistent as written: the harmonic potential cannot be periodic with the kick, so the discrete momentum lattice and the claimed mechanism are not justified; the non-resonant result is not new, and the resonant part is curve-fitting.","tokens_in":661,"tokens_out":819,"would_cite":false,"duration_ms":27577,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q12"],"pacs":["03.65.-w","05.45.Mt"],"model":"deepseek-v4-flash","headline":"In a PT-symmetric kicked harmonic oscillator, irrational frequency ratios produce a directed momentum current and ballistic energy growth, while rational ratios produce damped oscillations.","keywords":["PT-symmetric","kicked harmonic oscillator","non-Hermitian driving","ballistic diffusion","directed current","damped oscillation","Floquet dynamics","momentum lattice"],"falsifier":"Run the same split-step algorithm with a genuine parabolic potential on an unbounded line (no periodic boundary) and check whether $\\langle p \\rangle = Gt$ still holds; if the linear momentum growth disappears, the directed current is an artifact of the periodic-coordinate assumption.","tokens_in":11612,"feed_emoji":"⚛️","tokens_out":5355,"duration_ms":46566,"temperature":0.7,"pith_summary":"The paper asks what happens to a quantum kicked harmonic oscillator when the kicking potential is made PT-symmetric, with a real cosine plus an imaginary sine component. It claims the answer depends only on whether the ratio of the oscillator frequency to the kicking frequency is rational or irrational. For irrational ratios the system develops a linear growth of mean momentum and quadratic growth of energy, i.e. a directed current and ballistic diffusion; for rational ratios it responds like a damped oscillator, with momentum and energy oscillating as damped cosines of a common frequency. This matters because it connects non-Hermitian driving to transport and energy control in a simple, analytically guided model.","feed_headline":"A PT-symmetric kick turns harmonic motion into ballistic drift","feed_subtitle":"Irrational frequency ratios give a directed momentum current; rational ratios give damped oscillation.","key_machinery":"The central object is the Floquet operator $U = U_\\omega U_K$ acting on momentum eigenstates $|\\phi_m\\rangle$ with eigenvalues $m\\hbar_{\\text{eff}}$. The argument maps the non-resonant problem onto an effective lattice Hamiltonian (Eq. 9) in that momentum basis, where the non-Hermitian term $i\\lambda\\sin\\theta$ becomes asymmetric nearest-neighbour hopping $a_m^\\dagger a_{m+1} - a_{m+1}^\\dagger a_m$; the nonzero bias of this hopping is what produces the constant force and the ballistic current. The numerical propagation uses a split-step decomposition of $U_\\omega$ into alternating momentum and coordinate kicks.","core_discovery":"The central discovery is a dynamical phase dichotomy controlled by the frequency ratio. With an irrational ratio and nonzero imaginary kicking strength, the mean momentum grows as $\\langle p \\rangle = Gt$, the kinetic energy as $\\langle E_k \\rangle \\approx \\tfrac12 G^2 t^2$, the potential energy stays constant, and the total energy grows ballistically; the probability distribution remains a Gaussian moving at constant velocity, so the state behaves as a quasi-classical particle under a constant force. With a rational ratio and sufficiently strong imaginary kicking strength, $\\langle p \\rangle$, $\\langle E_k \\rangle$, $\\langle E_p \\rangle$, and $\\langle E \\rangle$ all oscillate as damped cosines with the same frequency $\\omega_c$, with kinetic and potential energy oscillating in antiphase, so the system acts as an effective damped harmonic oscillator created by resonant coupling between the non-Hermitian driving and the harmonic potential.","pith_inferences":["If the directed current is real, the same PT-symmetric kicking scheme could act as a ratchet for cold atoms in a harmonic trap, with the rational/irrational frequency ratio as a control switch between drift and ringing.","The mapping to nearest-neighbour hopping suggests the ballistic current should survive for other kicking potentials with a nonzero imaginary Fourier component, so the dichotomy may be generic among PT-symmetric Floquet drives.","The paper studies only one irrational and one rational ratio; a systematic scan of the frequency parameter would show whether the damped-oscillation frequency $\\omega_c$ follows a simple function of the rational ratio, a testable prediction the paper leaves implicit.","Because the state is renormalized after each kick, the reported observables describe the normalized dynamics; without that renormalization the norm growth would dominate, and the claimed damped oscillation could be an artifact of the renormalization rather than intrinsic physics."],"forward_implications":["For irrational frequency ratios, the system supplies a constant momentum kick per pulse, producing a directed current whose growth rate $G$ increases with the imaginary kicking strength $\\lambda$ and saturates near $2\\pi$.","For rational frequency ratios, momentum and total energy oscillate at a common frequency $\\omega_c$, with kinetic and potential energy exchanging energy in antiphase, so the system behaves as an effective damped harmonic oscillator.","The damped amplitude, phase shift, and saturation value of the oscillations can be tuned through the imaginary part of the kicking strength $\\lambda$.","The non-resonant regime reproduces, in the harmonic oscillator setting, the ballistic energy growth previously seen in the non-Hermitian kicked rotor, pointing to a common mechanism.","The behaviour is tied to nearest-neighbour hopping on the momentum lattice: the non-Hermitian term biases the hopping and generates the current."],"supporting_citations":[{"why":"Supplies the split-step numerical method for evolving the kicked harmonic oscillator.","marker":"[12]"},{"why":"Establishes directed momentum current and ballistic energy diffusion in the non-Hermitian kicked rotor, the effect the paper extends to the harmonic oscillator.","marker":"[34]"},{"why":"Provides the recent non-Hermitian kicked rotor results that the present work compares against for the non-resonant regime.","marker":"[35]"},{"why":"Provides the non-dimensionalization scheme that maps the physical parameters to the effective variables used in the model.","marker":"[36]"},{"why":"Gives the mapping method used to derive the effective lattice Hamiltonian from the kicked oscillator dynamics.","marker":"[40]"},{"why":"Supports the mapping to an effective bosonic lattice model with nearest-neighbour hopping.","marker":"[41]"}],"fun_headline_variants":["PT-symmetric kicks: irrational ratio gives ballistic drift","Kicked harmonic oscillator: frequency ratio decides drift or damp","Ballistic diffusion or damped oscillation? PT kick decides","PT-symmetric kicked oscillator: drift for irrational, damp for rational","One PT kick, two dynamical fates: ballistic or damped"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper relies on treating the harmonic potential as periodic with the same period as the kicking potential, so the coordinate is effectively taken on a circle; for a genuine unbounded parabolic trap this assumption fails, and the entire momentum-lattice picture rests on it.","fun_headline_variants_meta":{"raw":{"variants":["PT-symmetric kicks: irrational ratio gives ballistic drift","Kicked harmonic oscillator: frequency ratio decides drift or damp","Ballistic diffusion or damped oscillation? PT kick decides","PT-symmetric kicked oscillator: drift for irrational, damp for rational","One PT kick, two dynamical fates: ballistic or damped"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000416,"raw_usage":{"total_tokens":2102,"prompt_tokens":854,"completion_tokens":1248,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":470,"completion_tokens_details":{"reasoning_tokens":1164}},"tokens_in":470,"tokens_out":1248,"duration_ms":10352,"temperature":1.0,"reasoning_tokens":1164,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:05:51.433232+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same split-step algorithm with a genuine parabolic potential on an unbounded line (no periodic boundary) and check whether $\\langle p \\rangle = Gt$ still holds; if the linear momentum growth disappears, the directed current is an artifact of the periodic-coordinate assumption.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the split-step numerical method for evolving the kicked harmonic oscillator."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes directed momentum current and ballistic energy diffusion in the non-Hermitian kicked rotor, the effect the paper extends to the harmonic oscillator."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the recent non-Hermitian kicked rotor results that the present work compares against for the non-resonant regime."},{"cited_title":"Russomanno, M","cited_arxiv_id":null,"evidence_quote":"Gives the mapping method used to derive the effective lattice Hamiltonian from the kicked oscillator dynamics."}],"review_version":1}