Next Article in Journal
Approaches to Spherically Symmetric Solutions in f(T) Gravity
Next Article in Special Issue
Vacuum Polarization in a Zero-Width Potential: Self-Adjoint Extension
Previous Article in Journal
A Possible Resolution to Troubles of SU(2) Center Vortex Detection in Smooth Lattice Configurations
Previous Article in Special Issue
No-Slip Boundary Conditions for Electron Hydrodynamics and the Thermal Casimir Pressure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Actuation of Casimir Oscillators toward Chaos: Comparison of Topological Insulators and Metals

1
Department of Physics, Faculty of Physics and Chemistry, Alzahra University, Tehran 1993891167, Iran
2
Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
3
Department of Mechanical Engineering, University of Sistan and Baluchestan, Zahedan 981-35161, Iran
*
Author to whom correspondence should be addressed.
Universe 2021, 7(5), 123; https://doi.org/10.3390/universe7050123
Submission received: 13 April 2021 / Revised: 20 April 2021 / Accepted: 21 April 2021 / Published: 29 April 2021
(This article belongs to the Special Issue The Casimir Effect: From a Laboratory Table to the Universe)

Abstract

:
In the current study, we explore the sensitivity of the actuation dynamics of electromechanical systems on novel materials, e.g., Bi2Se3, which is a well-known 3D Topological Insulator (TI), and compare their response to metallic conductors, e.g., Au, that are currently used in devices. Bifurcation and phase portraits analysis in conservative systems suggest that the strong difference between the conduction states of Bi2Se3 and Au yields sufficiently weaker Casimir force to enhance stable operation. Furthermore, for nonconservative driven systems, the Melnikov function and Poincare portrait analysis probed the occurrence of chaotic behavior leading to increased risk for stiction. It was found that the presence of the TI enhanced stable operation against chaotic behavior over a significantly wider range of operation conditions in comparison to typical metallic conductors. Therefore, the use of TIs can allow sufficient surface conductance to apply electrostatic compensation of residual contact potentials and, at the same time, to yield sufficiently weak Casimir forces favoring long-term stable actuation dynamics against chaotic behavior.

1. Introduction

Nowadays, advancement in fabrication and, consequently, manufacturing techniques of micro/nanoelectromechanical systems (MEMS/NEMS) attract strong attention from the scientific and technology point of view in various sensor technologies, accelerometers, microswitches, etc. [1,2,3,4]. On the other hand, the omnipresent Casimir forces appear to become dominant by reduction of the scale of microdevices into the submicron range, because MEMS/NEMS have surface areas large enough but gaps small enough for this force to play a significant role. As the magnitude of the Casimir force increases, it can lead to the permanent adhesion of moving elements with adjacent surfaces, known as stiction, under certain conditions [4]. Strategies to reduce stiction have been studied in terms of a suitable choice of materials and for the development of knowledge in predicting stable device operation on a long-term basis for the general case of nonconservative driven MEMS/NEMS [5,6,7].
In fact, the Casimir force was predicted in 1948 [8] for perfectly reflecting parallel plates due to the perturbation of vacuum fluctuations of the electromagnetic (EM) field. Soon after, Lifshitz and coworkers considered the general case of dielectric plates by exploiting the fluctuation-dissipation theorem, which relates the dissipative properties of the plates (optical absorption by many microscopic dipoles) and the resulting EM fluctuations. In terms of the Lifshitz theory [9,10], the van der Waals and Casimir forces are the short- and long-range asymptotic limits, respectively, of the same force [1,2,3,4,8,9,10,11,12,13,14]. Nevertheless, the strong dependence of the Casimir force on the material optical properties can be utilized to tune the actuation of devices by proper choice of the interacting materials [5,6,7,15,16,17,18,19]. Several studies have shown that strong Casimir forces exist between components made of metals due to their high absorption of conduction electrons in the infrared range, while the less conductive materials can provide weaker Casimir forces and enhance the stability of microdevices suitable for operation in harsh environments [5,6,7,15,16,17,18,19].
Furthermore, due to these unique properties of their surface states, topological insulator (TI) materials can potentially open a new window of opportunity for MEMS engineering [20,21,22,23,24,25]. Indeed, there has recently been strong interest from both the theoretical and experimental point of view on TI materials, which provide a new quantum state of matter [20,21,22,23,24,25]. TIs show an insulating gap in the bulk but have gapless surface states that are protected topologically. Several recent studies have shown that a repulsive Casmir force can be achieved between two TI plates under certain conditions [23,24,25], which is promising for MEMS design and manufacturing. In these devices, the Casimir force is not only sufficiently strong to play an inevitable role, but the actuation dynamics can also abruptly change toward stiction due to the attractive forces. Hence, any conditions for the occurrence of a repulsive force are of paramount importance for stable motion. In addition, the attractive Casimir forces from TIs can be sufficiently weak to allow the stable operation of actuating devices with suitable compensation of contact potentials due to their conductive surface states [26], while stable operation cannot be achieved using good conductive materials, for example, Au, due to significantly stronger Casimir forces [6,7].
So far, several investigations have been conducted to answer the question of how the Casimir force changes with these novel materials under both trivial TI (with θ=0) and nontrivial TI (θ 0) conditions where θ defines the topological magnetoelectric polarizability (TMEP) [23,24,25]. However, to date, how significant the influence of TIs on stable device actuation and long-term performance remains unexplored, although it is well known that the occurrence of chaotic behavior is unavoidable by shrinking the size of devices.

