Pith. sign in

REVIEW 3 major objections 4 minor 41 references

A gate-tuned graphene drum turns one drive tone into integer harmonics, a dense phononic comb, then a reverse period-doubling exit.

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-11 14:33 UTC pith:7PV3F6DL

load-bearing objection Clean experimental sequence of integer HHG, dense combs, and reverse period-doubling under gate-tuned 1:2 IR in graphene; the two-mode attribution is plausible but only qualitatively backed. the 3 major comments →

arxiv 2607.04724 v1 pith:7PV3F6DL submitted 2026-07-06 cond-mat.mes-hall physics.app-ph

Tunable Nonlinear Landscapes in Graphene Nanoelectromechanical Systems

Ateesh K. Rathi , Rajan Singh , Javed A. Mondal , Arnab Sarkar , Ryan J.T. Nicholl , Kirill I. Bolotin , Saikat Ghosh This is my paper
classification cond-mat.mes-hall physics.app-ph
keywords graphene NEMSphononic frequency combhigh-harmonic generation1:2 internal resonanceperiod doublingquadratic nonlinearitygate-tunable resonatorDuffing coefficient
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.

This paper shows that a single suspended monolayer graphene drum, under a gate that breaks out-of-plane symmetry and locks the fundamental flexural mode into a 1:2 internal resonance with a higher mode, can be driven by one tone into three successive regimes: phase-locked integer high-harmonic generation, a dense phononic frequency comb, and a reverse period-doubling transition that doubles the comb spacing and returns energy to even-order lines. The gate sets both the modal frequency ratio and the strength of the quadratic nonlinearity that couples the modes, so the same device can be switched among frequency multiplier, broadband comb source, and chaotic generator simply by changing gate or drive. Measured spectra fix the membrane's quadratic and cubic nonlinear coefficients. A sympathetic reader cares because the work turns the usually avoided nonlinear landscape of 2D NEMS into a controllable resource for classical analogs of extreme nonlinear optics and for phononic signal processing on one gated chip.

Core claim

Under a gate-defined 1:2 internal resonance that strengthens quadratic intermodal coupling, a single AC drive on a monolayer graphene drum produces sequential phase-locked integer harmonics that fill into a dense, phase-locked phononic frequency comb and then undergo a reverse period-doubling transition (comb spacing doubles, line density halves, energy returns to even-order lines). The same spectra yield the membrane coefficients ζ ≈ 1.16×10^9 N/m² and β ≈ 1.01×10^20 N/m³.

What carries the argument

Gate-broken out-of-plane symmetry that produces a large effective quadratic nonlinearity ζ (or intermodal coupling α) and places the fundamental flexural mode on a 1:2 internal resonance with a higher mode near 2ω0; this coupling, together with the cubic Duffing term β, drives the deterministic multi-wave mixing that yields integer harmonics, comb formation, and the reverse period-doubling cascade.

Load-bearing premise

The claim rests on the idea that a two-mode model with gate-induced quadratic coupling plus one Duffing term is enough to explain the integer harmonics, comb, and reverse period-doubling, rather than multi-mode interactions or measurement artifacts.

What would settle it

At gate voltages that hold the 1:2 condition, raise drive through the reported thresholds while recording broadband spectra: absence of sequential integer harmonics with the claimed power-law scalings, or failure of the comb to reverse-period-double with energy returning to even-order lines, would falsify the central attribution to the quadratic intermodal coupling.

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

If this is right

  • One gated graphene drum can serve as a zero-quiescent-power RF frequency multiplier whose usable harmonic is selected by re-gating rather than redesign.
  • The same device supplies a broadband phononic frequency comb whose line density can be walked deliberately via reverse period doubling.
  • Controlled entry into and exit from the chaotic comb maps onto a hardware route for bifurcation annealing in mechanical Ising machines.
  • Tracking the harmonics offers a multi-spectral handle for mass sensing on an atomically thin membrane.
  • The extracted ζ and β fix quantitative bounds for modeling and engineering nonlinear 2D NEMS.

