Pith. sign in

REVIEW 3 major objections 4 minor 59 references

A single ringdown trace, read through the backbone relation, returns the conservative nonlinear frequency–amplitude expansion coefficients of a nanomechanical resonator with errors below 10%.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review

2026-08-04 16:45 UTC pith:3E6TPJ36

load-bearing objection A genuinely useful extension of ringdown backbone extraction to higher-order nonlinear coefficients, with convincing cross-method agreement for alpha1 and alpha2, but the alpha3 value rests on an unvalidated demodulator-bandwidth assumption that should be tested before the headline claim is trusted. the 3 major comments →

arxiv 2608.02008 v1 pith:3E6TPJ36 submitted 2026-08-03 cond-mat.mes-hall

Extracting higher-order nonlinearities in nanomechanical resonators using the backbone relation

classification cond-mat.mes-hall
keywords nanomechanical resonatorsbackbone curveringdownDuffing nonlinearityhigher-order nonlinearitiesinstantaneous frequencyfrequency driftlock-in demodulation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper seeks to establish that the instantaneous frequency of a freely decaying nanomechanical resonator, plotted against its decaying amplitude, directly encodes the conservative nonlinear coefficients alpha1, alpha2, alpha3 in the expansion of the nonlinear eigenfrequency. A single ringdown measurement captured by lock-in demodulation yields these coefficients with smaller uncertainty than conventional frequency-response fits, and it inherently absorbs small drifts of the bare eigenfrequency. The method explicitly includes symmetry-breaking (odd-order) nonlinearities, which earlier ringdown analyses omitted. The authors validate the approach by comparing backbone-derived coefficients with frequency-response fits on a high-Q silicon nitride string, finding agreement within a few percent and smaller fit errors from the backbone, especially for the lower-order terms.

Core claim

The paper establishes that the backbone relation, Eq. (2), truncated at increasing order, describes the instantaneous frequency observed during a ringdown. Fitting the amplitude–frequency data from a single ringdown with this relation returns the conservative nonlinear coefficients alpha1, alpha2, alpha3 that match the values obtained from fitting the driven frequency response, but with substantially smaller relative errors (1.7% for alpha1 and 8.3% for alpha2 vs 10.7% and 12.6% for response fits; alpha3 errors are comparable at about 4%). The backbone method also fits the detuning of the drive from the bare eigenfrequency as a nuisance parameter, which absorbs static frequency drift without

What carries the argument

The load-bearing object is the backbone relation, omega(A) = omega0 + M alpha1 A^2 omega0/2 + M^2 alpha2 A^4 omega0^2/4 + M^3 alpha3 A^6 omega0^3/8 + ..., the frequency–amplitude expansion of the nonlinear eigenfrequency of a weakly nonlinear mode (Eq. (2)). During ringdown, the amplitude decays exponentially while the instantaneous frequency is read as 2*pi*dphi/dt from the in-phase and quadrature outputs of a lock-in demodulator. Fitting this curve extracts the coefficients directly, with the fitted detuning delta-omega absorbing bare-eigenfrequency drifts. The action–angle derivation in the appendix provides the theoretical link: omega_I(I) = dH/dI expanded in powers of the action, conver

Load-bearing premise

The method assumes that the demodulated instantaneous frequency 2*pi*dphi/dt faithfully follows the true nonlinear eigenfrequency throughout the ringdown, which requires the lock-in bandwidth and sampling rate to track a frequency excursion of up to about 1.2 kHz without low-pass distortion — the paper states this requirement but never quantifies or independently validates it.

What would settle it

