Published: 15 February 2015

Electromechanical fluidic coupled vibrations for micro film considering the van der Waals force

Lizhong Xu1
Xiaorui Fu2
1, 2Mechanical Engineering Institute, Yanshan University, Qinhuangdao, 066004, China
Corresponding Author:
Lizhong Xu
Views 25
Reads 9
Downloads 1066

Abstract

In this paper, a multi-fields coupled dynamic equation of the micro film is presented which includes mechanical force, the van der Waals force, electrostatic force and fluidic damping force. The multi-fields coupled free vibration of the micro film is investigated. The multi-fields coupled forced response of the micro film near natural frequency is analyzed. The multi-fields coupled natural frequencies of the micro film and their changes along with system parameters are discussed. Amplitude frequency characteristics of the micro film under the coupled fields are given.

1. Introduction

Microelectromechanical Systems (MEMS) are the integrated systems consisting of microelectronics, microactuators and microsensors [1]. A typical MEMS consist of micro beams or micro films. An example is a small rectangular silicon film with sides in the order of mm and thickness of the order of microns, that deforms when subject to electric fields. Owing to its small size, significant forces and deformations can be obtained with the application of low voltages. Examples of devices that utilize vibrations of such films are synthetic microjets, and microspeakers, etc [2]. Some devices are proposed for generating electrical power from mechanical energy by resonance vibration of the micro beam or micro film [3-5].

Dynamics of the micro film is an important subject that should be developed. Rafael Nadal-Guardia used one-dimensional (1-D) lumped model of parallel film electrostatic transducers and studied its dynamic behavior [6]. Lu used a thin film model to interpret the thickness-dependent sum frequency generation spectra collected from the air/silica/poly (n-butyl methacrylate) /water system [7]. Akai experimentally investigated the vibration mode and transmission ultrasonic waves by pMUTs using epitaxial Pb(Zr,Ti)O3 thin films on the epitaxial γ-Al 2O3/Si substrates to radiate ultrasonic waves effectively [8]. Li presented an electrostatic film vibration voltage transducer [9]. Gualdino fabricated a thin-film silicon micro resonator by surface micromachining and studied vibration modes of micromechanical disk resonators made from hydrogenated amorphous silicon thin films [10]. Sorokin studied the dynamics of a circular film membrane with attached current-carrying conductors in zero gravity and calculated the spectrum of natural vibrations [11].

In a word, a number of the studies about dynamics of the micro film have been done. However, electrostatic excitation is the most frequently applied principle combining versatility and simple technology. In a MEMS device excited by the electrostatic force, the micro film is in an electromechanical-fluidic coupled field. If the clearance between the film and basement is small enough, the effects of the van der Waals force become obvious and could not be neglected. The dynamics problem of a micro film under electromechanical-fluidic-van der Waals force coupled fields has not been resolved yet.

In this paper, a multi-fields coupled dynamic equation of the micro film is presented which includes mechanical force field, the van der Waals force field, electrostatic field and fluidic damping force field. The multi-fields coupled free vibration and the amplitude frequency characteristics of the micro film are investigated. The multi-fields coupled natural frequencies of the micro film and their changes along with system parameters are discussed. A number of results are given which can be used in design and manufacture of the MEMS units.

2. Multi-fields coupled dynamics equations

Fig. 1 illustrates a micro film under multi-fields forces. Here u(x,y,t)=v-w(x,y,t), v is the initial clearance between the film and fixed plate, h is the thickness of the film, w is the transverse displacement of the film, a is the length of the film, b is the wideness of the film, q(x,y,t) is the transverse load per unit area on the film, U is the voltage between the film and fixed plate.

The dynamic equation of the micro film is:

1
ρh2wt2-F2w=qx,y,t,

where 2=2x2+2y2, ρ is the mass density of the film, t is the time, F is the tension in the film.

The electrostatic force per unit area is:

2
qe=FeA=ε0εr2Uv-w2,

where ε0 is permittivity constant of free space, εr is relative dielectric constant of the insulating layer.

Fig. 1Multi-fields coupled dynamic model of a micro film

Multi-fields coupled dynamic model of a micro film

The van der Waals force per unit area between the film and fixed plate is:

3
qr3=A6πv-w3,

where A is the Hamaker constant, A=10-19 J.

The damping force per unit area from air is:

