Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Phys Rev Lett. Author manuscript; available in PMC 2013 November 25.
Published in final edited form as:
Published online 2012 June 5.
PMCID: PMC3839317

Stress-Induced Variations in the Stiffness of Micro- and Nanocantilever Beams


The effect of surface stress on the stiffness of cantilever beams remains an outstanding problem in the physical sciences. While numerous experimental studies report significant stiffness change due to surface stress, theoretical predictions are unable to rigorously and quantitatively reconcile these observations. In this Letter, we present the first controlled measurements of stress-induced change in cantilever stiffness with commensurate theoretical quantification. Simultaneous measurements are also performed on equivalent clamped-clamped beams. All experimental results are quantitatively and accurately predicted using elasticity theory. We also present conclusive experimental evidence for invalidity of the longstanding and unphysical axial force model, which has been widely applied to interpret measurements using cantilever beams. Our findings will be of value in the development of micro- and nanoscale resonant mechanical sensors.

Physical or chemical adsorption onto the surface of micro- and nanoscale cantilevers is the basis for a rapidly growing field of biological sensing in the mechanical domain. Despite the remarkable achievements of the field [1,2], it is surprising that understanding the mechanism by which biomolecules affect the physics of these mechanical structures remains elusive [317]. The adsorption process can result in a multiplicity of concurrent effects including: mass loading of the device [4,5,18], enhancement of surface elasticity [12,13,17,19], increase in damping [20], and the imposition of surface stress [611,15]. Numerous studies over the past 35 years have claimed that surface stress can cause significant variations in device stiffness [314,17]. In contrast, other works ignore such effects, claiming that such stress loads do not affect the stiffness of cantilever devices, and argue for alternative interpretations. It is also striking that, within this extensive body of work [1,2], controlled, quantitative measurements of the effects of surface stress on cantilever stiffness with commensurate theoretical interpretation have yet to be reported.

The original theoretical model of Lagowski et al. [15] suggested that a net axial force is induced along the beam axis upon application of stress—a load similar to the case of doubly clamped beams derived within the framework of Euler-Bernoulli beam theory [21]. Lagowski et al. found that the axial force model was in quantitative and qualitative agreement with measurements. However, the axial force model as applied to cantilever beams has subsequently been shown to be in violation of Newton’s 3rd law; i.e., it does not satisfy the fundamental physical principle of force equilibrium [16,22,23], contrary to doubly clamped beams. In short, application of surface stress to a cantilever beam always induces stress of opposite sign within the beam material, resulting in zero net axial force along the beam. This phenomenon is contingent on the beam length greatly exceeding its width and thickness, a fundamental assumption of Euler-Bernoulli beam theory; see Refs. [16,22,23]. Hence, Euler-Bernoulli beam theory leads to the conclusion that cantilever beams should be insensitive to surface stress.

Though it was proven unphysical, numerous experimental studies have provided experimental data, seemingly consistent with the axial force model in micro- and nanoscale cantilever beams [69,11,15]. However, it is important to note that surface stress change in all such experimental studies is achieved using processes that may cause unspecified changes in the mechanical properties of the resonator [15,24]. It is thus not entirely evident whether these previous measurements are due solely to stress changes at the cantilever beam surface, or effects of an unspecified nature. This situation contrasts strongly to clamped-clamped beams whose stress-induced stiffness change can be systematically observed and theoretically predicted [21,25].

In this Letter we extend previous theoretical work, propose a new theoretical model and explain small but nonnegligible changes in the resonant frequencies of cantilever beams due to application of stress.We also present the first systematic measurements of both cantilever and doubly clamped beams using a robust and highly repeatable methodology, which shows remarkable agreement with theory. This resolves the above-mentioned long-standing debate and experimentally establishes the invalidity of the axial force model for cantilever beams.

In contrast to the case of doubly clamped beams, the application of surface stress to a long and thin cantilever beam does not generate a net axial force because the longitudinal displacement is not constrained. However, a cantilever generates a net in-plane stress in the immediate vicinity of its supporting clamp [22], which can affect its resonance frequency. Expressions for the relative frequency shifts Δf / fR of doubly clamped and cantilever beams due to the application of surface stress, and subsequent generation of net in-plane stress within the device, are shown in the second column of Table I; termed the stress effect. We emphasize that all results in Table I are derived for thin beams.