2. Optical Properties of Materials and Device Actuation

Prior to modeling the microdevice and its actuation dynamics, it is necessary to compute the Casimir force via Lifshitz theory (see Appendix A). For this purpose, the imaginary part of the frequency-dependent dielectric function ε ω is an essential input. Therefore, the optical properties of all samples, including Au, Bi2Se3, and Al2O3, were characterized with ellipsometry using the VUV-VASE (0.5–9.34 eV) and IR-VASE (0.03–0.5 eV) ellipsometers [17,26,27]. Bi2Se3 is a well-known 3D Topological insulator [26], and Au is used in devices due to its high conductivity (   ω p 2 /   ω τ | Au 1600 eV) [17,27]. The imaginary part ε’’(ω) of frequency-dependent dielectric response function ε(ω) of Bi2Se3 and Al2O3 and the corresponding functions at imaginary frequencies ε(iξ) are shown in Figure 1a,b. Since, in practice, the experimental data for the imaginary part ε’’ (ω) of the dielectric function ε(ω) cover only a limited frequency range, some extrapolations were made as shown in Appendix B.
For our actuation analysis, we considered a typical microswitch as shown in Figure 1c. The microswitch was constructed from 2 plates where the upper plate was suspended by a mechanical spring governed by the Hooke’s law while the lower plate was fixed [28]. The elastic restoring force F res = K d z of the spring with stiffness K counterbalanced the attractive Casimir force. It was assumed that both plates were coated with an optically bulk Au film (thickness: ~100 nm) or a 100 nm layer of Bi2Se3 deposited on a bulk sapphire (Al2O3) as described by the authors of [26]. For our calculation, we considered initial separations d between the actuating plates in the range 250 nm to 1 µm, and any thermal effects were not taken into account (assuming, basically, a temperature of 0 K). By considering all acting forces on the movable plate, the equation of motion for the microsystem can be written as:
M d 2 z dt 2 + ε M ω 0 Q dz dt = F res + F Cas + ε   F 0 cos ω t
where M is the mass of the moving plate and (M ω 0 /Q)(dz/dt) is the intrinsic energy dissipation of the actuating system with Q as its quality factor [29,30,31]. The frequency ω was assumed to have a value that is typical for many MEMS or atomic force microscope (AFM) cantilevers (typically ω = 300 kHz) [29,30,31]. The parameter ε was introduced to distinguish between the conservative frictionless autonomous operation of the actuating system ( ε = 0 ) and the nonconservative driven system by an external force ( ε = 1 ) in the presence of friction with a finite quality factor Q. Finally, we considered MEMS components with flat surfaces, because any nanoscale roughness has a significant influence at separations sufficiently below 100 nm [32]. In all cases, the lateral dimensions of the plates were considered to be L x = L y = 10   μ m .

3. Results and Discussion

To proceed with the analysis of the actuation dynamics, we introduced the bifurcation parameter δ Cas = F Cas m / kd , which is the ratio of the minimum Casimir force F Cas m = F Cas z = d to the maximum restoring force Kd . A small change of δ Cas can lead to an abrupt change in the actuation of the system [15]. Substitution of δ Cas yields for the equation of motion the more convenient form:
d 2 λ dT 2 + ε 1 Q d λ dT = 1 λ + δ Cas F Cas F Cas m + ε F 0 F res Max cos ω ω 0 T
with λ = z/d and T = ω 0 t .

3.1. Conservative Actuating System (ε = 0)