4
qc=cawt,

here, ca is the damping coefficient of the gas, it can be given as:

5
ca=μb3v-w3,

where μ is the gas viscosity, μ=1.86×10-5NSm-2.

Let ε=b/a, ca can be expressed as Taylor series as below:

6
ca=μb3v-w3=εμb2av311-wv3=εμb2av31+3wv+6(wv)2+.

Neglecting high order terms, and substituting Eq. (6) into Eq. (4), yields:

7
qc=cawt=εμb2av31+3wvwt.

Thus, the total force per unit area on the film is:

8
qx,y,t=qe+qr3+qc.

The displacement w of the micro film consists of a static component ws and a dynamic one Δw:

9
w=ws+Δw.

The load q(x,y,t) of static (qs) and dynamic (Δq) components as well:

10
qx,y,t=qs+Δq,
11
qs=U02ε0εr2v-ws2+A6πv-ws3.

From Δq=dqdwΔw, we know:

12
Δq=U02ε0εrv-ws3Δw+A2πv-ws4Δw-εμb2av+3wsv4Δwt-3εμb2aΔwv4Δwt.

Substituting Eq. (10) into Eq. (1), yields the following equations:

13
2ws=-qsF,
14
ρh2wt2-F2Δw=Δq.

3. Static displacement

Let:

15
qs=m=1n=1Cmnsinmπxasinnπyb,

here:

Cmn=4ab0a0bqssinmπxasinnπybdxdy=4qsmnπ21-cosmπ1-cosnπ.

Let:

16
ws=m=1n=1Amnsinmπxasinnπyb.

It should be noted that the Eqs. (15) and (16) only apply to analysis about the harmonic vibration and the weak nonlinear vibration of the film. For investigating propagating waves, transient processes, and chaotic dynamics of the film, a more complicated equation about the displacement distribution should be taken.

Substituting Eqs. (15) and (16) into Eq. (13), yields:

17
Amn=4q01-cosmπ1-cosnπFmnπ4m2a2+n2b2,

thus:

18
ws=m=1n=14qs1-cosmπ1-cosnπFmnπ4m2a2+n2b2sinmπxasinnπyb
=16qsπ4Fm=1,3,5...n=1,3,5...sinmπxasinnπybmnm2a2+n2b2.

The average static displacement of the film is:

19
ws¯=1ab0a0b16qsπ4Fm=1,3,5n=1,3,5sinmπxasinnπybmnm2a2+n2b2dxdy
=m=1,3,5...n=1,3,5...64qsπ6m2n2F(m2a2+n2b2).

From Eq. (10) we know the displacement ws is included in term qs. Here, we take the average displacement to replace the displacement ws, and then Eq. (19) is changed into a nonlinear equation about the average displacement. From the equation, the average displacement of the film can be obtained.

4. Free vibration

Substituting Eq. (12) into Eq. (14), yields:

20
ρh2Δwt2-F2Δw=U02ε0εv-w03Δw+A2πv-w04Δw
-εμb2av+3w0v4Δwt-3εμb2aΔwv4Δwt.

Let Δw=ϕ(x,y)q(t), and then substituting it into Eq. (20), yields:

21
q¨q+εHρhq˙q+εGρhϕq˙=F2ϕ+Eϕρhϕ,

where:

E=U02ε0εv-w03+A2πv-w04, H=μb2av+3w0v4, G=3μb2av4.

Letting Eq. (21) be equal to ω02, yields:

22
q¨+εHρhq˙+εGρhϕ-q˙+ω02q=0,
23
2Δϕx2+2Δϕy2+α2Δϕ=0,

where:

α2=E+ω02ρhF, ϕ-=1ab0a0bϕ(x,y)dxdy.

The mode function of the micro film is considered as:

24
Δϕmn=Asinmπxasinnπyb.

Substituting Eq. (24) into Eq. (23), yields:

25
Aπ2m2a2+n2b2-α2sinmπxasinnπyb=0.

From Eq. (25), we know:

26
π2m2a2+n2b2-α2=0.

From Eq. (26), the natural frequencies of the film can be given:

27
ωmn=1ρhFπ2m2a2+n2b2-U02ε0εv-ws3-A2πv-ws4.

Using Multi-scale method, ones assume:

28
q=q0T0,T1,T2+εq1T0,T1,T2+ε2q2T0,T1,T2.

