Data based identification and prediction of nonlinear and complex dynamical systems

We thank Dr. R. Yang (formerly at ASU), Dr. R.-Q. Su (formerly at ASU), and Mr. Zhesi Shen for their contributions to a number of original papers on which this Review is partly based. This work was supported by ARO under Grant No. W911NF-14-1-0504. W.-X. Wang was also supported by NSFC under Grants No. 61573064 and No. 61074116, as well as by the Fundamental Research Funds for the Central Universities, Beijing Nova Programme.Peer reviewedPostprin

Isomerization dynamics of a buckled nanobeam

We analyze the dynamics of a model of a nanobeam under compression. The model is a two mode truncation of the Euler-Bernoulli beam equation subject to compressive stress. We consider parameter regimes where the first mode is unstable and the second mode can be either stable or unstable, and the remaining modes (neglected) are always stable. Material parameters used correspond to silicon. The two mode model Hamiltonian is the sum of a (diagonal) kinetic energy term and a potential energy term. The form of the potential energy function suggests an analogy with isomerisation reactions in chemistry. We therefore study the dynamics of the buckled beam using the conceptual framework established for the theory of isomerisation reactions. When the second mode is stable the potential energy surface has an index one saddle and when the second mode is unstable the potential energy surface has an index two saddle and two index one saddles. Symmetry of the system allows us to construct a phase space dividing surface between the two "isomers" (buckled states). The energy range is sufficiently wide that we can treat the effects of the index one and index two saddles in a unified fashion. We have computed reactive fluxes, mean gap times and reactant phase space volumes for three stress values at several different energies. In all cases the phase space volume swept out by isomerizing trajectories is considerably less than the reactant density of states, proving that the dynamics is highly nonergodic. The associated gap time distributions consist of one or more `pulses' of trajectories. Computation of the reactive flux correlation function shows no sign of a plateau region; rather, the flux exhibits oscillatory decay, indicating that, for the 2-mode model in the physical regime considered, a rate constant for isomerization does not exist.Comment: 42 pages, 6 figure