Formulas determining the relative frequency shift of thin doubly clamped beams (top row) and cantilever beams (bottom row) when a normalized load, σ¯[equivalent](1ν)σsT/(Eh), is applied. Both stress (left column) ...

Application of a load to any elastic body changes its geometry due to elastic deformation. This geometric effect is typically ignored in the classical theory of linear elasticity. In the present context, application of surface stress induces a change in the beam length, width, thickness, and density, which alter the resonant frequency of both clamped-clamped and cantilever beams. Expressions for the relative frequency shift due to this geometric effect are shown in the third column of Table I. To account for the true geometry of the devices, finite element analysis is performed below. The complete effect is given by the sum of the stress and geometric effects.

The formulas in Table I indicate that the resonant frequencies of doubly clamped beams are more sensitive to surface stress changes than cantilever beams (considering typical devices dimensions, L > 10h). For clampedclamped beams, the stress effect listed in Table I dominates, whereas for cantilever beams the situation is more complex: the stress effect is dominant for thin structures (h [double less-than sign] b), with geometric effects prevailing for thicker devices. Importantly, if the geometric effect dominates, the relative frequency shift Δf / fR does not change with length. Alternatively, if the stress effect prevails then Δf / fR strongly depends on length.

To examine the predictions of the above theoretical model, we perform systematic measurements on nanoscale cantilever and doubly clamped beams of identical geometries. Our measurements are achieved by fabricating multilayered mechanical bars out of piezoelectric material, utilizing recent advances in nanoelectromechanical systems (NEMS) fabrication techniques [26]. Thus, stress changes in our devices are internally produced piezoelectrically, which allows for reproducible and controllable variation. It should be noted that in our experiment the stress is applied to a layer of finite thickness, rather than to the surface of the devices. To account for any differences that this may induce, we conduct rigorous numerical finite element simulations of both devices that includes their full structure and load (see [27]). Doubly clamped beams with identical dimensions provide a benchmark for the measurements on cantilevers.

Cantilever and doubly clamped nanoscale beams are fabricated on the same chip from a 320 nm 4-layer stack (comprised of: 20 nm aluminum nitride (AlN), 100 nm molybdenum (Mo), 100 nm AlN, and 100 nm Mo), using a process described elsewhere [26]. A standard wafer curvature measurement yields built-in stresses in all four layers close to zero. SEM micrographs of the devices are shown in Fig. 1. All beams possess an identical width b = 900 nm and total thickness htot = 320 nm, whereas their lengths L are 6, 8, and 10 µm. The absence of curvature in cantilever beams indicates the absence of a stress gradient in the structural layers. Note that the only difference between cantilevers and doubly clamped beams is that the boundary condition at one end is changed— material and geometric properties are otherwise identical. An electric field is produced in the active piezoelectric layer by applying a dc voltage between the top and bottom molybdenum layers. This induces an axial stress along the beam axis via its inherent piezoelectric properties. Ensuring the active piezoelectric layer is offset from the neutral axis of the beam enables the fundamental flexural mode of the device to be actuated through application of an additional ac voltage. Measurements of the resulting normal deflection are performed in vacuum using optical interferometry (see [27]). Resonance measurements are shown in Figs. 1(a) and 1(b) for doubly clamped beams and cantilevers, respectively. We monitor the resonant frequency while changing dc voltage by means of a computer-controlled phase locked loop. Controllable variations in the stiffness of both types of devices due to stress changes are observed.