The conservative system is discussed first in order to study the equilibrium points of the dynamical system. In fact, the locus of equilibrium points is obtained by the condition F total = F res + F Cas = 0 , since the presence of periodic solution describes regions with sufficient restoring force to prevent stiction of the plates. In this case, Equation (2) yields the condition λ = 1 δ Cas F Cas z / F Cas m d for λ . Figure 2 shows plots of δ Cas vs. λ for different initial distances illustrating how different interaction materials generate considerable change on the stability conditions of the system. Notably, there is a significant difference between the bifurcation curves around the maximum where one approaches critically unstable behavior. In Figure 2, the dashed lines show the stable area where the restoring force F res is strong enough to ensure stable periodic motion. The solid lines indicate unstable regions where the moving plate undergoes instable motion leading to permanent lockdown. When δ Cas < δ Cas MAX , two equilibrium points exist. The equilibrium point closer to λ = 1 (dashed line) is a stable center point, and the one closer to λ = 0 (solid line) is the unstable saddle point. By increasing δ Cas or equivalently weakening the restoring force ( δ Cas ~ 1 / k ), the distance between the equilibrium points decreases until δ Cas reaches the maximum saddle point δ Cas MAX . Clearly, when δ Cas ~ δ Cas MAX for the Au-Au system leads to loss of its stability, it is still δ Cas < δ Cas MAX for the TI-TI system, ensuring the presence of two equilibrium points and enhancing the possibility for stable motion.
Besides the bifurcation diagrams in Figure 2, the Poincare portraits in Figure 3 directly depict the strong sensitivity of actuation dynamics on its initial conditions and demonstrate how the size of the stable area is strongly dependent on the optical properties. In Figure 3, the enclosed area by the homoclinic orbit corresponds to periodic solutions and stable actuation. Otherwise, stiction would occur within one period. The phase portraits in Figure 3 also demonstrate that a system with stronger Casimir attractive force has significantly less phase space available for stable motion as is the case when comparing the distinct behavior of the Au-Au and TI-TI systems.
Because the extrapolation of the dielectric function for Au in the low frequency range was performed by the Drude model, Figure 4 shows the influence of extrapolation in the low frequency range by considering the effect of both the Drude (D) and Plasma (P) models. Due to the strong difference between the conduction states of Au and TI, the stable area for actuation of both the Au-Au and Au-TI systems is smaller than the TI-TI system for both extrapolation models. For the Au-Au system, the use of the Drude or Plasma models to extrapolate at low frequencies induces significant changes in the phase space available for stable actuation. However, for the Au-TI and the TI-TI systems, the extrapolation effect by the D/P models is reduced, e.g., in Au-TI systems, or fully eliminated, e.g., in TI-TI systems. For the latter, no extrapolation is performed for the optical data of the TI due to the absence of any measurable absorption in the far-infrared range (probably due to the insensitivity of the ellipsometry measurements to probe the conductive surface layer of the TI). In this plot, we considered δCas = 0.027 for both the Au-Au and Au-TI systems. According to the bifurcation diagram, by increasing the separation distance d for actuation (where the spring is not stretched during fabrication of the device) and therefore decreasing the value of F Cas m , the chosen value for δCas moves to the maximum of the bifurcation curve (when d changes from 250 nm to 1000 nm) where the sensitivity on the optical properties is most pronounced. As a result, the ratio of the phase space area for stable actuation in Figure 4, SAA/STT and SAT/STT, decreases with increasing separation distance for both extrapolation models at low frequencies beyond the far-infrared range where optical data are no longer available.

3.2. Nonconservative Driven System (ε = 1)