Where Pith is reading between the lines

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

  • If the reverse period-doubling exit is as sharp as reported, the device could function as a controllable chaos-to-periodic switch for phononic logic without needing separate resonators.
  • The gate-tunable quadratic landscape may generalize to other 2D membranes, suggesting a materials-agnostic design rule for on-chip high-harmonic generation.
  • Simultaneous multi-harmonic output under a single drive could simplify multi-frequency metrology or sensing architectures that currently require multiple sources.

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 manuscript reports that a gate-tuned monolayer graphene drum, driven into a 1:2 internal resonance between the fundamental flexural mode and a higher mode near 2ω₀, produces phase-locked integer high harmonics under a single drive tone. At higher drive the harmonics fill into a dense phononic frequency comb; further increase yields a reverse period-doubling transition in which comb spacing doubles, line density halves, and energy returns to even-order lines. Spectral maps versus drive amplitude and frequency (Figs. 2–3, S.2–S.3), power-law scalings A₂∝V_ac², A₃∝V_ac³, A₄∝V_ac⁸ in the integer-harmonic window, and supporting two-mode harmonic-balance plus Runge–Kutta simulations are used to attribute the sequence to gate-broken quadratic intermodal coupling. The same spectra are inverted for the membrane quadratic and cubic coefficients ζ ≈ 1.16×10⁹ N/m² and β ≈ 1.01×10²⁰ N/m³.

Significance. If the dynamical sequence and its attribution hold, the work supplies a concrete, gate-selectable route from integer high-harmonic generation through dense phononic combs to a reverse period-doubling exit on a single atomically thin resonator. That combination is of clear interest for zero-quiescent-power RF multiplication, multi-spectral sensing, and hardware bifurcation annealing. The experimental maps and power-law onsets are a genuine advance over earlier graphene-comb reports that emphasized Neimark–Sacker quasi-periodicity. The supporting analytic reduction (SI harmonic balance under 1:2 IR) and the explicit extraction formulas for ζ and β are useful and falsifiable once placed on a quantitative footing.

major comments (3)
  1. [Supporting Information, Extraction of the Nonlinear Coefficients] SI “Extraction of the Nonlinear Coefficients,” Eqs. (S.22) and (S.27): ζ and β are obtained from single-mode force-balance formulas that treat the second and third harmonics as driven solely by ζ x² and β x³ on an isolated fundamental. The main text and the preceding SI section, however, establish that under the same gate conditions the device sits on a strong 1:2 internal resonance, so A₂ is resonantly enhanced by the higher mode (Eqs. S.9, S.16; β_eff diverges as ω₂ o 2ω_d). Using the off-resonance single-mode expressions therefore undercuts the quantitative claim that the measured spectra yield the membrane nonlinear landscape. Either re-derive the extraction under the two-mode harmonic-balance solution already given in the SI, or restrict the extraction to gate voltages demonstrably off the 1:2 condition and show consistency.
  2. [Supporting Information, Numerical Simulations] SI Numerical Simulations, Eqs. (S.1)–(S.2)/(S.6)–(S.7) and Figs. S.4–S.5: the Runge–Kutta maps are stated to “reproduce the experimental progression,” yet no measured parameters (extracted ζ or α from modal overlap, γ_i, m_i, or the observed Δ_c and reverse-doubling thresholds) are inserted, and no experimental thresholds or comb spacings are overlaid. Without at least one quantitative match of onset voltages, Δ_c, or the reverse-doubling point, the two-mode model remains an illustration rather than a validated attribution of the reverse period-doubling cascade. A minimal fix is to fix α, β from the data (or from the corrected extraction above) and show that the simulated thresholds and Δ_c lie within experimental uncertainty.
  3. [Discussion; Figure 1C] Discussion and Fig. 1C: graphene drums support a dense spectrum of flexural modes. The reverse period-doubling cascade and the dense comb could in principle involve additional near-resonant modes or higher-order resonances not captured by the two-mode truncation. The manuscript should either (i) show that neighboring modes remain detuned and weakly driven throughout the mapped windows, or (ii) acknowledge multi-mode participation as an open alternative and state what additional measurement would rule it out. This is load-bearing for the claim that the observed sequence is controlled by the gate-tuned quadratic intermodal coupling alone.
