{"id":"6a1970e5-08ce-4933-b655-00fa31557b9c","arxiv_id":"2411.10105","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"A time-delayed feedback loop can sustain or suppress parametric autoresonance depending on whether the delay strength exceeds a threshold kth = γω / sin(ωβ̄).","lead":"The paper studies how a time delay in a feedback loop affects a nonlinear oscillator that automatically locks to a slowly changing driving frequency. It finds a critical delay strength above which the oscillation amplitude keeps growing, and tests this with computer simulations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Chirp term in slow-flow Eq. (11) mixes fast and slow times without an ordering for ζ, so the quasi-steady fixed point and the threshold (18) lack a valid asymptotic derivation.","rationale":"The reader's conditional verdict is well placed. I focused on the chirp term in Eq. (11) because it undermines the derivation before the stability analysis even starts; the trace/determinant question is secondary since a quick check at k = kth shows det(J) > 0 for the reported parameters, so the trace-only condition may be numerically adequate there even though it is not rigorously justified. The ζτ0 issue is more fundamental: Eq. (11) mixes the fast time τ0 into a supposedly slow equation without specifying the order of ζ. If ζ is O(ε), the term behaves as (ζ/ε)τ1 and the system is non-autonomous on the slow scale; if ζ is O(ε²), the term should not appear at leading order. Either way, the 'quasi-steady fixed point' in Section II-B is not a fixed point of a consistently derived slow flow, so the amplitude law (13) and threshold (18) lack a valid asymptotic basis. The paper itself acknowledges in Section III and Fig. 6 that the k–β interaction is 'more complex' than predicted, which is consistent with this concern. The paper provides useful numerics and a falsifiable formula, but the analytical derivation needs to be redone with explicit chirp ordering and a full stability analysis; therefore the reader's conditional verdict should stand.","tokens_in":9061,"tokens_out":11549,"duration_ms":112277,"concrete_test":"Independently re-derive the slow-flow phase equation from Eq. (2) with an explicit ordering for the chirp, e.g. ζ = ε^p ζ1 for p = 1, 2, 3, and require that the solvability conditions contain only slow-time variables. Check whether the ζτ0 term in Eq. (11) survives as a leading-order slow variable or is a residual fast-scale term; if the latter, recompute the quasi-steady fixed point and the resulting threshold, and compare with Eq. (18). Additionally, simulate the full DDE (2) for the parameters of Fig. 3 and compare the instantaneous phase difference ψ(t) with the solution of Eqs. (10)–(11); a systematic drift on the fast time scale would confirm the slow-flow reduction is inconsistent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central threshold (18) rests on the slow-flow equations (10)–(11), and Eq. (11) contains the term ζτ0, where τ0 is the fast time of the multiple-scales expansion. For a legitimate slow-flow reduction, the right-hand sides must depend on slow variables only; a term linear in the fast time cannot be balanced by D1η(τ1), an O(1) function of the slow time, unless the chirp rate is given a specific small ordering that the paper never states. The authors declare h, Γ, Λ, g to be O(ε) but say nothing about ζ = μ/ω²; if ζ is O(ε) then ζτ0 = (ζ/ε)τ1 is a non-autonomous forcing of order τ1, and if ζ is O(ε²) the term should be an O(ε) correction, not the leading-order term appearing in Eq. (11). In either case the 'quasi-steady fixed point' used to obtain the amplitude law (13) and the threshold (18) is not a true fixed point of a derived slow-flow system, because ψ keeps changing on the fast scale through ζτ0. Since Fig. 1 and Eq. (13) both derive from this fixed-point solution, the analytical claim that kth = γω/sin(ωβ̄) marks the onset of autoresonance is not backed by a consistent asymptotic derivation. The numerical agreement may be limited to the very small ζ used (ζ = 4×10^-5), and does not by itself validate the method.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a parametrically excited Duffing-type oscillator with constant negative delayed feedback and a linearly chirped parametric drive. Using a two-timescale perturbation expansion, the authors derive a slow-flow system, a quasi-steady amplitude law (Eq. 13), and a critical delay strength kth = γω / sin(ωβ̄) (Eq. 18), claiming that autoresonant amplitude growth occurs only above this threshold. Numerical simulations with exponential envelope fits and parameter scans in k, β̄, and h are presented as corroboration. The paper concludes that time delay can serve as a control mechanism for parametric autoresonance.","tokens_in":9365,"tokens_out":10967,"duration_ms":115554,"significance":"The problem is timely: time-delayed control of autoresonance is much less studied than externally forced autoresonance, and the paper gives an explicit, falsifiable threshold formula that is empirically testable and appears to organize the numerical data in the two main cases examined. The numerical study is more thorough than a single-trajectory check, using envelope fits, maxima-minima diagrams, and two-parameter scans. The main analytical derivation, however, has a serious asymptotic inconsistency involving the fast-time variable in the slow-flow equation, and the stability analysis checks only one half of the necessary linear-stability conditions. If the asymptotic reduction can be repaired and the determinant condition verified, the threshold result would be a useful contribution; in its current form the analytical claim is not rigorously established.","major_comments":[{"comment":"The slow-flow equation for ψ contains the term ζτ0, where τ0 is the fast time variable of the multiple-scales expansion defined in Eq. (4). In a standard two-time reduction, the slow-flow right-hand sides may depend only on A, ψ, and the slow time τ1; a term linear in τ0 cannot be absorbed into a slow-flow fixed point. Indeed, because ψ is defined as η + ζτ0²/2, the derivative of ζτ0²/2 with respect to the independent slow variable τ1 is zero, while if one instead treats τ0 as the physical time, the term ζτ0 produces a non-uniformity on the long time scale. The paper never states the ordering of ζ = μ/ω² relative to ε; if ζ is O(ε), then ζτ0 is not a small uniformly slow correction, and the 'quasi-steady fixed point' D1A = D1ψ = 0 used for Eqs. (12)-(13) and (18) is not a true fixed point of a derived autonomous slow-flow system. A consistent derivation would introduce a slow phase variable that absorbs the chirp before constructing the slow flow. As it stands, the derivation of Eq. (13) and Eq. (18) is not asymptotically justified, so the central threshold claim lacks a rigorous analytical basis, regardless of the numerical agreement.","section":"Section II.A, Eqs. (9)-(11)"},{"comment":"The threshold is derived from Tr(J) = 0 alone, but linear stability of the quasi-steady branch also requires the determinant condition det(J) > 0 (in the notation of Eq. (16)). Using Eq. (10) at the fixed point, the paper's expression (17) for det(J) reduces to det(J) = (3/8)Λ h A² cos(2ψ). Thus det(J) > 0 imposes cos(2ψ) > 0, a condition that is neither derived nor checked in the numerical parameter range. Without this verification, the point k = kth is not shown to be the actual stability boundary of the quasi-steady branch; the correspondence between the sign change of the fitted exponent b and Eq. (18) is therefore empirical rather than established by the linear-stability calculation. In particular, the claimed h-independence of kth follows only from the trace condition and does not address the full stability boundary.","section":"Section II.B, Eqs. (16)-(18)"},{"comment":"The numerical validation uses the sign of the exponent b in an exponential envelope fit x = a e^{bt} as the criterion for autoresonance, but the analytical amplitude law (13) predicts algebraic growth proportional to t^{1/2}, not exponential growth. The fit exponent is therefore not a direct test of the predicted amplitude law, and no quantitative overlay of the numerical envelope on Eq. (13) is provided. Moreover, the text reports early autoresonance onset for k values below the predicted threshold, e.g., 'around k ≈ 0.015' and 'k ≈ 0.005 in the h = 0.002 case' in the Fig. 3 discussion; these pre-threshold events are acknowledged but not reconciled with the claim that kth marks a sharp transition. The paper should either justify the exponential fit as a valid proxy or compare the numerical envelope directly with Eq. (13).","section":"Section III, Figs. 1-4"}],"minor_comments":[{"comment":"The layout of the determinant in Eq. (15) is ambiguous: the bottom-right entry appears to combine -h/2 sin(2ψ) and -s, making it unclear whether the characteristic matrix is J - sI or J + sI. Please present the characteristic equation explicitly with a standard matrix form.","section":"Section II.B, Eq. (15)"},{"comment":"The sentence describing pre-threshold autoresonance, 'there is a set of k values, around k ≈ 0.015, and the value k ≈ 0.005 in the h = 0.002 case', is unclear; please specify which panel and which parameter scan these values refer to.","section":"Section III, Fig. 3 discussion"},{"comment":"The paper uses 'stable autoresonance' to describe amplitude growth above the threshold, whereas the slow-flow fixed point itself would be unstable in the usual Lyapunov sense when the effective damping is negative. Please define what is meant by stability in this context.","section":"Throughout"},{"comment":"Reference [19] duplicates Reference [7]; please remove the duplicate or cite distinct works.","section":"References"},{"comment":"Equation (13) contains the combination (ω² - ω₀²)/ε, which is dimensionally and asymptotically unclear; please explain the ordering and the precise meaning of ε in that expression.","section":"Section II.A, Eq. (13)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a modest but legitimate extension of parametric autoresonance to include time-delayed feedback. The new piece is the threshold formula kth = γω / sin(ωβ̄), which is derived by a standard multiple-scales calculation and then checked numerically. The numerics do show that increasing k past this value changes the asymptotic behavior from decaying to growing, so the control handle is real.\n\nWhat the paper does well: it identifies a gap in the literature (delay effects on parametric autoresonance, as opposed to external autoresonance), produces an analytic expression for the threshold, and then does a systematic numerical study in the k–β̄ plane. The agreement between the predicted and numerical thresholds (0.02 vs 0.0201, 0.0066 vs 0.0068) is genuinely good.\n\nThe soft spots are three. First, the slow-flow equation (11) contains the term ζτ0, where τ0 is the fast time. The paper never states the ordering of ζ relative to the O(ε) parameters. If ζ ~ ε, then ζτ0 = (ζ/ε) τ1, which is a legitimate slowly growing term; but the paper should say so and treat it as such instead of leaving a τ0 in the equation. As written, the derivation looks inconsistent even if it can be fixed. Second, the stability threshold is found from Tr(J)=0 alone, with no check that det(J)>0 on that curve. A Hopf boundary requires both; if det(J) is negative there, the correct transition could be different. This is probably a minor omission, but it should be verified. Third, the numerical results themselves contain patches of autoresonance below threshold (the authors admit this around k ≈ 0.015 and k ≈ 0.005, and in the β̄ scans). That is not fatal, but it directly contradicts the abstract's clean claim, so the abstract and conclusions need to be softened. Also, the Fig. 2 caption has a typo (β̄ = 0.0092 instead of the value used elsewhere), which makes the figure hard to trust until corrected.\n\nThis paper is for an audience that cares about control of autoresonance in mechanical/electrical systems; it doesn't break new mathematical ground but does give a new control parameter. I think it deserves peer review—the derivation issues are fixable and the numerical work is substantial—but it needs a careful revision, not a quick accept.","headline":"A plausible incremental result with a useful threshold formula, but the paper needs to clean up the slow-flow derivation and soften its threshold claims.","tokens_in":9893,"tokens_out":9160,"would_cite":false,"duration_ms":86863,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34C15","34K20","70K28"],"pacs":[],"model":"deepseek-v4-flash","headline":"A delayed feedback switch controls whether parametric autoresonance grows or fades, with a sharp critical delay strength $k_{\\mathrm{th}} = \\gamma\\omega/\\sin(\\omega\\bar{\\beta})$.","keywords":["autoresonance","parametric autoresonance","time-delayed feedback","multiple-scale analysis","stability threshold","chirped forcing","Duffing oscillator","amplitude growth"],"falsifier":"Numerically integrate the original delayed equation with $\\gamma=0.002$, $\\omega=0.5$, $\\alpha=0.03$, $h=0.002$, and $\\bar{\\beta}=0.1$, scanning $k$ through $0.02$: if the fitted envelope exponent $b$ does not change sign near $k\\approx0.0201$, or if amplitude growth appears for $k$ well below that value, the trace-condition threshold is not the actual onset. A second check is to evaluate $\\det(J)$ along the $\\mathrm{Tr}(J)=0$ curve to see whether the fixed point loses stability before the predicted threshold.","tokens_in":8836,"feed_emoji":"📈","tokens_out":6509,"duration_ms":65225,"temperature":0.7,"pith_summary":"The paper establishes that a constant negative time delay in the feedback of a parametrically chirped nonlinear oscillator acts as a control switch for autoresonance. Using multiple-scale perturbation theory, it derives a critical delay strength $k_{\\mathrm{th}} = \\gamma\\omega / \\sin(\\omega\\bar{\\beta})$: for $k>k_{\\mathrm{th}}$ the phase-locked amplitude grows, while for $k<k_{\\mathrm{th}}$ it decays. This matters because time delay is unavoidable in real feedback loops, and the result turns that nuisance into a tunable control parameter for resonance growth. Numerical integration of the original delayed equation confirms the threshold: for $\\bar{\\beta}=0.1$ the predicted $k_{\\mathrm{th}}=0.02$ matches the observed onset near $0.0201$, and for $\\bar{\\beta}=0.3$ the predicted $0.0066$ matches about $0.0068$.","feed_headline":"Delay strength flips autoresonance from decay to growth","feed_subtitle":"A critical threshold $k_{\\mathrm{th}}=\\gamma\\omega/\\sin(\\omega\\bar{\\beta})$ separates growing from fading oscillations; numerics confirm it.","key_machinery":"The central object is the slow-flow system obtained by a two-time-scale expansion of the delayed parametric oscillator: two first-order equations, $D_1 A = F_1(A,\\psi)$ and $D_1\\psi = F_2(A,\\psi)$, for the amplitude $A$ and the phase difference $\\psi$. The quasi-steady fixed point of these equations gives an analytical amplitude formula, and the stability of that fixed point is analyzed through the Jacobian matrix $J$; the condition $\\mathrm{Tr}(J)=0$, combined with the fixed-point condition $D_1 A=0$, produces the critical delay strength $k_{\\mathrm{th}} = \\gamma\\omega / \\sin(\\omega\\bar{\\beta})$. This trace condition is the load-bearing step that turns the delay parameters $(k,\\bar{\\beta})$ into a predicted control threshold.","core_discovery":"The paper argues that for a parametrically chirped Duffing-type oscillator with negative delayed feedback, the slow-flow amplitude and phase equations admit a quasi-steady fixed point whose stability is controlled by the delay strength. Setting the trace of the linearized Jacobian to zero at this fixed point yields the critical value $k_{\\mathrm{th}} = \\gamma\\omega / \\sin(\\omega\\bar{\\beta})$, below which the fixed point is unstable and the oscillation envelope decays, and above which autoresonant amplitude growth is sustained. The threshold is independent of the drive amplitude $h$ and of the nonlinearity coefficient $\\alpha$, and it implies an inverse relation $\\bar{\\beta}_{\\mathrm{th}} = (1/\\omega)\\arcsin(\\gamma\\omega/k)$ between delay time and delay strength. Numerical simulations confirm the predicted onset: the fitted envelope exponent $b$ changes from negative to positive at $k \\approx 0.0201$ for $\\bar{\\beta}=0.1$ and at $k \\approx 0.0068$ for $\\bar{\\beta}=0.3$, in good agreement with the analytical formula.","pith_inferences":["Beyond the paper's own results, the stability threshold is derived from the trace condition alone; verifying whether $\\det(J)>0$ also holds along the $\\mathrm{Tr}(J)=0$ boundary could reveal a stricter onset or a Hopf bifurcation that the single-threshold formula does not capture.","The slow-flow phase equation contains the term $\\zeta\\tau_0$, which mixes the fast time $\\tau_0$ into the slow dynamics; if this step is not rigorously justified, the amplitude formula and threshold may need modification at larger chirp rates.","The striped regions in the $(k,\\bar{\\beta})$ plane where autoresonance appears before the analytically predicted threshold suggest additional resonance tongues; mapping those tongues is a natural extension that could test the limits of the trace-condition prediction.","Because $k_{\\mathrm{th}}$ depends only on damping, frequency, and delay time, one could design a gain-scheduled delay controller that switches resonance growth on and off; testing this at several damping values would directly probe the formula's validity."],"forward_implications":["If the threshold claim is correct, the delay strength $k$ can serve as an experimentally adjustable control for autoresonance, complementing the usual control through the chirp rate $\\mu$.","The threshold is independent of the drive amplitude $h$, so the same onset value $k_{\\mathrm{th}}$ should be observed for different forcing strengths, as the numerical results in Fig. 3 indicate.","For a fixed delay strength, the required time delay is $\\bar{\\beta}_{\\mathrm{th}} = (1/\\omega)\\arcsin(\\gamma\\omega/k)$, giving a reciprocal design relation between delay time and delay strength.","The sign of the exponential envelope fit parameter $b$ provides a direct observable criterion for the onset of autoresonance: positive $b$ means growing amplitude, negative $b$ means decaying amplitude.","In mechanical and electrical systems where sensing and actuation introduce lag, the delay can be tuned to sustain or suppress resonance rather than being an uncontrolled disturbance."],"supporting_citations":[{"why":"Supplies the baseline parametric autoresonance model that this paper extends to include time-delayed feedback.","marker":"[32]"},{"why":"Provides a physical realization of parametric autoresonance in Faraday waves, a setting where delay control could be applied.","marker":"[33]"},{"why":"Supplies the capture-into-parametric-autoresonance mechanism that the phase-locking assumption relies on.","marker":"[34]"},{"why":"Extends parametric autoresonance to nonlinear wave equations, framing the broader model class.","marker":"[35]"},{"why":"Gives the time-delayed state-feedback resonance framework for Duffing oscillators whose perturbation methods are adapted here.","marker":"[36]"},{"why":"Provides theory and numerics for time-delayed feedback in Duffing oscillators, supporting the slow-flow stability analysis.","marker":"[38]"},{"why":"Documents delay-induced resonance in the time-delayed Duffing oscillator, motivating the use of delay as a control parameter.","marker":"[39]"},{"why":"The most direct antecedent, studying time delay in an autoresonant chain and providing the comparison point for this parametric single-oscillator extension.","marker":"[42]"}],"fun_headline_variants":["Critical delay strength turns decay into growth","Time delay threshold governs sustained amplitude","Delay feedback: a threshold for amplitude growth","Autoresonance growth requires sufficient delay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The threshold formula rests on treating the slow-flow phase equation, which contains the fast time $\\tau_0$ inside the term $\\zeta\\tau_0$, as a legitimate slow system, and on reading the stability boundary from the trace condition alone.","fun_headline_variants_meta":{"raw":{"variants":["Critical delay strength turns decay into growth","Time delay threshold governs sustained amplitude","Delay feedback: a threshold for amplitude growth","Autoresonance growth requires sufficient delay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000274,"raw_usage":{"total_tokens":1611,"prompt_tokens":889,"completion_tokens":722,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":670}},"tokens_in":505,"tokens_out":722,"duration_ms":8372,"temperature":1.0,"reasoning_tokens":670,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:58:30.317024+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the original delayed equation with $\\gamma=0.002$, $\\omega=0.5$, $\\alpha=0.03$, $h=0.002$, and $\\bar{\\beta}=0.1$, scanning $k$ through $0.02$: if the fitted envelope exponent $b$ does not change sign near $k\\approx0.0201$, or if amplitude growth appears for $k$ well below that value, the trace-condition threshold is not the actual onset. A second check is to evaluate $\\det(J)$ along the $\\mathrm{Tr}(J)=0$ curve to see whether the fixed point loses stability before the predicted threshold.","supporting_citations":[{"cited_title":"Khain and B","cited_arxiv_id":null,"evidence_quote":"Supplies the baseline parametric autoresonance model that this paper extends to include time-delayed feedback."},{"cited_title":"Assaf and B","cited_arxiv_id":null,"evidence_quote":"Provides a physical realization of parametric autoresonance in Faraday waves, a setting where delay control could be applied."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the capture-into-parametric-autoresonance mechanism that the phase-locking assumption relies on."},{"cited_title":"Friedland and A","cited_arxiv_id":null,"evidence_quote":"Extends parametric autoresonance to nonlinear wave equations, framing the broader model class."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the time-delayed state-feedback resonance framework for Duffing oscillators whose perturbation methods are adapted here."},{"cited_title":"Jeevarathinam, S","cited_arxiv_id":null,"evidence_quote":"Provides theory and numerics for time-delayed feedback in Duffing oscillators, supporting the slow-flow stability analysis."},{"cited_title":"Cantis´ an, M","cited_arxiv_id":null,"evidence_quote":"Documents delay-induced resonance in the time-delayed Duffing oscillator, motivating the use of delay as a control parameter."},{"cited_title":"Chac´ on, F","cited_arxiv_id":null,"evidence_quote":"The most direct antecedent, studying time delay in an autoresonant chain and providing the comparison point for this parametric single-oscillator extension."}],"review_version":1}