Record the same ringdown at decreasing lock-in demodulator bandwidths (or with a different detection chain such as an optical interferometer with a frequency counter) and check whether the fitted alpha1, alpha2, alpha3 remain constant. If the coefficients shift when the bandwidth approaches the observed detuning range (up to 1208 Hz), or when the sampling rate is reduced, the backbone extraction is distorted and the reported accuracy is not intrinsic to the method.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • Nonlinear characterization reduces to a single ringdown trace (about 0.2 s), much faster than a high-Q frequency response sweep, reducing exposure to drift.
  • The fitted detuning parameter absorbs static eigenfrequency drift, so the extracted alpha_i are stable even when the bare frequency shifts by several linewidths between measurements.
  • The method naturally includes symmetry-breaking (odd-order) nonlinearities, unlike earlier ringdown-based extraction schemes.
  • For a drive where the response deviates by more than 5% from the model at the current truncation order, the next alpha_i is included; the paper demonstrates this for alpha1, alpha2, alpha3.
  • The phase-space trajectory's sense-of-rotation reversal at small detunings provides a visual consistency check linking the observed amplitude to the backbone crossing frequency.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Because the extracted coefficients are only as trustworthy as the instantaneous-frequency readout, comparing backbone fits against an independent frequency reference (e.g., a phase-locked loop or a heterodyne counter) at the largest detuning would separate genuine nonlinear frequency shifts from demodulator artifacts.
  • The same single-trace approach could be extended to dissipative nonlinearities by letting the damping rate itself depend on amplitude; the exponential-decay check in the paper (App. B) provides a ready diagnostic for when such an extension becomes necessary.
  • The observed rotation reversal in phase space could be turned into a calibration point: the amplitude at which the trajectory reverses should coincide with the amplitude at which the backbone crosses the reference frequency, giving an in-situ cross-check on the amplitude calibration.
  • For high-order coefficients like alpha3, information is concentrated near the top of the ringdown; a weighted fit or a shorter, higher-power initial segment could further reduce their uncertainty.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper presents a method for extracting conservative nonlinear frequency-amplitude coefficients α1, α2, α3 of a nanomechanical resonator from a single ringdown measurement. By recording in-phase and quadrature components with a lock-in amplifier, the authors compute the instantaneous frequency from the time derivative of the demodulated phase (Eq. (6)) and fit the backbone relation (Eq. (2)) with a simultaneous detuning parameter to absorb eigenfrequency drift. The method is demonstrated on a SiN string at three drive powers and benchmarked against fits of the driven response using Eq. (A13). Reported backbone coefficients are α1 = 4.097×10^25 kg^-1 m^-2, α2 = 8.4162×10^46 s kg^-2 m^-4, α3 = 5.783×10^68 s^2 kg^-3 m^-6, compared with response-fit values 4.170×10^25, 8.402×10^46, and 6.176×10^68, respectively. The paper also reports an action-ringdown validation and a phase-space criterion for the ringdown start time.

Significance. If the method is sound, it offers a practical and faster alternative to frequency-response fits for characterizing higher-order conservative nonlinearities, including symmetry-breaking terms, with reduced sensitivity to low-frequency drift. The paper's strengths include the use of two independent measurement channels, a thermomechanical calibration, a direct check that the action ringdown is exponential, and explicit inclusion of detuning in the backbone fit. The cross-method agreement for α1 and α2 is strong and supports the framework. The main caveat is that the central equation relies on unvalidated measurement-chain assumptions, so the absolute accuracy and the α3 value rest on a premise that needs direct testing.

major comments (3)
  1. [Section II, Eq. (6)] The identification ω_inst = 2π dφ/dt = ω(A) − ω_d assumes the lock-in demodulator can track the instantaneous frequency without phase distortion over the full ringdown, including detunings up to 1208 Hz at −26 dBm. The text only says the bandwidth 'must be set large enough' and gives no bandwidth value, filter order, sampling rate, or validation against a known chirp. Any frequency-dependent phase lag or amplitude roll-off will bias all α_i, most strongly α3, whose weight lies in the earliest high-detuning points (App. E). The frequency-response comparison does not clear this because those measurements are quasi-static. Please quantify the demodulator settings and provide a direct phase-fidelity test.
  2. [Table II / App. E] The reported errors are fit standard errors and omit propagation from the fixed lower-order coefficients and from the calibration chain (a, c, M, Γ). For α2 and α3, the backbone fits use α1 (and α2) as fixed inputs, so their uncertainty is not included. The α3 backbone value differs from the response value by 6.8%, while the individual fit errors are about 4.3% and 3.5%; a proper covariance or total-error budget is needed to justify the statement that the coefficients are accurate to 'below 10%'.
  3. [Sec. III and App. B] The method is presented as requiring no prior knowledge of the nonlinearities, but the extraction of α2 and α3 is done with lower-order coefficients fixed from previous fits. In addition, the conclusions are based on single ringdown traces per drive power. Repeating a ringdown several times would allow a direct assessment of run-to-run variability and would strengthen the precision claim beyond the internal fit errors.
