Generic Rotating-Frame-Based Approach to Chaos Generation in Nonlinear Micro- and Nanoelectromechanical System Resonators

Samer Houri, Motoki Asano, Hiroshi Yamaguchi, Natsue Yoshimura, Yasuharu Koike, and Ludovico Minati

Phys. Rev. Lett. 125, 174301 – Published 23 October 2020

Abstract

This Letter provides a low-power method for chaos generation that is generally applicable to nonlinear micro- and nanoelectromechanical systems (MNEMS) resonators. The approach taken is independent of the material, scale, design, and actuation of the device in question; it simply assumes a good quality factor and a Duffing type nonlinearity, features that are commonplace to MNEMS resonators. The approach models the rotating-frame dynamics to analytically constrain the parameter space required for chaos generation. By leveraging these common properties of MNEMS devices, a period-doubling route to chaos is generated using smaller forcing than typically reported in the literature.

Received 26 February 2020
Revised 26 July 2020
Accepted 9 September 2020

DOI:https://doi.org/10.1103/PhysRevLett.125.174301

Physics Subject Headings (PhySH)

Chaos

Chaotic systems Micromechanical devices Nanomechanical devices

Time series analysis

Nonlinear DynamicsCondensed Matter, Materials & Applied Physics

Authors & Affiliations

Samer Houri ^*, Motoki Asano, and Hiroshi Yamaguchi

NTT Basic Research Laboratories, NTT Corporation, 3-1 Morinosato-Wakamiya, Atsugi-shi, Kanagawa 243-0198, Japan

Natsue Yoshimura, Yasuharu Koike, and Ludovico Minati ^†

Institute of Innovative Research, Tokyo Institute of Technology, Yokohama 226-8503, Japan

^*Samer.Houri.dg@hco.ntt.co.jp
^†minati.l.aa@m.titech.ac.jp

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand

Issue

Vol. 125, Iss. 17 — 23 October 2020

Reuse & Permissions

Access Options

Author publication services for translation and copyediting assistance advertisement

Images

Figure 1
(a) Bistability map plotted as a function of dimensionless force and detuning, showing the region of bistability for a lossless driven Duffing resonator (gray area), and for a low-loss ( $Q = 1000$ ) Duffing resonator (area between the dashed blue lines). (b) Amplitude (arbitrary units) versus detuning response of a lossless Duffing taken for $F_{1} \sqrt{α} = 10^{- 4}$ . The corresponding phase-space plot (also in arbitrary units) for a detuning of $δ = 2.5 \times 10^{- 3}$ is shown in (c). The stable fixed points and the saddle point are shown as black and green dots, respectively, and the black traces correspond to the homoclinic orbits. Small amplitude libration orbits around the high-amplitude branch (blue) and low-amplitude branch (red) are shown. $α$ is the Duffing parameter.
Reuse & Permissions
Figure 2
(a) Experimental lock-in amplifier data for the frequency response obtained by a single tone sweep showing the linear (black trace $100 m V_{PP}$ ) and the Duffing regimes ( $3 V_{PP}$ ); the latter shows bistability upon performing a forward (blue trace) and a backward sweep (red trace). $H$ denotes the relative amplitude response, expressed in mV per V drive. (b) Scatter plot of periodically sampled libration oscillation under the effect of a two-tone excitation, with one fixed tone (indicated by $F_{1}$ having $ω_{1} / 2 π = 1566.5 kHz$ and $F_{1} = 3 V_{PP}$ ) and one swept tone ( $F_{2} = 2.1 V_{PP}$ , $1558 kHz < ω_{2} / 2 π < 1576 kHz$ ), shown for the lower (red) and higher (blue) amplitude branches. Period 1, period 2, and chaotic oscillations are detected and shown in (c)–(e), respectively, for both the high (blue) and low (red) amplitude branches. The black dots correspond to the experimentally-obtained fixed points, and the crosses correspond to the calculated saddle point.
Reuse & Permissions
Figure 3
Numerical simulations obtained for $δ = 4.2 \times 10^{- 3}$ (i.e., $ω_{1} / 2 π = 1566.5 kHz$ ), and $F_{1} \sqrt{α} = 1.5 \times 10^{- 4}$ (i.e., $F_{1} = 3 V_{PP}$ ). (a) Basins of attraction under the effect of $F_{2} = 0.6 F_{1}$ , $Ω = 3.275 \times 10^{- 3}$ for the high-branch, and $F_{2} = 0.5 F_{1}$ , $Ω = - 2.7 \times 10^{- 3}$ for the low-branch showing that only one of the two basins is disturbed depending on the value of $Ω$ . The white dots and crosses mark the location of the stable and saddle points in the unperturbed system, respectively. (b) Two-dimensional maps (both panels) showing the location of autocorrelation peak as a function of detuning and forcing ( $Ω$ , $F_{2}$ ). Values greater than 1 (red and bright areas) indicate period-doubling bifurcations and chaos. The experimentally obtained bifurcation areas are equally shown (delineated by the dashed lines). The area bounded by the analytical model, i.e., $F_{2 C}$ , is shown as the solid yellow lines. The solid vertical lines indicate where in the ( $F_{2}$ , $Ω$ ) parameter space the basins in (a) are located.
Reuse & Permissions
Figure 4
Comparison between numerical (dots) and analytical (solid lines) values of $F_{2 \min}$ (b) and $Ω_{\min}$ (c) required for $P 2$ and chaos for both the low amplitude branch (red) and high amplitude branch (blue) of the Duffing resonator shown in (a). For both numerical and analytical data, $P 2$ only appears in the range where bistability exists (indicated by the dashed vertical gray lines, and denoted $A$ and $B$ , respectively). As the high amplitude branch approaches the saddle-node bifurcation its libration motion becomes more unstable, hence the data points do not reach the saddle-node bifurcation.
Reuse & Permissions

Physical Review Letters