{"id":"3a1493e2-447a-4fa9-bea0-c733dbfa6762","arxiv_id":"2412.07547","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A numerical study of a delayed Duffing driver and a periodically forced Duffing response identifies a resonance at response frequency about 2, caused by the interplay of the coupling and the external forcing.","lead":"This paper simulates a delayed Duffing oscillator driving a non-delayed Duffing oscillator that is also forced by a periodic signal, and reports a new resonance near frequency 2, which it calls coupling-forcing resonance. The interest is in distinguishing this effect from two previously known coupling-induced resonance effects in the same system.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed coupling-forcing resonance may be an amplitude-dependent nonlinear resonance of the response Duffing oscillator whose restoring force is statically modified by the coupling term, not a genuinely new interaction mechanism.","rationale":"The reader's weakest-assumption diagnosis is essentially correct: the comparison C = 0 versus C = 1.66 does not control for the coupling-induced modification of the response oscillator's effective restoring force. I partially disagree with the wording of the reader's concern, because the relevant shift is not simply the linear natural frequency sqrt(1 + C) = sqrt(2.66) approximately 1.63; the observed peak near omega2 = 2 is better explained by the amplitude-dependent backbone curve of the Duffing oscillator, omega^2 = (1 + C) + (3/4)A^2, which for C = 1.66 places resonance at omega2 = 2 for A approximately 1.3. This sharpens rather than removes the concern. The paper's attribution of the new resonance to an interaction between the driver system and the response forcing would be refuted if a surrogate with a constant driver reproduces the Fig. 5(b) peak; it would be supported if the surrogate fails to do so. The paper has no quantitative resonance criterion beyond visual comparison and no error or convergence analysis, but those are secondary. The central observation of high amplitudes near omega2 = 2 is plausible; the interpretation as a distinct mechanism is not fully established. Therefore the conditional verdict is appropriate, and no change to the reader's verdict is needed.","tokens_in":7179,"tokens_out":4427,"duration_ms":46013,"concrete_test":"For the tau = 1 case, replace the driver by its steady-state fixed point X (solving -gamma X - alpha X + alpha X^3 = 0, which gives X = sqrt(1.5) for the branch reached from the stated history), and integrate the surrogate equation x2'' + mu x2' - x2 + x2^3 = C(X - x2) + f cos(omega2 t) with C = 1.66, mu = 0.01 over the same f-omega2 grid used in Fig. 5(b). If this constant-driver surrogate reproduces the high-amplitude region localized near omega2 = 2, the claimed interaction mechanism is not needed to explain the observations. If the peak shifts or disappears, the interaction interpretation is supported. As a robustness check, overlay the numerically computed autonomous backbone curve, obtained by setting f = 0 and recording the response amplitude A as a function of omega2, on the surrogate amplitude map.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the high-amplitude band near omega2 = 2 arises from a dynamical interaction between the driver forcing C(x1 - x2) and the response forcing f cos(omega2 t), rather than from a parameter shift of the response oscillator's own resonance. Rewriting Eq. (2) as x2'' + mu x2' + (-1 - C)x2 + x2^3 = C x1 + f cos(omega2 t) shows that the coupling contributes both a forcing term, C x1, and a linear stiffness modification, -C x2. For tau = 1, the paper states that the driver converges to a fixed point, so C x1 is a static constant. The remaining equation is a standard damped Duffing oscillator with linear coefficient -(1 + C) and cubic coefficient +1, whose backbone relation is omega^2 = (1 + C) + (3/4)A^2 for amplitude A. With C = 1.66, this predicts a resonant band near omega2 = 2 for A approximately 1.3, close to the observed oscillation amplitudes. Thus the comparison between C = 0 and C = 1.66 does not isolate a new interaction mechanism; it conflates the coupling with a change in the response oscillator's nonlinear resonance frequency. The manuscript never computes the autonomous frequency-amplitude curve of Eq. (2) with x1 replaced by its steady-state value, nor does it vary C to check whether the resonance frequency tracks the predicted backbone rather than remaining pinned at 2.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies a unidirectionally coupled pair of Duffing oscillators in which a time-delayed Duffing driver acts on a non-delayed Duffing response through the coupling term C(x1 - x2), while the response is also driven by an external periodic force f cos(omega2 t). The authors compare the C = 0 and C = 1.66 cases, observe a high-amplitude band near omega2 approximately 2 in the response, and attribute it to a new 'coupling-forcing resonance' arising from the interaction between the coupling mechanism and the response forcing. They further claim that this resonance can coexist with the transmitted resonance and the coupling-induced resonance in the same parameter set, and they present f-omega2 and C-f scans for tau = 1 and tau = 2.","tokens_in":7504,"tokens_out":5108,"duration_ms":49367,"significance":"If the interpretation is correct, the paper would add a new member to the family of coupling-related resonances and would provide a useful mapping of where the three resonance phenomena overlap. The numerical approach is standard and includes appropriate exploratory tools: adaptive-step DDE integration, C = 0 control runs, FFT-based frequency panels, and parameter scans over f-omega2 and C-f. These are strengths of the paper. However, the central discriminating test, namely whether the peak near omega2 = 2 is a genuine dynamical interaction or simply the ordinary resonance of a Duffing oscillator whose stiffness has been shifted by the C term, is not performed. The current evidence does not exclude the latter, so the main claim is not yet established.","major_comments":[{"comment":"Rewriting Eq. (2) as x2'' + mu x2' + (-1 - C) x2 + x2^3 = C x1 + f cos(omega2 t) shows that the coupling strength C enters the response equation both as a forcing term C x1 and as a linear stiffness modification -C x2. For tau = 1 the driver converges to a fixed point, so C x1 is a static constant and the response is a standard damped Duffing oscillator with linear coefficient -(1 + C). Its backbone relation is omega^2 = (1 + C) + (3/4) A^2; with C = 1.66 and A approximately 1.3 this predicts resonance near omega2 = 2, exactly the band reported in Figs. 3 and 5(b). The C = 0 versus C = 1.66 comparison therefore does not isolate an interaction mechanism; it conflates the putative new resonance with a parameter shift of the response oscillator's own nonlinear resonance. The authors should compute the autonomous frequency-amplitude curve of Eq. (2) with x1 at its steady-state value, and/or vary C while scanning omega2, to determine whether the observed peak tracks the backbone rather than remaining pinned at omega2 = 2.","section":"II, Eq. (2), and III A"},{"comment":"The identification of the resonance is based on visual inspection of amplitude maps: no quantitative threshold defines 'high oscillation amplitudes,' no error bars or convergence checks are reported, and the color-scale conventions are not specified. Because the central claim is that the coupling-forcing resonance is localized at a specific frequency and is absent when either perturbation is removed, the paper should provide an objective criterion, for example an amplitude exceeding the C = 0 baseline by a stated factor, and should demonstrate that the high-amplitude band is robust to integration tolerances and to variations of initial conditions.","section":"III A, Figs. 3 and 5"},{"comment":"The C-f scans in Fig. 7 are performed at two fixed omega2 values and do not include the omega2-C plane, so they cannot test whether the resonance frequency varies with C as predicted by the backbone relation discussed in the first comment. The statement that C = 1.66 'consistently corresponds to elevated oscillation amplitudes' is not quantified; in the displayed panels the black line appears to pass through both high- and low-amplitude regions. The authors should report a quantitative measure, such as the amplitude at the forced resonance as a function of C, or a thresholded region in the C-omega2 plane.","section":"III B, Fig. 7"}],"minor_comments":[{"comment":"There is a typo in the sentence introducing the coupling mechanism: 'the the positions' should read 'the positions.'","section":"II"},{"comment":"The sentence 'in Fig. 2, we show the oscillation amplitudes of the Duffing oscillator without coupling for omega2 = 0.5' is inconsistent with Fig. 2 being an f-omega2 map; either the figure or the sentence should be corrected.","section":"III A"},{"comment":"The notation omega2C for the FFT-based oscillation frequency is undefined and could be misread as a product; please define it explicitly.","section":"Figures 5(e) and 5(f)"},{"comment":"The Conclusions state that the resonance appears at 'omega2 = 2' without qualification, while the figures show a band of high amplitudes; the statement should be phrased as approximate.","section":"IV"},{"comment":"The manuscript does not report the numerical tolerance used in Matlab's ddesd; adding this would strengthen reproducibility.","section":"II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a natural continuation of the authors' earlier work [24,27], but the novelty hinges on excluding the simplest alternative explanation, namely that the coupling term changes the linear stiffness of the response oscillator and thereby shifts its ordinary Duffing resonance. In my view, the paper requires the additional backbone and parameter-shift analysis before the central claim can be evaluated. If that analysis shows the peak tracks the backbone, the paper would need to be substantially reframed; if it shows the opposite, the paper would be publishable after adding the quantitative resonance criterion requested in the major comments."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: the paper is clearly written and the numerical observation is real, but the central interpretation is not established. The new 'coupling-forcing resonance' near omega2 = 2 can be explained, at least for the tau = 1 case, as the ordinary amplitude-dependent resonance of the response oscillator after the coupling term has changed its potential. Their control comparison, C = 0 versus C = 1.66, changes both the forcing and the stiffness; it does not isolate an interaction.\n\nThe math behind this: for tau = 1 the driver converges to a fixed point, so C x1 is static. Equation (2) becomes x2'' + mu x2' - (1+C)x2 + x2^3 = const + f cos(omega2 t). That is a standard driven Duffing oscillator with a C-dependent equilibrium and natural frequency. The peak near omega2 = 2 is then plausibly the intrawell resonance of this modified oscillator, not a new mechanism. The paper never computes the backbone curve, nor does it vary C to see whether the peak tracks the predicted frequency. Fig. 7 varies C only at two fixed drive frequencies, which does not settle the question. The stress-test note's specific backbone formula has a sign error, but the structural point holds.\n\nWhat is genuinely good: the paper is honest about the difficulty of separating the three resonance regimes, the FFT panels are a sensible diagnostic, and the simulations are standard. For this model, the observation at omega2 = 2 is new. The authors also explicitly flag that mixed cases are hard to disambiguate, which is commendable.\n\nThe main soft spot is the missing control for the parameter-shift explanation; this is load-bearing. A secondary issue is that the resonance criterion is purely visual--no quantitative peak definition, no convergence checks. Those are minor in comparison.\n\nWho gets value from this? Specialists in delayed and coupled oscillators, and anyone interested in how parameter changes can masquerade as new mechanisms. It deserves a serious referee, but the referee should require a backbone calculation and a control with x1 held at its steady-state value. Without that, I would not cite 'coupling-forcing resonance' as a distinct phenomenon.\n\nMy recommendation: send it to peer review, conditional on the authors addressing the natural-frequency-shift alternative. The fix is well-defined and the paper is otherwise serviceable.\n\nBest","headline":"The claimed coupling-forcing resonance is likely a shifted natural-frequency resonance of the response oscillator, not a distinct interaction mechanism as presented.","tokens_in":8008,"tokens_out":7432,"would_cite":false,"duration_ms":75600,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34C15"],"pacs":[],"model":"deepseek-v4-flash","headline":"A new resonance, the coupling-forcing resonance, emerges from the interaction between the time-delayed driver's coupling and the response system's periodic forcing, localized near forcing frequency $\\omega_2 \\approx 2$.","keywords":["Duffing oscillator","time-delayed system","coupling-forcing resonance","transmitted resonance","coupling-induced resonance","continuous control coupling","periodic forcing","nonlinear resonance"],"falsifier":"Compute the response oscillator's linearized natural frequency from the forced equation as a function of the coupling constant $C$: with $\\alpha=-1$, the effective restoring coefficient is $C+\\alpha$, so the linear frequency candidate is $\\sqrt{C-1}$. If, as $C$ is varied, the high-amplitude peak in the forcing-frequency plane tracks $\\sqrt{C-1}$ instead of remaining pinned at $\\omega_2\\simeq2$, the coupling-forcing resonance is ordinary resonance of a coupling-shifted oscillator rather than a distinct interaction effect.","tokens_in":6993,"feed_emoji":"📈","tokens_out":16020,"duration_ms":134073,"temperature":0.7,"pith_summary":"The paper claims that a non-delayed Duffing oscillator (a nonlinear oscillator with a cubic restoring term), driven both by a time-delayed Duffing oscillator through a unidirectional coupling and by its own periodic forcing, develops a distinct resonance, the coupling-forcing resonance, when the two perturbations act together. The resonance is localized near a response-forcing frequency of $\\omega_2 \\approx 2$, and it appears at coupling strengths and forcing amplitudes for which neither system alone shows large oscillations. This matters because it gives a third, separable route to large-amplitude oscillations in coupled oscillator systems, alongside the transmitted resonance and the coupling-induced resonance, and it tells a researcher where to look in frequency and coupling strength to observe or avoid such amplification. The paper also shows that the three resonance mechanisms can coexist and overlap in the same parameter set, which complicates attributing an observed amplitude increase to any one mechanism.","feed_headline":"A coupling-forcing resonance appears at frequency 2","feed_subtitle":"A delayed driver plus a periodic response forcing creates large oscillations neither can produce alone.","key_machinery":"The load-bearing object is the two-oscillator model: a time-delayed Duffing oscillator (a nonlinear oscillator with a cubic restoring term) as the driver, and a non-delayed Duffing oscillator as the response, coupled unidirectionally through the continuous-control term $C(x_1-x_2)$, with the response also forced by $f\\cos\\omega_2 t$. The mechanism under study is the interaction between that coupling term and the response's external forcing; the paper isolates it by comparing $C=0$ with $C=1.66$ and by checking the driver's amplitude to rule out transmission. The named phenomenon, the coupling-forcing resonance, is identified as the high-amplitude region around $\\omega_2\\simeq2$ that appears only when both the coupling and the forcing are active.","core_discovery":"The central discovery is that the coupling term $C(x_1-x_2)$ and the external forcing $f\\cos\\omega_2 t$ on the response system cooperate to produce a resonance that is absent when either perturbation acts alone. For $C=1.66$ and with the driver unforced, the response oscillator shows a region of high-amplitude oscillations around $\\omega_2\\simeq2$ for sufficiently large $f$, whereas the $C=0$ response does not and the driver's own amplitude shows no corresponding peak. The paper names this phenomenon the coupling-forcing resonance and demonstrates in the $F$--$\\omega_2$ plane that it can coexist with the transmitted resonance (when the driver is forced by its own periodic signal) and with the coupling-induced resonance (when the response forcing is absent), sometimes overlapping so strongly that the three cannot be cleanly separated. In the coupling-constant plane, the high-amplitude region is more pronounced at $\\omega_2=2$ than at $\\omega_2=0.5$, and $C=1.66$ is shown to lie in a band of elevated amplitudes across the two delay regions studied.","pith_inferences":["Inference: the same interaction mechanism should appear in other driver-response pairs, such as a delayed linear oscillator driving a different nonlinear oscillator, whenever the coupling signal and the response forcing share a frequency near the response's effective natural frequency.","Inference: the overlap of the three resonances could be disentangled by phase-resolved or frequency-resolved measurements, because the coupling-forcing resonance should leave a spectral signature in the response that is absent in the driver, whereas the transmitted resonance should carry the driver's forcing frequency.","Inference: for applications where large oscillations are harmful, these results suggest avoiding response forcing frequencies near $\\omega_2\\simeq2$ at moderate coupling strengths; for energy harvesting or weak-signal detection, operating there deliberately would exploit the amplification."],"forward_implications":["A measured large-amplitude peak near $\\omega_2\\simeq2$ in this coupled system is a signature of the coupling-forcing resonance, not of the driver's own dynamics or the uncoupled response.","Tuning the coupling constant $C$ and the response forcing amplitude $f$ gives practical control of the amplification: $C=1.66$ sits in a high-amplitude band, and the effect is stronger at $\\omega_2=2$ than at $\\omega_2=0.5$.","Because the resonance appears for delay values in two different driver regimes, the phenomenon persists whether the driver settles to a fixed point or a limit cycle.","In parameter regions where transmitted and coupling-induced resonances also occur, attributing an amplitude increase to one mechanism requires checking the driver's amplitude and the response's frequency content, since the three can overlap."],"supporting_citations":[{"why":"Supplies the baseline coupling-induced resonance and fixes the working value C=1.66 as a coupling strength that already enhances the response without external forcing.","marker":"[24]"},{"why":"Defines the transmitted resonance that the new phenomenon is compared against and shown to coexist with.","marker":"[27]"},{"why":"Establishes the delay-induced bifurcation regions of the time-delayed Duffing oscillator used to choose the two delay values studied.","marker":"[21]"},{"why":"Documents delay-induced resonance in the delayed Duffing oscillator, the background for the driver's oscillation regimes.","marker":"[22]"},{"why":"Introduces the continuous-control coupling approach underlying the C(x1-x2) term.","marker":"[25]"},{"why":"Provides the synchronization-by-continuous-control formulation used in the unidirectional coupling.","marker":"[26]"}],"fun_headline_variants":["At frequency 2, coupling and forcing spark a resonance","New resonance when delayed driver and periodic drive combine","Coupled driver plus periodic forcing yields resonance at ω=2","Two drives together create a resonance neither can produce alone"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the high-amplitude region near $\\omega_2=2$ is a genuinely new interaction effect, rather than the coupling term simply shifting the response oscillator's natural frequency so that the forcing at $\\omega_2=2$ becomes ordinary resonance; the paper never checks that alternative.","fun_headline_variants_meta":{"raw":{"variants":["At frequency 2, coupling and forcing spark a resonance","New resonance when delayed driver and periodic drive combine","Coupled driver plus periodic forcing yields resonance at ω=2","Two drives together create a resonance neither can produce alone"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000863,"raw_usage":{"total_tokens":3726,"prompt_tokens":911,"completion_tokens":2815,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":527,"completion_tokens_details":{"reasoning_tokens":2750}},"tokens_in":527,"tokens_out":2815,"duration_ms":18956,"temperature":1.0,"reasoning_tokens":2750,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:43:36.613363+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the response oscillator's linearized natural frequency from the forced equation as a function of the coupling constant $C$: with $\\alpha=-1$, the effective restoring coefficient is $C+\\alpha$, so the linear frequency candidate is $\\sqrt{C-1}$. If, as $C$ is varied, the high-amplitude peak in the forcing-frequency plane tracks $\\sqrt{C-1}$ instead of remaining pinned at $\\omega_2\\simeq2$, the coupling-forcing resonance is ordinary resonance of a coupling-shifted oscillator rather than a distinct interaction effect.","supporting_citations":[{"cited_title":"Vibrational resonance","cited_arxiv_id":null,"evidence_quote":"Supplies the baseline coupling-induced resonance and fixes the working value C=1.66 as a coupling strength that already enhances the response without external forcing."},{"cited_title":"Bogdanov–Takens resonance in time- 15 delayed systems","cited_arxiv_id":null,"evidence_quote":"Defines the transmitted resonance that the new phenomenon is compared against and shown to coexist with."},{"cited_title":"Stochastic resonance","cited_arxiv_id":null,"evidence_quote":"Establishes the delay-induced bifurcation regions of the time-delayed Duffing oscillator used to choose the two delay values studied."},{"cited_title":"Stochastic Resonance","cited_arxiv_id":null,"evidence_quote":"Documents delay-induced resonance in the delayed Duffing oscillator, the background for the driver's oscillation regimes."},{"cited_title":"Delay-induced resonance in the time- delayed duffing oscillator","cited_arxiv_id":null,"evidence_quote":"Introduces the continuous-control coupling approach underlying the C(x1-x2) term."},{"cited_title":"Delay-induced resonance sup- presses damping-induced unpredictability","cited_arxiv_id":null,"evidence_quote":"Provides the synchronization-by-continuous-control formulation used in the unidirectional coupling."}],"review_version":1}