minor comments (4)
  1. [Table II] The numbers in the first column appear garbled (e.g., '14.170·10^25' and '22.096·10^25'); I assume these are 4.170 and 2.096. Please correct the formatting.
  2. [App. B, Fig. 4] The action axis is labeled 'kg·V^2·s^-1', which is an unusual unit combination. Please clarify the scaling and how the phase-space area is converted to action.
  3. [Sec. III] The 5% threshold for adding the next nonlinear order is empirical. A formal model-selection criterion or a discussion of how the threshold affects the extracted coefficients would make the procedure more robust.
  4. [App. A] Equations (A6) and (A13) are attributed to prior work [42,51]. Since these are load-bearing for the response fits, a brief derivation or at least a statement of assumptions would help the reader assess the comparison.

Circularity Check

0 steps flagged

No significant circularity: the nonlinear coefficients are fit outputs, and the cross-method agreement is validation rather than a reduction to the model inputs.

full rationale

The paper's central claim is that the backbone curve from a ringdown, fitted with Eq. (2) (equivalently Eq. A9), yields accurate values of the conservative nonlinear coefficients α1, α2, α3. These coefficients are least-squares fit parameters, not predictions derived from first principles. The benchmark is agreement with independent frequency-response fits; such cross-method agreement is a validation, not a circular reduction. Eq. (6), ω_inst(A) = ω(A) − ω_d, is a definition of the demodulated instantaneous frequency relative to the drive, not a disguised imposition of the fit result. The fitted detuning δω in the backbone fit absorbs frequency drifts, but this is an explicit fitting parameter and is not claimed to be an independent prediction. The theoretical expressions in Appendix A, including Eqs. (A6) and (A13), are attributed in part to the authors' prior work [42,51]; however, these formulas are also supported by textbook references [49,50,55], are standard results, and are not the load-bearing element of the extraction—the α_i are fit directly from the measured backbone without mapping back to the potential coefficients γ_i. The manuscript contains self-citations but none is load-bearing in the sense of forcing the conclusion: the extracted coefficients are not outputs of those cited equations alone, and the cited results are externally grounded. The concern about unquantified demodulator bandwidth is a potential systematic-error/correctness issue, not a circularity: it challenges the validity of Eq. (6) under experimental conditions, but does not make the argument self-referential. Accordingly, the derivation chain is self-contained with respect to circularity.

Axiom & Free-Parameter Ledger

5 free parameters · 7 axioms · 0 invented entities

The framework rests on standard action-angle perturbation theory plus three domain assumptions (single-mode isolation, linear damping, near-sinusoidal motion), two of which are experimentally checked. The unquantified premise is demodulator fidelity at large detunings. All extracted coefficients are least-squares fits, so the paper delivers a calibrated measurement method rather than a first-principles derivation of nonlinear coefficients. No new physical entities are postulated.

free parameters (5)
  • Empirical 5% threshold for adding the next nonlinear order = 0.05 relative deviation of omega_max
    Sec. III: the rule for when to include alpha2 or alpha3 is admitted to be empirical ('This threshold is empirical, but applicable for our case') and its influence on the extracted coefficients is not analyzed.
  • Nonlinear coefficients alpha1, alpha2, alpha3 (fit outputs) = alpha1 ~ 4.10e25 kg^-1 m^-2; alpha2 ~ 8.42e46 s kg^-2 m^-4; alpha3 ~ 5.9e68 s^2 kg^-3 m^-6
    These are the quantities the method extracts by least-squares fits of Eq. (2) to data. They are the measured results, not assumptions, but their absolute SI values depend on calibration constants.
  • Displacement calibration a = 3.5181e-5 m/V
    App. C: fitted to the thermomechanical noise spectrum assuming T = 293 K and M = 2.11e-15 kg. All SI alpha_i values scale as a^{-2i}, so calibration error propagates into every reported coefficient.
  • Drive calibration combination c = (b/a)^2 = 0.2645
    App. C: slope of squared response amplitude versus squared drive voltage across linear-response curves; used to derive b. It anchors the absolute force scale but does not enter the extracted alphas directly.
  • Detuning delta_omega in backbone fits = 135 Hz (-35 dBm), 459 Hz (-30 dBm), 1208 Hz (-26 dBm)
    Free parameter added to absorb bare-eigenfrequency drift. Its small stated errors support the drift-robustness claim, but in principle it trades against the alpha coefficients in the fit.
