pith. sign in

arxiv: 2602.22513 · v2 · submitted 2026-02-26 · ⚛️ physics.class-ph · nlin.AO· physics.app-ph

Nonlinear stabilization of chiral modes in space-time modulated parametric oscillators

Scott Lambert , Elise Jaremko , Jayson Paulose This is my paper

Pith reviewed 2026-05-15 19:33 UTC · model grok-4.3

classification ⚛️ physics.class-ph nlin.AOphysics.app-ph
keywords parametric oscillatorschiral modesspace-time modulationnonlinear stabilizationaveraged equationscubic nonlinearitynonreciprocal amplificationelastic resonators
0
0 comments X

The pith

Cubic nonlinearity arrests growth of chiral modes in space-time modulated parametric oscillators, yielding steady finite-amplitude motion.

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

In a trio of coupled parametric oscillators whose modulation phases are chosen to amplify a circulating chiral mode in the linear regime, a cubic nonlinearity halts the exponential growth. The resulting motion reaches a steady finite amplitude while preserving the expected chirality. Space-time symmetry of the modulation reduces the full three-oscillator dynamics to a single averaged equation that predicts the nonlinear trajectories, steady amplitudes, and time scales. These chiral states have finite basins of attraction and remain reachable over wide ranges of initial conditions and parameters. Finite-element simulations of elastic-plate resonators reproduce the same features, confirming that the reduced description applies to realistic continuum systems.

Core claim

We establish the existence of nonlinear chiral steady states in a trio of coupled parametric oscillators with modulation phases chosen to selectively amplify a circulating mode in the linearized system. We find that a cubic nonlinearity arrests exponential growth of the amplified mode, producing a steady finite-amplitude motion that retains the expected chirality. By exploiting space-time symmetry, we reduce the dynamics to a single averaged equation that quantitatively predicts nonlinear trajectories, steady-state amplitudes, and characteristic time scales. The chiral steady states possess finite basins of attraction and are accessible from wide ranges of initial conditions and system par

What carries the argument

Space-time symmetry of the chosen modulation phases, which reduces the three-oscillator system to one averaged equation that incorporates the cubic nonlinearity and predicts finite-amplitude chiral motion.

If this is right

  • Chiral steady states with finite basins of attraction become stable and reachable from broad ranges of initial conditions and parameters.
  • Linear features such as chirality and directional amplification persist into strongly nonlinear regimes.
  • The reduced averaged equation supplies explicit predictions for amplitudes and characteristic time scales once nonlinearity is significant.
  • Finite-element models of continuum elastic resonators reproduce the same chiral steady states, showing the reduction applies beyond discrete oscillators.

Where Pith is reading between the lines

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

  • The same symmetry-reduction approach could be tested in optical or acoustic parametric resonators to achieve nonreciprocal routing at larger scales.
  • If the modulation phases are altered while preserving the space-time symmetry, the averaged equation should still close and predict analogous chiral states.
  • Adding higher-order nonlinearities would likely shift the steady amplitude but leave the chirality intact if the symmetry reduction remains valid.

Load-bearing premise

The space-time symmetry of the chosen modulation phases permits an exact reduction of the three-oscillator system to a single averaged equation whose predictions remain quantitatively accurate once the cubic nonlinearity becomes order-one.

What would settle it

If measured steady-state amplitudes or oscillation trajectories in simulations or experiments deviate substantially from the predictions of the single averaged equation when the cubic term is order-one, the quantitative accuracy of the symmetry reduction would be falsified.

Figures

Figures reproduced from arXiv: 2602.22513 by Elise Jaremko, Jayson Paulose, Scott Lambert.

