Multistability of the Vibrating System of a Micro Resonator

Multiple attractors and their fractal basins of attraction can lead to the loss of global stability and integrity of Micro Electro Mechanical Systems (MEMS). In this paper, multistability of a class of electrostatic bilateral capacitive micro-resonator is researched in detail. First, the dynamical model is established and made dimensionless. Second, via the perturbating method and the numerical description of basins of attraction, the multiple periodic motions under primary resonance are discussed. It is found that the variation of AC voltage can induce safe jump of the micro resonator. In addition, with the increase of the amplitude of AC voltage, hidden attractors and chaos appear. The results may have some potential value in the design of MEMS devices.


Introduction
Multistability, i.e., the coexistence of multiple attractors, is a common dynamical phenomenon in MEMS/NEMS [1,2]. Based on it, there are many applications such as MEMS-based memory [3] and switches [4]. In addition, considering the loss of global stability that multistability may trigger, there are some devices that should avoid the appearance of multiple attractors in their vibrating systems, such as filters [5], microvalves [6], and micro-relays [7]. As one of the fastest developing MEMS products [8], electrostatic micro-resonators should assure that the resonators undergo periodic vibration whose amplitudes vary continuously with the driven voltages. However, in practical applications of electrostatic micro-resonators [9], there are many complex dynamic behaviors such as multistability [10,11], quasi-periodic motion [12] and periodic-n motion [13], chaos [14,15], and pull-in instability.
It is of great significance to study the multistability and necessary conditions for inducing it either for avoiding this phenomenon or making use of it. Thus, multistability of vibrating systems of micro resonators has been studied experimentally and numerically during these decades [16]. Via experiments, Mohammadreza investigated the dynamic response of an electrostatic micro-actuator in the vicinity of the primary resonance and the parametric one [17]. Siewe et al. [18] studied the vibration of a double-side MEMS resonator numerically and found the variation of the driven voltage could induce the coexistence of chaos and quasi-periodic motions. Shang et al. [19] found the coexisting chaos and dynamical pull-in in the vibrating system of a single-side electrostatic micro sensor. Haghighi et al. [20] found the coexisting periodic-n motion and the chaotic motion of micromechanical resonators with electrostatic forces on both sides, and then discussed the global bifurcation of its vibrating system by approximately expressing its homoclinic orbits as the ones of a typical duffing equation. When amplifying signals of a nanomechanical duffing resonator, Almog et al. [21] found that multistability was an interesting dynamical phenomenon of nonlinear systems and could be explored for many applications. Gusso et al. [22] studied chaos of a typical micro/nanoelectromechanical beam resonator with two-sided electrodes experimentally and observed multiple attractors in a significant region of the relevant parameter space, involving periodic and chaotic attractors.
By applying cell-mapping method to depict the basins of attraction for all the attractors, they also found that the basin boundaries were fractal under certain conditions of the excitations, indicating that the attractors are strongly intermingled. Liu et al. [23] applied the method of multiple scales (MMS) to analyze the multiple periodic motions induced by the local bifurcation, and used the Melnikov method to predict necessary conditions for chaos and its control. The corresponding numerical results were also presented by the basins of attraction and spectrum diagrams. Angelo et al. [24] investigated the effect of the linear and nonlinear stiffness terms and damping coefficients on dynamical behaviors of a microelectromechanical resonator and controlled the chaotic motion by forcing it into an orbit obtained analytically via the harmonic balance method. However, most study concentrated on describing or observing the phenomenon itself rather than studying its mechanism, which is still not that clear yet.
To this end, we consider a typical electrostatic driven bilateral capacitive microresonator and study the possible multistability and its mechanism in its vibrating system. The paper is organized as follows. In Section 2, the dynamical model is constructed and made dimensionless. In Sections 3 and 4, two different cases for coexisting multiple periodic attractors, fractal basins of attraction, and other complex attractors of the systems are discussed both theoretically and numerically. In Section 5, the conclusions are presented.