Substituting Eq. (28) into Eq. (22), and equating coefficients of like powers of ε in both sides, ones obtain:

29a
D02q0+ω02q0=0,
29b
D02q1+ω02q1=-2D0D1q0-A1+Bq0D0q0,
29c
D02q2+ω02q2=-2D0D1q1-D12q0-2D0D2q0-A1+Bq0D0q1
-A1+Bq0D1q0-BD0q0q1,

where:

Dn=Tn, A1=Hρh=μb2av+3w0ρhv4, B=Gϕρh.

The solution of Eq. (29a) is:

30
q0=KT1,T2eiω0T0+K-T1,T2e-iω0T0.

Substituting Eq. (30) into Eq. (29b), yields:

31
D02q1+ω02q1=-2D1iω0K-iω0A1Keiω0T0-Biω0e2iω0T0K2+M,

where M is conjugate complex of other terms.

In order to remove secular item, let:

32
-2D1iω0K-iω0A1K=0.

Thus:

33
D1K=-A1K2.

Substituting Eq. (33) into Eq. (31), yields:

34
q1=iBK2e2iω0T03ω0+M.

Substituting Eqs. (30) and (34) into Eq. (29c), yields:

35
D02q2+ω02q2=-2iω0D2K+A24K+13B2K2K-eiω0T0
-16ABK2e2iω0T0+B2K3e3iω0T0+12ABKK-+M.

In order to remove secular item, let:

36
-2iω0D2K+A124K+13B2K2K-=0.

Thus:

37
D2K=-iA12K8ω0-iB2K2K-6ω0.

Substituting Eq. (37) into Eq. (35), yields:

38
q2=A1BK2e2iω0T018ω02-B2K3e3iω0T08ω02+A1BKK-2ω02+M.

We know:

39
dKdt=D0K+εD1K+ε2D2K.

Substituting Eqs. (33) and (37) into Eq. (39), yields:

40
dKdt=-εA1K2-iε2A12K8ω0-iε2B2K2K-6ω0.

Letting Kt=1/2k(t)eiθt and substituting it into Eq. (40), yields:

41
k'=-εA12k,θ'=-ε2A128ω0-ε2B2k224ω0.

From Eq. (41), it is known:

42
k=k0e-12εA1t,θ=-ε2A12t8ω0+εB2k02eiεA1t24ω0A1+θ0,

where k0 and θ0 is determined by initial conditions.

Substituting Eq. (42) into Kt=1/2k(t)eiθt, and then substituting it into Eqs. (30), (34), and (38), yields:

43a
q0=kcosω0t+θ,
43b
q1=-Bk26ω0sin2ω0t+2θ,
43c
q2=ABk236ω02cos2ω0t+2θ-B2k332ω02cos3ω0t+3θ+ABk24ω02.

Substituting Eq. (43) into Eq. (28), yields:

44
qt=kcosφ-εBk26ω0sin2φ+ε2A1Bk236ω02cos2φ-ε2B2k332ω02cos3φ+A1Bk2ε24ω02,

where φ=ω0t+θ.

5. On linear forced response near natural frequency

Let the disturbing voltage be ΔU=Ecos(ωet), which is subjected between the micro film and the base plate. The exciting electric field force will be caused. Considering weak nonlinear influence of the above exciting electric field force, the dynamic force per unit area is changed from Eq. (12) into following form:

45
Δq=U02ε0εv-w03Δw+A2πv-ws4Δw+εU0ε0εrv-ws2Ecosωet
-εμb2av+3w0v4Δwt-3εμb2av4ΔwΔwt.

Substituting Eq. (45) into Eq. (14), yields:

46
q¨+εA1+Bqq˙-εCcosωet+ω02q=0,

where:

A1=μb2av+3w0ρhv4, B=3μb2aϕ-ρhv4, C=U0ε0ερhv-w02ϕ-E,

ωe is exciting frequency of the voltage.

Using Multi-scale method, from Eq. (46), ones can obtain the magnitude and frequency relationship of the film vibrations as below:

47
12εA1k2+-ωe2+ω022ωe-ε2A128ωe-ε2B2k224ωe-ωe2-ω0228ωe32K2
=εω02C8ωe3-5εC8ωe2+ε2A1C8ωe22.

6. Results and discussions