FIG. 1
(color online). Resonant response of piezoelectric beams. a, SEM micrograph of the doubly clamped beams used for the experiments. On top of the micrograph, we show resonant responses of each of the beams, yielding resonant frequencies of 38.3 MHz (length ...

We initially study the effect of an applied stress on the fundamental resonant frequencies of doubly clamped beams. Formulas connecting resonance frequency shifts to surface stress are presented in Table I (see Supplemental Material [27] for corresponding expressions for piezoelectric loads). Measurement results of doubly clamped beams are presented in Fig. 2, which clearly exhibit a linear variation in resonant frequency with voltage (stress) in the piezoelectric layer, as predicted theoretically. No variation of quality factor with bias voltage is observed. This provides strong validation for the robustness of the experimental methodology and the ability to tune the stiffness of doubly clamped beams. Use of the resonant frequency shift expressions enables the axial force (or the corresponding surface stress, as commonly reported) to be calibrated against the applied voltage (see [27]), which we determine to be 2:3 N/m per V . In agreement with theory (see Table I), the relative frequency shift shown in Fig. 2(a) depends strongly on the beam length, scaling asL2 (fR is the original resonant frequency), whereas the absolute shift in resonant frequency Δf is found to be independent of the beam length [see Fig. 2(b)].

FIG. 2
(color online). Frequency shift results for doubly clamped beams. a, Relative frequency shift Δf / fR for the three doubly clamped beams in Fig. 1(a), showing the dependence of Δf / fR on length (as predicted by theory). b, Absolute frequency ...

In Fig. 3 we present complementary results for the resonant frequency shift in cantilever beams. These nanomechanical devices have identical dimensions and material properties to the doubly clamped beams in Fig. 2. Linear tuning of the resonant frequency with the applied piezoelectric layer voltage is also observed. However, these frequency shifts are 2 orders of magnitude smaller than those exhibited by doubly clamped beams, even though the piezoelectric stress loads are identical—this constitutes direct experimental evidence of the invalidity of the unphysical axial force model for cantilever beams; see Eq. (S9). Not only the magnitude but also the scaling behavior of the resonant frequency shift Δf with the beam length is inconsistent with the axial force model (see [27]).

FIG. 3
(color online). Frequency shift results for cantilever beams. a, Relative frequency shift Δf / fR for the three cantilever beams in Fig. 1(b), showing independence of length, which is consistent with the theoretical model. Both experimental (lighter ...

Notably, measurements of Δf / fR for cantilever beams display independence on the cantilever length and are identical in sign to those obtained for doubly clamped beams; see Figs. 2(a) and 3(a). This observation is completely consistent with the geometric effect listed in Table I, which is expected to dominate since these devices are relatively thick (b=h ~ 3). Importantly, the formulas given in Table I are derived under the assumption of a thin beam (h [double less-than sign] b), where the stress load is applied to an infinitesimal layer at the surface. To account for the true device geometry and properties, as noted above, we calculate the combined contribution from geometric and stress-induced effects using full 3D finite element simulations (see [27]). The agreement between predictions from these simulations (shown as scatter plots in Figs. 2 and and3)3) and the experimentally measured data is within 15%. The apparent slight asymmetry between positive and negative voltages is within the experimental error (see [27]). These results constitute the first repeatable measurements and theoretical quantification of stress-induced changes in the frequency shifts (stiffness) of cantilever beams. They also provide compelling experimental evidence for the invalidity of the axial force model, which has been widely applied to the interpretation of measurements made using cantilever beams.

It remains to understand why several previously reported cantilever measurements display results that are inconsistent with the present findings [69,11,14,15]. These studies report stress-induced stiffness changes that appear to be in agreement with the unphysical axial force model that erroneously predicts effects much larger in magnitude than the results reported here. These previous, now anomalous, measurements were typically performed using surface adsorption or modification processes. While these processes are known to induce surface stress, additional unspecified and uncontrolled effects are also possible. There is certainly the distinct possibility of over-layer formation, since material is adsorbed to [12,13,17,19] or removed from the beam surface [15]. Some studies have shown excellent agreement between classical composite beam theory and such measurements [12,13,17,19], while other measurements remain unexplained [69,14,15]. Our findings strongly suggest that previous reports of stressinduced changes in cantilever beams originate from other uncontrolled surface phenomena.

Miniaturization of NEMS technologies is key to enhanced sensitivity and ultrafast measurements. Our study establishes that as cantilever thickness is reduced, a transduction mechanism arises that is different to the one observed in the present devices. Specifically, we predict that cantilever devices made of ultrathin materials (such as graphene) allow for gigantic tunability in their resonance properties. Such materials would exhibit an additional inverse-squared thickness dependence, enabling very strong modification of their stiffness. This could ultimately permit the development of UHF bandpass cantilever filters, with broad control of their frequency range for applications in biosensing, telecommunications and medical technologies.

Supplementary Material

Supplementary Info


We would like to thank X. L. Feng, J. Xiang, and M. J. Lachut for useful suggestions and discussions. We would also like to thank E. Defaÿ, G. Le Rhun, and C. Marcoux from CEA-LETI for providing us with Aluminum Nitride material. L. G.V. acknowledges financial support from the European Commission (PIOF-GA-2008-220682) and Professor A. Boisen. J. E. S. acknowledges support from the Australian Research Council grants scheme.


1. Arlett JL, Myers EB, Roukes ML. Nature Nanotech. 2011;6:203. [PMC free article] [PubMed]
2. Boisen A, Dohn S, Sylvest Keller S, Schmid S, Tenje M. Rep. Prog. Phys. 2011;74 036101.
3. Lavrik NV, Sepaniak MJ, Datskos PG. Rev. Sci. Instrum. 2004;75:2229.
4. Yang YT, Callegari C, Feng XL, Ekinci KL, Roukes ML. Nano Lett. 2006;6:583. [PubMed]
5. Naik AK, Hanay MS, Hiebert WK, Feng XL, Roukes ML. Nature Nanotech. 2009;4:445. [PMC free article] [PubMed]
6. Chen GY, Thundat T, Wachter EA, Warmack RJ. J. Appl. Phys. 1995;77:3618.
7. Cherian S, Thundat T. Appl. Phys. Lett. 2002;80:2219.
8. Lee JH, Kim TS, Yoon KH. Appl. Phys. Lett. 2004;84:3187.
9. McFarland AW, Poggi MA, Doyle MJ, Bottomley LA, Colton JS. Appl. Phys. Lett. 2005;87 053505.
10. Dorignac J, Kalinowski A, Erramilli S, Mohanty P. Phys. Rev. Lett. 2006;96:186105. [PubMed]
11. Hwang KS, Eom K, Lee JH, Chun DW, Cha BH, Yoon DS, Kim TS, Park JH. Appl. Phys. Lett. 2006;89:173905.
12. Tamayo J, Ramos D, Mertens J, Calleja M. Appl. Phys. Lett. 2006;89:224104.
13. Ramos D, Tamayo J, Mertens J, Calleja M, Zaballos A. J. Appl. Phys. 2006;100:106105.
14. Gavan KB, Westra HJR, van der Drift EWJM, Venstra WJ, van der Zant HSJ. Appl. Phys. Lett. 2009;94:233108.
15. Lagowski J, Gatos HC, Sproles ES. Appl. Phys. Lett. 1975;26:493.
16. Lu P, Lee HP, Lu C, O’Shea SJ. Phys. Rev. B. 2005;72 085405.
17. Ramos D, Tamayo J, Mertens J, Calleja M, Villanueva LG, Zaballos A. Nanotechnology. 2008;19 035503. [PubMed]
18. Li M, Tang HX, Roukes ML. Nature Nanotech. 2007;2:114. [PubMed]
19. Gupta AK, Nair PR, Akin D, Ladisch MR, Broyles S, Alam MA, Bashir R. Proc. Natl. Acad. Sci. U.S.A. 2006;103 13 362. [PubMed]
20. Yang YT, Callegari C, Feng XL, Roukes ML. Nano Lett. 2011;11:1753. [PMC free article] [PubMed]
21. Timoshenko S, Goodier JN. Theory of Elasticity. New York: McGraw-Hill; 1969.
22. Lachut MJ, Sader JE. Phys. Rev. Lett. 2007;99:206102. [PubMed]
23. Gurtin ME, Markenscoff X, Thurston RN. Appl. Phys. Lett. 1976;29:529.
24. Pini V, Tamayo J, Gil-Santos E, Ramos D, Kosaka P, Tong H-D, van Rijn C, Calleja M. ACS Nano. 2011;5:4269. [PubMed]
25. Villanueva LG, Karabalin RB, Matheny MH, Kenig E, Cross MC, Roukes ML. Nano Lett. 2011;11:5054. [PubMed]
26. Karabalin RB, Matheny MH, Feng XL, Defaÿ E, Le Rhun G, Marcoux C, Hentz S, Andreucci P, Roukes ML. Appl. Phys. Lett. 2009;95:103111.
27. See Supplemental Material at for details.