For the more realistic case of externally actuated devices, we performed investigations of the existence of chaotic behavior in dissipative MEMS driven by an external periodic force F o   cos ω   t [33]. Chaotic behavior could occur if the separatrix (homoclinic orbit) of the conservative system splits, and this behavior can be studied by the so-called Melnikov function and Poincare portrait analysis [33,34]. If we define the homoclinic solution of the conservative system as φ hom C T , then the Melnikov function for the oscillating system ( ε = 1 ) is given by [28,29]:
M T 0 = 1 Q + d φ hom C T dT 2 dT + τ 0   τ res MAX   + d φ hom C T dT cos ω ω 0 T + T 0 dT
The separatrix splits if the Melnikov function has simple zeros such that M T 0 = 0 and M T 0 0 . If M T 0 have no zeros. Thus, the device motion will not be chaotic. The conditions of no simple zeros where M T 0 = 0 and M T 0 = 0 provide the threshold condition for possible chaotic motion [33,34]. If we define:
μ hom c = + d φ hom C T dT 2 dT   and   β ω = H   Re   F d φ hom C T dT ,
then the threshold condition for chaotic motion α = β ω / μ hom c with α = 1 / Q F 0   / F res MAX 1 obtains the form:
α = γ ω 0 d F 0 = H   Re   F d φ hom C T dT / + d φ hom C T dT 2 dT
with γ =   M ω o / Q , and H denoting the Hilbert transform [33,34].
Figure 5 shows the threshold Melnikov curves α = γ ω 0   d / F 0 vs. the driving frequency ratio ω / ω o . For relatively large values of α (above the curve), the dissipation dominates the driving force ( α ~ γ / F 0 ), leading to regular motion, which asymptotically approaches the stable periodic orbit of the conservative system. However, for parameter values below the curve, the splitting of the separatrix could lead to chaotic motion. Clearly, for the Au-Au system with the highest conductivity of interacting materials, which lead also to stronger Casimir forces, chaotic motion is more likely to occur as it is manifested by the larger area below the threshold curves. However, for the TI-TI system, this area is decreased significantly, indicating higher possibilities for stable driven actuation. Hence, Figure 5 clearly demonstrates the strong dependence of the region below the threshold curve on the optical properties of the interacting materials for all values of initial distance d prior to initiation of its actuation.
Furthermore, the Poincare portraits in Figure 6 and Figure 7 illustrate the sensitive dependence of the chaotic motion on the optical properties for the Au and TI systems for both different initial distances d and threshold values of the parameter α . When the occurrence of chaotic motion is possible with the decreasing value of α , there is a region of initial conditions where the distinction between qualitatively different solutions is unclear. If we compare Figure 5 with Figure 3 where chaotic motion does not occur, it is revealed that, for chaotic motion, there is no simple smooth boundary between the stable (red regions) and unstable solutions (blue regions). With decreasing the value of the parameter α , the possibility of occurrence of chaotic motion is enhanced as it is manifested by the dramatic shrinkage of the stable solutions (elliptic-like red area) in Figure 7 for the Au-Au system. This is also in strong contrast with the TI-TI system where stable operation is much more efficiently ensured even for very low values of the parameter α within the area favoring unstable motion.
From Figure 6 and Figure 7, it is evident how the chaotic motion of Casimir oscillators can impose a considerable risk of stiction after several periods of oscillation, limiting the long-term prediction of the device operation. This undesirable malfunction occurs more prominently for the more conductive systems since the Casimir force is stronger. On the other hand, for the actuating system where TIs are involved, the measurement of the residual contact potentials Vo, as performed by the authors of [26] for the Au-Bi2Se3 system leading to potential values Vo ≈ 200 mV for TI thicknesses in the range 10–100 nm, allows subsequent voltage compensation (for TI thickness above 30 nm [26]). At the same time, the TI (e.g., Bi2Se3) yields weaker Casimir forces that allow stable operation over a wider area of initial conditions in phase space. Indeed, according to the authors of [6,7], the stable area for operation of actuating devices undergoes further shrinkage with additional electrostatic forces. This effect is very crucial for operation under condition with Melnikov parameter α < 1 (corresponding to strong driving force F o and weak dissipation γ) where, for the Au-Au system, any phase space for stable operation in Figure 7 is absent. Moreover, for the TI-TI system, it has still significant size.

4. Conclusions

In conclusion, we investigated the sensitivity of actuation dynamics of electromechanical systems on novel materials, for example, Bi2Se3, a well-known 3D Topological Insulator (TI), and compared their response to good conductors, e.g., Au. Analysis in conservative systems using bifurcation and phase portraits suggests that the strong difference between the conduction states of Bi2Se3 and Au yield sufficiently weaker Casimir force to enhance stable operation. Furthermore, for nonconservative driven systems, the Melnikov function and Poincare portraits analysis probed the occurrence of chaotic behavior leading to increased risk for stiction. It was found that the presence of the TI enhances stable operation against chaotic behavior over a significantly wider range of operation conditions in comparison to typical metallic conductors. Therefore, the use of TIs can allow sufficient surface conductance to apply electrostatic compensation of any contact potentials [23] and, at the same, to yield sufficient weak Casimir forces to allow long-term prediction of actuation dynamics against chaotic behavior. This effect will be very crucial for operation under Melnikov parameters α < 1 (relatively strong driving force and weak dissipation) where, for good conductors, any phase space for stable operation is absent. For the systems where TIs are involved, operation could have still a significant available range of initial conditions to perform stable actuation on longer-term basis.

Author Contributions

All authors contributed to the research reported in this article. Conceptualization, F.T., Z.B., and G.P.; methodology, F.T., Z.B., and G.P.; software, F.T., Z.B., and M.S., validation, F.T., Z.B., M.S.; formal analysis, F.T., G.P.; investigation, F.T., Z.B., M.S., G.P.; writing—original draft preparation, F.T., and G.P.; writing—review and editing, F.T., Z.B., M.S., and G.P.; funding acquisition, F.T., and G.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partly funded by the Netherlands Organization for Scientific Research (NWO) under grant number 16PR3236.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The computer-generated data used in this article are available by request to the authors.

Acknowledgments

G.P. acknowledges support from the Netherlands Organization for Scientific Research (NWO) under grant number 16PR3236. F.T. acknowledges support from the Department of Physics at Alzahra University. We would like to acknowledge useful discussions with D.T. Yiman for the analysis of the optical data for the Bi2Se3 system.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Lifshitz Theory for Stratified Media