Dynamical Model
We choose to study a class of bilateral micro resonator whose simplified diagram is shown in Figure 1. The driven forces on the resonator are electrostatic ones between the moving electrode and the fixed electrode [25]. The driven voltage in Figure 1 is the combination of alternate current (AC) and direct current (DC) actuation. In the figure, x is the vertical displacement of the moving electrode at moment t, d the initial gap width between the moving electrode and each fixed one, V b the DC bias voltage, V AC sinΩt the AC voltage where V AC is the amplitude and Ω the frequency. Suppose that the amplitude of the AC voltage V AC is much lower than the bias DC voltage V b , i.e., V AC V b . According to the Second Law of Newton, the vibrating system of the moving electrode can be expressed as a nonlinear system as follows: where m represents the effective lumped mass of the moving electrode, k 1 its linear mechanical stiffness, k 2 its cubic nonlinear stiffness, c the damping coefficient, C 0 the initial capacitance of the parallel-plate structure. with two-sided electrodes experimentally and observed multiple attractors in a significant region of the relevant parameter space, involving periodic and chaotic attractors. By applying cell-mapping method to depict the basins of attraction for all the attractors, they also found that the basin boundaries were fractal under certain conditions of the excitations, indicating that the attractors are strongly intermingled. Liu et al. [23] applied the method of multiple scales (MMS) to analyze the multiple periodic motions induced by the local bifurcation, and used the Melnikov method to predict necessary conditions for chaos and its control. The corresponding numerical results were also presented by the basins of attraction and spectrum diagrams. Angelo et al. [24] investigated the effect of the linear and nonlinear stiffness terms and damping coefficients on dynamical behaviors of a microelectromechanical resonator and controlled the chaotic motion by forcing it into an orbit obtained analytically via the harmonic balance method. However, most study concentrated on describing or observing the phenomenon itself rather than studying its mechanism, which is still not that clear yet.
To this end, we consider a typical electrostatic driven bilateral capacitive micro-resonator and study the possible multistability and its mechanism in its vibrating system. The paper is organized as follows. In Section 2, the dynamical model is constructed and made dimensionless. In Sections 3 and 4, two different cases for coexisting multiple periodic attractors, fractal basins of attraction, and other complex attractors of the systems are discussed both theoretically and numerically. In Section 5, the conclusions are presented.

Dynamical Model
We choose to study a class of bilateral micro resonator whose simplified diagram is shown in Figure 1. The driven forces on the resonator are electrostatic ones between the moving electrode and the fixed electrode [25]. The driven voltage in Figure 1 is the combination of alternate current (AC) and direct current (DC) actuation. In the figure, x is the vertical displacement of the moving electrode at moment t, d the initial gap width between the moving electrode and each fixed one, Vb the DC bias voltage, VACsinΩt the AC voltage where VAC is the amplitude and Ω the frequency. Suppose that the amplitude of the AC voltage VAC is much lower than the bias DC voltage Vb, i.e., AC b V V  . According to the Second Law of Newton, the vibrating system of the moving electrode can be expressed as a nonlinear system as follows: where m represents the effective lumped mass of the moving electrode, k1 its linear mechanical stiffness, k2 its cubic nonlinear stiffness, c the damping coefficient, C0 the initial capacitance of the parallel-plate structure. Introducing the following dimensionless variables Introducing the following dimensionless variables which is a dimensionless system. Since in the original system (1), the viscous damping coefficient of air c is very tiny, and V AC V b , the parameters µ and γ in (3) will be both small and can be considered as perturbed parameters. Thus, considering µ = 0 and γ = 0 in Equation (3), one has the unperturbed system that can be expressed as below: Letting the right side of Equation (4) be zero, one can determine equilibria of the dimensionless system (3). Equation (5) is a Hamilton system with the Hamiltonian and the function of potential energy (P.E.) Concerning Equation (4), the number of the equilibria, and the shapes and positions of the possible potential wells of the unperturbed system (4) depend on the parameters α and β. The same as in [20], the values of the parameters in the system (1) are given by: m = 5 × 12 −12 kg, c = 5 × 12 −8 kg/s, k 1 = 5 µN/µm, k 2 = 15 µN/µm 3 , d = 2µm, C 0 = 1.875 × 10 −18 mF. (7) Accordingly, in system (4), α = 12. Different equilibria and potential energy diagrams of the unperturbed system under different values of the parameter β can be seen in Figure 2. It shows that there are three P.E. poles when β = 0.211, five P.E. poles when β increases to 0.338, and only one P.E. pole when β increases to 0.6. Under different values of β, the potential wells and unperturbed orbits are shown in Figure 3. When β = 0.211, there are three equilibria (two non-trivial equilibria are saddles and the origin is a center) as well as one well surrounded by heteroclinic orbits (see Figure 3a). As β increases to 0.338, there will be five equilibria among which two non-trivial equilibria S 1 (−0.196339,0) and S 2 (0.196339,0) are centers of the two wells surrounded by homoclinic orbits; the other three equilibria are unstable. When β = 0.6, no wells or nontrivial equilibria of the unperturbed system (4) exist. The P.E. poles in Figure 2 correspond to the fixed points shown in Figure 3. Therefore, according to Equations (2) and (7), when the structural parameters are fixed, the number of centers will depend on the value of DC bias voltage V b : when the DC bias voltage is very low, there will be a center of the system (4) as well as a stable point attractor of the system (3) without AC voltage. Under a higher DC bias voltage, there may be two centers of the system (4). As is well known, periodic vibration can often be attributed to the perturbation of the centers. Since the number and the location of the centers in Figure 3a,b are totally different, the mechanism for the possible multiple periodic attractors of the vibrating system of the micro-resonator can be different as well. Therefore, in Sections 3 and 4, we discuss the different mechanism of multi-stability for these two different cases, i.e., the only center (the origin) and the two non-trivial centers, respectively.

