Equilibrium forms branching of a nanolayers system

A. V. Lukin1 , I. A. Popov2 , D. Yu. Skubov3 , L. V. Shtukin4

1, 2, 3, 4SPbPU, Saint-Petersburg, Russia

3, 4IPME RAS, Saint-Petersburg, Russia

3Corresponding author

Vibroengineering PROCEDIA, Vol. 8, 2016, p. 517-521.
Received 7 September 2016; accepted 12 September 2016; published 7 October 2016

Copyright © 2016 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF Acknowledgements References
Cite this article
Views 15
Reads 7
Downloads 973

In this work we analyze equilibrium forms branching of three parallel conductive micro/nanolayers which can be used as a sensor of internal or external mechanical and electrical parameters such as conductivity, tension force or gap between layers. Alteration of physical parameters can lead to symmetry breaking, i.e. asymmetrical deflection of outermost layers and nonzero deflection of middle, initially neutral layer. By slow current intensity variation in one of outermost layers, initial symmetry of the system can be achieved, thus determining the parameters alteration which had occurred.

Keywords: conductive micro/nanolayers, equilibrium forms branching, MEMS, NEMS, symmetry breaking.

1. Introduction

In recent years, great interest of physicists, biologists and electrical engineers aroused the development of micro- and nanotechnology due to possibility of sensors production capable for nano- and microscale measurements of physical and biological parameters [1-6] such as molecular weight, quantum state, properties of biochemical reactions and others. Nanomechanical sensors have supplied sensitivity to extremely small forces and masses of order zeptonewtons (1 zN = 10-21 N) and zeptograms correspondingly. Such sensitivity is caused by achievement of Mega and Giga-hertz frequency range of MEMS internal dynamic characteristics. For example, the silicone element 1 uM in length, 0.12 uM in width and with thickness of 0.075 uM has a natural frequency 1.03 GHz a Q-factor of order 1014 [1]. Not less widespread are micro sensors based on “pull-in” effect, which gave an opportunity to identify microchanges of external or internal physical characteristics [7].

In this paper we study the equilibrium forms branching of three parallel micro/nanolayers at direct current (Fig. 1). This type of structure is really feasible at the nanoscale [2]. This micro/nano structure can be used as a sensor of external and internal mechanical and electrical parameters or as a nanoresonator for measuring the mass of adhering nanoparticles [8, 9]. It is important to note that in the case of micro/nano schemes non-linear electromechanical interaction in them must be taken into account, as it is the basis of the mechanical or electrical characteristics identification.

Fig. 1. Scheme of a system of nanolayers

Scheme of a system of nanolayers

Physical parameters alteration can lead to symmetry breaking in the system, i.e. asymmetrical deflection of outermost layers and nonzero deflection of middle, initially neutral layer. By slow current intensity variation in one of outermost layers, initial symmetry of the system can be achieved, thus determining the parameters alteration which had occurred.

We will estimate graphene nanostrings (one-dimensional layers) parameters [2], at which significant interaction between layers is possible: layer length l= 1000 nm, width δ= 10 nm, thickness b= 0.3 nm. Thus, the resonant frequency of the vibrating string is 30 MHz with Q ~ 100 [6]. Using the expression for the natural frequency ω- of a stretched string, we obtain the required tension force T=ρVSl2ω-2/π2= 1,5∙10-12 N. Introducing the parameter λ2=μ0I2/2πTh2 which characterizes the force interaction between layers, and assuming that its value is of order of the gap h~ 10 nm, we obtain the value of current I0=hλ/lπT/μ0= 40 mkA.

2. Mathematical formulation

Introducing the notation of direct currents in layers in case of their equality I1=I2=I3=I=const, we obtain the approximate expression for distributed attracting forces between layers excluding edge effects as q=μ0I22πrr2, where r – vector of distance between the layers which is directed normally to them [12]. Not taking into account the interaction between outermost layers, we get equilibrium equations in the form:

T u 1 ' ' + q 12 = 0 ,             T u 2 ' ' + q 21 + q 23 = 0 ,           T u 3 ' ' + q 32 = 0 .

where T is tension force, qij – attractive force acting from j conductor to i, ui – vertical deflections in layers; positive values of all forces and deflections correspond to vertical downwards direction.

In general case of different currents in layers attractive forces can be written in the form:

q i j = μ 0 I i I j 2 π 1 h + u j - u i ,           q j i = - q i j .

where h is an initial gap between layers. We introduce dimensionless coordinates and displacements s=x/l, νi=ui/h and main physical parameter λ2=μ0I2l2/2πTh2. Thus we get system of equilibrium equations in non-dimensional form:

v 1 ' ' + α 1 λ 2 1 + v 2 - v 1 = 0 ,         v 2 ' ' - α 1 λ 2 1 + v 2 - v 1 + α 2 λ 2 1 + v 3 - v 2 = 0 ,         v 3 ' ' - α 2 λ 2 1 + v 3 - v 2 = 0 ,

where the parameters α1, α2 characterize the deviation from symmetry, which can be caused, for example, by change of currents.

Assuming symmetrical deformed state relatively to central section of three layers (i.e., section normal to layers at s= 0.5), we come to a conclusion that in every section through layers length the sum of first derivatives of deflections is equal to zero. Therefore, due to zero boundary conditions, sum of deflections of three layers is equal to zero through the length: v1+v2+v3=0.

Introducing deflections differences z1=v1-v2, z2=v2-v3 and substracting Eq. (2) from Eq. (1) and Eq. (3) from Eq. (2), we obtain a system of two nonlinear Eq. (4):

z 1 ' ' + 2 α 1 λ 2 1 - z 1 - α 2 λ 2 1 - z 2 = 0 ,           z 2 ' ' + 2 α 2 λ 2 1 - z 2 - α 1 λ 2 1 - z 1 = 0 ,

with zero boundary conditions. Deflections of three layers in symmetrical case can be thus written in the form:

v 1 = 2 z 1 + z 2 3 ,           v 2 = z 2 - z 1 3 ,           v 3 = - 2 z 2 + z 1 3 .

The achieved boundary value problem (Eq. (4)) in case of equal currents (α1=α2) allows symmetrical solution with non-deformed middle layer and outermost layers equally deflected one towards another: v2=0, v1=-v3. Then we have z1=z2=z and two equations in Eq. (4) transform to one: z''+λ2/1-z=0 with zero boundary conditions z0=0, z1=0. The above-named problem is equivalent to already considered problem of two nanostrings system equilibria at direct current [12]. Its bifurcation diagram contains two equilibrium forms (stable one with smaller amplitude and non-stable with bigger amplitude) at λ<λ* 1.53, which disappear at λ=λ*, what physically corresponds to coalescence of layers.

System of equations Eq. (4) allows for energy integral which is symmetrical relatively to z1, z2:

( z 1 ' ) 2 + ( z 2 ' ) 2 + z 1 ' z 2 ' - 3 λ 2 ln 1 - z 1 - 3 λ 2 ln 1 - z 2 = h = c o n s t .

