Relating Elastic Modulus To Indentation Response Using Atomic Force Microscopy
M. R. VanLandingham1, S. H. McKnight2, G. R. Palmese1, R. F. Eduljee1, J. W. Gillespie, Jr.1, and R. L. McCullough1
1Center For Composite Materials, University
of Delaware, Newark, DE 19716
2Army Research Laboratory -- Materials
Directorate, Wilmington, DE 19898
Atomic force microscopy (AFM) has become an extremely useful tool in materials science applications. In addition to imaging surfaces with nanometer resolution, the AFM can also be used to determine surface roughness [1], probe local changes in friction [2], measure surface forces [1, 3], and assess changes in local elasticity over a sample surface. Recently much attention has been focused on the interpretation of AFM force curves for characterizing tip-sample interactions [1, 3-7]. Through appropriate analysis of force curves, a wealth of information can be obtained regarding the mechanical, chemical, and adhesive properties of the surface. To date, most studies using force mode AFM [3-6] have concentrated on the attractive and transition regions of the force curve. Very little effort has been dedicated to the examination of the contact portion (or so-called repulsive region) of the force curve. This region contains valuable information regarding the nanoscale mechanical response of the sample [7]. In this paper, a technique is developed that (1) relates the elastic modulus of the sample to the sample response measured using AFM force curves; and (2) illuminates the importance of the relative stiffnesses of the cantilever probe and the sample to the determination of elastic modulus.
During the force mode of AFM, the probe tip is first lowered into contact with the sample, then indented into the surface, and finally lifted off of the sample surface. Concurrently, a laser beam is reflected off the top of the probe and onto a segmented photodiode, thus producing a measurement of the probe tip deflection. A plot of this tip deflection signal as a function of the vertical displacement of the piezo actuator is produced. Such a plot, termed a "force curve," is shown in Fig. 1. As the piezo displacement in the z-direction increases, the tip deflection voltage (Vtd) remains constant until the probe tip makes contact with the surface (A to B in Fig. 1). Just before tip-sample contact is made, the probe tip can be pulled down to the surface by attractive forces, causing a small decrease in Vtd (B to C in Fig. 1). Decreasing the piezo height further causes the cantilever to deflect in the opposite direction, resulting in an increased Vtd signal (C to D in Fig. 1). During unloading (D to E in Fig. 1), the piezo retracts, reducing the cantilever deflection until the tip separates from the sample surface. Often, however, the tip adheres to the surface, causing a further decrease in Vtd until the tip jumps out of contact (E to F in Fig. 1).

Fig. 1 – AFM force curve of a polyurethane sample (E = 0.05 GPa) measured using a 1N/m probe.
To extract the nanoscale mechanical response of the material from the contact portion of the resulting force curve, the tip-sample interaction is modeled as two springs in series, as shown in Fig. 2. After contact is made between the probe tip and the sample surface, piezo displacement results in both probe tip deflection and sample indentation, the amounts of which depend on the relative stiffnesses of the sample and the cantilever probe. The relationship between the displacement of the piezoelectric actuator (Dzp), the displacement due to tip deflection (Dzt), and the indentation displacement (Dzi) is simply
|
|
|
The tip displacement can be directly related to the applied force, P, by
|
|
|
where kc is the spring constant of the cantilever probe. The accuracy of the calculated forces is generally on the order of 1 to 10 nN, depending on cantilever stiffness, because the noise on the tip displacement detector limits depth resolution to 0.1 nm [1].