axioms (7)
  • domain assumption The driven mode is isolated and described by a single-coordinate Hamiltonian with anharmonic potential U(q) of Eq. (1)
    Sec. I and App. A: justified by tuning the DC voltage to avoid hybridization with other modes; not verified by a direct intermodal measurement.
  • domain assumption Linear, amplitude-independent damping with Gamma/omega0 << 1, giving exponential action decay (Eq. A7)
    App. B verifies exponential decay of action and amplitude; nonlinear damping is stated to be negligible in the explored amplitude range.
  • domain assumption Weak nonlinearity, so vibrations remain near-sinusoidal and A is approximately sqrt(2I/M omega0) (Eq. A8)
    App. A: 'to the lowest order in the amplitude' the action-amplitude mapping is linear; checked indirectly through exponential action decay (App. B), with higher overtones neglected.
  • standard math omega_I(I) admits the Taylor expansion of Eq. (A4) truncated at the fitted order
    Standard action-angle perturbation theory from textbooks [35, 36]; the truncation order is selected by the empirical 5% rule rather than by a convergence argument.
  • standard math The mappings of Eq. (A6) relating alpha1, alpha2 to potential coefficients gamma_i are correct
    Canonical perturbation results referenced to the authors' own PRX 2022 paper [51]; the Helmholtz-Duffing subset also appears in [49, 50].
  • standard math The driven-response amplitude formula Eq. (A13) is valid
    Rotating-wave/averaged response of a weakly nonlinear damped oscillator, taken from [42, 51, 55].
  • domain assumption The demodulated phase derivative equals the true instantaneous nonlinear eigenfrequency throughout the ringdown
    Sec. II: the demodulator bandwidth 'must be set large enough to capture the complete energy of the resonator as it rings down the backbone curve', but the bandwidth and sampling rate are never quantified or verified against an independent reference.

reviewed 2026-08-04 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Extracting higher-order nonlinearities in nanomechanical resonators using the backbone relation." pith.science (2026). https://pith.science/paper/3E6TPJ36

@misc{pith2026260802008,
  author       = {Pith},
  title        = {Pith review of: Extracting higher-order nonlinearities in nanomechanical resonators using the backbone relation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3E6TPJ36}},
  note         = {Machine review of arXiv:2608.02008}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Nanomechanical resonators are a powerful platform for studying nonlinear dynamics with high sensitivity and precision. We explore the nonlinear response of a high-Q nanomechanical string resonator in and beyond the Duffing regime and introduce a robust framework for accurately extracting its conservative nonlinearities. The method is based on the backbone curve obtained from ringdown measurements, making it inherently resilient to small frequency fluctuations while explicitly accounting for both symmetry-breaking and non-symmetry-breaking nonlinearities. To validate the approach, we perform complementary ringdown and frequency-response measurements on the nanostring resonator and benchmark the backbone-based extraction against established frequency-response techniques. The comparison confirms the accuracy of the proposed framework and demonstrates its advantages over conventional methods for nonlinear characterization.

Figures