The Casimir force F Cas d between two plates which are covered with a thin layer of Bi2Se3 on a bulk of Al2O3 (Figure 1) in Equation (2) is given by [8,9,10]:
F Cas d = ћ 2 π 2   0 d ξ   ν = TE , TM   0 dk   k   k 0 exp 2 k 0 d R ν + R ν 1 1   .  
ξ are the imaginary frequencies. The reflection coefficients between the plates are given by [4,9,10]:
R TM + i ξ , k   = r TM 0 ,   + B +   r TM + A , + B   exp 2   k   + B   a 1 +   r TM 0 ,   + B +   r TM + A , + B   exp 2   k   + B   a  
R TE + i ξ , k   = r TE 0 ,   + B +   r TE + A , + B   exp 2   k   + B   a 1 +   r TE 0 ,   + B +   r TE + A , + B   exp 2   k   + B   a ,
where ‘a’ is thickness of Bi2Se3, and the index of +B (-B) and +A (-A) define regions that correspond to Bi2Se3 and Al2O3 for a plate at the right (left) side, respectively. The vacuum between the plates is referred by the index 0. The reflection coefficients on the various boundaries are given by   r TE n ,   m = k n   k m / k n +   k m and r TM n ,   m = ε m   k n   ε n   k m / ε m   k n +   ε n   k m for the transverse electric (TE) and magnetic (TM) field polarizations, respectively. k i = ε i   i ξ   ξ 2 / c 2 +   k 2 i = 0 , 1 , 2   represents the out-off plane wave vector in the gap between the interacting plates ( k 0 ) and in each of the interacting plates ( k i = 1 , 2 ). k is the in-plane wave vector.
Furthermore, ε i ξ is the dielectric function evaluated at imaginary frequencies, which is considered as a vital input to calculate the Casimir force between real materials using the Lifshitz theory. Applying the Kramers–Kronig relation, ε(iξ) is given by [9,10]:
ε i ξ = 1 + 2 π 0 ω ε ω ω 2 +   ξ 2   d ω .

Appendix B. Dielectric Function of Materials with Extrapolations

For the calculation of the integral in Equation (A4), the measured data for the imaginary part of the frequency dependent dielectric function ε ω are required. In this study, for the Au sample, the experimental data of the imaginary part of the dielectric function cover only a limited range of frequencies ω 1   = 0.03   ev <   ω   <   ω 2   = 8.9   ev . Therefore, for the low optical frequencies ( ω < ω 1 ), we extrapolated using the Drude model [9,10,17,35]:
ε L ω = ω p 2   ω τ ω   ω 2 +   ω τ 2
where ω p is the plasma frequency and ω τ is the relaxation frequency. For the high optical frequencies ( ω > ω 2 ), we extrapolated using the expression [9,10,17,35]:
ε H ω = A ω 3
Using Equations (A4)–(A6), the function ε i ξ in terms of the Drude model is given by [17,35]:
ε i ξ D = 1 + 2 π + ω 1 ω 2 ω ε exp ω ω 2 + ξ 2   d ω + Δ L ε i ξ + Δ H ε i ξ
with
Δ L ε i ξ = 2 π 0 ω 1 ω ε L ω ω 2 + ξ 2   d ω = 2 ω p 2   ω τ π ξ 2 ω τ 2 arctan ω 1 ω τ ω τ arctan ω 1 ξ ξ
and
Δ H ε i ξ = 2 π ω 2 ω ε H ω ω 2 +   ξ 2   d ω = 2 ω 2 3 ε ( ω 2 ) π ξ 2 1 ω 2 π 2   arctan ω 2 ξ ξ
Finally, for the plasma model, one must replace the term Δ L ε i ξ in Equation (A8) with ω p 2 / ξ 2 . Therefore, for the plasma model, ε i ξ is given by the simpler expression:
ε i ξ P = 1 + 2 π ω 1 ω 2 ω ε exp ω ω 2 + ξ 2 d ω + ω p 2 ξ 2 + Δ H ε i ξ .
For the case of Bi2Se3 and Al2O3, there is no any extrapolation because they do not have any measurable Drude tail indicating absorption for the imaginary part at low frequencies.