Fig. 2 – Schematic and spring model of tip-sample interaction in contact portion of force curve.
For a sample which is infinitely stiff with respect to the probe, no indentation will occur and Dzp = Dzt. From simple beam theory, the angle change (Dqt) and the tip displacement (Dzt) are related by
|
|
|
where Lc is the length of the cantilever probe. This angle change is directly related to the change in Vtd by system conversion and amplification factors which can be lumped into a constant, Cq, such that
|
|
|
Combining Equations 3 and 4, the slope of the force curve, S, is shown to have an upper limit, S*, given by
|
|
|
For a sample that deforms due to the force applied by the probe, S will be reduced from S* by a reduction factor, f (i.e., f = S / S*), such that larger values of f correspond to stiffer samples. Also, the amount of indentation, Dzi, at different points along the force curve can be calculated using Equation 1. The applied force, P, can also be calculated using Equation 2 if kc is known. Using contact mechanics for the indentation of an elastic half space [8, 9], the indentation displacement, Dzi, can be related to the applied force, P, through the relation
|
|
|
where x is a constant which depends on probe tip geometry, E and n are the sample elastic modulus and Poisson's ratio, respectively, and m is a power law exponent characteristic of the indentation behavior. For behavior characteristic of a flat cylindrical punch, m = 1, while for behavior characteristic of spheres and paraboloids of revolution, m = 1.5. Normally, m is determined from a curve fit of P as a function of Dzi for a material of known E and n.
For the experimental investigations, a Digital Instruments D3000 scanning probe microscope was used to produce force curves from the interaction of two separate silicon cantilever probes and a number of polymer samples. The spring constants of the two probes were estimated to be 1 and 60 N/m respectively. The polymer samples consisted of polyurethanes, visible-light-cured acrylated epoxies, thermally-cured epoxies, and polyetheretherketone (PEEK). The elastic moduli of these polymer samples were determined independently using dynamic mechanical analysis and range from 0.02-3.0 GPa. In order to provide a calibration reference, titanium (E = 110 GPa [10]) was used as an "infinitely stiff" material for the determination of Cq from Equation 4. For each sample, a minimum of 10 individual force curves was obtained. Each force curve was then analyzed according to the methodology outlined previously.
The results obtained show promising correlation between the generated force curve data and the sample moduli. The slope reduction factor, f, was observed to increase as the stiffness of the sample increased. For each cantilever probe, however, the f values leveled off above a certain modulus value (approximately 0.5 GPa for the 60 N/m probe and 0.1 GPa for the 1 N/m probe). The response of polymers with moduli greater than 0.5 and 0.1 GPa, respectively for the 60 N/m and 1 N/m probes, approaches that of the "infinitely stiff" titanium sample. Therefore, above these modulus values, samples of differing moduli become indistinguishable.
The use of the reduction factor (f), however, does not take into account the actual indentation process. Using Equations 1 through 4, the loads applied to the sample surface and the corresponding indentation displacements can be determined. The ratio of the maximum indentation displacement relative to a maximum probe tip displacement of 100 nm versus sample modulus is plotted in Fig. 3. This ratio is shown to be quite sensitive to changes in modulus up to 0.5 and 0.1 GPa for the 60 N/m and 1 N/m probes, respectively. When the amount of probe tip displacement exceeds the indentation displacement, Dzi / Dzt becomes insensitive to modulus changes. In other words, the force applied by the probe to the sample surface is insufficient to produce adequate deformation of the stiffer samples. Because the 60 N/m probe applies higher forces to the sample than does the 1 N/m probe for a given tip displacement, it is sensitive up to higher modulus values and thus can be used to evaluate higher modulus samples. Also shown in Fig. 3 are curves generated using the spring model shown schematically in Fig. 2 and outlined in Equations 1, 2, and 6. Input values typical of experimentally measured responses (m = 1.2 and x/(1-n2) = 80) were used to calculate values of Dzi corresponding to Dzt = 100 nm for probe spring constants of 1, 60, and 500 N/m. The model results correspond well to the experimental data and indicate that extremely stiff probes (kc = 500 N/m) will be required to probe modulus changes in polymers commonly used in engineering applications. Both the experimental results and the results of the spring model indicated that cantilever selection should be based on the expected range of moduli.

Fig. 3 – Ratio of maximum indentation displacement, Dzimax, to a maximum tip displacement, Dztmax, of 100 nm plotted as a function of sample modulus, E.
Future work will address the use of stiffer cantilevers to probe higher modulus (E > 1.0 GPa) polymeric materials. While this technique shows promise, several factors must be addressed to extend the capabilities of the AFM as a probe of local property variations. One such area to be studied will be the stochastic aspects of the indentation process. Material heterogeneity and other stochastic effects can be significant when probing nanoscale properties and must be addressed. A second concern is the effect of time-dependent behavior on the measured response. However, changing several independent force mode parameters should allow for a more detailed characterization of the relative amounts of viscoelastic and plastic behavior involved. Finally, techniques currently under development should allow for more accurate measurements of the spring constants of cantilever microbeams so that the applied forces can be determined more precisely. With a firm understanding of these effects, a more quantitative evaluation of the indentation mechanics can be made. Calibration curves can then be generated which relate the nanoscale mechanical response of the sample to values of modulus for various probe-polymer combinations. These calibration curves will allow for the study of modulus changes in heterogeneous regions typically found in polymer blends, adhesive systems, coatings, and in the vicinity of fibers in composite materials.
References
1. S. M. Hues, R. J. Colton, E. Meyer, and H.-J. Guntherodt, MRS Bull. XVIII (1993) 41.
2. J. A. Ruan and B. Bhushan, Trans. ASME J. Tribol. 116 (1994) 378.
3. N. A. Burnham, R. J. Colton, and H. M. Pollock, Nanotechnol. 4 (1993) 64.
4. J. P. Aime, Z. Elkaakour, C. Odin, T. Bouhacina, D. Michel, J. Curely, and A. Dautant, J. Appl. Phys. 76 (1994) 754.
5. H. A. Mizes, K.-G. Loh, R. J. D. Miller, S. K. Ahuja, and E. F. Grabowski, Appl. Phys. Lett. 59 (1991) 2901.
6. S. P. Jarvis, T. P. Weihs, A. Oral, and J. B. Pethica, in Proceedings of Materials Research Society Symposium, 308, San Francisco, April 1993, edited by P. H. Townsend, T. P. Weihs, J. E. Sanchez, Jr., and P. Borgesen (Materials Research Society, Pittsburgh, PA) p. 127.
7. S. M. Hues, C. F. Draper, and R. J. Colton, J. Vac. Sci. Technol. B 12 (1994) 2211.
8. I. N. Sneddon, Int. J. Engng. Sci. 3 (1965) 47.
9. W. C. Oliver and G. M. Pharr, J. Mater. Res. 7 (1992) 1564.
10. J. F. Shackelford "Introduction to Materials Science for Engineers" (MacMillen, New York, 1985).
Copyright 1997 Chapman and Hall, all rights reserved.