Above equations are utilized for the free vibration analysis of the micro film. The parameters of the numerical example are shown in Table 1. Comparison between natural frequencies of the film is shown in Tables 2-4 and 5-7. ω0 is the natural frequency of the film without considering the van der Waals force, ω is the natural frequency of the film considering the van der Waals force, δ is the relative error between them. Table 2-4 gives natural frequencies of the film for various clearances between film and fixed plate (here, F= 1 N, U0= 2 V). Table 5-7 gives natural frequencies of the film for various film tensions (here, v= 2 μm, U0= 1 V).

Table 1Parameters of the micro film

a (mm)
b (mm)
h (μm)
v (μm)
ρ (kg∙m-3)
εr
2
1
5
2
2.7×103
1

Table 2Natural frequencies of the film vibration (v= 2 μm)

Order
ω0 (rad/s)
ω (rad/s)
δ (%)
Mode(1,1)
23874.67
18421.61
22.84
Mode(2,2)
57532.60
55505.97
3.52
Mode(3,3)
88713.02
87465.27
1.41
Mode(4,4)
119436.18
118520.63
0.77

Table 3Natural frequencies of the film vibration (v= 2.5 μm)

Order
ω0 (rad/s)
ω (rad/s)
δ (%)
Mode(1,1)
27202.94
25558.12
6.05
Mode(2,2)
58991.52
58264.69
1.23
Mode(3,3)
89682.77
89241.45
0.49
Mode(4,4)
120124.94
119837.40
0.24

Table 4Natural frequencies of the film vibration (v= 3 μm)

Order
ω0 (rad/s)
ω (rad/s)
δ (%)
Mode(1,1)
28512.15
27826.97
2.40
Mode(2,2)
59581.88
59294.99
0.48
Mode(3,3)
90088.85
89917.51
0.19
Mode(4,4)
120445.01
120341.70
0.09

Table 5Natural frequencies of the film vibration (F= 0.2 N)

Order
ω0 (rad/s)
ω (rad/s)
δ (%)
Mode(1,1)
12680.74
12631.14
0.39
Mode(2,2)
26629.17
26605.58
0.09
Mode(3,3)
40286.04
40270.45
0.04
Mode(4,4)
53873.54
53861.88
0.02

Table 6Natural frequencies of the film vibration (F= 0.6 N)

Order
ω0 (rad/s)
ω (rad/s)
δ (%)
Mode(1,1)
22964.35
22955.57
0.038
Mode(2,2)
46607.89
46603.56
0.009
Mode(3,3)
70098.87
70095.99
0.004
Mode(4,4)
93552.29
93550.14
0.002

Table 7Natural frequencies of the film vibration (F= 1 N)

Order
ω0 (rad/s)
ω (rad/s)
δ (%)
Mode(1,1)
29884.63
29880.19
0.015
Mode(2,2)
60288.04
60285.84
0.004
Mode(3,3)
90575.43
90573.97
0.002
Mode(4,4)
120834.08
120832.9
0.001

From Tables 2-4 and 5-7, the following observations are worth noting:

1) Without considering the van der Waals force, the natural frequencies of the film vibration are larger than those of considering the van der Waals force. It is because the van der Waals force system is equivalent to a soft spring system. The van der Waals force can cause decrease of the natural frequencies of the micro film.

2) The deviation between the natural frequencies of the film with and without the van der Waals force decreases with increasing the order number of the mode. For mode(1,1) and v= 2 μm, the relative error between the natural frequencies is 22.84 %. For mode(4,4) and v= 2 μm, the relative error between the natural frequencies is 0.77 %. It shows that influence of the van der Waals force on the natural frequencies decreases with increasing the order number of the mode.

3) As initial clearance between the micro film and the base plate increases, the deviation between the natural frequencies of the film with and without the van der Waals force decreases. For mode(1,1) and v= 2 μm, the relative error between the natural frequencies is 22.84 %. For mode(1,1) and v= 3 μm, the relative error between the natural frequencies is 2.4 %.

4) As the film tension increases, the natural frequency of the film increases significantly. For mode(1,1), when the film tension grows from F= 0.2 N to F= 1 N, the natural frequency of the film increases by 136 %. For mode(4,4), when the film tension grows from F= 0.2 N to F= 1 N, the natural frequency of the film increases by 124 %. It shows that the film tension has obvious effects on the natural frequency of the film for various modes. As the film tension increases, the deviation between the natural frequencies of the film with and without the van der Waals force decreases slightly.