Multiple Periodic Attractors in the Neighborhood of the Origin
Considering the case where the DC bias voltage is low, and the periodic vibration of the microstructure is induced by the perturbation of the only center (see Figure 3a, where Vb = 3 V), one may use the Method of Multiple Scales (MMS) to analyze the periodic solutions in the neighborhood of the origin. Expanding the fractional terms of the dimensionless system (3) as Taylor series in the neighborhood of u = 0, and neglecting the higherorder-than-three terms of u, one has: As mentioned in Section 2, the values of the parameters μ and γ in the above system are small; one can introduce a small parameter ε satisfying 0 1 ε <  , and can rescale the two parameters in the system (8) as: 2 2 , .
Then Equation (8) becomes To apply MMS, one may rescale some terms in the system (10) that and

Multiple Periodic Attractors in the Neighborhood of the Origin
Considering the case where the DC bias voltage is low, and the periodic vibration of the microstructure is induced by the perturbation of the only center (see Figure 3a, where Vb = 3 V), one may use the Method of Multiple Scales (MMS) to analyze the periodic solutions in the neighborhood of the origin. Expanding the fractional terms of the dimensionless system (3) as Taylor series in the neighborhood of u = 0, and neglecting the higherorder-than-three terms of u, one has: As mentioned in Section 2, the values of the parameters μ and γ in the above system are small; one can introduce a small parameter ε satisfying 0 1 ε <  , and can rescale the two parameters in the system (8) as: 2 2 , .
Then Equation (8) becomes To apply MMS, one may rescale some terms in the system (10) that and

