{"id":"8ac78e82-79c9-4cac-8813-07bda566a1b0","arxiv_id":"2501.11889","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A driven asymmetric double-well oscillator shows rare bursts of irregular mixed-mode oscillations, classified as extreme events, with velocity spikes as a proposed precursor.","lead":"Researchers report a new kind of rare event in a driven Helmholtz-Duffing oscillator: bursts of irregular mixed-mode oscillations caused by rare hopping between two asymmetric potential wells. The study uses simulations to classify these bursts and proposes the oscillator velocity as a warning sign before they occur.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fixed-step RKF45 at dt=0.01 is never validated; in a double-well system near crisis, numerical noise could manufacture the rare barrier crossings that define the reported bursts.","rationale":"Good-faith reading: the paper's central contribution is the discovery of rare irregular MMO bursts in the HD oscillator and their characterization; the strongest claim also includes a velocity lead indicator. I initially considered attacking the velocity claim directly, since y=dx/dt and a large x excursion must be preceded by an increase in y; a quantitative predictor with false-positive rates would be needed. That concern is real and the paper's Fig. 7 is only visual. However, the more fundamental load-bearing assumption is that the rare bursts are physical. All statistical and mechanistic conclusions inherit from the numerical trajectories. In a chaotic double-well system near a crisis, fixed-step integration without error control can add effective stochastic forcing; because the events are rare, even a small spurious escape rate can dominate. The presence of inconsistent parameter values in captions reinforces the need for reproducible numerical validation. The reader's weakest assumption already flagged fixed-step RKF45, so agreement is 'agree.' The requested check is cheap and would resolve the concern; until it is run, CONDITIONAL remains the appropriate verdict.","tokens_in":12108,"tokens_out":6312,"duration_ms":70853,"concrete_test":"Recompute the central case at c2=91.5464884 with an adaptive Dormand-Prince integrator (tol=1e-8, 1e-10, 1e-12) and with fixed steps h=0.005 and h=0.002 over an identical long interval; record burst count, inter-burst intervals, and threshold exceedance rate. If the irregular MMO bursts disappear or their rate changes by more than a factor of 2 between settings, the reported extreme events are numerical artifacts. Optionally, add explicit deterministic noise of amplitude comparable to the local truncation error to see whether barrier crossings are noise-triggered.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The numerical integration of Eq. (2) uses RKF45 with a fixed step h=0.01 and no error control or convergence test. The claimed phenomenon is precisely a rare transition between the two wells across the unstable fixed point at (0,0), occurring only in a narrow window just below the crisis at c2=91.5464885. Close to a crisis, the chaotic attractor is extremely sensitive, and deterministic integration error acts like a small noise that can induce spurious escapes over the barrier; no evidence is given that the bursts persist under tighter integration. The quoted positive largest Lyapunov exponent 0.0396177 is not backed by the method, integration time, or convergence. The manuscript also contains inconsistent c2 values in figure captions (e.g., Fig. 4 uses c2=5.6494884 whereas the text uses c2=91.5464884) and provides no code or data, so it is not possible to confirm that the numerics explored the intended parameter set. Since every downstream conclusion - POT classification, GEV fit, mechanism, and the velocity predictor - is based on these trajectories, an unvalidated fixed-step integrator is the single most load-bearing uncertainty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a new class of extreme events in the driven Helmholtz-Duffing oscillator, termed extreme irregular mixed-mode oscillatory bursts. The authors argue that, for a narrow interval of the quadratic nonlinearity coefficient c2 near a crisis, the trajectory intermittently escapes the deeper potential well and performs irregular back-and-forth hopping between the two wells, producing bursts that are classified as extreme events by a peak-over-threshold (POT) criterion. They further report that the threshold exceedances are well described by a generalized extreme value (GEV) distribution with a negative shape parameter, that the bursts occur during the rising phase of the external drive, and that the velocity variable provides a reliable lead indicator of an imminent burst. The paper is primarily numerical: it integrates the equation of motion with a fixed-step RKF45 scheme and analyzes the resulting time series, phase portraits, and histograms.","tokens_in":12406,"tokens_out":4107,"duration_ms":44050,"significance":"If the numerical results are robust, the paper identifies a phenomenologically distinct type of extreme event: rather than isolated large spikes, the extreme event is an entire burst of irregular mixed-mode oscillations, and the mechanism is rare chaotic inter-well hopping in an asymmetric double-well potential. The observation that the oscillator returns to periodic mixed-mode oscillations inside the deeper well between bursts is a notable qualitative feature that distinguishes this system from previously reported extreme-event mechanisms. The paper also makes a concrete, testable claim that velocity can serve as a lead indicator. However, the significance is currently limited by the absence of numerical convergence tests, inconsistent reporting of the central parameter c2 in figure captions, an arbitrary and unvalidated threshold, and a statistically under-supported GEV fit. The velocity-predictor claim is supported only by a visual inspection of a single time series. These issues are load-bearing because every conclusion in the paper derives from the same unverified trajectories.","major_comments":[{"comment":"The fixed-step RKF45 integration with h=0.01 is not accompanied by any error-control or convergence test. Since the claimed extreme events are rare crossings of the unstable fixed point at (0,0) in a narrow parameter window just below crisis, deterministic integration error can act as a small perturbation that may induce or suppress these crossings. Please provide a convergence study (e.g., halving h, comparing with an adaptive solver, or reporting event counts as a function of tolerance) and state the total integration time and transient discarded. Without such evidence, the central observation could be a numerical artifact.","section":"Section II, Eq. (2)"},{"comment":"The parameter value at which extreme events occur is reported inconsistently. The text states that extreme events appear at c2=91.5464884, but the captions of Fig. 4 and Fig. 5 state c2=5.6494884, and the caption of Fig. 3 states c2=91.564884. Because the claimed phenomenon occurs in an extremely narrow c2 interval of width 1e-7, this inconsistency prevents the reader from knowing which trajectories were actually computed. Please correct all captions and verify that the displayed data correspond to the stated parameter value.","section":"Sections II and III, Figs. 2, 4, 5"},{"comment":"The POT threshold with n=5 is arbitrary, and no sensitivity analysis is provided; the number of exceedances and all subsequent statistics depend on this choice. Moreover, the GEV density in Eq. (4) is fitted to peak-over-threshold exceedances, but the standard asymptotic model for POT exceedances is the generalized Pareto distribution, not the GEV, and no goodness-of-fit test (e.g., Kolmogorov-Smirnov, Anderson-Darling, or QQ plot) is reported. Please justify the threshold choice, scan over n, and test the distributional fit; otherwise the claim of a Weibull-type distribution in Table I is unsupported.","section":"Section II, Eq. (3) and Fig. 3, Table I"},{"comment":"The velocity lead indicator is identified by visual inspection of the same time series in which the bursts occur. No quantitative rule is given for what constitutes a velocity 'shoot-up,' no prediction horizon is defined, and no false-alarm or miss rates are reported. The claim that velocity is a reliable lead indicator therefore requires a defined detection algorithm and validation on held-out data or on multiple independent realizations; the current evidence is anecdotal.","section":"Section IV, Fig. 7"},{"comment":"The largest Lyapunov exponent of 0.0396177 is quoted without describing the algorithm, integration time, or convergence. Given the extreme sensitivity of the dynamics near crisis, this value cannot be verified from the information provided. Please report the method and parameters used to compute the Lyapunov exponent, or remove the quantitative claim.","section":"Section III"}],"minor_comments":[{"comment":"There is a typo in 'Hemholtz-Duffing' (should be 'Helmholtz-Duffing').","section":"Introduction"},{"comment":"The mechanism name 'Pommeau-Maneville intermittency' should be 'Pomeau-Manneville intermittency.'","section":"Introduction"},{"comment":"The caption says 'for c2 = 91.5464883 in subplots (g-h)', but the third row contains three columns; the histogram panel is presumably (i). Please correct the subplot references.","section":"Fig. 2 caption"},{"comment":"The statement 'n ∈ R \\ {0} and n>1' is redundant and inconsistent with 'n can take any value except 0 and ±1'; please clarify the admissible range of n.","section":"Section II, Eq. (3)"},{"comment":"The fitting was performed using 'the MATLAB Distribution Fitter App,' which is not reproducible. Please provide the fitting procedure, likelihood optimization details, or a script/data file.","section":"Section II, Table I"},{"comment":"The figure caption states c2=91.564884, which differs from the text value 91.5464884; this should be corrected along with the other parameter inconsistencies.","section":"Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"To the editor: This is a numerically driven paper and the authors do not provide code or data. Given the very narrow parameter window and the rare-event nature of the claim, I would strongly encourage the journal to ask for a reproducibility statement or data-sharing as part of the revision. The parameter inconsistencies in the figure captions are not merely cosmetic; they undermine confidence in the reported results until corrected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper reports a genuinely new subtype of extreme event: rare, irregular mixed-mode oscillatory bursts caused by chaotic inter-well hopping in an asymmetric Helmholtz-Duffing oscillator. That observation looks real, and the distinction from earlier work (where extreme events are isolated large spikes rather than MMO bursts) is legitimate and clearly explained. The mechanism section is the strongest part: the authors show how increasing c2 deepens the left well and makes inter-well excursions rarer, which is a clean, testable narrative.\n\nThe soft spots are real but not fatal to the core observation. The fixed-step RKF45 at h=0.01 is the most load-bearing issue because the phenomenon is defined by rare escapes near a crisis; a convergence test (or a second integrator with error control) is essential to rule out numerical noise manufacturing the escapes. The inconsistent c2 values in the figure captions (Figs. 4 and 5 say 5.6494884 while the text uses 91.5464884) are sloppy and must be fixed, but the main results use the latter value consistently. The POT threshold n=5 is arbitrary, and the GEV fit is reported without a goodness-of-fit test; that is common in this literature but worth a sentence. The velocity-as-precursor claim is the weakest: it is identified from the same time series in which the bursts appear, with no false-alarm or miss rates, so \"reliable lead indicator\" overstates what is currently a visually motivated observation.\n\nI am less worried than the stress-test note about the fixed-step integrator manufacturing the bursts: the phase portraits and the monotonic transition with c2 are consistent with a real dynamical effect. But the concern is legitimate, and the authors can answer it with a two-line convergence check. No code or data is provided, which is frustrating for a numerical paper.\n\nWho is this for? Researchers working on extreme events in oscillators, especially the MMO and rare-hopping crowd. It is a solid contribution to that niche, not a breakthrough. I would send it to peer review—the phenomenon deserves scrutiny and the paper is honest about its scope—but I would ask for the numerical validation and the precursor analysis to be tightened before publication.","headline":"Genuinely new extreme-event subtype with a clear mechanism, but the numerical validation and the velocity-precursor claim need real work before the details can be trusted.","tokens_in":12858,"tokens_out":1740,"would_cite":false,"duration_ms":17770,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34C15","37D45","37M05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper reports a new kind of extreme event—irregular mixed-mode oscillatory bursts—in a driven Helmholtz-Duffing oscillator with an asymmetric double-well potential, and argues that the system velocity is a reliable lead indicator of…","keywords":["extreme events","mixed-mode oscillations","Helmholtz-Duffing oscillator","peak-over-threshold method","generalized extreme value distribution","inter-well hopping","velocity lead indicator","double-well potential"],"falsifier":"Recompute the same parameter values with an adaptive Runge-Kutta solver at a tolerance much smaller than the fixed step of 0.01, and check whether the irregular mixed-mode bursts and their preceding velocity spikes still occur at $c_2=91.5464884$; if they vanish or the velocity lead disappears, the reported events and indicator are artifacts of the fixed-step integration.","tokens_in":11921,"feed_emoji":"⚡","tokens_out":6923,"duration_ms":61663,"temperature":0.7,"pith_summary":"The paper reports a new type of extreme event: extreme irregular mixed-mode oscillatory bursts in a driven, asymmetric double-well Helmholtz-Duffing oscillator. It argues that rare chaotic hopping between the two potential wells generates short bursts of irregular mixed-mode oscillations embedded in otherwise periodic mixed-mode oscillations, and that these bursts qualify as extreme events by the peak-over-threshold criterion. It further shows the burst statistics fit a generalized extreme value distribution of Weibull type, and that each burst is preceded by a sharp upward deviation in the system velocity, making velocity a reliable lead indicator. If correct, this adds a new mechanism and a measurable precursor for extreme events in a widely used oscillator model.","feed_headline":"Velocity spikes herald extreme bursts in a driven oscillator","feed_subtitle":"Chaotic hopping between two wells creates the bursts; speed surges give early warning.","key_machinery":"The central object is the asymmetric double-well potential $V(x)=c_1 x^2/2 + c_2 x^3/3 + c_3 x^4/4$ with negative $c_1$ and positive $c_2$ and $c_3$, where the quadratic coefficient $c_2$ tunes the depth asymmetry between the wells. The mechanism that carries the argument is rare chaotic inter-well hopping: the trajectory spends long stretches in the deeper well in periodic mixed-mode motion, and occasionally crosses the unstable fixed point at $x=0$ to reach the shallower well, briefly showing irregular mixed-mode oscillations before returning. The velocity variable $y=\\dot{x}$ is the proposed lead indicator because it deviates from its periodic burst pattern just before the position variable $x$ produces an extreme mixed-mode burst.","core_discovery":"For the driven Helmholtz-Duffing oscillator with $\\mu=0.1$, $c_1=-1.0$, $c_3=534.53$, $F=0.5$, and $\\omega=0.42$, the paper finds that choosing the quadratic nonlinearity coefficient near $c_2=91.5464884$ places the system in a regime where the trajectory is mostly confined to the deeper left well, exhibiting periodic mixed-mode oscillations with one large and eight small peaks. Rarely, the trajectory chaotically crosses the unstable fixed point at the origin and hops to the shallower right well, producing an irregular mixed-mode burst whose peaks cross the peak-over-threshold threshold; these bursts are identified as extreme events. The hop occurs during the rising phase of the external drive, and the burst statistics are fitted to the generalized extreme value distribution with a negative shape parameter, indicating a Weibull-type tail. The paper further states that every burst is consistently preceded by a sudden upshot in the velocity variable, which it proposes as a reliable lead indicator of the event.","pith_inferences":["The velocity-precursor claim is presented qualitatively; a natural extension is to define a quantitative threshold on $y$, measure detection rates and false-alarm rates over long simulations, and compare against random-chance prediction.","The same asymmetry-induced inter-well bursting mechanism might appear in other double-well systems with a periodic mixed-mode baseline, such as buckled beam snap-through or energy harvesters, where a deeper well similarly suppresses rare escapes.","Because the generalized extreme value fit has a negative shape parameter (bounded tail), one could test how the tail and the extreme-event rate vary as $c_2$ approaches the transition, providing a parameter sweep the paper does not report."],"forward_implications":["Extreme events in nonlinear oscillators need not appear only as isolated large-amplitude spikes in a chaotic time series; they can also appear as intermittent bursts of irregular mixed-mode oscillations embedded in periodic mixed-mode oscillations.","The depth asymmetry of the double well, controlled by $c_2$, acts as a tunable switch: increasing $c_2$ deepens the left well and monotonically shortens the time spent in the right well until inter-well hopping becomes impossible, while decreasing $c_2$ makes the bursts so frequent that they no longer qualify as extreme events.","The rising phase of the external drive gates the inter-well transition, so the phase of the drive is part of the event mechanism and could be used in conjunction with the velocity indicator.","The velocity variable can be tracked as a single measurable early-warning signal for the bursts, without needing to know the full phase-space state in advance."],"supporting_citations":[{"why":"Supplies the peak-over-threshold criterion used to classify extreme events and the distinction between rare and extreme.","marker":"[2]"},{"why":"Frames extreme-event prediction as a grand challenge, motivating the search for reliable indicators.","marker":"[25]"},{"why":"Provides the forced Liénard system where interior crisis and Pomeau-Manneville intermittency generate extreme events, the baseline mechanism this paper contrasts with.","marker":"[32]"},{"why":"Reports mixed-mode oscillations and extreme events in a fractional-order Bonhoeffer-van der Pol oscillator, the closest prior study combining MMO with extreme events.","marker":"[39]"},{"why":"Describes extreme bursting oscillations via pulse-shaped explosion, a prior burst-type extreme event whose structure is compared to the new irregular MMO bursts.","marker":"[52]"},{"why":"Reviews the Helmholtz-Duffing oscillator and its applications, establishing the model's practical relevance.","marker":"[61]"},{"why":"Reports extreme events in a Liénard system with an asymmetric double-well potential, the prior asymmetry-based mechanism the authors differentiate from their inter-well hopping bursts.","marker":"[62]"}],"fun_headline_variants":["Velocity surge predicts chaotic well-hopping bursts","Speed spikes forecast extreme oscillator bursts","Velocity spike flags imminent rare bursts","Velocity spike is early warning for extreme bursts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the fixed-step numerical integration at step size 0.01 reproduces the rare bursts as genuine dynamics rather than numerical artifacts, and that the peak-over-threshold choice of $n=5$ correctly marks the boundary between ordinary oscillations and extreme events.","fun_headline_variants_meta":{"raw":{"variants":["Velocity surge predicts chaotic well-hopping bursts","Speed spikes forecast extreme oscillator bursts","Velocity spike flags imminent rare bursts","Velocity spike is early warning for extreme bursts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000281,"raw_usage":{"total_tokens":1643,"prompt_tokens":904,"completion_tokens":739,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":688}},"tokens_in":520,"tokens_out":739,"duration_ms":6704,"temperature":1.0,"reasoning_tokens":688,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:44:56.803978+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the same parameter values with an adaptive Runge-Kutta solver at a tolerance much smaller than the fixed step of 0.01, and check whether the irregular mixed-mode bursts and their preceding velocity spikes still occur at $c_2=91.5464884$; if they vanish or the velocity lead disappears, the reported events and indicator are artifacts of the fixed-step integration.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the peak-over-threshold criterion used to classify extreme events and the distinction between rare and extreme."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Frames extreme-event prediction as a grand challenge, motivating the search for reliable indicators."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the forced Liénard system where interior crisis and Pomeau-Manneville intermittency generate extreme events, the baseline mechanism this paper contrasts with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports mixed-mode oscillations and extreme events in a fractional-order Bonhoeffer-van der Pol oscillator, the closest prior study combining MMO with extreme events."},{"cited_title":"Sudharsan, A","cited_arxiv_id":null,"evidence_quote":"Describes extreme bursting oscillations via pulse-shaped explosion, a prior burst-type extreme event whose structure is compared to the new irregular MMO bursts."},{"cited_title":"Roy and S","cited_arxiv_id":null,"evidence_quote":"Reviews the Helmholtz-Duffing oscillator and its applications, establishing the model's practical relevance."},{"cited_title":"Roy and S","cited_arxiv_id":null,"evidence_quote":"Reports extreme events in a Liénard system with an asymmetric double-well potential, the prior asymmetry-based mechanism the authors differentiate from their inter-well hopping bursts."}],"review_version":1}