In a word, the van der Waals force should be considered for small clearance between the film and fixed plate, and small film tension.

When the exciting frequency is near the natural frequency of the micro film, changes of the vibration amplitudes with respect to exciting frequency under different modes and parameters are investigated. The results are shown in Figs. 2-5. Fig. 2 gives amplitude frequency characteristics of micro film with and without van der Waals force. Fig. 3 gives amplitude frequency characteristics of micro film for different mechanical parameters. Fig. 4 gives amplitude frequency characteristics of micro film for different fluidic viscosity. Fig. 5 gives amplitude frequency characteristics of micro film for different bias voltage. They show:

1) For mode(1,1), when van der Waals force is considered, the natural frequency of the film vibration is smaller than that without van der Waals force, but its vibrating amplitude is larger than that without van der Waals force (here, vibrating displacement peak increases by 25 % under resonance).

2) For mode(3,3), when van der Waals force is considered, both natural frequency and vibrating amplitude of the film vibration is near to ones without van der Waals force. It shows that the effects of the van der Waals force on both natural frequency and vibrating amplitude of the film become weak for high order modes.

Fig. 2Amplitude frequency characteristics under van der Waals force

Amplitude frequency characteristics under van der Waals force

a) Mode(1,1) with van der Waals force

Amplitude frequency characteristics under van der Waals force

b) Mode(3,3) with van der Waals force

Amplitude frequency characteristics under van der Waals force

c) Mode(1,1) without van der Waals force

Amplitude frequency characteristics under van der Waals force

d) Mode(3,3) without van der Waals force

3) When the micro film length increases, the natural frequency of the film vibration drops and its vibrating amplitude grows. For high order modes, effects of the micro film length on the natural frequency become weak, but effects of the micro film length on its vibrating amplitude are still obvious.

4) For mode(1,1), when nonlinear parameter ε increases, the natural frequency and its vibrating amplitude of the film vibration drops. The resonance peak becomes smaller and smaller with increasing parameter ε. The resonant peak of the first mode will vanish with further increasing nonlinear parameter ε. For mode(3,3), when nonlinear parameter ε increases, the natural frequency and the resonance peak of the film vibration does not change nearly. When the exciting frequency deviates from the natural frequency, the vibrating amplitude of the film with a larger parameter ε is larger than that with a small parameter ε.

5) When the micro film thickness increases, the natural frequency of the film vibration drops and its vibrating amplitude grows significantly. For high order modes, effects of the micro film thickness on the natural frequency and its vibrating amplitude are more obvious than those for low order modes.

6) As the static voltage between the film and fixed plate increases, the natural frequency of the film vibration drops and its vibrating amplitude grows significantly. For low order modes, effects of the static voltage between the film and fixed plate on the natural frequency and its vibrating amplitude are more obvious than those for high order modes.

Fig. 3Amplitude frequency characteristics for different mechanical parameters

Amplitude frequency characteristics for different mechanical parameters

a) Mode(1,1) and a changes

Amplitude frequency characteristics for different mechanical parameters

b) Mode(3,3) and a changes

Amplitude frequency characteristics for different mechanical parameters

c) Mode(1,1) and ε changes

Amplitude frequency characteristics for different mechanical parameters

d) Mode(3,3) and ε changes

Amplitude frequency characteristics for different mechanical parameters

e) Mode(1,1) and h changes

Amplitude frequency characteristics for different mechanical parameters

f) Mode(3,3) and h changes

Fig. 4Amplitude frequency characteristics for different fluidic viscosity

Amplitude frequency characteristics for different fluidic viscosity

a) Mode(1,1) and μ changes

Amplitude frequency characteristics for different fluidic viscosity

b) Mode(3,3) and μ changes

Fig. 5Amplitude frequency characteristics for different bias voltage

Amplitude frequency characteristics for different bias voltage

a) Mode(1,1) and U0 changes

Amplitude frequency characteristics for different bias voltage

b) Mode(3,3) and U0 changes

7) As the gas viscosity between the micro film and the fixed plate decreases, the vibrating amplitudes of the micro film decrease significantly. It has not nearly effects on the natural frequency of the film vibration for the high and low order modes.