minor comments (4)
  1. [Figure 2D] Figure 2D and SI Fig. S.2B: the V_ac^8 scaling quoted for A₄ is unusually steep; a short note on the cascade order that produces it (or a fit residual) would help the reader judge whether the power is robust or an effective local slope.
  2. [Eq. (1); SI Eqs. (S.5)–(S.7)] Main text Eq. (1) writes the potential with coefficients K, ζ, β while the SI equations of motion use m, α, β; a single consistent notation table (or explicit conversion α ↔ ζ ∫ φ₁² φ₂ dA) would remove ambiguity.
  3. [Period-Doubling Cascade and Comb Spacing] The Feigenbaum cascade is identified by the sequence Δ_c o Δ_c/2 o Δ_c/4 o chaos, but the text correctly notes that resolution is insufficient for δ. Soften the language from “matches the universal period-doubling route” to “consistent with” to avoid over-claim.
  4. [Abstract; Figure 3B] Abstract and Introduction list “chaotic generator” among the device functions; the data show a cascade into broadband spectra labeled “Chaos,” but no Lyapunov or correlation-dimension diagnostic is given. A brief caveat that chaos is inferred from the spectral hierarchy would be appropriate.

Circularity Check

0 steps flagged

No circular derivation: experimental spectra under independently measured 1:2 IR yield coefficients via standard force balance; model is interpretive, not self-forcing.

full rationale

The paper's load-bearing claims are direct experimental observations (integer harmonics with measured power-law onsets, dense combs with spacing Δc, reverse period-doubling to 2Δc) obtained by gating the measured modal dispersion (Fig. 1C) onto ω2 ≈ 2ω0 and increasing Vac. The 1:2 condition is set by independent linear-regime frequency tracking, not by the nonlinear coefficients later extracted. ζ and β are obtained post-hoc from measured |A2|/|A1|^{2} and |A3|/|A1|^{3} via the textbook single-mode harmonic-balance formulas (SI Eqs. S.22, S.27); they are not inputs that are then re-sold as predictions of the same spectra. The two-mode equations (S.1–S.2 / S.6–S.7) and qualitative Runge–Kutta maps serve only as consistency checks that reproduce the sequence of regimes; no fitted α, eta or thresholds are used to “predict” the data that defined them. Self-citations to the group’s prior graphene NEMS work are background and do not supply uniqueness theorems or ansatzes that force the present spectral features. The derivation chain is therefore self-contained against external benchmarks and contains no self-definitional, fitted-as-prediction, or load-bearing self-citation reductions.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The central claims rest on continuum membrane mechanics with gate-induced symmetry breaking, a two-mode 1:2 IR reduction, and standard harmonic-balance extraction of nonlinear coefficients from measured amplitudes. Free parameters are the extracted ζ and β (and effective modal parameters implicit in the fits). No new particles or forces are postulated; reverse period-doubling is an observed dynamical regime, not an invented entity.

free parameters (3)
  • quadratic coefficient ζ = 1.16e9 N/m^2
    Extracted from measured |A2| and |A1| via SI Eq. (S.22); value 1.16×10^9 N/m² is data-dependent and enters the claim that the gate drives a strong χ⁽²⁾-like regime.
  • cubic (Duffing) coefficient β = 1.01e20 N/m^3
    Extracted from measured |A3| and |A1| via SI Eq. (S.27); value 1.01×10^20 N/m³ is data-dependent and used to quantify the nonlinear landscape.
  • intermodal quadratic coupling α (and modal masses/dampings in simulations)
    α ≡ ζ ∫ φ1² φ2 dA and γi, βi, mi in the RK maps (SI Eqs. S.1–S.2) are chosen/effective parameters that reproduce qualitative spectral maps; not independently measured for every figure.
axioms (4)
  • domain assumption Gate-induced out-of-plane displacement breaks membrane symmetry and generates an effective quadratic force term ζX² (potential term ⅓ζX³) that couples modes at 1:2 resonance.
    Stated in Introduction/Experimental Setup and SI ‘Symmetry breaking and the intermodal coupling’; standard for buckled/gated membranes but load-bearing for attributing HHG to quadratic three-wave mixing.
  • domain assumption Two-mode truncation (fundamental + higher mode near 2ω0) with coupling terms 2α x1 x2 and α x1² plus Duffing on the driven mode captures the observed harmonic, comb, and period-doubling spectra.
    SI Eqs. (S.6)–(S.7) and numerical maps; multi-mode content is not exhaustively excluded experimentally.
  • standard math Single-term harmonic balance and leading-order cascade scalings relate measured An to ζ and β (and explain A2∝V², A3∝V³, A4∝V⁸).
    SI harmonic-balance section and main-text power-law fits; standard nonlinear-oscillator technique.
  • domain assumption Spectral line-spacing halvings Δc → Δc/2 → Δc/4 followed by broadband fill constitute a Feigenbaum period-doubling route to chaos, and the later doubling of spacing is reverse period-doubling.
    Main text Period-Doubling Cascade section; Feigenbaum constant not measured; identification is spectral only.