Multiple Periodic Attractors in the Neighborhood of the Origin
Considering the case where the DC bias voltage is low, and the periodic vibration of the microstructure is induced by the perturbation of the only center (see Figure 3a, where V b = 3 V), one may use the Method of Multiple Scales (MMS) to analyze the periodic solutions in the neighborhood of the origin. Expanding the fractional terms of the dimensionless system (3) as Taylor series in the neighborhood of u = 0, and neglecting the higher-order-than-three terms of u, one has: ..
As mentioned in Section 2, the values of the parameters µ and γ in the above system are small; one can introduce a small parameter ε satisfying 0 < ε 1, and can re-scale the two parameters in the system (8) as: Then Equation (8) becomes where To apply MMS, one may rescale some terms in the system (10) that and Comparing the coefficients of ε 1 , ε 2 , and ε 3 in the system (10), respectively, one obtains that and To solve Equation (14), one can assume that where Substituting Equations (17) and (18) into Equation (15), and eliminating the secular terms of Equation (15), one will have: Solving Equation (15), one may assume: Substituting Equation (20) into Equation (16), and eliminating secular terms of Equation (16), one will obtain: Substituting Equations (19) and (21) into Equation (22), and expressing it by the original dimensionless parameters of Equation (3), one has: where According to Equation (18), it is obvious that the amplitude of the periodic solution a is the function of the time scale T 1 where T 1 is a one-order term of ε; thus, one can assume the amplitude of the solution u of Equation (10) a as: Letting . a = 0, and . θ = 0, one can obtain: Eliminating the triangulation function of Equation (26), one can get: According to Equations (17), (18) and (25), the periodic solution can be expressed as: To determine the stability of the periodic solutions, one can get the corresponding characteristic equation of the periodic solution based on Equation (23). It shows that the periodic solution will lose its stability when its amplitude a satisfies: Based on Equations (26)-(29), the variation of the amplitude of the periodic solutions of the system (3) and their stability with AC voltage is shown in Figure 4 where the frequency and the amplitude of AC voltage are considered as the control parameters in Figure 4a and Figure 4b, respectively. In Figure 4a, where V AC = 0.01 V, when ω is lower than 0.45, there is only one periodic attractor in the system (3) whose amplitude changes continuously with the increase in ω. Comparatively, when ω ranges from 0.46 to 0.69, the global dynamical behaviors of the system (3) will change to bistable periodic attractors, which can be attributed to Hopf bifurcation. As ω continues to increase from 0.7, the periodic attractor with the higher amplitude will disappear, only the periodic attractor with the lower amplitude will exist, and its amplitude will decrease continuously with the increase of ω. Similarly, the change in global dynamical behaviors in Figure 4b also shows that in a certain range of V AC , there will be two periodic attractors coexisting, which can be due to Hopf bifurcation of the system (3). The accuracy of the theoretical prediction in Figure 4 are verified by the numerical results.
Eliminating the triangulation function of Equation (26), one can get: According to Equations (17), (18) and (25), the periodic solution can be expressed as: To determine the stability of the periodic solutions, one can get the corresponding characteristic equation of the periodic solution based on Equation (23). It shows that the periodic solution will lose its stability when its amplitude a  satisfies: Based on Equations (26)-(29), the variation of the amplitude of the periodic solutions of the system (3) and their stability with AC voltage is shown in Figure 4 where the frequency and the amplitude of AC voltage are considered as the control parameters in Figure 4a and Figure 4b, respectively. In Figure 4a, where VAC = 0.01 V, when ω is lower than 0.45, there is only one periodic attractor in the system (3) whose amplitude changes continuously with the increase in ω. Comparatively, when ω ranges from 0.46 to 0.69, the global dynamical behaviors of the system (3) will change to bistable periodic attractors, which can be attributed to Hopf bifurcation. As ω continues to increase from 0.7, the periodic attractor with the higher amplitude will disappear, only the periodic attractor with the lower amplitude will exist, and its amplitude will decrease continuously with the increase of ω. Similarly, the change in global dynamical behaviors in Figure 4b also shows that in a certain range of VAC, there will be two periodic attractors coexisting, which can be due to Hopf bifurcation of the system (3). The accuracy of the theoretical prediction in Figure 4 are verified by the numerical results.      (3), different initial conditions may lead to different periodic attractors. Accordingly, it is necessary to classify the basins of attraction for the two different periodic attractors. Here, the 4th order Runge-Kutta approach and the cell-mapping method are applied to depict the basins of attraction of the system (3). The time step is taken as 1/10 2 of the period of excitation. To investigate the long-term dynamical behaviors, it is supposed that an initial condition will be safe if the vibration in this initial condition keeps satisfying |u(T)| < 1 within 10 5 excited circles; otherwise, the micro resonator will undergo pull in [20]. The union of all initial conditions leading to the same periodic motion will be the basin of attraction for that attractor which will surely be marked in the same color in the initial plane. The basins of attraction of system (3) are drawn in sufficiently large ranges for the initial position and velocity of the proof mass defined as |u(0)| < 1 and | . u(0)| < 1.5 by generating an 200 × 100 array of initial points. The change of attractors and the area and nature of their basins of attraction with frequency ω is shown in Figure 5 where the amplitude of AC voltage V AC is fixed as 0.01 V.
proof mass defined as |u(0)| < 1 and |u (0)| < 1.5 by generating an 200 × 100 array of initial points. The change of attractors and the area and nature of their basins of attraction with frequency ω is shown in Figure 5 where the amplitude of AC voltage VAC is fixed as 0.01 V.
According to Figure 5, with the increase in parameter ω, the number of attractors and the boundary of basins of attraction will both change. When ω = 0.45 in system (3), there will only be one periodic attractor whose basin of attraction is comparatively bigger with a smooth boundary (see Figure 5a1,a2). However, with a small increase of ω, i.e., ω = 0.45, the global dynamics are totally different, as shown in Figure 5b1,b2 where two periodic attractors coexist, whose basins of attraction mix each other and are both fractal. It means that the dynamical behavior of system (3) is highly sensitive to initial conditions. In other words, system (3) may undergo a safe jump. A similar phenomenon can be seen in Figure  5c1,c2 under a higher ω. As ω increases, the basin of attraction of the periodic attractor with the higher amplitude becomes small (see the red regions in Figure 5b2,c2 and d2). Specifically, in Figure 5d2, the regions of attraction are almost blue, and there is very little area of basin of attraction for that periodic attractor. As ω = 0.70 (see Figure 5e1,e2), the periodic attractor with the higher amplitude disappears, and there is only the other attractor whose basin of attraction is almost the same as that in Figure 4b, showing that when the frequency ω increases enough, the periodic attractor with the lower amplitude replaces the initial one.

