{"id":"abc2c8ff-82ef-4f1a-872b-74372abaad2c","arxiv_id":"2607.27242","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A micromechanical parametric oscillator with tunable nonlinear friction exhibits two coexisting pairs of period-two states whose onset follows a swallowtail catastrophe.","lead":"This paper shows that a single micromechanical oscillator can host two pairs of stable vibration states at the same time when its friction is engineered through a two-phonon coupling. That expands the standard picture of parametric resonance and offers a platform for studying rare catastrophe structures and multistable dynamics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The swallowtail catastrophe claim rests on an asserted 4D-to-1D normal-form reduction in Appendix C whose coefficients are never computed; without verifying the vanishing quartic term, the catastrophe class is not established.","rationale":"The reader's weakest assumption precisely identifies the unvalidated normal-form reduction in Appendix C, and I concur that this is the most load-bearing concern in the paper. The central experimental observation—two coexisting pairs of period-two states—is supported by direct data and by a full model whose parameters, though incompletely reported, are described as fixed by ringdown measurements rather than fit to the multistability data. However, the additional claim that the boundary of multistability is a swallowtail catastrophe is a topological classification that depends entirely on the reduction of the four-dimensional slow dynamics to a one-dimensional normal form. Because the coefficients are never computed and the quartic term is not shown to vanish, the A4 classification is not established. This does not undermine the core multistability result, but it does mean the paper's most theoretical claim is not yet quantitatively demonstrated. The proposed concrete test—a numerical center-manifold reduction followed by comparison of the predicted and full-model bifurcation sets—would settle the matter. I agree with the reader's CONDITIONAL verdict: the paper is solid but needs this validation and improved reproducibility. I would not change the verdict, hence UNCHANGED.","tokens_in":16166,"tokens_out":11058,"duration_ms":101085,"concrete_test":"Perform a numerical center-manifold reduction of equations (B1)-(B2) near the purported swallowtail point: (1) continue the full model to locate the codim-3 point where two cusp lines meet in (eps, h, fp) space; (2) compute the Jacobian and identify the unique soft mode (zero eigenvalue); (3) project the Taylor expansion of the vector field onto the center manifold to fifth order; (4) verify that the coefficient of z^4 vanishes (within numerical tolerance) and the coefficient of z^5 is positive; (5) derive the effective a, b, c as functions of eps, fp, hd and compare the predicted bifurcation curves with a direct numerical continuation of the full model, including the scaling of the distance between cusp ridges as fp approaches the swallowtail value. If the quartic term is nonzero or the predicted curves disagree, the swallowtail claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim includes that the onset of multistability is 'governed by a swallowtail catastrophe.' This rests on Appendix C, where the four-dimensional slow dynamics (B1)-(B2) is reduced to a one-dimensional effective potential V_eff = z^5/5 + a z^3/3 + b z^2/2 + c z. However, the coefficients a, b, c are never computed from the full equations, and the reduction itself is not validated. In particular, the soft-mode eigenvector is not identified, the projection of the vector field onto the center manifold is not performed, and the key requirement for an A4 (swallowtail) singularity—that the coefficient of the quartic term z^4 vanishes after reduction—is not checked. If the quartic term survives, the relevant singularity is the butterfly (A5), whose bifurcation set differs from the swallowtail. Moreover, the claim that a, b, c are linear combinations of the deviations of eps, fp, and hd is only true for a universal unfolding; the actual mapping requires computing parameter derivatives of the center-manifold reduction, which is not supplied. Without this, the 'swallowtail catastrophe' is an interpretation imposed on observed bifurcation topology rather than a derived consequence of the model. The experimental multistability itself is supported by direct measurements in Figs. 4(c-d) and by the full-model calculations, but the theoretical identification of the catastrophe class is under-supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports experiments on a micromechanical resonator whose fundamental mode is coupled to a faster-decaying auxiliary mode by a sideband pump, producing nonlinear friction. After characterizing the nonmonotonic amplitude-dependent damping via ringdown measurements, the authors apply parametric modulation and map the stationary response. They observe an isolated branch (isola) for moderate drive and, for another sideband-pump amplitude, a second pair of stable period-two states coexisting with the conventional pair and with the stable zero-amplitude state. Theoretical response curves and bifurcation lines are obtained from a two-mode model, Eqs. (1)-(2)/(B1)-(B2), with parameters stated to be fixed by ringdown data. The paper interprets the merging of cusp points in the three-parameter space (detuning, parametric-drive amplitude, sideband-pump amplitude) as a swallowtail catastrophe and proposes a one-dimensional effective potential of the form z^5/5 + a z^3/3 + b z^2/2 + c z.","tokens_in":16573,"tokens_out":4557,"duration_ms":51406,"significance":"If the claims hold, the paper establishes a qualitatively new phenomenon for a single parametric oscillator: coexistence of two distinct pairs of stable period-two states, induced by controlled nonlinear friction. This extends the canonical single-pair picture and is directly relevant to nanomechanics, Ising machines, and cat-qubit implementations. The experimental multistability is supported by direct measurements in Figs. 4(c)-(d) and by full-model calculations, and the use of two devices strengthens the result. However, the identification of the multistability boundary as a swallowtail catastrophe is presently asserted rather than derived in Appendix C, and the quantitative comparison lacks reported parameter values and measurement uncertainties. These issues are fixable but are load-bearing for the paper's central headline claim.","major_comments":[{"comment":"The swallowtail classification is asserted, not derived. The text states that one must linearize around v_j^sw to identify the soft mode, then writes z = Σ c_i w_i and says that 'after rescaling' V_eff has the swallowtail form. The center-manifold projection is never performed: the soft-mode eigenvector c_i is not given, the coefficient of z^4 is not shown to vanish, and the mapping from (ϵ, f_p, h_d) to (a,b,c) is never computed. The A4 classification therefore rests on an assumed universal unfolding. Please supply this reduction or verify it numerically from Eqs. (B1)-(B2). The experimental multistability by itself does not establish the catastrophe class.","section":"Appendix C, Eqs. (C1)-(C3)"},{"comment":"Measured bifurcation points are plotted without uncertainty intervals, and the text reports no repeated-sweep statistics. Claims of 'quantitatively map[ping]' the bifurcation structure (Abstract; Sec. 4) need at least representative error bars on ϵ, h, and the extracted bifurcation frequencies, as well as a statement of how measurement noise propagates to the displayed theory lines.","section":"Figs. 2(c), 3(c), 4(a)-(d)"},{"comment":"The theory lines use nonlinear coefficients γ1, γ2, γ (and hence Λ11, Λ22, Λ12) and effective masses, but Table 1 lists only ω1,2 and Γ1,2. The numerical values of these parameters and their extraction procedure (presumably from ringdown data) are not stated. Without them the comparison in Figs. 2-4 is not reproducible and the parameter-fixing protocol cannot be checked.","section":"Secs. 2.1-2.2, Eqs. (1)-(4)"}],"minor_comments":[{"comment":"The fourth equation reads 'ẏ1 = F4'; this should be 'ẏ2 = F4'.","section":"Appendix C, Eq. (C1)"},{"comment":"Reference 20 appears as 'arXiv:2062.06559v1 (2026)'; this identifier is not in a valid arXiv format and should be corrected or replaced.","section":"Ref. 20"},{"comment":"The normal-form potential is written as V = x^5 + a x^3 + b x^2 + c x in Sec. 3.3 but as V_eff = z^5/5 + a z^3/3 + b z^2/2 + c z in Eq. (C3). The coefficient conventions should be unified or the rescaling explained.","section":"Sec. 3.3 vs. Appendix C"},{"comment":"The table formatting is garbled; the device labels and units should be typeset cleanly.","section":"Table 1"},{"comment":"The text says the merging at critical point C is 'directly measured,' but Fig. 2(d) is described as a numerical calculation. Please clarify which quantities are experimental and which are simulated.","section":"Sec. 3.1 and Fig. 2(d)"}],"recommendation":"major_revision","confidential_remarks":"The experimental result — two coexisting pairs of period-two states — appears solid and important. My main concern is that the swallowtail catastrophe identification is asserted in Appendix C rather than derived; the authors should be asked to provide the actual center-manifold reduction or a numerical verification of the normal form. They should also report the nonlinear model parameters and experimental uncertainties. I recommend major revision rather than rejection because the central experimental phenomenon is directly evidenced and the theoretical gap is local, though essential to the headline claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the most useful thing to know: the central experimental claim is credible. Two distinct pairs of stable period-two states coexisting in a single parametric oscillator is directly visible in Figs. 4(c)-(d), and the theory lines are not fit to that data—they come from equations of motion with parameters extracted from ringdown measurements. Two devices give the same picture. That part is solid.\n\nWhat is genuinely new here is the experimental realization of multistable parametric resonance beyond the standard single pair, and the demonstration that the induced nonlinear friction mechanism—reported by the same group in 2018—leads to this richer bifurcation structure. The isola branch and the merging of b3 and b4 are nice additions. If the experimental result survives scrutiny, it has consequences for Ising machines and cat-qubit schemes where people assume a single pair of period-two states.\n\nNow the soft spots. The largest is the theoretical identification of the swallowtail catastrophe. Appendix C states that the four-dimensional slow dynamics reduces to a one-dimensional effective potential z^5/5 + a z^3/3 + b z^2/2 + c z, with a,b,c linear in the parameter deviations. But the coefficients are never computed from equations (B1)-(B2), the soft mode is not explicitly identified, and, crucially, the vanishing of the quartic term is never checked. If that term does not vanish, the singularity is a butterfly (A5), not a swallowtail (A4). So the 'governed by a swallowtail catastrophe' claim is presently an interpretation imposed on the observed bifurcation topology, not a derived result. The experimental observation of two cusp points merging as fp is varied is consistent with a swallowtail, but it is not by itself a proof of the normal form.\n\nOther weaknesses are more conventional: there are no error bars on the measured bifurcation points, the nonlinear coupling coefficients and effective masses are not reported, and there is no raw data or code deposit. The data-availability statement just says data are available on request. This is below the standard for a paper making this level of claim.\n\nNet assessment: the experimental multistability is a real result and the paper deserves a serious referee. The normal-form argument needs to be fixed or softened, and the reproducibility gaps need to be addressed. I would send it to peer review, with a clear request to substantiate the catastrophe classification and provide the missing parameters and error estimates.","headline":"The two-pair period-two multistability is real and well supported; the swallowtail catastrophe label is asserted rather than derived, and the paper needs to fix that gap plus report its parameters and error bars.","tokens_in":17023,"tokens_out":3453,"would_cite":true,"duration_ms":34605,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that a single micromechanical parametric oscillator, when its nonlinear friction is produced by resonant coupling to a faster-decaying mode, can host two distinct pairs of stable period-two oscillations at once, and that th","keywords":["parametric resonance","nonlinear friction","period-two states","multistability","swallowtail catastrophe","bifurcation","micromechanical oscillator","two-phonon coupling"],"falsifier":"Compute the bifurcation set directly from equations (B1)–(B2) and check whether the two cusp ridges actually meet at a single swallowtail point when ε, h and fp are varied; alternatively, measure the stationary probability distribution near the expected swallowtail point: the theory predicts a specific non-Gaussian scaling (log of the distribution proportional to the sixth power of the soft-mode coordinate) and the simultaneous coalescence of four states, so the absence of that scaling or of the predicted merging would disprove the claim.","tokens_in":16027,"feed_emoji":"🦋","tokens_out":5969,"duration_ms":57338,"temperature":0.7,"pith_summary":"The paper establishes that, contrary to the usual expectation of a single pair of phase-opposite period-two states, a micromechanical parametric oscillator can simultaneously support two distinct pairs of stable period-two oscillations. This occurs when the oscillator's nonlinear friction is implemented by a resonant drive that couples the main mode to a faster-decaying auxiliary mode, because the induced friction is non-monotonic in amplitude: it first grows, then falls as amplitude-dependent frequency shifts break the resonance condition. The authors show experimentally and with a two-mode model that the resulting bifurcation diagram contains an isolated branch and, when a third control parameter is varied, a swallowtail catastrophe in which two cusps merge and two pairs of stable states appear. Establishing this matters because multiple coexisting period-two states would extend the canonical picture of parametric resonance and could enable multistate encoding or richer fluctuation physics in parametric-oscillator platforms.","feed_headline":"Two pairs of period-two states appear in one parametric oscillator","feed_subtitle":"Nonmonotonic nonlinear friction turns the bifurcation diagram into a swallowtail catastrophe, enabling multistate parametric oscillators.","key_machinery":"The central object is the drive-induced two-phonon coupling between the low-frequency plate mode (mode 1) and a high-frequency, faster-decaying beam mode (mode 2), produced by a sideband pump at ω2 − 2ω1. This realizes controllable nonlinear friction: energy leaves mode 1 two quanta at a time through mode 2. The crucial feature is that the friction is non-monotonic in the vibration amplitude: at small amplitudes it grows quadratically, but once the conservative nonlinearities shift the mode frequencies by an amount comparable to Γ2, the resonance condition breaks and the friction falls back toward its zero-amplitude value. That non-monotonicity creates an isolated branch (isola) of period-tw","core_discovery":"On the authors' terms: with conventional nonlinear friction, a parametrically driven oscillator has one pair of stable period-two states and one stable zero-amplitude state; the new claim is that controlled nonlinear friction — two-phonon loss to a faster-decaying mode induced by a sideband pump — changes the bifurcation structure so that two pairs of stable period-two states can coexist, together with two pairs of unstable period-two states and the zero-amplitude state, in a triangular region of the (detuning, parametric-drive) plane. The onset of this multistability is governed by a swallowtail catastrophe: as the sideband-pump amplitude is varied, the two cusp points of the bifurcation di","pith_inferences":["Editorial inference: the non-monotonic friction mechanism is generic to any driven two-phonon-resonant system, so superconducting cat-qubit setups with engineered two-photon loss should also exhibit multiple period-two manifolds if the auxiliary mode's frequency shifts with amplitude.","Editorial inference: the appendices leave the coefficients a, b, c of the swallowtail normal form undetermined; explicitly computing them from equations (B1)-(B2) would convert the structural identification into a quantitative prediction, e.g., of the exact triangular region's boundaries.","Editorial inference: near the swallowtail, four stationary states coalesce, so switching rates between the two stable period-two pairs and the zero state should display nontrivial scaling with ε, h, fp; measuring these escape rates would test the normal form beyond static bifurcation topology.","Editorial inference: the existence of two stable period-two pairs suggests a natural multistate memory element, but the isolated branch is only reachable by the two-step excitation protocol; a single-parameter read/write scheme would need to be engineered."],"forward_implications":["A single parametric oscillator can host nine coexisting stationary states: two pairs of stable period-two oscillations, two pairs of unstable period-two oscillations, and the stable zero-amplitude state.","The multistability boundary is a swallowtail catastrophe, so reaching it requires tuning three independent control parameters; the system provides a controlled experimental map of a codimension-three catastrophe.","Near the critical point where five states merge, the stationary distribution is strongly non-Gaussian, with the logarithm of the distribution behaving as a sixth power of the soft-mode coordinate rather than a parabola.","The non-monotonic friction produces an isolated response branch that cannot be reached by sweeping the drive frequency alone; the paper shows a two-step protocol (adjusting the parametric drive strength) that accesses it.","The results establish electromechanical oscillators as a platform for quantitatively studying catastrophe theory and multistable nonequilibrium dynamics."],"fun_headline_variants":["Parametric oscillator hosts two stable period-two pairs via nonlinear friction","Swallowtail catastrophe enables multistable parametric oscillators","Nonlinear friction adds a second period-two state pair in one oscillator","Two-phonon loss creates coexisting period-two states in micromechanical oscillator","Micromechanical oscillator shows two period-two pairs with swallowtail bifurcation"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The claim that the onset is a swallowtail catastrophe rests on the reduction in Appendix C, where the four-dimensional slow-amplitude dynamics is replaced by the one-dimensional normal form V(z)=z^5/5 + a z^3/3 + b z^2/2 + c z; the coefficients a, b, c are asserted, not computed from the underlying equations, so if higher-order terms break that reduction the swallowtail identification becomes an interpretation rather than a derived result.","fun_headline_variants_meta":{"raw":{"variants":["Parametric oscillator hosts two stable period-two pairs via nonlinear friction","Swallowtail catastrophe enables multistable parametric oscillators","Nonlinear friction adds a second period-two state pair in one oscillator","Two-phonon loss creates coexisting period-two states in micromechanical oscillator","Micromechanical oscillator shows two period-two pairs with swallowtail bifurcation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1223,"prompt_tokens":739,"completion_tokens":484,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":483,"completion_tokens_details":{"reasoning_tokens":392}},"tokens_in":483,"tokens_out":484,"duration_ms":4353,"temperature":1.0,"reasoning_tokens":392,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T23:16:14.739353+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the bifurcation set directly from equations (B1)–(B2) and check whether the two cusp ridges actually meet at a single swallowtail point when ε, h and fp are varied; alternatively, measure the stationary probability distribution near the expected swallowtail point: the theory predicts a specific non-Gaussian scaling (log of the distribution proportional to the sixth power of the soft-mode coordinate) and the simultaneous coalescence of four states, so the absence of that scaling or of the predicted merging would disprove the claim.","supporting_citations":[],"review_version":1}