pith-pipeline@v1.1.0-grok45 · 17886 in / 3679 out tokens · 31403 ms · 2026-07-11T14:33:08.932110+00:00 · methodology

0 comments
read the original abstract

Nonlinear nanomechanical resonators give convenient solid-state access to classical analogs of extreme nonlinear optics and to phononic signal processing. Here we report integer high-harmonic generation and phononic frequency combs in a suspended monolayer graphene drum. A gate voltage breaks the out-of-plane symmetry of the membrane and tunes its fundamental flexural mode onto a 1:2 internal resonance with a higher mode at twice the frequency, where the quadratic coupling between the two modes becomes large. A single drive tone then generates phase-locked integer harmonics in sequence, and at larger drive these fill in to a dense frequency comb. Raising the drive further, we find a reverse period-doubling transition: the comb spacing doubles, the line density halves, and energy flows back into the even-order comb lines. The measured spectra yield the quadratic ($\zeta$) and cubic ($\beta$) nonlinear coefficients of the membrane. These results show how the tunable nonlinear landscape of graphene supports distinct dynamical regimes on demand, allowing a single gated device to act in turn as a frequency multiplier, a broadband comb source, and a chaotic generator.

Figures

Figures reproduced from arXiv: 2607.04724 by Arnab Sarkar, Ateesh K. Rathi, Javed A. Mondal, Kirill I. Bolotin, Rajan Singh, Ryan J.T. Nicholl, Saikat Ghosh.

Figure 1
Figure 1. Figure 1: Device, effective potential, and gate-tunable 1:2 internal resonance. (A)(a) Scanning electron micrograph of the device. (b) Capacitive geometry: an AC voltage Vac drives the membrane and a DC gate voltage Vdc sets the static tension. (B) Effective potential V (X) of the fundamental mode, showing the competing quadratic and cubic nonlinear terms. (C) Dispersion of the mechanical resonances versus Vdc. (D) … view at source ↗
Figure 2
Figure 2. Figure 2: Integer high-harmonic generation and amplitude scaling. (A) Fundamental-mode spectrum with increasing drive amplitude at Vdc = 172 V; the spectral dip deepens as the nonlinear mode coupling strengthens. (B) Spectrum at fixed drive amplitude, showing the mode splitting from internal resonance. The black vertical line marks the drive frequency used in the following measurements. (C) Broadband spectrum under … view at source ↗
Figure 3
Figure 3. Figure 3: From integer harmonics through frequency combs to chaos. (A) Spectral map versus drive amplitude at fixed ωd/2π = 4.8338 MHz and Vdc = 158 V. Three regimes are marked: (i) integer-harmonic generation (Vac ≤ 6.8 V); (ii) dense comb (6.8 ≤ Vac ≤ 8.8 V); (iii) reverse period-doubled comb with doubled spacing (8.8 ≤ Vac ≤ 9.4 V). (B) Spectral map versus drive frequency at fixed drive amplitude, showing the evo… 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