In a word, the van der Waals force has obvious effects on both natural frequency and vibrating amplitude of the film for the mode(1,1); the static voltage between the film and fixed plate has obvious effects on the natural frequency for the mode(1,1) and the vibrating amplitude of the film for all of the modes; the gas viscosity only has obvious effects on the vibrating amplitude of the film for all of the modes; among mechanical parameters, the micro film thickness has obvious effects on both natural frequency and vibrating amplitude of the film for all of the modes, and the micro film length has obvious effects on both natural frequency and vibrating amplitude of the film for the mode(1,1).

7. Conclusions

In this paper, a multi-fields coupled dynamic equation of the micro film is presented. The multi-fields coupled free vibration of the micro film is investigated. The multi-fields coupled forced response of the micro film near natural frequency is analyzed. The multi-fields coupled natural frequencies of the micro film and their changes along with system parameters are discussed. Amplitude frequency characteristics of the micro film and its changes along with system parameters are given. A number of useful results are given. The results show:

1) The van der Waals force should be considered for small clearance between the film and fixed plate, small film tension, and low order modes.

2) The bias voltage between the film and fixed plate has obvious effects on the natural frequency for the mode(1,1) and the vibrating amplitude of the film for all of the modes.

References

  • Thielicke E., Obermeier E. Microactuators and their technologies. Mechatronics, Vol. 10, 2000, p. 431-455.
  • Bao Z. P., Subrata M. Electrostatic BEM for MEMS with thin conducting plates and shells. Engineering Analysis with Boundary Elements, Vol. 28, 2004, p. 1427-1435.
  • El-hami M., Glynne-Jones P., White N. M., Hill M., Beeby S., James E., Brown A. D.,Ross J. N. Design and fabrication of a new vibration-based electromechanical power generator. Sensors and Actuators, A: Physical, Vol. 92, Issue 1-3, 2001, p. 335-342.
  • Wang P. H., Dai X. H., Fang D. M. Design, fabrication and performance of a new vibration-based electromagnetic micro power generator. Microelectronics Journal, Vol. 38, Issue 12, 2007, p. 1175-1180.
  • Kuehne I., Frey A., Marinkovic D. Power MEMS-A capacitive vibration-to-electrical energy converter with built-in voltage. Sensors and Actuators, A: Physical, Vol. 142, Issue 1, 2008, p. 263-269.
  • Rafael N. G., Anna M. B., Alfons D. AC transfer function of electrostatic capacitive sensors based on the 1-D equivalent model: application to silicon microphones. Journal of Microelectromechanical Systems, Vol. 12, Issue 6, 2003, p. 972-978.
  • Lu X. L., Clarke M. L., Li D. A sum frequency generation vibrational study of the interference effect in poly (n-butyl methacrylate) thin films sandwiched between silica and water. Journal of Physical Chemistry C, Vol. 115, Issue 28, p. 13759-13767.
  • Akai D., Ozaki K., Numata Y. Vibration analysis and transmission characteristics of piezoelectric micromachined ultrasonic transducers using epitaxial Pb(Zr,Ti)O3 thin films on γ-Al2O3/Si substrate. Japanese Journal of Applied Physics, Vol. 51, Issue 11, 2011.
  • Li J. B., Xu Q. F. An electrostatic film vibration voltage transducer. Applied Mechanics and Materials. 2013, p. 321-324, 1443-1448, 2013.
  • Gualdino A., Chu V., Conde J. P. Study of the out-of-plane vibrational modes in thin-film amorphous silicon micromechanical disk resonators. Journal of Applied Physics, Vol. 113, Issue 17, 2013.
  • Sorokin V. M., Yashchenko A. K. Vibrations of a frameless film membrane stabilized by the Ampère force in zero gravity. Journal of Applied Mechanics and Technical Physics, Vol. 54, Issue 6, 2013, p. 885-893.

About this article

Received
14 July 2014
Accepted
18 September 2014
Published
15 February 2015
SUBJECTS
Flow induced structural vibrations
Keywords
multi-fields coupled
micro film
MEMS
free vibration
amplitude frequency characteristics
Acknowledgements

This project is supported by Key Basic Research Foundation in Hebei Province of China (13961701D).

Conflict of interest

Authors declare that there is no conflict of interests regarding the publication of this article.