References

  1. Rodriguez, A.W.; Capasso, F.; Johnson, S.G. The Casimir effect in microstructured geometries. Nat. Photonics 2011, 5, 211. [Google Scholar] [CrossRef]
  2. Capasso, F.; Munday, J.N.; Iannuzzi, D.; Chan, H.B. Casimir Forces and Quantum Electrodynamical Torques: Physics and Nanomechanics. IEEE J. Sel. Top. Quant. Electron. 2007, 13, 400. [Google Scholar] [CrossRef]
  3. Ball, P. Fundamental physics: Feel the force. Nature 2007, 447, 77. [Google Scholar] [CrossRef]
  4. Bordag, M.; Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. Advances in the Casimir Effect; Oxford University Press: New York, NY, USA, 2009. [Google Scholar]
  5. Tajik, F.; Sedighi, M.; Palasantzas, G. Sensitivity on materials optical properties of single beam torsional Casimir actuation. J. Appl. Phys. 2017, 121, 174302. [Google Scholar] [CrossRef] [Green Version]
  6. Tajik, F.; Sedighi, M.; Khorrami, M.; Masoudi, A.A.; Palasantzas, G. Chaotic behavior in Casimir oscillators: A case study for phase-change materials. Phys. Rev. E 2017, 96, 042215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Tajik, F.; Sedighi, M.; Khorrami, M.; Masoudi, A.A.; Waalkens, H.; Palasantzas, G. Dependence of chaotic behavior on optical properties and electrostatic effects in double-beam torsional Casimir actuation. Phys. Rev. E 2018, 98, 02210. [Google Scholar] [CrossRef] [Green Version]
  8. Casimir, H.B.G. Zero Point Energy Effects on Quantum Electrodynamics. Proc. K. Ned. Akad. Wet. 1948, 51, 793. [Google Scholar]
  9. Lifshitz, E.M. The Theory of Molecular Attractive Forces between Solids. Sov. Phys. JETP 1956, 2, 73. [Google Scholar]
  10. Dzyaloshinskii, I.E.; Lifshitz, E.M.; Pitaevskii, L.P. General theory of van der waals forces. Sov. Phys. Uspekhi 1961, 4, 153. [Google Scholar] [CrossRef]
  11. Decca, R.S.; L’opez, D.; Fischbach, E.; Klimchitskaya, G.L.; Krause, D.E.; Mostepanenko, V.M. Novel constraints on light elementary particles and extra-dimensional physics from the Casimir effect. Eur. Phys. J. C 2007, 51, 963. [Google Scholar] [CrossRef] [Green Version]
  12. Decca, R.S.; López, D.; Fischbach, E.; Klimchitskaya, G.L.; Krause, D.E.; Mostepanenko, V.M. Tests of new physics from precise measurements of the Casimir pressure between two gold-coated plates. Phys. Rev. D 2007, 75, 077101. [Google Scholar] [CrossRef] [Green Version]
  13. Craighead, H.G. Nanoelectromechanical systems. Science 2000, 290, 1532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bochobza-Degani, O.; Nemirovsky, Y. Modeling the pull-in parameters of electrostatic actuators with a novel lumped two degrees of freedom pull-in model. Sens. Actuators A 2002, 97–98, 659. [Google Scholar]
  15. Sedighi, M.; Broer, W.H.; Palasantzas, G.; Kooi, B.J. Sensitivity of micromechanical actuation on amorphous to crystalline phase transformations under the influence of Casimir forces. Phys. Rev. B 2013, 88, 165423. [Google Scholar] [CrossRef] [Green Version]
  16. Sedighi, M.; Palasantzas, G. Casimir and hydrodynamic force influence on microelectromechanical system actuation in ambient conditions. Appl. Phys. Lett. 2014, 104, 074108. [Google Scholar] [CrossRef]
  17. Sedighi, M.; Svetovoy, V.B.; Broer, W.H.; Palasantzas, G. Casimir forces from conductive silicon carbide surfaces. Phys. Rev. B 2014, 89, 195440. [Google Scholar] [CrossRef]
  18. Tajik, F.; Sedighi, M.; Babamahdi, Z.; Masoudi, A.A.; Waalkense, H.; Palasantzas, G. Dependence of non-equilibrium Casimir forces on material optical properties towards chaotic motion during device actuation. Chaos 2019, 29, 093126. [Google Scholar] [CrossRef] [PubMed]
  19. Tajik, F.; Sedighi, M.; Babamahdi, Z.; Masoudi, A.A.; Waalkense, H.; Palasantzas, G. Sensitivity of non-equilibrium Casimir forces on low frequency optical properties towards chaotic motion of microsystems. Chaos 2020, 30, 023108. [Google Scholar] [CrossRef]
  20. Hasan, M.Z.; Kane, C.L. Colloquium: Topological insulators. Rev. Mod. Phys. 2010, 82, 3045. [Google Scholar] [CrossRef] [Green Version]
  21. Moore, J.E. The birth of topological insulators. Nature 2010, 464, 194. [Google Scholar] [CrossRef]
  22. Hsieh, D.; Qian, D.; Wray, L.; Xia, Y.; Hor, Y.; Cava, R.J.; Hasan, M. A Topological Dirac insulator in a quantum spin Hall phase. Nature 2008, 452, 970. [Google Scholar]
  23. Grushin, A.G.; Corteijo, A. Tunable Casimir Repulsion with Three-Dimensional Topological Insulators. Phys. Rev. Lett. 2011, 106, 020403. [Google Scholar] [CrossRef]
  24. Grushin, A.G.; Rodriguez-Lopez, P.; Corteijo, A. Effect of finite temperature and uniaxial anisotropy on the Casimir effect with three-dimensional topological insulators. Phys. Rev. B 2011, 84, 045119. [Google Scholar] [CrossRef] [Green Version]
  25. Martinez, J.C.; Jalil, M.B.A. Tuning the Casimir force via modification of interface properties of three-dimensional topological insulators. J. Appl. Phys. 2013, 113, 204302. [Google Scholar] [CrossRef]
  26. Babamahdi, Z.; Svetovoy, V.B.; Yimam, D.T.; Kooi, B.J.; Banerjee, T.; Moon, J.; Enache, S.; Oh, M.; Stöhr, M.; Palasantzas, G. Casimir and electrostatic forces from Bi2Se3 thin films of varying thickness. Phys. Rev. B Lett. 2021, 103, L161102. [Google Scholar]
  27. Svetovoy, V.B.; van Zwol, P.J.; Palasantzas, G.; de Hosson, J.T.M. Optical properties of gold films and the Casimir force. Phys. Rev. B 2008, 77, 035439. [Google Scholar] [CrossRef] [Green Version]
  28. Pelesko, J.A.; Bernstein, D.H. Modeling MEMS and NEMS; Chapman & Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
  29. Garcıa, R.; Perez, R. Dynamic atomic force microscopy methods. Surf. Sci. Rep. 2002, 47, 197. [Google Scholar] [CrossRef]
  30. Li, M.; Tang, H.X.; Roukes, M.L. Ultra-sensitive NEMS-based cantilevers for sensing, scanned probe and very high-frequency applications. Nat. Nanotechnol. 2007, 2, 114. [Google Scholar] [CrossRef] [PubMed]
  31. Mamin, H.J.; Rugar, D. Sub-attonewton force detection at millikelvin temperatures. Appl. Phys. Lett. 2011, 79, 3358. [Google Scholar] [CrossRef]
  32. Torricelli, G.; van Zwol, P.J.; Shpak, O.; Palasantzas, G.; Svetovoy, V.B.; Binns, C.; Kooi, B.J.; Jost, P.; Wuttig, M. Casimir force contrast between Amorphous and Crystalline Phases of AIST. Adv. Funct. Mater. 2012, 22, 3729. [Google Scholar] [CrossRef]
  33. Broer, W.; Waalkens, H.; Svetovoy, V.B.; Knoester, J.; Palasantzas, G. Nonlinear actuation dynamics of driven Casimir oscillators with rough surfaces. Phys. Rev. Appl. 2015, 4, 054016. [Google Scholar] [CrossRef] [Green Version]
  34. Hirsch, M.W.; Smale, S.; Devaney, R.L. Differential Equations, Dynamical Systems, and an Introduction to Chaos; Elsevier Academic Press: San Diego, CA, USA, 2004. [Google Scholar]
  35. Sedighi, M.; Svetovoy, V.B.; Palasantzas, G. Casimir force measurements from silicon carbide surfaces. Phys. Rev. B 2016, 93, 085434. [Google Scholar] [CrossRef]