Multistability in the Neighborhood of Non-Trivial Equilibria
In this section, the case that the DC voltage is higher, and the periodic vibration of the microstructure is induced by the perturbation of the two non-trivial centers (see Figure  3b) is considered; thus, we set β = 0.338, i.e., According to Figure 5, with the increase in parameter ω, the number of attractors and the boundary of basins of attraction will both change. When ω = 0.45 in system (3), there will only be one periodic attractor whose basin of attraction is comparatively bigger with a smooth boundary (see Figure 5a1,a2). However, with a small increase of ω, i.e., ω = 0.45, the global dynamics are totally different, as shown in Figure 5b1,b2 where two periodic attractors coexist, whose basins of attraction mix each other and are both fractal. It means that the dynamical behavior of system (3) is highly sensitive to initial conditions. In other words, system (3) may undergo a safe jump. A similar phenomenon can be seen in Figure 5c1,c2 under a higher ω. As ω increases, the basin of attraction of the periodic attractor with the higher amplitude becomes small (see the red regions in Figure 5b2,c2,d2). Specifically, in Figure 5d2, the regions of attraction are almost blue, and there is very little area of basin of attraction for that periodic attractor. As ω = 0.70 (see Figure 5e1,e2), the periodic attractor with the higher amplitude disappears, and there is only the other attractor whose basin of attraction is almost the same as that in Figure 4b, showing that when the frequency ω increases enough, the periodic attractor with the lower amplitude replaces the initial one.