Figure 1
Figure 1. Figure 1: FIG. 1. The “trimer” model and its normal modes. Each col [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Floquet spectra of the trimer at various parameter [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Numerical solution of the full equations of mo [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a): Phase portraits of the reduced averaged equa [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Rescaled trimer steady-state amplitude from numer [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Phase space radius [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Chiral steady states in a continuum model of mechanical parametric oscillators. (a) Trimer of plate resonators arranged [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 1
Figure 1. Figure 1: However, COMSOL does not decompose the degen￾erate subspace into the clockwise and counterclockwise￾rotating modes, but rather into arbitrary linear combi￾nations of them with the degeneracy lifted by the mesh. To recreate the dynamics of the discrete model in the elastic system, we modulated the in-plane prestresses within each of the three quadrilateral domains (indexed by j) as σj (t) = σ0  1 + δ ′ cos… view at source ↗
read the original abstract

Phase control of parametric modulation in coupled oscillator networks enables sculpting of dynamical states with desired spatiotemporal symmetries. Symmetry-aware Floquet analysis successfully predicts such states in linear systems, but whether their symmetry properties persist under nonlinearity remains largely unexplored. Here, we establish the existence of nonlinear chiral steady states in a trio of coupled parametric oscillators with modulation phases chosen to selectively amplify a circulating mode in the linearized system. We find that a cubic nonlinearity arrests exponential growth of the amplified mode, producing a steady finite-amplitude motion that retains the expected chirality. By exploiting space-time symmetry, we reduce the dynamics to a single averaged equation that quantitatively predicts nonlinear trajectories, steady-state amplitudes, and characteristic time scales. The chiral steady states possess finite basins of attraction and are accessible from wide ranges of initial conditions and system parameters. Finite-element simulations of elastic plate resonators quantitatively reproduce these features, establishing the relevance of the reduced model to realistic continuum systems. Our results demonstrate that desirable properties of linear time-modulated systems, such as chirality and directional amplification, persist into strongly nonlinear regimes, opening pathways to robust nonreciprocal signal routing and amplification in parametrically driven platforms.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper claims that in a system of three coupled parametric oscillators with space-time modulation phases chosen to amplify a circulating chiral mode, a cubic nonlinearity arrests exponential growth and produces a finite-amplitude steady state that retains chirality. Exploiting the space-time symmetry, the dynamics reduce to a single averaged equation whose predictions for nonlinear trajectories, steady-state amplitudes, and timescales are asserted to be quantitatively accurate and to match finite-element simulations of elastic plate resonators, showing that linear symmetry properties persist into the nonlinear regime.

Significance. If the quantitative accuracy of the symmetry-based reduction holds beyond the perturbative regime, the work would demonstrate that desirable features such as chirality and directional amplification from linear Floquet analysis can be preserved under strong nonlinearity. This has potential implications for nonreciprocal signal routing in parametrically driven mechanical systems. The exact reduction enabled by the imposed symmetry and the direct comparison to continuum simulations are notable strengths that support the central claim.

major comments (2)
  1. [Derivation of the reduced averaged equation] The reduction to a single averaged equation is presented as following directly from the space-time symmetry, yet the claim that this equation remains quantitatively predictive once the cubic nonlinearity reaches order-one strength requires explicit justification. Standard averaging methods assume small parameters; the manuscript should provide an error bound, residual-term analysis, or regime of validity (in the section deriving the averaged equation) to address possible higher-harmonic corrections or residual couplings that could appear when the nonlinearity is no longer perturbative.
  2. [Finite-element validation and results] The abstract and results sections assert that the reduced model quantitatively reproduces trajectories, amplitudes, and timescales in finite-element simulations, but without reported error metrics (e.g., relative L2 discrepancies), the specific parameter values employed in the comparisons, or discussion of any adjustments, the strength of this agreement cannot be fully assessed.
minor comments (2)
  1. [Figures] Figure captions should explicitly list the numerical values of the cubic nonlinearity coefficient, modulation depth, and coupling strengths used in each simulation panel to facilitate reproducibility.
  2. [Notation and equations] Clarify the notation for the modulation phases across the text, equations, and figures to avoid potential ambiguity in the symmetry arguments.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive suggestions. The comments highlight opportunities to strengthen the justification of the symmetry reduction and to make the validation more quantitative. We address each point below and will incorporate the requested clarifications and metrics into the revised manuscript.

read point-by-point responses

  1. Referee: [Derivation of the reduced averaged equation] The reduction to a single averaged equation is presented as following directly from the space-time symmetry, yet the claim that this equation remains quantitatively predictive once the cubic nonlinearity reaches order-one strength requires explicit justification. Standard averaging methods assume small parameters; the manuscript should provide an error bound, residual-term analysis, or regime of validity (in the section deriving the averaged equation) to address possible higher-harmonic corrections or residual couplings that could appear when the nonlinearity is no longer perturbative.

    Authors: The space-time symmetry of the modulation phases enforces an exact decoupling: the two counter-circulating modes remain identically zero when the system is initialized in the chiral subspace, even for finite-amplitude cubic nonlinearity. Consequently the reduction to a single complex amplitude equation is exact (not perturbative) before any averaging is applied. Averaging is used only to obtain the slow envelope dynamics of this isolated mode. In the revision we will add a dedicated subsection that (i) derives the exact symmetry reduction without invoking small-parameter assumptions, (ii) states the residual higher-harmonic terms that are neglected by averaging, and (iii) supplies numerical checks of those residuals over the range of nonlinearity strengths explored in the paper, thereby delineating the regime of quantitative validity. revision: yes

  2. Referee: [Finite-element validation and results] The abstract and results sections assert that the reduced model quantitatively reproduces trajectories, amplitudes, and timescales in finite-element simulations, but without reported error metrics (e.g., relative L2 discrepancies), the specific parameter values employed in the comparisons, or discussion of any adjustments, the strength of this agreement cannot be fully assessed.

    Authors: We agree that explicit error metrics and parameter transparency are required. In the revised manuscript we will (i) report relative L2 discrepancies between the reduced-model trajectories and the finite-element time series for both displacement and velocity fields, (ii) list the precise numerical values of modulation depth, detuning, cubic coefficient, and damping used in each comparison, and (iii) describe any minor mesh-convergence adjustments that were made to ensure the continuum model lies within the regime where the three-mode truncation remains accurate. revision: yes

Circularity Check

0 steps flagged

Symmetry reduction and cubic model term are independent of target predictions

full rationale

The paper imposes space-time modulation phases that enforce symmetry, then derives an exact reduction of the three-oscillator system to one averaged equation. The cubic nonlinearity is added as a phenomenological term, not fitted to the steady-state amplitudes or trajectories. No self-citations appear in the provided text as load-bearing for the reduction or quantitative claims. The predictions follow directly from solving the reduced equation under the stated symmetry, making the derivation self-contained rather than circular.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the averaging reduction under space-time symmetry and the assumption that cubic nonlinearity dominates saturation without introducing new instabilities or breaking chirality.

free parameters (1)
  • cubic nonlinearity coefficient
    The strength of the stabilizing cubic term is a model parameter whose value determines the saturated amplitude; its specific numerical value is not given in the abstract.
axioms (1)
  • domain assumption Space-time symmetry of the modulation permits exact reduction to a single averaged equation
    Invoked to collapse the three-oscillator dynamics and obtain quantitative predictions for amplitude and time scales.

pith-pipeline@v0.9.0 · 5504 in / 1369 out tokens · 43171 ms · 2026-05-15T19:33:38.990918+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

78 extracted references · 78 canonical work pages

  1. [1]

    Computing trajectories We numerically integrated the trimer equations of mo- tion using the solve ivp function in SciPy [61]. We used the default RK4 integrator with a relative tolerance of 10−6 and an absolute tolerance of 10 −9 when comput- ing damped trajectories, and relative and absolute toler- ances of 10 −7 and 10 −10, respectively, when computing ...

    work page

  2. [2]

    First, differ- entiate Eq

    Derivation of the general averaged complex amplitude equation(7) Here we derive the general equation (7). First, differ- entiate Eq. (6) to obtain ¨xj =iω ˙Ajeiωt −ω 2Ajeiωt + c.c. (B1) and insert this result into Eq. (3) to obtain h (1−ω 2)Ajeiωt +iω ˙Ajeiωt + c.c. i +h j = 0 (B2) Next, add the two coordinate change equations (5)–(6) to find 1 2 xj −i ˙x...

    work page

  3. [3]

    5) Here, we derive the expressions Eq

    Analysis of conserved quantity (Fig. 5) Here, we derive the expressions Eq. (C17) and Eq. (C23) for the period in terms of the complete elliptic integral of the first kind, starting from their respective in- tegral expressions Eq. (C16) and Eq. (C22). We use the following convention for the complete elliptic integral of the first kind, after [40]: K(m) = ...

    work page

  4. [4]

    5) To extract the period of the amplitude oscillations from numerical trimer trajectories (Fig

    Determining the period of a numerical integration (Fig. 5) To extract the period of the amplitude oscillations from numerical trimer trajectories (Fig. 5, dots), we computed the upper envelope of the trajectory, used a Butterworth filter to remove high-frequency features, and computed the time between the adjacent local maxima of the fil- tered envelope. ...

    work page

  5. [5]

    closed” and those that do not will be called “open

    Parameter sweeps in Sec. III D For the integrations where the system parameters were varied (Figures 6 and 7), we varied the system parame- tersk,δ,ϵ, andγone at a time around a set of baseline parameters. When each parameter was varied, the oth- ers were kept at their baseline values, with the exception that whenkwas varied, the modulation frequencyγwas ...

    work page arXiv 2000

  6. [6]

    Background: the trimer fixed point in 6D phase space In order to understand what this means, we need to describe what the constant-amplitude state correspond- ing to the fixed point (32) represents in the full six- dimensional phase space of the unaveraged system (1). The averaging procedure described in section III cap- tures this state as single point i...

    work page

  7. [7]

    basin of attraction,

    Defining the radius of the basin of attraction We define thebasin of attraction of the fixed point am- plitudeas the set of all points (x 1, x2, x3, v1, v2, v3) that flow to oscillatory steady states having the amplituder ∗ given in Eq. (32). We will abuse notation and refer to this set as the “basin of attraction,” although strictly speak- ing it is not ...

    work page

  8. [8]

    In other words, the indicator function is 0 if the phase- space sphere of radiuspis entirely within the basin of attraction, and 1 if any portion of it lies outside the basin

    Finding the basin radius as a function ofk We begin by defining theindicator functionf I(p) to be 0 if any point at distancepfromy ∗ converges to a state with the fixed-point amplituder ∗, and 1 otherwise. In other words, the indicator function is 0 if the phase- space sphere of radiuspis entirely within the basin of attraction, and 1 if any portion of it...

    work page 2000

  9. [9]

    We first process these three traces to produce a single amplitude trace

    Extract a single amplitude trace The finite element method (FEM) traces are produced by recording the vertical displacement of a particular 23 point on each plate. We first process these three traces to produce a single amplitude trace. The procedure is as follows:

    work page

  10. [10]

    Numerically compute the local maxima of each con- tinuum trace to obtain three amplitude traces

    work page

  11. [11]

    Aggregate the three numerical amplitude traces to create a single trace, ordered by time

    work page

  12. [12]

    Generate a piecewise-linear interpolant of the ag- gregated traces to establish a constant time incre- ment

    work page

  13. [13]

    We use the savgol filter function of signal module within the Python library Scipy

    Filter the interpolant with a Savitsky-Golay filter. We use the savgol filter function of signal module within the Python library Scipy. This produces the FEM amplitude trace. We use a second-order filter with a window width corresponding to ten times the period 2π/ω 0 of the trimer fundamental mode. For Fig. 9,ω 0 = 52.607 MHz, and the dimensionless freq...

    work page

  14. [14]

    The frequency scale for the FEM trace is the frequencyω0 of the fundamental mode, which can be obtained directly as an output fromCOMSOL

    Rescale time and length units Next, we rescale the time and amplitude coordinates of the FEM amplitude traces to make them dimensionless. The frequency scale for the FEM trace is the frequencyω0 of the fundamental mode, which can be obtained directly as an output fromCOMSOL. We rescale the time points from the FEM trace to be dimensionless by multiplying ...

    work page

  15. [15]

    Perform fits in different dynamic regions We use information from the different dynamical re- gions to fit the dimensionless parameters of the discrete model. We first process the FEM amplitude trace in the ringdown region—the region where the trace amplitude undergoes damped oscillations about the fixed point—by subtracting out the fixed point value, so ...

    work page

  16. [16]

    9(c)), we numerically integrate Eq

    Integrate the amplitude equation To obtain the time-evolution of the chiral mode am- plitude (red curve in Fig. 9(c)), we numerically integrate Eq. (10). But first, we need to translate the FEM system initial conditions into initial conditions for the reduced averaged equations. Note that since the ansatz used to obtain the reduced averaged equations from...

    work page

  17. [17]

    LetR ∗0 be the numerically-computed steady-state am- plitude of the FEM trace for a set of FEM simulation parameters (α′ 0, δ′ 0, ϵ′ 0, ω′ 0)

    Steady-state amplitude forδ ′ = 0.2 Once the parameters are found, we can estimate the effect of changes in the parameters on changes in the FEM steady-state amplitude. LetR ∗0 be the numerically-computed steady-state am- plitude of the FEM trace for a set of FEM simulation parameters (α′ 0, δ′ 0, ϵ′ 0, ω′ 0). Hereδ ′ 0 corresponds to theδ ′ in Eq. (I1),ω...

    work page

  18. [18]

    W. W. Mumford, Some notes on the history of parametric transducers, Proceedings of the IRE48, 848 (2007)

    work page 2007

  19. [19]

    Smith, R

    A. Smith, R. Sandell, J. Burch, and A. Silver, Low noise microwave parametric amplifier, IEEE Transactions on Magnetics21, 1022 (1985)

    work page 1985

  20. [20]

    Aumentado, Superconducting parametric amplifiers: The state of the art in josephson parametric amplifiers, IEEE Microwave magazine21, 45 (2020)

    J. Aumentado, Superconducting parametric amplifiers: The state of the art in josephson parametric amplifiers, IEEE Microwave magazine21, 45 (2020)

    work page 2020

  21. [21]

    L.-A. Wu, M. Xiao, and H. Kimble, Squeezed states of light from an optical parametric oscillator, Journal of the Optical Society of America B4, 1465 (1987)

    work page 1987

  22. [22]

    Rugar and P

    D. Rugar and P. Gr¨ utter, Mechanical parametric am- plification and thermomechanical noise squeezing, Phys- ical Review Letters67, 699 (1991), publisher: American Physical Society

    work page 1991

  23. [23]

    Louisell,Coupled Mode and Parametric Electronics (John Wiley & Sons, New York, 1960)

    W. Louisell,Coupled Mode and Parametric Electronics (John Wiley & Sons, New York, 1960)

    work page 1960

  24. [24]

    Eichler and O

    A. Eichler and O. Zilberberg,Classical and Quantum Parametric Phenomena(Oxford University Press, 2023)

    work page 2023

  25. [25]

    Leuch, L

    A. Leuch, L. Papariello, O. Zilberberg, C. L. Degen, R. Chitra, and A. Eichler, Parametric symmetry break- ing in a nonlinear resonator, Physical review letters117, 214101 (2016)

    work page 2016

  26. [26]

    J. F. Rhoads and S. W. Shaw, The impact of nonlinear- ity on degenerate parametric amplifiers, Applied Physics Letters96(2010)

    work page 2010

  27. [27]

    Papariello, O

    L. Papariello, O. Zilberberg, A. Eichler, and R. Chi- tra, Ultrasensitive hysteretic force sensing with paramet- ric nonlinear oscillators, Physical Review E94, 022201 (2016)

    work page 2016

  28. [28]

    Z. Wang, A. Marandi, K. Wen, R. L. Byer, and Y. Ya- mamoto, Coherent ising machine based on degenerate optical parametric oscillators, Physical Review A88, 063853 (2013)

    work page 2013

  29. [29]

    Marandi, Z

    A. Marandi, Z. Wang, K. Takata, R. L. Byer, and Y. Ya- mamoto, Network of time-multiplexed optical parametric oscillators as a coherent Ising machine, Nature Photonics 8, 937 (2014)

    work page 2014

  30. [30]

    P. L. McMahon, A. Marandi, Y. Haribara, R. Hamerly, C. Langrock, S. Tamate, T. Inagaki, H. Takesue, S. Utsunomiya, K. Aihara, R. L. Byer, M. M. Fejer, H. Mabuchi, and Y. Yamamoto, A fully programmable 100-spin coherent Ising machine with all-to-all connec- tions, Science354, 614 (2016)

    work page 2016

  31. [31]

    T. L. Heugel, O. Zilberberg, C. Marty, R. Chitra, and A. Eichler, Ising machines with strong bilinear coupling, Physical Review Research4, 013149 (2022)

    work page 2022

  32. [32]

    Calvanese Strinati, L

    M. Calvanese Strinati, L. Bello, A. Pe’er, and E. G. Dalla Torre, Theory of coupled parametric oscillators be- yond coupled ising spins, Physical Review A100, 023835 (2019)

    work page 2019

  33. [33]

    Calvanese Strinati, L

    M. Calvanese Strinati, L. Bello, E. G. Dalla Torre, and A. Pe’er, Can nonlinear parametric oscillators solve ran- dom ising models?, Physical Review Letters126, 143901 (2021)

    work page 2021

  34. [34]

    T. L. Heugel, M. Oscity, A. Eichler, O. Zilberberg, and R. Chitra, Classical Many-Body Time Crystals, Physical Review Letters123, 124301 (2019), publisher: American Physical Society

    work page 2019

  35. [35]

    N. Y. Yao, C. Nayak, L. Balents, and M. P. Zaletel, Classical discrete time crystals, Nature Physics16, 438 (2020)

    work page 2020

  36. [36]

    Yi-Thomas and J

    S. Yi-Thomas and J. D. Sau, Theory for dissipative time crystals in coupled parametric oscillators, Physical Re- view Letters133, 266601 (2024)

    work page 2024

  37. [37]

    Lifshitz and M

    R. Lifshitz and M. Cross, Response of parametrically driven nonlinear coupled oscillators with application to micromechanical and nanomechanical resonator arrays, Physical review B67, 134302 (2003)

    work page 2003

  38. [38]

    Lifshitz and M

    R. Lifshitz and M. C. Cross, Nonlinear dynamics of nanomechanical and micromechanical resonators, Re- views of nonlinear dynamics and complexity1, 1 (2008)

    work page 2008

  39. [39]

    Margiani, O

    G. Margiani, O. Ameye, O. Zilberberg, and A. Eichler, Three strongly coupled kerr parametric oscillators form- ing a boltzmann machine, Physical Review Letters135, 097201 (2025)

    work page 2025

  40. [40]

    Ameye, A

    O. Ameye, A. Eichler, and O. Zilberberg, Parametric in- stability landscape of coupled Kerr parametric oscilla- tors, Physical Review Research7, 033204 (2025)

    work page 2025

  41. [41]

    Bello, M

    L. Bello, M. Calvanese Strinati, E. G. Dalla Torre, and A. Pe’er, Persistent Coherent Beating in Coupled Para- metric Oscillators, Physical Review Letters123, 083901 (2019)

    work page 2019

  42. [42]

    C. A. Downing, D. Zueco, and L. Mart´ ın-Moreno, Chiral Current Circulation and PT Symmetry in a Trimer of Os- cillators, ACS Photonics 10.1021/acsphotonics.0c01208 (2020)

    work page doi:10.1021/acsphotonics.0c01208 2020

  43. [43]

    Melkani and J

    A. Melkani and J. Paulose, Space-time symmetry and nonreciprocal parametric resonance in mechanical sys- tems, Physical Review E110, 015003 (2024)

    work page 2024

  44. [44]

    Kruss and J

    N. Kruss and J. Paulose, Nondispersive One-Way Signal Amplification in Sonic Metamaterials, Physical Review Applied17, 024020 (2022)

    work page 2022

  45. [45]

    Melkani and J

    A. Melkani and J. Paulose, Space-time Floquet opera- tor: Non-reciprocity and fractional topology of space- time crystals (2025), arXiv:2510.16562 [cond-mat]

    work page arXiv 2025

  46. [46]

    Zangeneh-Nejad and R

    F. Zangeneh-Nejad and R. Fleury, Active times for acous- tic metamaterials, Reviews in Physics4, 100031 (2019)

    work page 2019

  47. [47]

    Ventsel and T

    E. Ventsel and T. Krauthammer,Thin Plates and Shells: Theory, Analysis, and Applications(Taylor & Francis Group, 2001)

    work page 2001

  48. [48]

    Karki and J

    P. Karki and J. Paulose, Stopping and Reversing Sound via Dynamic Dispersion Tuning in a Phononic Metama- terial, Physical Review Applied15, 034083 (2021), pub- lisher: American Physical Society

    work page 2021

  49. [49]

    Yakubovich, V

    V. Yakubovich, V. Starzhinskii, and G. Schmidt,On Linear Differential Equations with Periodic Coefficients, Vol. 1 (Wiley, London, 1976) pp. 222–222

    work page 1976

  50. [50]

    L. D. Landau and E. M. Lifshitz,Mechanics, 2nd ed., Vol. 1 (Butterworth Heinemann, 1960)

    work page 1960

  51. [51]

    S. H. Strogatz,Nonlinear dynamics and Chaos: with ap- plications to physics, biology, chemistry, and engineering, Studies in nonlinearity (Addison-Wesley Pub, Reading, Mass, 1994)

    work page 1994

  52. [52]

    J. A. Sanders, F. Verhulst, J. A. Murdock, and J. Mur- dock,Averaging methods in nonlinear dynamical systems, 2nd ed., Applied mathematical sciences No. 59 (Springer, New York, NY, 2007)

    work page 2007

  53. [53]

    Guckenheimer, Dynamics of the Van der Pol equa- tion, IEEE Transactions on Circuits and Systems27, 983 (1980)

    J. Guckenheimer, Dynamics of the Van der Pol equa- tion, IEEE Transactions on Circuits and Systems27, 983 (1980)

    work page 1980

  54. [54]

    Kovacic, R

    I. Kovacic, R. Rand, and S. Mohamed Sah, Mathieu’s 26 equation and its generalizations: overview of stability charts and their features, Applied Mechanics Reviews70, 020802 (2018)

    work page 2018

  55. [55]

    We add an overall (−) sign to the phases because the am- plified mode travels in the direction opposite the modu- lation wave [26]

    work page

  56. [56]

    When ∆γ̸= 0, this quantity is no longer conserved and there are no longer closed orbits, but in the extended phase space each time slice of the flow field is identical— up to a shift in the phase offsetc—to the resonance flow field by Eq. (12)

    work page

  57. [57]

    Abramowitz and I

    M. Abramowitz and I. A. Stegun,Handbook of mathe- matical functions with formulas, graphs, and mathemati- cal tables, Vol. 55 (US Government printing office, 1948)

    work page 1948

  58. [58]

    Mahboob and H

    I. Mahboob and H. Yamaguchi, Bit storage and bit flip operations in an electromechanical oscillator, Nature Nanotechnology3, 275 (2008)

    work page 2008

  59. [59]

    Mahboob, M

    I. Mahboob, M. Mounaix, K. Nishiguchi, A. Fujiwara, and H. Yamaguchi, A multimode electromechanical para- metric resonator array, Scientific Reports4, 4448 (2014)

    work page 2014

  60. [60]

    J. P. Mathew, R. N. Patel, A. Borah, R. Vijay, and M. M. Deshmukh, Dynamical strong coupling and parametric amplification of mechanical modes of graphene drums, Nature Nanotechnology11, 747 (2016)

    work page 2016

  61. [61]

    M. J. Seitner, M. Abdi, A. Ridolfo, M. J. Hartmann, and E. M. Weig, Parametric Oscillation, Frequency Mix- ing, and Injection Locking of Strongly Coupled Nanome- chanical Resonator Modes, Physical Review Letters118, 254301 (2017)

    work page 2017

  62. [62]

    J. M. Miller, D. D. Shin, H.-K. Kwon, S. W. Shaw, and T. W. Kenny, Phase Control of Self-Excited Parametric Resonators, Physical Review Applied12, 044053 (2019)

    work page 2019

  63. [63]

    Mestre, S

    L. Mestre, S. Singh, G. Margiani, L. Catalini, A. Eichler, and V. Dumont, Network of parametrically driven silicon nitride mechanical membranes, Physical Review Applied 24, 044072 (2025)

    work page 2025

  64. [64]

    Prasad, N

    P. Prasad, N. Arora, and A. K. Naik, Parametric am- plification in MoS 2 drum resonator, Nanoscale9, 18299 (2017)

    work page 2017

  65. [65]

    Naserbakht, A

    S. Naserbakht, A. Naesby, and A. Dantan, Stress- Controlled Frequency Tuning and Parametric Amplifica- tion of the Vibrations of Coupled Nanomembranes, Ap- plied Sciences9, 4845 (2019)

    work page 2019

  66. [66]

    Z.-J. Su, Y. Ying, X.-X. Song, Z.-Z. Zhang, Q.-H. Zhang, G. Cao, H.-O. Li, G.-C. Guo, and G.-P. Guo, Tunable parametric amplification of a graphene nanomechanical resonator in the nonlinear regime, Nanotechnology32, 155203 (2021)

    work page 2021

  67. [67]

    A. M. v. d. Zande, R. A. Barton, J. S. Alden, C. S. Ruiz- Vargas, W. S. Whitney, P. H. Q. Pham, J. Park, J. M. Parpia, H. G. Craighead, and P. L. McEuen, Large-scale arrays of single-layer graphene resonators,Nano Letters, Nano Letters10, 4869 (2010)

    work page 2010

  68. [68]

    J. Cha, K. W. Kim, and C. Daraio, Experimental realiza- tion of on-chip topological nanoelectromechanical meta- materials, Nature564, 229 (2018)

    work page 2018

  69. [69]

    Carter, U

    B. Carter, U. F. Hernandez, D. J. Miller, A. Blaikie, V. R. Horowitz, and B. J. Alem´ an, Coupled Nanomechanical Graphene Resonators: A Promising Platform for Scalable NEMS Networks, Micromachines14, 2103 (2023)

    work page 2023

  70. [70]

    Cha and C

    J. Cha and C. Daraio, Electrical tuning of elastic wave propagation in nanomechanical lattices at mhz frequen- cies, Nature Nanotechnology13, 1016 (2018)

    work page 2018

  71. [71]

    T. Mei, J. Lee, Y. Xu, and P. X.-L. Feng, Frequency tuning of graphene nanoelectromechanical resonators via electrostatic gating, Micromachines9, 312 (2018)

    work page 2018

  72. [72]

    Blaikie, D

    A. Blaikie, D. Miller, and B. J. Alem´ an, A fast and sensitive room-temperature graphene nanomechanical bolometer, Nature Communications10, 1 (2019)

    work page 2019

  73. [73]

    Apffel and R

    B. Apffel and R. Fleury, Experimental observation of topological transition in linear and nonlinear paramet- ric oscillators, Physical Review E109, 054204 (2024)

    work page 2024

  74. [74]

    Inagaki, Y

    T. Inagaki, Y. Haribara, K. Igarashi, T. Sonobe, S. Ta- mate, T. Honjo, A. Marandi, P. L. McMahon, T. Umeki, K. Enbutsu, O. Tadanaga, H. Takenouchi, K. Aihara, K.- i. Kawarabayashi, K. Inoue, S. Utsunomiya, and H. Take- sue, A coherent Ising machine for 2000-node optimization problems, Science354, 603 (2016)

    work page 2000

  75. [75]

    K. H. Matlack, M. Serra-Garcia, A. Palermo, S. D. Hu- ber, and C. Daraio, Designing perturbative metamate- rials from discrete models, Nature Materials17, 323 (2018)

    work page 2018

  76. [76]

    Fleury, D

    R. Fleury, D. L. Sounas, C. F. Sieck, M. R. Haberman, and A. Al` u, Sound Isolation and Giant Linear Nonre- ciprocity in a Compact Acoustic Circulator, Science343, 516 (2014)

    work page 2014

  77. [77]

    Nassar, B

    H. Nassar, B. Yousefzadeh, R. Fleury, M. Ruzzene, A. Al` u, C. Daraio, A. N. Norris, G. Huang, and M. R. Haberman, Nonreciprocity in acoustic and elastic mate- rials, Nature Reviews Materials , 1 (2020)

    work page 2020

  78. [78]

    Virtanen, R

    P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Pe- terson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙I. Po- lat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henr...

    work page 2020