Figures reproduced from arXiv: 2608.02008 by Daniel K. J. Bone{\ss}, Eva M. Weig, Maria Kallergi, Maximilian Seitner, Wolfgang Belzig.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Scanning electron micrograph of a doubly-clamped silicon nitride string resonator (green) and flanking electrodes (gold). [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Ringdown measurements in and beyond the Duffing regime. Turquoise and black dots show the response amplitude before and after [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Frequency response at a drive power of [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Action decay measured at the upper bifurcation point for [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: It was obtained by averaging over 1, 000 traces since the signal is weak and buried in noise. A Lorentzian fit with only the calibration parameter 𝑎 as fit parameter in Eq. (C4) and assuming 𝑇 = 293 K, 𝑀 = 2.11 · 10−15 kg, 𝜔0/(2𝜋) = 6.519 MHz and 2Γ/(2𝜋) = 21 Hz is shown in [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Measured thermomechanical spectrum of the fundamental [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. (a) Ringdown measurement reproduced from Fig. 2(a) of the [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

59 extracted references

  1. [1]

    The system is best described using so-called action- angle variables defined by [35] 𝐼=(2𝜋) −1 ∮ 𝑝𝑑𝑞and𝜓=𝜕 𝐼 ∫ 𝑝𝑑𝑞.(A3) Anisolatedmodeperformsperiodicvibrationswith𝐼=const

    Isolated Mode If the mode is isolated, i.e.,𝐹d =Γ=0, the energy is con- served. The system is best described using so-called action- angle variables defined by [35] 𝐼=(2𝜋) −1 ∮ 𝑝𝑑𝑞and𝜓=𝜕 𝐼 ∫ 𝑝𝑑𝑞.(A3) Anisolatedmodeperformsperiodicvibrationswith𝐼=const.. The phase linearly accumulates in time,𝜓=𝜔 I(𝐼)𝑡. The vibration frequency is𝜔I(𝐼)=d𝐸/d𝐼, where𝐸is the m...

  2. [2]

    Then the energy is no longer conserved, and similarly the action becomes time dependent

    Damped Mode Wenowstudythedynamicsofthemodeundertheinfluence of dissipation. Then the energy is no longer conserved, and similarly the action becomes time dependent. Importantly, in the regime whereΓ/𝜔 0 ≪1, the action𝐼only decays slowly, while the phase𝜓performs fast oscillations with the frequency𝜔 I(𝐼)beingapproximatelyconstantoveroneperiod. Averaging o...

  3. [3]

    In action-angle variables Eq

    Driven Mode We now consider the case where the oscillator is driven at frequency𝜔 d with an amplitude𝐹d. In action-angle variables Eq. (A1) takes the form ¤𝜓=𝜔 I(𝐼)− 𝜕𝑞 𝜕𝐼(𝐹d cos(𝜔 d𝑡)−2Γ𝑝),(A10) ¤𝐼= 𝜕𝑞 𝜕𝜓(𝐹d cos(𝜔 d𝑡)−2Γ𝑝).(A11) Thedrivenmodeisexpectedtoperformvibrationsatthedrive frequencywithaconstantamplitude𝐴 0. WemaketheAnsatz 𝑞=𝐴 0 cos(𝜓(𝑡)),(A12) ...

  4. [4]

    Bachtold, J

    A. Bachtold, J. Moser, and M. I. Dykman, Mesoscopic physics ofnanomechanicalsystems,Rev.Mod.Phys.94,045005(2022)

  5. [5]

    Rep.534, 89 (2014)

    M.ImbodenandP.Mohanty,Dissipationinnanoelectromechan- ical systems, Phys. Rep.534, 89 (2014)

  6. [6]

    Sementilli, E

    L. Sementilli, E. Romero, and W. P. Bowen, Nanomechani- cal dissipation and strain engineering, Adv. Funct. Mater.32, 2105247 (2021)

  7. [7]

    N. J. Engelsen, A. Beccari, and T. J. Kippenberg, Ultrahigh- quality-factor micro- and nanomechanical resonators using dis- sipation dilution, Nat. Nanotechnol.19, 725 (2024)

  8. [8]

    Lifshitz and M

    R. Lifshitz and M. C. Cross, Nonlinear dynamics of nanome- chanicalandmicromechanicalresonators,inReviews of Nonlin- ear Dynamics and Complexity,editedbyH.G.Schuster(Wiley- VCH, Weinheim, 2008) pp. 1–52

  9. [9]

    J.S.AldridgeandA.N.Cleland,Noise-enabledprecisionmea- surements of a Duffing nanomechanical resonator, Phys. Rev. Lett.94, 156403 (2005)

  10. [10]

    Q. P. Unterreithmeier, T. Faust, and J. P. Kotthaus, Nonlinear switching dynamics in a nanomechanical resonator, Phys. Rev. B81, 241405 (2010)

  11. [11]

    H. B. Chan and C. Stambaugh, Activation barrier scaling and crossoverfornoise-inducedswitchinginmicromechanicalpara- metric oscillators, Phys. Rev. Lett.99, 060601 (2007)

  12. [12]

    H. B. Chan, M. I. Dykman, and C. Stambaugh, Paths of fluctu- ation induced switching, Phys. Rev. Lett.100, 130602 (2008)

  13. [13]

    W. J. Venstra, H. J. R. Westra, and H. S. J. van der Zant, Stochastic switching of cantilever motion, Nat. Commun.4, 2624 (2013)

  14. [14]

    Chowdhury, S

    A. Chowdhury, S. Barbay, M. G. Clerc, I. Robert-Philip, and R. Braive, Phase stochastic resonance in a forced nanoelec- tromechanical membrane, Phys. Rev. Lett.119, 234101 (2017)

  15. [15]

    J. S. Huber, G. Rastelli, M. J. Seitner, J. Kölbl, W. Belzig, M. I. Dykman, and E. M. Weig, Spectral evidence of squeezing of a weaklydampeddrivennanomechanicalmode,Phys.Rev.X10, 021066 (2020)

  16. [16]

    J. S. Ochs, M. Seitner, M. I. Dykman, and E. M. Weig, Ampli- fication and spectral evidence of squeezing in the response of a strongly driven nanoresonator to a probe field, Phys. Rev. A 103, 013506 (2021)

  17. [17]

    Almog, S

    R. Almog, S. Zaitsev, O. Shtempluck, and E. Buks, Noise squeezing in a nanomechanical Duffing resonator, Phys. Rev. Lett.98, 078103 (2007)

  18. [18]

    J. Suh, M. D. LaHaye, P. M. Echternach, K. C. Schwab, and M. L. Roukes, Parametric amplification and back-action noise squeezing by a qubit-coupled nanoresonator, Nano Lett.10, 3990 (2010)

  19. [19]

    Mahboob, E

    I. Mahboob, E. Flurin, K. Nishiguchi, A. Fujiwara, and H. Ya- maguchi, Enhanced force sensitivity and noise squeezing in an electromechanical resonator coupled to a nanotransistor, Appl. Phys. Lett.97, 253105 (2010)

  20. [20]

    Rugar and P

    D. Rugar and P. Grütter, Mechanical parametric amplification andthermomechanicalnoisesqueezing,Phys.Rev.Lett.67,699 (1991)

  21. [21]

    F.Yang,M.Fu,B.Bosnjak,R.H.Blick,Y.Jiang,andE.Scheer, Mechanically modulated sideband and squeezing effects of membrane resonators, Phys. Rev. Lett.127, 184301 (2021)

  22. [22]

    R. B. Karabalin, R. Lifshitz, M. C. Cross, M. H. Matheny, S. C. Masmanidis, and M. L. Roukes, Signal amplification by sensitive control of bifurcation topology, Phys. Rev. Lett.106, 094102 (2011)

  23. [23]

    N.KacemandS.Hentz,Bifurcationtopologytuningofamixed behavior in nonlinear micromechanical resonators, Appl. Phys. Lett.95, 183104 (2009)

  24. [24]

    P. M. Polunin, Y. Yang, M. I. Dykman, T. W. Kenny, and S. W. Shaw, Characterization of MEMS resonator nonlinearities us- ing the ringdown response, J. Microelectromech. Syst.25, 297 (2016)

  25. [25]

    Huang, S

    L. Huang, S. M. Soskin, I. A. Khovanov, R. Mannella, K. Nin- ios, and H. B. Chan, Frequency stabilization and noise-induced spectralnarrowinginresonatorswithzerodispersion,Nat.Com- mun.10, 3930 (2019)

  26. [26]

    J.M.L.Miller,A.Gomez-Franco,D.D.Shin,H.-K.Kwon,and T.W.Kenny,Amplitudestabilizationofmicromechanicaloscil- lators using engineered nonlinearity, Phys. Rev. Res.3, 033268 (2021)

  27. [27]

    Kaisar, J

    T. Kaisar, J. Lee, D. Li, S. W. Shaw, and P. X.-L. Feng, Non- linear stiffness and nonlinear damping in atomically thin MoS2 nanomechanical resonators, Nano Lett.22, 9831 (2022)

  28. [28]

    H.J.R.Westra,M.Poot,H.S.J.vanderZant,andW.J.Venstra, Nonlinear modal interactions in clamped-clamped mechanical resonators, Phys. Rev. Lett.105, 117205 (2010)

  29. [29]

    M. H. Matheny, L. G. Villanueva, R. B. Karabalin, J. E. Sader, andM.L.Roukes,Nonlinearmode-couplinginnanomechanical systems, Nano Lett.13, 1622 (2013)

  30. [30]

    Antonio, D

    D. Antonio, D. H. Zanette, and D. López, Frequency stabiliza- tioninnonlinearmicromechanicaloscillators,Nat.Commun.3, 806 (2012)

  31. [31]

    Eichler, M

    A. Eichler, M. del Álamo Ruiz, J. A. Plaza, and A. Bachtold, Strong coupling between mechanical modes in a nanotube res- onator, Phys. Rev. Lett.109, 025503 (2012)

  32. [32]

    C. Chen, D. H. Zanette, D. A. Czaplewski, S. Shaw, and D. López, Direct observation of coherent energy transfer in nonlinearmicromechanicaloscillators,Nat.Commun.8,15523 (2017)

  33. [33]

    Y. Yan, X. Dong, L. Huang, K. Moskovtsev, and H. B. Chan, Energy transfer into period-tripled states in coupled electrome- chanical modes at internal resonance, Phys. Rev. X12, 031003 (2022)

  34. [34]

    M. Wang, D. J. Perez-Morelo, D. Lopez, and V. A. Aksyuk, Persistent nonlinear phase-locking and nonmonotonic energy dissipation in micromechanical resonators, Phys. Rev. X12, 041025 (2022)

  35. [35]

    D.A.Czaplewski,C.Chen,D.Lopez,O.Shoshani,A.M.Eriks- son,S.Strachan,andS.W.Shaw,Bifurcationgeneratedmechan- ical frequency comb, Phys. Rev. Lett.121, 244302 (2018)

  36. [36]

    Eichler, J

    A. Eichler, J. Moser, J. Chaste, M. Zdrojek, I. Wilson-Rae, and A.Bachtold,Nonlineardampinginmechanicalresonatorsmade fromcarbonnanotubesandgraphene,Nat.Nanotechnol.6,339 (2011)

  37. [37]

    Zaitsev, O

    S. Zaitsev, O. Shtempluck, E. Buks, and O. Gottlieb, Nonlinear damping in a micromechanical oscillator, Nonlinear Dyn.67, 859 (2011)

  38. [38]

    D.Landauand E.M

    L. D.Landauand E.M. Lifshitz,Mechanics,3rd ed.,Courseof Theoretical Physics, Vol. 1 (Butterworth-Heinemann Limited, Oxford, 1976) Chap. §49

  39. [39]

    (Addison-Wesley, San Francisco, 2002) Chap

    H.Goldstein,C.P.Poole,andJ.L.Safko,Classical Mechanics, 3rd ed. (Addison-Wesley, San Francisco, 2002) Chap. 10.6

  40. [40]

    Duffing,Erzwungene Schwingungen bei veränderlicher Eigenfrequenz und ihre technische Bedeutung, Vol

    G. Duffing,Erzwungene Schwingungen bei veränderlicher Eigenfrequenz und ihre technische Bedeutung, Vol. Heft 41/42 (F. Vieweg, Braunschweig, 1918) 134 Seiten. 13

  41. [41]

    Samanta, N

    C. Samanta, N. Arora, and A. K. Naik, Tuning of geometric nonlinearity in ultrathin nanoelectromechanical systems, Appl. Phys. Lett.113, 113101 (2018)

  42. [42]

    Benedettini and G

    F. Benedettini and G. Rega, Non-linear dynamics of an elastic cable under planar excitation, Int. J. Non-Linear Mech.22, 497 (1987)

  43. [43]

    Eichler, J

    A. Eichler, J. Moser, M. I. Dykman, and A. Bachtold, Sym- metry breaking in a mechanical resonator made from a carbon nanotube, Nat. Commun.4, 2843 (2013)

  44. [44]

    Nabholz, L

    U. Nabholz, L. Lamprecht, J. E. Mehner, A. Zimmermann, and P.Degenfeld-Schonburg,Parametricamplificationofbroadband vibrational energy harvesters for energy-autonomous sensors enabled by field-induced striction, Mech. Syst. Signal Process. 139, 106642 (2020)

  45. [45]

    J. S. Ochs, G. Rastelli, M. Seitner, M. I. Dykman, and E. M. Weig,Resonantnonlinearresponseofananomechanicalsystem with broken symmetry, Phys. Rev. B104, 155434 (2021)

  46. [46]

    Keskekler, H

    A. Keskekler, H. Arjmandi-Tash, P. G. Steeneken, and F. Ali- jani,Symmetry-breaking-inducedfrequencycombsingraphene resonators, Nano Lett.22, 6048 (2022)

  47. [47]

    A.Keşkekler,V.Bos,A.M.Aragón,P.G.Steeneken,andF.Al- ijani, Multimode nonlinear dynamics of graphene resonators, Phys. Rev. Appl.20, 064020 (2023)

  48. [48]

    J. M. Londoño, S. A. Neild, and J. E. Cooper, Identification of backbone curves of nonlinear systems from resonance decay responses, J. Sound Vib.348, 224 (2015)

  49. [49]

    Q. P. Unterreithmeier, E. M. Weig, and J. P. Kotthaus, Univer- sal transduction scheme for nanomechanical systems based on dielectric forces, Nature458, 1001 (2009)

  50. [50]

    Faust, P

    T. Faust, P. Krenn, S. Manus, J. P. Kotthaus, and E. M. Weig, Microwave cavity-enhanced transduction for plug and playnanomechanicsatroomtemperature,Nat.Commun.3,728 (2012)

  51. [51]

    J.Rieger,T.Faust,M.J.Seitner,J.P.Kotthaus,andE.M.Weig, Frequency and Q factor control of nanomechanical resonators, Appl. Phys. Lett.101, 103110 (2012)

  52. [52]

    A. H. Nayfeh and D. T. Mook, Nonlinear oscillations (John Wiley & Sons, 2024) pp. 165, 197

  53. [53]

    I.KovacicandM.J.Brennan,TheDuffingequation: Nonlinear oscillators and their behaviour (John Wiley & Sons, 2011) pp. 97,281

  54. [54]

    J. S. Ochs, D. K. J. Boneß, G. Rastelli, M. Seitner, W. Belzig, M. I. Dykman, and E. M. Weig, Frequency comb from a sin- gle driven nonlinear nanomechanical mode, Phys. Rev. X12, 041019 (2022)

  55. [55]

    Antoni, K

    T. Antoni, K. Makles, R. Braive, T. Briant, P.-F. Cohadon, I. Sagnes, I. Robert-Philip, and A. Heidmann, Nonlinear me- chanics with suspended nanomembranes, EPL100, 68005 (2012)

  56. [56]

    Catalini, Y

    L. Catalini, Y. Tsaturyan, and A. Schliesser, Soft-clamped phononicdimersformechanicalsensingandtransduction,Phys. Rev. Appl.14, 014041 (2020)

  57. [57]

    Güttinger, A

    J. Güttinger, A. Noury, P. Weber, A. M. Eriksson, C. Lagoin, J. Moser, C. Eichler, A. Wallraff, A. Isacsson, and A. Bach- told, Energy-dependent path of dissipation in nanomechanical resonators, Nat. Nanotechnol.12, 631 (2017)

  58. [58]

    10 (CRC Press, 1961) p

    N.N.BogoliubovandI.Mitropolski,Asymptotic Methods in the Theory of Non-Linear Oscillations, Vol. 10 (CRC Press, 1961) p. 239

  59. [59]

    B. D. Hauer, C. Doolin, K. S. D. Beach, and J. P. Davis, A gen- eralprocedureforthermomechanicalcalibrationofnano/micro- mechanical resonators, Ann. Phys.339, 181 (2013)

This paper was first reviewed by deepseek-v4-flash on August 4, 2026.