Multistability in the Neighborhood of Non-Trivial Equilibria
In this section, the case that the DC voltage is higher, and the periodic vibration of the microstructure is induced by the perturbation of the two non-trivial centers (see Figure 3b) is considered; thus, we set β = 0.338, i.e., V b = 3.8V. In addition, we consider the effect of AC voltage on the global dynamics of the system (3). To begin with, setting where u c is the abscissa of the right center (see S 1 in Figure 3), rescaling he two parameters µ and γ in the system (3) by Equation (9), expanding the fractional terms of the dimensionless system (3) as a Taylor series in the neighborhood of the non-trivial equilibria and ignoring the higher-order-than-cubic terms ofû, the system (3) becomes .. wherê To apply the Method of Multiple Scale in Equation (32), one can assume in this equation that: Comparing the coefficients of ε 1 , ε 2 and ε 3 , one has: and One can set the solution of Equation (34) as: Substituting the equation above into Equation (36) and eliminating its secular terms, one can have: Now setting and substituting Equations (38) The periodic solution of the system (3) satisfies . b = 0, and . ϕ = 0, i.e., The periodic solution can be expressed analytically as: According to the characteristic solutions of Equation (42), it shows that the theoretical periodic solution expressed by Equation (44) will become unstable if: Based on Equations (43)-(45), the evolution of the periodic solutions of system (3) with the amplitude of AC voltage when ω = 0.6 is shown in Figure 6. Obviously, when V AC increases from 0, the two non-trivial equilibria lose their stability; instead, there are two periodic attractors coexisting. The amplitudes of the two periodic attractors increases with the amplitude of AC voltage. The coexistence of multiple periodic attractors can be attributed to the disturbance of the bistable non-trivial equilibria of the system (3) when V AC = 0 V.
Based on Equations (43)-(45), the evolution of the periodic solutions of system (3) with the amplitude of AC voltage when ω = 0.6 is shown in Figure 6. Obviously, when VAC increases from 0, the two non-trivial equilibria lose their stability; instead, there are two periodic attractors coexisting. The amplitudes of the two periodic attractors increases with the amplitude of AC voltage. The coexistence of multiple periodic attractors can be attributed to the disturbance of the bistable non-trivial equilibria of the system (3) when VAC = 0 V. In Figure 6, when VAC varies from 0 to 0.055 V, the numerical simulation is in great agreement with the theoretical solution. However, when VAC exceeds 0.056 V, the theoretical prediction of the periodic attractor in the neighborhood of the right non-trivial equilibria is not that accurate, which may be due to the limitation of the Method of Multiple Scale. It will then be essential for us to apply numerical simulation to investigate the evolution of the attractors with the change in AC voltage. The basic settings for the simulation, such as the time step and initial plane, are the same as that in Section 3. The change of the attractors and the area and nature of their basins of attraction with VAC are shown in Figure 7, where ω = 0.6. The evolution of global dynamics of system (3) with the increase in VAC can be separated into the following five stages.
Firstly, when VAC = 0 V, there are two point attractors coexisting whose basins of attraction are fractal and trigger each other (see Figure 7a1,a2). According to Figure 7a2, in a small neighborhood of each point attractor, the attractor of system (3) is locally stable. Otherwise, a small disturbance of initial conditions will lead to a different point attractor, meaning that it is easy to induce a safe jump.
Secondly, when VAC increases from 0 to 0.01 V (see Figure 7b1-d2), the number of the periodic attractors increases with VAC. At VAC = 0.005 V, the two point attractors become two periodic attractors; apart from these two periodic attractors predicted theoretically, a new periodic attractor appears suddenly, marked by the yellow curve in Figure 7b1, and its basin of attraction is discrete (see the yellow regions in Figure 7b2). It shows that the new periodic attractor is a hidden attractor [26]. When VAC increases to 0.006 V, another hidden attractor appears, which is almost symmetric to the former one (see the blue curve of Figure 7c1 and the blue regions of Figure 7c2). When VAC = 0.01 V, there are five periodic attractors coexisting, as shown in Figure 7d1. A new periodic attractor appears (see the green curve in Figure 7d1), whose amplitude is much bigger than the other ones.
Thirdly, as VAC increases from 0.01 V to 0.116 V, the number of attractors will decrease. Comparing Figure 7e1 with Figure 7d1, it is obvious that when VAC increases to 0.02 V, the yellow periodic attractor disappears whose basin of attraction is eroded by that In Figure 6, when V AC varies from 0 to 0.055 V, the numerical simulation is in great agreement with the theoretical solution. However, when V AC exceeds 0.056 V, the theoretical prediction of the periodic attractor in the neighborhood of the right non-trivial equilibria is not that accurate, which may be due to the limitation of the Method of Multiple Scale. It will then be essential for us to apply numerical simulation to investigate the evolution of the attractors with the change in AC voltage. The basic settings for the simulation, such as the time step and initial plane, are the same as that in Section 3. The change of the attractors and the area and nature of their basins of attraction with V AC are shown in Figure 7, where ω = 0.6. The evolution of global dynamics of system (3) with the increase in V AC can be separated into the following five stages.
Firstly, when V AC = 0 V, there are two point attractors coexisting whose basins of attraction are fractal and trigger each other (see Figure 7a1,a2). According to Figure 7a2, in a small neighborhood of each point attractor, the attractor of system (3) is locally stable. Otherwise, a small disturbance of initial conditions will lead to a different point attractor, meaning that it is easy to induce a safe jump.
Secondly, when V AC increases from 0 to 0.01 V (see Figure 7b1-d2), the number of the periodic attractors increases with V AC . At V AC = 0.005 V, the two point attractors become two periodic attractors; apart from these two periodic attractors predicted theoretically, a new periodic attractor appears suddenly, marked by the yellow curve in Figure 7b1, and its basin of attraction is discrete (see the yellow regions in Figure 7b2). It shows that the new periodic attractor is a hidden attractor [26]. When V AC increases to 0.006 V, another hidden attractor appears, which is almost symmetric to the former one (see the blue curve of Figure 7c1 and the blue regions of Figure 7c2). When V AC = 0.01 V, there are five periodic attractors coexisting, as shown in Figure 7d1. A new periodic attractor appears (see the green curve in Figure 7d1), whose amplitude is much bigger than the other ones.
Thirdly, as V AC increases from 0.01 V to 0.116 V, the number of attractors will decrease. Comparing Figure 7e1 with Figure 7d1, it is obvious that when V AC increases to 0.02 V, the yellow periodic attractor disappears whose basin of attraction is eroded by that of the green attractor; thus, the basin of attraction of the green attractor can be much bigger in Figure 7e2 than in Figure 7d2. When V AC continues to increase, the other three periodic attractors, i.e., the blue attractor, the red one, and the black one, disappear successively (see Figure 7f1,h1,j1) whose basins of attraction are aggressed by the basin of attraction of the green attractor, as shown in Figure 7e1-j2. Till V AC becomes 0.116 V, there will be a single periodic attractor left whose basin of attraction is not fractal but with a smooth boundary (see Figure 7j1,j2).
Besides, when V AC increases to 0.128 V, there will be a new complex attractor coexisting with the former green periodic attractor. It is a period-3 attractor (see the purple curve in Figure 7k1) whose basin of attraction is fractal and eroded to the basin of attraction of the periodic attractor (see Figure 7k2). It follows that a small change of initial conditions possibly shifts the dynamical behavior of the system (3) from a periodic motion to a period-3 motion, which is another type of safe jump.
Finally, as V AC continues to increase, another type of complex dynamical behavior is induced. According to the phase map, Poincare map, and frequency spectrum in Figure 8a-c, there is only a chaotic attractor when V AC = 0.28 V, and the boundary of its basin of attraction is not fractal (see Figure 8d). attractors, i.e., the blue attractor, the red one, and the black one, disappear successively (see Figure 7f1,h1,j1) whose basins of attraction are aggressed by the basin of attraction of the green attractor, as shown in Figure 7e1-j2. Till VAC becomes 0.116 V, there will be a single periodic attractor left whose basin of attraction is not fractal but with a smooth boundary (see Figure 7j1,j2).
Besides, when VAC increases to 0.128 V, there will be a new complex attractor coexisting with the former green periodic attractor. It is a period-3 attractor (see the purple curve in Figure 7k1) whose basin of attraction is fractal and eroded to the basin of attraction of the periodic attractor (see Figure 7k2). It follows that a small change of initial conditions possibly shifts the dynamical behavior of the system (3) from a periodic motion to a period-3 motion, which is another type of safe jump.
Finally, as VAC continues to increase, another type of complex dynamical behavior is induced. According to the phase map, Poincare map, and frequency spectrum in Figure  8a-c, there is only a chaotic attractor when VAC = 0.28 V, and the boundary of its basin of attraction is not fractal (see Figure 8d).