41 extracted references · 12 canonical work pages

  1. [1]

    Lifshitz and M

    R. Lifshitz and M. C. Cross, Nonlinear dynamics of nanomechanical and micromechanical resonators, in Reviews of Nonlinear Dynamics and Complexity (John Wiley & Sons, Ltd, 2008) Chap. 1, pp. 1–52

    2008

  2. [2]

    Bachtold, J

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

    2022

  3. [3]

    Tadokoro, H

    Y. Tadokoro, H. Tanaka, and M. I. Dykman, Driven nonlinear nanomechanical resonators as digital signal detectors, Scientific Reports8, 10.1038/s41598-018-29572-7 (2018)

    work page doi:10.1038/s41598-018-29572-7 2018

  4. [4]

    Maillet, X

    O. Maillet, X. Zhou, R. R. Gazizulin, R. Ilic, J. M. Parpia, O. Bourgeois, A. D. Fefferman, and E. Collin, Measuring frequency fluctuations in nonlinear nanomechanical resonators, ACS Nano12, 5753–5760 (2018)

    2018

  5. [5]

    Weber, J

    P. Weber, J. G¨ uttinger, A. Noury, J. Vergara-Cruz, and A. Bachtold, Force sensitivity of multilayer graphene optomechan- ical devices, Nature Communications7, 12496 (2016)

    2016

  6. [6]

    Manzaneque, M

    T. Manzaneque, M. K. Ghatkesar, F. Alijani, M. Xu, R. A. Norte, and P. G. Steeneken, Resolution limits of resonant sensors, Phys. Rev. Appl.19, 054074 (2023). 7

    2023

  7. [7]

    Davidovikj, F

    D. Davidovikj, F. Alijani, S. J. Cartamil-Bueno, H. S. J. van der Zant, M. Amabili, and P. G. Steeneken, Nonlinear dynamic characterization of two-dimensional materials, Nature Communications8, 10.1038/s41467-017-01351-4 (2017)

    doi:10.1038/s41467-017-01351-4 2017

  8. [8]

    X. Fan, C. He, J. Ding, Q. Gao, H. Ma, M. C. Lemme, and W. Zhang, Graphene mems and nems, Microsystems & Nanoengineering10, 10.1038/s41378-024-00791-5 (2024)

    work page doi:10.1038/s41378-024-00791-5 2024

  9. [9]

    Chen, Graphene nanoelectromechanical resonators and oscillators, Ph.D

    C. Chen, Graphene nanoelectromechanical resonators and oscillators, Ph.D. thesis, Columbia University (2013)

    2013

  10. [10]

    S. O. Erbil, U. Hatipoglu, C. Yanik, M. Ghavami, A. B. Ari, M. Yuksel, and M. S. Hanay, Full electrostatic control of nanomechanical buckling, Phys. Rev. Lett.124, 046101 (2020)

    2020

  11. [11]

    Singh, R

    R. Singh, R. J. Nicholl, K. I. Bolotin, and S. Ghosh, Motion transduction with thermo-mechanically squeezed graphene resonator modes, Nano Letters18, 6719–6724 (2018)

    2018

  12. [12]

    Mahboob, K

    I. Mahboob, K. Nishiguchi, A. Fujiwara, and H. Yamaguchi, Phonon lasing in an electromechanical resonator, Phys. Rev. Lett.110, 127202 (2013)

    2013

  13. [13]

    Ganesan, C

    A. Ganesan, C. Do, and A. Seshia, Phononic frequency comb via intrinsic three-wave mixing, Phys. Rev. Lett.118, 033903 (2017)

    2017

  14. [14]

    Singh, A

    R. Singh, A. Sarkar, C. Guria, R. J. Nicholl, S. Chakraborty, K. I. Bolotin, and S. Ghosh, Giant tunable mechanical nonlinearity in graphene–silicon nitride hybrid resonator, Nano Letters20, 4659–4666 (2020)

    2020

  15. [15]

    G. Ye, R. Sun, J. Zhao, and F. Ma, Magnetostrictive mechanical frequency combs, Nature Communications16, 10.1038/s41467-025-64878-x (2025)

    work page doi:10.1038/s41467-025-64878-x 2025

  16. [16]

    Y. He, T. Kuang, X. Han, Z. Feng, X. Chen, W. Xiong, S. Jin, Z. Tan, Q. Zhang, H. Luo, H. Jing, and G. Xiao, Coherent acoustic frequency comb via floquet engineering of optical tweezer phonon lasers, Science Advances11, eadv9984 (2025), https://www.science.org/doi/pdf/10.1126/sciadv.adv9984

    work page doi:10.1126/sciadv.adv9984 2025

  17. [17]

    Sarkar, Anurag, J

    A. Sarkar, Anurag, J. A. Mondal, R. Singh, A. A. Makki, A. K. Rathi, R. J. Nicholl, S. Chakraborty, K. I. Bolotin, and S. Ghosh, Observation of tunable discrete-time-crystalline phases, Phys. Rev. Appl.25, L021001 (2026)

    2026

  18. [18]

    T. J. Kippenberg, A. L. Gaeta, M. Lipson, and M. L. Gorodetsky, Dissipative kerr solitons in optical microresonators, Science361, 10.1126/science.aan8083 (2018)

    doi:10.1126/science.aan8083 2018

  19. [19]

    Antonio, D

    D. Antonio, D. H. Zanette, and D. L´ opez, Frequency stabilization in nonlinear micromechanical oscillators, Nature Com- munications3, 10.1038/ncomms1813 (2012)

    work page doi:10.1038/ncomms1813 2012

  20. [20]

    Houri, D

    S. Houri, D. Hatanaka, M. Asano, and H. Yamaguchi, Demonstration of multiple internal resonances in a microelectrome- chanical self-sustained oscillator, Phys. Rev. Appl.13, 014049 (2020)

    2020

  21. [21]

    S. M. E. H. Yousuf, J. Lee, S. W. Shaw, and P. X.-L. Feng, Phononic frequency combs in atomically thin nanoelectrome- chanical resonators via 1:1 and 2:1 internal resonances, Journal of Microelectromechanical Systems32, 335 (2023)

    2023

  22. [22]

    Mouharrar, S

    H. Mouharrar, S. Rahmanian, R. Abdelrahman, Y. S. Shama, M. Akbari, S. Basrour, K. Musselman, D. Mu˜ noz-Rojas, M. Yavuz, and E. Abdel-Rahman, Generation of soliton frequency combs in nems, Nano Letters24, 10834–10841 (2024)

    2024

  23. [23]

    C. Chen, D. H. Zanette, D. A. Czaplewski, S. Shaw, and D. L´ opez, Direct observation of coherent energy transfer in nonlinear micromechanical oscillators, Nature Communications8, 10.1038/ncomms15523 (2017)

    work page doi:10.1038/ncomms15523 2017

  24. [24]

    Zhang, H

    H. Zhang, H. Li, J. Sun, S. Kirkbride, G. Teng, Z. Liu, D. Chen, M. Parajuli, M. Pandit, G. Sobreviela, C. Zhao, W. Yuan, H. Chang, and A. A. Seshia, Coherent energy transfer in coupled nonlinear microelectromechanical resonators, Nature Communications16, 10.1038/s41467-025-59292-2 (2025)

    work page doi:10.1038/s41467-025-59292-2 2025

  25. [25]

    J. Wu, S. Zang, P. Song, W. Zhang, and L. Shao, Giant energy exchange rate in mode-coupled resonators enables super- continuum mechanical frequency combs, Microsystems & Nanoengineering12, 10.1038/s41378-026-01168-6 (2026)

    work page doi:10.1038/s41378-026-01168-6 2026

  26. [26]

    D. A. Czaplewski, C. Chen, D. Lopez, O. Shoshani, A. M. Eriksson, S. Strachan, and S. W. Shaw, Bifurcation generated mechanical frequency comb, Phys. Rev. Lett.121, 244302 (2018)

    2018

  27. [27]

    Asadi, J

    K. Asadi, J. Yeom, and H. Cho, Strong internal resonance in a nonlinear, asymmetric microbeam resonator, Microsystems & Nanoengineering7, 10.1038/s41378-020-00230-1 (2021)

    work page doi:10.1038/s41378-020-00230-1 2021

  28. [28]

    Ke¸ skekler, H

    A. Ke¸ skekler, H. Arjmandi-Tash, P. G. Steeneken, and F. Alijani, Symmetry-breaking-induced frequency combs in graphene resonators, Nano Letters22, 6048–6054 (2022)

    2022

  29. [29]

    J. Sun, S. Yu, H. Zhang, D. Chen, X. Zhou, C. Zhao, D. D. Gerrard, R. Kwon, G. Vukasin, D. Xiao, T. W. Kenny, X. Wu, and A. Seshia, Generation and evolution of phononic frequency combs via coherent energy transfer between mechanical modes, Phys. Rev. Appl.19, 014031 (2023)

    2023

  30. [30]

    Brabec and F

    T. Brabec and F. Krausz, Intense few-cycle laser fields: Frontiers of nonlinear optics, Rev. Mod. Phys.72, 545 (2000)

    2000

  31. [31]

    E. J. Takahashi, P. Lan, O. D. M¨ ucke, Y. Nabekawa, and K. Midorikawa, Attosecond nonlinear optics using gigawatt-scale isolated attosecond pulses, Nature Communications4, 10.1038/ncomms3691 (2013)

    work page doi:10.1038/ncomms3691 2013

  32. [32]

    P. S. Linsay, Period doubling and chaotic behavior in a driven anharmonic oscillator, Phys. Rev. Lett.47, 1349 (1981)

    1981

  33. [33]

    Mahboob, R

    I. Mahboob, R. Dupuy, K. Nishiguchi, A. Fujiwara, and H. Yamaguchi, Hopf and period-doubling bifurcations in an electromechanical resonator, Applied Physics Letters109, 10.1063/1.4960735 (2016)

    work page doi:10.1063/1.4960735 2016

  34. [34]

    Eichler, J

    A. Eichler, J. Moser, M. Dykman, and A. Bachtold, Symmetry breaking in a mechanical resonator made from a carbon nanotube, Nature Communications4, 10.1038/ncomms3843 (2013)

    work page doi:10.1038/ncomms3843 2013

  35. [35]

    J. S. Ochs, G. Rastelli, M. Seitner, M. I. Dykman, and E. M. Weig, Resonant nonlinear response of a nanomechanical system with broken symmetry, Phys. Rev. B104, 155434 (2021)

    2021

  36. [36]

    M. J. Feigenbaum, Quantitative universality for a class of nonlinear transformations, Journal of Statistical Physics19, 25–52 (1978)

    1978

  37. [37]

    J. Wu, P. Song, S. Zang, Z. Mao, W. Zhang, and L. Shao, Self-injection locked and phase offset-free micromechanical frequency combs, Phys. Rev. Lett.134, 107201 (2025)

    2025

  38. [38]

    H. Wang, D. Nezich, J. Kong, and T. Palacios, Graphene frequency multipliers, IEEE Electron Device Letters30, 547–549 (2009). 8

    2009

  39. [39]

    T. L. Heugel, O. Zilberberg, C. Marty, R. Chitra, and A. Eichler, Ising machines with strong bilinear coupling, Phys. Rev. Res.4, 013149 (2022). XXX Supporting Information: Tunable Nonlinear Landscapes in Graphene Nanoelectromechanical Systems Ateesh K. Rathi, 1 Rajan Singh, 1 Javed A. Mondal, 1 Arnab Sarkar,1 Ryan J.T. Nicholl, 2 Kirill I. Bolotin, 3 and...

    2022

  40. [40]

    (S.8) in Eq

    Balancing the cos 2(ωdt+φ) components and dropping the off-resonant damping term, k2 −4m 2ω2 d X2 + 1 2 αX2 1 = 0, so X2 = αX2 1 2m2 4ω2 d −ω 2 2 , ω 2 2 = k2 m2 .(S.9) For the fundamental, Eq. (S.8) in Eq. (S.6) with the resonant parts of the nonlinear terms, x3 1 =X 3 1 h 3 4 cos(ωdt+φ) + 1 4 cos 3(ωdt+φ) i ,(S.10) x1x2 =X 1X2 cos(ωdt+φ) cos 2(ωdt+φ) = ...

  41. [41]

    15 Cubic Nonlinearityβ For the cubic term, the equation of motion is m¨x+mγ˙x+mω 2 0x+βx 3 =Fcos(ω dt).(S.23) A cubic nonlinearity drives a response at 3ω d

    Balancing terms at 2ω d, A2 = ζA2 1 m 4ω2 d −ω 2 0 +i2γω d ,(S.20) and taking the magnitude, |A2|= |ζ| |A1|2 m q 4ω2 d −ω 2 0 2 + 2γωd 2 .(S.21) The quadratic coefficient is therefore ζ= m|A 2| |A1|2 q 4ω2 d −ω 2 0 2 + 2γωd 2 .(S.22) Substituting the measured amplitudes and parameters givesζ= 1.16×10 9 N/m2. 15 Cubic Nonlinearityβ For the cubic term, the ...