We present a method for characterizing ultrathin films using sensitivity-enhanced atomic force acoustic microscopy, where a concentrated-mass cantilever having a flat tip was used as a sensitive oscillator. Evaluation was aimed at 6-nm-thick and 10-nm-thick diamond-like carbon (DLC) films deposited, using different methods, on a hard disk for the effective Young's modulus defined as E/(1 - ν2), where E is the Young's modulus, and ν is the Poisson's ratio. The resonant frequency of the cantilever was affected not only by the film's elasticity but also by the substrate even at an indentation depth of about 0.6 nm. The substrate effect was removed by employing a theoretical formula on the indentation of a layered half-space, together with a hard disk without DLC coating. The moduli of the 6-nm-thick and 10-nm-thick DLC films were 392 and 345 GPa, respectively. The error analysis showed the standard deviation less than 5% in the moduli.
Keywords:Atomic force acoustic microscopy; Thin film; Elastic modulus; Diamond-like carbon; Concentrated-mass cantilever
The protective coating for hard disks, namely a diamond-like carbon (DLC) film, is now targeted for thickness less than 3 nm because of the reduced spacing between the magnetic layer and the read/write head . The mechanical properties become very important for reliability of the devices. The chemical structure of DLC significantly depends on the deposition process and influences the mechanical properties such as elasticity and hardness. Especially the Young's modulus E drastically varies with a content of sp3-bonds, which form three-dimensional interlinks in the amorphous network of carbons (E ≈ 100–800 GPa) [2-4]. Therefore, the modulus is useful to identify the chemical structure of films.
Various approaches for the determination of the elastic properties of thin films have been previously used, including nanoindentation , laser spectroscopic methods , and removed substrate methods . However, it is still a challenging problem to evaluate ultrathin films like DLC films with thickness less than 10 nm.
Atomic force acoustic microscopy (AFAM)  is a promising method, which belongs to a family of dynamic techniques of atomic force microscope (AFM) such as micro-deformation microscopy  and ultrasonic atomic force microscopy . AFAM measures the resonant frequency f of an AFM cantilever whose sensor tip is in contact with a sample oscillated by a piezoelectric device. If an appropriate order of the vibration mode is selected, f varies with the contact stiffness k*, namely the interactive force gradient between a tip and a sample. The effective Young's modulus of a sample, defined as (Es: the Young's modulus, νs: the Poisson's ratio), is evaluated using contact mechanics relating k* to .
Characterization of a 50-nm-thick Ni film deposited on a Si substrate was demonstrated in AFAM, where f was observed without the substrate effects . In regard to DLC thin films, only relative evaluation was performed . These studies required a blunt tip with a radius of about 200 nm and a stiff cantilever of spring constant kc ≈ 50 N/m to realize reproducible measurements. However, the requirement reduced the spatial resolution and the sensitivity in detection of the contact force.
When attempting to analyze difficult samples like a DLC film with thickness less than 10 nm, higher performance of AFAM is required on the detection of k* and the spatial resolution. We previously proposed a concentrated-mass (CM) cantilever as a way of enhancing the sensitivity in k*-detection without trade-offs . A CM cantilever assures the maximum sensitivity for any sample material. Also, a flat tip with ductile-metal coating, keeps a stable contact area of a radius less than 5 nm and drastically simplifies the relation between k* and .
The method we previously developed, termed sensitivity-enhanced AFAM , is extended in this letter to the determination of the elastic modulus of ultrathin films. The demonstration was carried out for DLC films with thickness of 6 and 10 nm, deposited on a hard disk. A curve relating f to was determined from multiple measurements on reference samples. The uncertainty was discussed by error analysis. In the evaluation of the DLC-coated samples, the substrate effect was taken into account by using an analytical model for indentation of a layered half-space .
Experimental Procedure and Theory
CM Cantilever and Apparatus
The experimental procedure is described elsewhere in detail [11,12]. We will briefly explain it here. The main body of a CM cantilever was a rectangular cantilever made of single-crystalline silicon (μMasch Co. Ltd., kc = 0.65 N/m, fundamental resonant frequency 40.9 kHz). The silicon tip had an apex radius of about 10 nm and was coated with a 25-nm-thick Pt/Ti film. The coated tip was plastically deformed on a flat diamond surface under a contact load of 2 μN to give it a flat-ended shape. This plastic deformation also induced a work-hardening of the coating, which would prolong the lifetime of the coated tip . For the concentrated mass, a tungsten (W) particle of 35 × 33 × 20 μm in size was micro-machined from a W sheet of 20 μm thick by focused ion beam (FIB). The particle's mass was about 445 ng, which corresponds to a mass ratio of 10.9, namely the ratio of the particle's mass to the silicon-cantilever's mass. The particle was attached adhesively to the free end of the cantilever by micromanipulation. Figure 1 shows a scanning electron micrograph of the CM cantilever. The main difference from the previous works [11,12] was in the use of the micro-machined particle instead of a deoxidized random particle for the concentrated mass. Another difference was in the process that a flat tip was formed from a virgin tip, not from a tip wasted after several tens of scans for imaging.
Figure 1. The concentrated-mass (CM) cantilever. The CM was micro-machined from a 20-μm-thick tungsten film using a focused ion beam (FIB).
An atomic force microscope (SII Co. Ltd., SPI3700-SPA270) was used in so-called contact mode for observing the contact resonance spectra. The amplitude of cantilever vibration was acquired with a lock-in amplifier through a heterodyne down-converter. A piezoelectric device placed beneath a sample was used for the oscillation. The time-averaged cantilever deflection signal, which corresponds to the contact force Fe, was maintained through a built-in feedback circuit, where the electronic circuit is not subjected to sinusoidal signals at ultrasonic frequencies. The resonant frequency was measured at five to ten different locations on a sample to confirm reproducibility. All experiments were carried out at a temperature of 20–25°C and relative humidity of about 40–50%.
Reference and DLC Samples
The reference samples and the elastic moduli are listed in Table 1. We employed a sapphire (0001) wafer in addition to silicon wafers and a diamond (100) used in the previous work . These values were deduced from the crystal moduli determined by ultrasonic velocity techniques for bulk samples (see appendix for sapphire).
Table 1. Elastic moduli of reference samples
DLC films of 6 nm thick and 10 nm thick were deposited on a substrate by sputtering a carbon target in Ar gas and by plasma-assisted chemical vapor deposition (CVD), respectively. The film thickness was estimated based on the deposition time. The substrate was a hard disk, which consisted of metallic multi-layers for magnetic record and a glass substrate, namely (50-nm-thick Co-Cr-alloy layer)/(70-nm-thick Ti-alloy layer)/(0.6-mm-thick glass substrate). Also, the substrate without DLC coating was tested for the elastic modulus.
Theory for Evaluation of Thin Films
The resonant frequency (f) of a CM cantilever increases with the contact stiffness (k*) in accordance with the spring-mass model, namely k*/kc = (f/f0)2, where f0 is the fundamental resonant frequency in the absence of a sample. A flat tip maintains a constant contact area independent of the adhesion force and the contact force. This also ensures constant k*. The theoretical formula k* = 2aE* for a flat-ended punch  is applicable, where a is the radius of the contact area. E* is the effective Young's modulus of the contact region, defined as . is the effective Young's modulus of a tip. These equations give the formula relating f to :
where is a factor proportional to the contact radius. Both and A can be determined from the f measurements for reference samples.
Analytical models on indentation of a layered half-space for a circular punch proved the validity of the following empirical formula :
where and are the effective Young's moduli of a film and a substrate, respectively. The coefficient γ is a function of a/t, where t is the film thickness. The numerical result on a relation of γ and a/t was graphically shown in reference . Note that the symbol a in reference  is defined as the square root of the contact area, which differs from the definition of a (the radius of the contact area) in this letter, and then γ multiplied by π1/2 equals the symbol α in reference . Examples of the numerical result are indicated with circles in Figure 2. The numerical data can be well fitted by the following formula.
where c0 = 0.4684, c2 = 0.009968, c3 = 1.004, n0 = 0.4910, n1 = 1.736, and n2 = 6.607 are the coefficients determined by a nonlinear least-square fit.
Errors in for a sample come from uncertainties in the predetermination of A and and in the f measurement for the sample, which are represented by the standard deviations σA, σEt, and σf, respectively. The standard deviation σEs of can be estimated by the error propagation on Eq. 1: , where and . In term of the covariance between A and , the correlation coefficient is set to -1, the validity of which was confirmed in the fitting of Eq. 1. Assuming negligible errors in γ and a/t, the standard deviation σEfilm of is estimated by the error propagation on Eq. 2: , where , , and σEsub is the standard deviation of .
Results and Discussion
Effective Young's Modulus of a Flat Tip and the Contact Radius
The CM cantilever in free space measured f0 = 9.917 kHz for the fundamental resonant frequency. Figure 3 shows spectra for the reference, Si (100) wafer. The resonant frequency seems to become independent of the contact force (Fe) when increasing Fe. This reflects the constant contact area observed in the case of the flat tip.
Figure 3. The spectra of the CM cantilever vibration in contact with silicon (100) when increasing the contact force (Fe).
To measure f, we set Fe to be a value as small as possible, at which the resonant peak was clear and settled in frequency. The value depended on the sample material. The resonance frequencies for Si (100), Si (111), Al2O3 (0001), and diamond (100) were f = 199.3 ± 1.3 kHz, 218.6 ± 1.9 kHz, 254.5 ± 1.1 kHz, and 281.0 ± 1.1 kHz, where Fe is set to 300, 400, 500, and 700 nN, respectively. The errors show the 95% confidence regions (±2σ). The excellent reproducibility was attained in the measurements for 5–10 different positions on each reference surface. Figure 4 shows examples of spectra for the reference samples.
Figure 4. The spectra for reference samples [Si (100), Si (111), Al2O3 (0001), and diamond (100)] and hard disk samples [6 nm-DLC (Sputter), 10 nm-DLC (CVD) and Co-Cr alloy (hard disk without DLC)].
Fitting Eq. 1 to the relationship between the resonant frequencies measured for reference and the effective Young's moduli listed in Table 1, we determined and , which are hard to measure or estimate directly. Figure 5 shows the least-squares fit obtained for the reference samples, which yielded A = 0.2496 ± 0.0061 (± 2σ) m/kg and = 184.6 ± 8.8 (± 2σ) GPa. The errors for A and correlate, and the error's sign is taken opposite to each other.
Figure 5. The theoretical curve, which relates the resonant frequency to the effective Young's modulus of a sample, fitted to the experimental data (ο) for the reference. The error bars and the broken curves indicate the 95% confidence regions, namely twice the standard deviations.
Use of the values of A, kc, and f0 produced a reasonable contact radius a = 1.7 nm. Also, the value of is comparable to the averaged value for bulk platinum (196 GPa) and bulk titanium (129 GPa), but close to the value for platinum differently from the previous work ( = 152.3 GPa) . This would be on account of the contact area smaller than that in the previous work (a = 4.4 nm) . In the contact deformation of the present tip, the contribution of surface layer (Pt) would dominate rather than the insert layer (Ti).
The square of the correlation coefficient (r2 = 0.9987) of the fit confirms the validity of the theory on a CM cantilever with a flat tip. The error bar for each data point and the broken curves in Figure 5 indicate the 95% confidence regions.
Evaluation of DLC Thin Films
The samples coated with the 6-nm-thick DLC film (Sputter) and the 10-nm-thick DLC film (CVD) measured f = 240.4 ± 1.6 (± 2σ) kHz and f = 239.6 ± 0.5 (± 2σ) kHz, where Fe is set to 600 and 800 nN, respectively. The DLC coating shifted the resonant frequency to higher than that of the sample without DLC coating [f = 229.8 ± 1.6 (± 2σ) kHz (Fe = 500 nN)]. Also, the values of f for the two DLC films were alike despite the different thickness. This does not mean that the resonance is free from the substrate effects.
The effective Young's modulus of a sample was determined from the curve in Figure 5 to be = 310.5 ± 11.4 GPa, 305.2 ± 3.5 GPa, and 247.8 ± 8.2 GPa for the hard disks with 6-nm-thick DLC (Sputter), 10-nm-thick DLC (CVD), and without DLC coating, respectively. The errors are in the 95% confidence regions. The last one corresponds to . The value of was similar to the modulus of Co-Cr alloys (230–280 GPa) [15,16].
Substituting the values of and into Eq. 2, we obtained the effective Young's modulus of a film ( ), where γ was calculated using a = 1.7 nm and t = 6 nm or 10 nm. The moduli were = 391.8 ± 34.7 (± 2σ) GPa and 345.1 ± 8.5 (± 2σ) GPa for the 6-nm DLC (Sputter) and the 10-nm DLC (CVD), respectively. The presence of substrate effects was clear in that the values of for the 6-nm-film-coated and 10-nm-film-coated samples were 20 and 10% less than the corresponding values of , respectively. The values of were within the range of values reported for several DLC films, from 100 to 800 GPa [2-4]. Also, a good precision of 2σ < 10% was attained.
An error in a/t, which was neglected in the present evaluation, also causes uncertainty of the results. A postulated error of 20% in a/t results in a relatively small error of about 5 and 2.5% in for the DLC films of 6 nm thick (a/t = 0.283) and 10 nm thick (a/t = 0.17), respectively. The resulting error increases with a/t. Therefore, the contact radius (a) should be minimized.
The indentation depth δs, namely the total displacement δ (= Fe/k*) minus the tip deformation, can be estimated by taking account of the contribution of a sample, , in the contact stiffness. The estimate was = 0.57 nm and 0.77 nm for the 6-nm-DLC and 10-nm-DLC samples, respectively. These indentation depths are 10% or less of the film thickness. The substrate effect should be carefully considered even when AFAM is applied. The present method provides the AFAM method of determining the elastic modulus for ultrathin films, eliminating the influence of a substrate. The sensitivity-enhanced AFAM proved to be sensitive enough for the determination of the ultrathin film elasticity and to have the excellent repeatability and reliability.
The 6-nm-thick and 10-nm-thick diamond-like carbon (DLC) films on a hard disk were evaluated for the effective Young's modulus using sensitivity-enhanced atomic force acoustic microscopy. The modulus was determined from the measurements of the contact resonant frequency with the aid of a theory on indentation of a layered half-space. The moduli of the 6-nm DLC and 10-nm DLC were 391.8 ± 34.7 GPa and 345.1 ± 8.5 GPa, respectively, which reflected difference in the method of film deposition. The errors, the 95% confidence regions (± 2σ), show that this method gives a precise estimate of the effective Young's modulus.
Calculations of the modulus are alike for cubic crystals like diamond and silicon  and for trigonal crystals like sapphire (Al2O3). A trigonal crystal, however, has more constants, c11, c12, c13, c14, c33, and c44, than a cubic crystal. They relate stresses σi to strains εi (i = 1 - 6) as follows:
where the subscripts adopt an abbreviated notation (e.g., σ1 = σ11, and σ4 = σ23) . All components are referred to Cartesian coordinates xi (i = 1 - 3), where the x1 and x3 axes are taken along the a1-axis  and the c-axis , respectively. Calculation of the effective Young's modulus for the C-plane (0001) requires the Young's modulus E in the direction of the c-axis. Eliminating ε1, ε3, and ε4 in Eq. 4 under the condition of uniaxial stressing σi = 0 (i ≠ 3), we can obtain the following formula from E = σ3/ε3.
While the Poisson's ratio ν depends on the direction in which lateral strain is measured, we let ε1 (= ε2) represent lateral strain for simplicity. This simplifies the calculation of ν:
We extend our appreciation to Mr. Fukunori Izumida of Iwate Industrial Research Institute and Mr. Shunsuke Goto of LIGHTOM, Inc., for their assistance in micro-machining the concentrated mass, to Dr. Jun Ariake of Akita Institute of Advanced Technology for providing the DLC-coated hard disk samples, and to Mr. Yasuyuki Hasebe of MAHK Co., Ltd., for his generosity in proffering the sapphire wafer. This work was supported by the Japan Society for the Promotion of Science under a Grant-in-Aid for Science Research (S) 18106003, (A) 20246028, and (B) 20360049.
IEEE Trans Mag. 2000, 36:67. Publisher Full Text
Surf Coat Technol. 1998, 98:1097. Publisher Full Text
Appl Phys Lett. 1999, 75:1893. Publisher Full Text
Diam Rel Mater. 2001, 10:2069. Publisher Full Text
Appl Phys Lett. 1996, 68:779. Publisher Full Text
Rev Sci Instrum. 1996, 67:3281. Publisher Full Text
Appl Phys Lett. 1993, 62:829. Publisher Full Text
Jpn J Appl Phys. 1996, 35:3787. Publisher Full Text
Nanotechnology. 2005, 16:703. Publisher Full Text
Thin Solid Films. 2001, 392:75. Publisher Full Text
JSME Int J A. 2002, 45:567. Publisher Full Text
Nanotechnology. 2005, 16:542. Publisher Full Text
Int J Solids Struct. 1987, 23:1657. Publisher Full Text
Int J Eng Sci. 1965, 3:47. Publisher Full Text
Biomed Mater Eng. 1994, 4:231. PubMed Abstract