Conclusions
In this paper, a typical electrostatic bilateral micro-resonator is considered. The theory of local bifurcation and numerical approaches are applied to analyze the global dynamics of the vibrating system of the micro resonator. The main conclusions are presented

Conclusions
In this paper, a typical electrostatic bilateral micro-resonator is considered. The theory of local bifurcation and numerical approaches are applied to analyze the global dynamics of the vibrating system of the micro resonator. The main conclusions are presented as follows: (1) DC bias voltage has some effect on the dynamics of the micro resonator. Without AC voltage, when the DC bias voltage is low, there will be only one stable point attractor in its vibrating system; when the DC bias voltage increases, there may be bistable point attractors. (2) In the case of a low bias DC voltage, multiple periodic attractors and the corresponding safe jump occur due to Hopf bifurcation when varying the frequency or amplitude of AC voltage in certain ranges. (3) Under a higher bias DC voltage that can induce bistable point attractors, when increasing the value of the amplitude of AC voltage, there will be multiple periodic attractors attributed to the loss of stability of the two non-trivial point attractors; apart from this, there will be some other complex dynamical behaviors of the microresonator vibrating system, such as safe jump, hidden attractors, period-n attractor, and chaos.
Our results provide some theoretical reference in avoiding complex dynamics of micro resonators, thus having some potential values in the design of micro sensors. The hidden attractors are depicted numerically, but their mechanism is still not that clear, which will be discussed in our future study.