Searching partial solution z1=z2=z, the energy integral is transformed to (z')2-2λ2ln1-z=h=const, which correspond to the problem of two nanostrings.

3. Results and discussion

The obtained boundary value problem (Eq. (3)) is solved numerically, using the bvp4c Matlab subprogram. In Fig. 2 there are shown calculated symmetrical and non-symmetrical nanolayers equilibrium forms at different values of parameter λ.

Fig. 2. Forms of equilibrium

Forms of equilibrium

Analytical solution is achieved by Galerkin method solving system of Eq. (4) in assumption of deflection symmetry. Projection conditions of Galerkin method are obtained by searching the solution in the form zi=xisinπs, i= 1, 2 and integrating equations in Eq. (4), multiplied by zi, along layers length:

π 2 2 2 x 1 + x 2 - 3 α 1 λ 2 π x 1 π 1 1 - x 1 2 - 1 + 2 arcsin x 1 1 - x 1 2 = 0 ,
π 2 2 2 x 2 + x 1 - 3 α 2 λ 2 π x 2 π 1 1 - x 2 2 - 1 + 2 arcsin x 2 1 - x 2 2 = 0 .

Graphical solution of Eq. (7) for various values of λ is shown in Fig. 3 (grey and dark lines keep with first and second equation in Eq. (7) correspondingly).

Fig. 3. Graphical solution of Eq. (7) for various values of λ

Graphical solution of Eq. (7) for various values of λ

Numerically calculated bifurcation diagram of general system of equations (Eq. (3)) in case of equal currents is shown in Fig. 4, where deflection norms (maximum of absolute values of deflections along layer length) of one outermost (v1) and middle (v3) layers are presented in dependence of parameter λ. The notations v1, v3 and v5 for top, middle and bottom layer result from rewriting ODE system (Eq. (3)) in normal form. Solution is obtained using AUTO 07P, a software for continuation and bifurcation problems [13].

Fig. 4. Bifurcation diagram

Bifurcation diagram

For λ<λ**1.34 there are three qualitatively different equilibrium forms: one stable form with non-deflected middle layer and symmetrically deflected outermost layers; and two non-stable forms with deflection of middle layer and smaller or bigger deflection of outermost layers (for example, for λ0v1 0.5, v3 0.5, v5 0 and vice versa for v1 and v5 – see Fig. 4, branches 1, 2, 3 for v1 and branches 1, 2 for v3). For λ**<λ<λ* 1.53 non-deflected state of middle layer becomes non-stable and aforenamed non-symmetrical forms don’t exist. For λ>λ* there are no equilibrium forms in the system, what physically corresponds to coalescence of layers. It is remarkable that for λ<λ* there also exist equilibria with outermost layers (v1 or v5) deflection norm higher than 0.5 (see Fig. 4, branches 1,4-7 for v1 and branches 1, 3, 4 for v3). Among these forms exist non-symmetrical ones relatively to central section of layers, i.e. symmetry breaking occurs (see Fig. 4, branches 4, 5 for v1 and branches 3 for v3). Naturally, these forms could be found only by solving general boundary value problem (Eq. (3)), as the system Eq. (7) already contains assumption about central symmetry of solutions. It is notable that analytical solution by Galerkin method shown good agreement with numerical results.

4. Conclusions

In this work the branching theory was applied for the analysis of new type of nano- micro-electromechanical system. The obtained structure of branching can be used for determination of alteration of nano- and micro-scale physical parameters of the system.


Work is performed with assistance of a Grant No. RFBR 14-01-00845.


  1. Eom Kilho, Park Harold S., Yoon Dae Sung, Kwon Taeyun Nanomechanical resonators and their applications in biological/chemical detection. Nanomechanics principles. Physics Reports, Vol. 503, 2011, p. 115-163. [Search CrossRef]
  2. He X. Q., Kitipornchai S., Liew K. M. Resonance analysis of multi-layered graphene sheets used as nanoscale resonators. Nanotechnology, Vol. 16, 2005, p. 2086-2091. [Search CrossRef]
  3. Lui Yilun, Xu Zhipping, Zheng Quanshui The interlayer shear effect on grapheme multilayer resonators. Journal of Mechanics and Physics of Solids, Vol. 59, 2011, p. 1613-1622. [Search CrossRef]
  4. Natsuki Toxhiaki, Shi Jin-Xing, Ni Qing-Quing Vibration analysis of nanomechamical mass sensors using double-layered grapheme sheets resonators. Journal of Applied Physics, Vol. 114, 2013, p. 904307. [Search CrossRef]
  5. van der Zande А.М., et al. Large-scale arrays of single-layer graphene resonators. Nano Letters, Vol. 10, 2010, p. 4869-4873. [Search CrossRef]
  6. Chen C., Hone J. Graphene nanoelectromechanical systems. Proceedings of the IEEE. Vol. 101, Issue 7, 2013, p. 1766-1779. [Search CrossRef]
  7. Zhang Wen-Ming, Yan Han, Peng Zhi-Ke, Meng Guang Electrostatic pull-in instability in MEMS/NEMS: a review. Sensors and Actuators A: Physical, Vol. 214, 2014, p. 187-218. [Search CrossRef]
  8. Morozov N. F., Berinsky I. E., Indeitsev D. A., Skubov D. Yu., Shtukin L. V. Differential Graphene Resonator. Doklady Physics, Vol. 59, Issue 7, 2014, p. 37-40, (in Russian). [Search CrossRef]
  9. Shtukin L. V., Berinsky I. E., Indeitsev D. A., Morozov N. F., Skubov D. Yu. Electromechanic models of nanoresonators. Physical Mesomechanics, RAS, Vol. 19, Issue 1, 2016, p. 24-30, (in Russian). [Search CrossRef]
  10. Bolotin V. V. The Dynamic Stability of Elastic Systems. Holden-Day, San Francisco, 1964. [Search CrossRef]
  11. Keller Joseph B., Antman Stuart Bifurcation Theory and Nonlinear Eigenvalue Problems. W. A. Benjamin, New York, 1969. [Search CrossRef]
  12. Skubov D. Yu., Khodzhaev K. Sh. Non-linear Electromechanics. Springer, 2008. [Search CrossRef]
  13. Doedel E. J., Oldeman B. E. AUTO-07P: Continuation and Bifurcation Software for Ordinary Differential Equations. Concordia University, Montréal, Canada, 2009. [Search CrossRef]
JVE Journals is rebranding to Extrica

Inspired by innovations from the previous century and the rapid growth during the last years, we are improving for excellence in your publishing experience

Read to know more
JVE Journals