Figure 1. (a) Dielectric functions at imaginary frequency for Au, Al2O3, and Bi2Se3. (b) Imaginary part of the frequency-dependent dielectric functions for Au, Al2O3, and Bi2Se3 obtained by ellipsometry. (c) Schematic of the MEMS made of Al2O3 covered with an Au or Bi2Se3 layer.
Figure 1. (a) Dielectric functions at imaginary frequency for Au, Al2O3, and Bi2Se3. (b) Imaginary part of the frequency-dependent dielectric functions for Au, Al2O3, and Bi2Se3 obtained by ellipsometry. (c) Schematic of the MEMS made of Al2O3 covered with an Au or Bi2Se3 layer.
Universe 07 00123 g001
Figure 2. Bifurcation diagrams   δ Cas vs. λ (=z/d) for different interacting materials. The initial distance d between the plates is (a) 250 nm, (b) 500 nm, and (c) 1000 nm. The dashed and solid lines represent the unstable and stable equilibrium points, respectively.
Figure 2. Bifurcation diagrams   δ Cas vs. λ (=z/d) for different interacting materials. The initial distance d between the plates is (a) 250 nm, (b) 500 nm, and (c) 1000 nm. The dashed and solid lines represent the unstable and stable equilibrium points, respectively.
Universe 07 00123 g002
Figure 3. Contour plots of the transient times to stiction in the phase plane d λ / dt vs. λ for systems made of several pairs of interacting materials. The initial distance between the plates is indicated in the plots. Here, we considered δ Cas = 0.065 for d = 250 nm and δ Cas = 0.031 for d = 1000 nm. For the calculations, we used 150 × 150 initial conditions ( λ ,   d λ / dt ). The red (elliptical in shape central area) region surrounded by the homoclinic orbit contains the initial conditions that lead to stable oscillations.
Figure 3. Contour plots of the transient times to stiction in the phase plane d λ / dt vs. λ for systems made of several pairs of interacting materials. The initial distance between the plates is indicated in the plots. Here, we considered δ Cas = 0.065 for d = 250 nm and δ Cas = 0.031 for d = 1000 nm. For the calculations, we used 150 × 150 initial conditions ( λ ,   d λ / dt ). The red (elliptical in shape central area) region surrounded by the homoclinic orbit contains the initial conditions that lead to stable oscillations.
Universe 07 00123 g003
Figure 4. Ratio of the stable area of actuation between different systems (Au-Au and Au-TI) and the TI-TI system vs. initial separation distance for the extrapolations at low frequencies for Au via the Drude (D) and Plasma (P) models. For the calculations, we considered δ Cas = 0.027 .
Figure 4. Ratio of the stable area of actuation between different systems (Au-Au and Au-TI) and the TI-TI system vs. initial separation distance for the extrapolations at low frequencies for Au via the Drude (D) and Plasma (P) models. For the calculations, we considered δ Cas = 0.027 .
Universe 07 00123 g004
Figure 5. Threshold curve α   = γ ω 0   d / F 0 vs. driving frequency ω / ω o (with ω o the natural frequency of the system). The area below the curve defines the condition that can possibly lead to chaotic motion. The values of δ Cas and the initial distances d considered here were (a) δ Cas = 0.065 and d = 250 nm and (b) δ Cas = 0.031 and d = 1   μ m .
Figure 5. Threshold curve α   = γ ω 0   d / F 0 vs. driving frequency ω / ω o (with ω o the natural frequency of the system). The area below the curve defines the condition that can possibly lead to chaotic motion. The values of δ Cas and the initial distances d considered here were (a) δ Cas = 0.065 and d = 250 nm and (b) δ Cas = 0.031 and d = 1   μ m .
Universe 07 00123 g005
Figure 6. Contour plot of the transient times to stiction in the phase plane d λ / dt vs. λ for Scheme 250. nm, α = 0.12 , ω / ω o = 0.5, and δ Cas = 0.065 . For the calculations, we used 150 × 150 initial conditions ( λ ,   d λ / dt ). The stable actuation (red) region increased significantly from Au-Au to TI-TI systems.
Figure 6. Contour plot of the transient times to stiction in the phase plane d λ / dt vs. λ for Scheme 250. nm, α = 0.12 , ω / ω o = 0.5, and δ Cas = 0.065 . For the calculations, we used 150 × 150 initial conditions ( λ ,   d λ / dt ). The stable actuation (red) region increased significantly from Au-Au to TI-TI systems.
Universe 07 00123 g006
Figure 7. Contour plot of the transient times to stiction in the phase plane d λ / dt vs. λ for systems made of different pairs of materials. The initial distance be 1000 nm, ω / ω o = 0.5, and δ Cas = 0.031 . The parameter α is indicated in the plots for both columns. For the calculations, we used 150 × 150 initial conditions ( λ ,   d λ / dt ).
Figure 7. Contour plot of the transient times to stiction in the phase plane d λ / dt vs. λ for systems made of different pairs of materials. The initial distance be 1000 nm, ω / ω o = 0.5, and δ Cas = 0.031 . The parameter α is indicated in the plots for both columns. For the calculations, we used 150 × 150 initial conditions ( λ ,   d λ / dt ).
Universe 07 00123 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tajik, F.; Babamahdi, Z.; Sedighi, M.; Palasantzas, G. Nonlinear Actuation of Casimir Oscillators toward Chaos: Comparison of Topological Insulators and Metals. Universe 2021, 7, 123. https://doi.org/10.3390/universe7050123

AMA Style

Tajik F, Babamahdi Z, Sedighi M, Palasantzas G. Nonlinear Actuation of Casimir Oscillators toward Chaos: Comparison of Topological Insulators and Metals. Universe. 2021; 7(5):123. https://doi.org/10.3390/universe7050123

Chicago/Turabian Style

Tajik, Fatemeh, Zahra Babamahdi, Mehdi Sedighi, and George Palasantzas. 2021. "Nonlinear Actuation of Casimir Oscillators toward Chaos: Comparison of Topological Insulators and Metals" Universe 7, no. 5: 123. https://doi.org/10.3390/universe7050123

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop