TECHNICAL FIELD

[0001] The scope includes a method for measuring the bulk mechanical properties of a material.

CROSS REFERENCE

[0001] This application claims the benefit of U.S. Provisional Application No. 62/084,032, filed December 8, 2015, the subject matter of which is expressly incorporated by reference.

BACKGROUND ART

[0002] Scanning probe microscopy (SPM) is a branch of microscopy that forms images of surfaces of a sample or specimen using a physical probe that scans the sample. An image of the surface is obtained by mechanically moving the probe in a raster (line by line) scan of the sample and recording the interaction between the probe and the sample surface as a function of the position of the probe relative to the surface of the sample. An atomic force microscope (AFM) is a type of SPM that employs a cantilever with an atomically sharp tip at the end as the probe. Referring now to Fig. 1, a prior art AFM is illustrated. The AFM of Fig. 1 comprises a controller 12, a microscale cantilever 14 with a sharp tip (probe) at one end that is used to scan the surface of the sample 18, a laser 20, a photo-diode 22 (while the term photo-diode is used herein, an array of photo-diodes is typically used), and an image processor 24. When the probe tip is brought into proximity of the sample surface, forces between the tip and the sample cause a deflection of the cantilever, and this deflection may be measured using the laser and photo-diode. A laser spot (illustrated by line 28) is reflected (illustrated by line 30) by the top of the cantilever onto the photo-diode, and the resulting electrical signal (illustrated by line 32) is sent to the controller 12. The deflection of the cantilever causes an amplitude change in the reflected laser spot (and in the resulting electrical signal), such that the amplitude at any specific time corresponds to the surface contour of the sample at a specific location on the sample. As this scanning is performed over the entire surface of the sample, the resulting amplitude data corresponds to the surface contour of the entire sample. The location and amplitude data is provided to the image processor (illustrated by line 34), such that the image processor is able to create an image of the surface of the sample.

[0003] There are two primary modes in which the AFM is operated: static and dynamic modes. In static modes, the cantilever is brought into contact with the sample, and its deflection directly represents the topography of the probed sample. In dynamic modes, the cantilever is driven at user- specified frequencies and amplitudes by a piezo stack at its base and the amplitude and phase of the resulting oscillation is directly measured. Topography and, in principle, mechanical properties of the sample can be determined using these methods. The mechanical properties such as the (reduced) Young's modulus are of particular modern interest, and efforts to develop quantitative techniques measure them using the AFM has comprised much of the research in atomic force microscopy. All current methods assume the sample surface is quasi-static, and most drive the cantilever at a single frequency and directly measure an effective stiffness of the sample that must be further resolved. The effective stiffness is a function of the contact radius between the cantilever tip and the sample. A separate experiment is typically conducted to characterize the tip geometry, which is typically assumed to be hemi-spherical, using a calibrated sample. Using this information, the contact area in the experiment of interest can be estimated, and the reduced Young's modulus can be extracted from the effective stiffness. Ignoring sample dynamics and assuming the sample is quasi-static precludes measuring the full set of viscoelastic parameters of a material and is also what forces a separate experiment to be conducted to characterize the tip geometry. It also leads to the erroneous conclusion that the sample stiffness is dependent on the driving frequency.

SUMMARY

[0004] Previous work has been done to obtain the effective mechanical properties via harmonic amplitude spectrum atomic force microscopy (HAS -AFM) using generic peripheral equipment in the patent application 62/084,032. HAS-AFM drives the cantilever at a plurality of frequencies to measure effective damping and mass in addition to effective stiffness. These effective parameters also involve the contact area, but additionally include corrections to the usual quasi-static equations employed in most quantitative AFM techniques. In the proposed method herein, we leverage the HAS-AFM apparatus to obtain the cantilever response at a plurality of frequencies and fit a novel model of dynamic contact mechanics to determine (reduced) Young's modulus and bulk damping parameters. [0005] Harmonic amplitude spectrum atomic force microscopy (HAS-AFM) is employed using a field programmable gate array (FPGA) peripheral device to send and receive signals from the signal access module (SAM) of a generic atomic force microscope (AFM) system. The FPGA delivers signals at a plurality of frequencies to the SAM to excite the AFM cantilever. The cantilever interacts with a material sample and generates signals at a plurality of frequencies and their harmonics. The response of the cantilever at the plurality of frequencies and their harmonics is received by and stored on the FPGA. Logic, detailed herein, is provided to determine the bulk mechanical properties of the sample from the recorded cantilever response. The logic uniquely utilizes sample dynamics in its measurement.

BRIEF DESCRIPTION OF DRAWINGS

[0006] Fig.1 illustrates a prior art atomic force microscope.

[0007] Fig.2 illustrates the peripheral equipment and software the field programmable gate array (FPGA) delivery vector.

[0008] Fig.3 illustrates the implementation of the harmonic amplitude spectrum atomic force microscopy (HAS-AFM) technology by interfacing with the AFM through the signal access module (SAM) with the FPGA delivery vector.

[0009] Fig.4 illustrates an example of a pressure distribution at the upper boundary of a sample.

[0010] Fig.5 illustrates the method of images employed to obtain the model of dynamic contact mechanics.

[0011] Fig.6 depicts a chart illustrating the predicted fractional change in the magnitude of the sample Green function as a function of factor multiplied by cantilever fundamental frequency.

DISCLOSURE OF INVENTION

[0012] The method to measure bulk mechanical properties of probed samples is enumerated as follows

A. Record cantilever response when in contact with a sample is obtained at a plurality of frequencies via HAS-AFM. This step is comprised of 1. Connecting a field programmable gate array (FPGA) to a signal access module (SAM) for an atomic force microscope

2. Driving the cantilever in free space at a plurality of frequencies and recording the response using logic on the FPGA

3. Characterizing the cantilever based on the recorded response

4. Engaging the sample surface with the cantilever tip

5. Driving at the plurality of frequencies while engaged and recording the cantilever response using logic on the FPGA

B. Utilize a novel model of dynamic contact mechanics to solve for bulk mechanical parameters of a sample using the cantilever response by

1. Choosing the parameters in the model that make the predicted motion of the model optimally fit the recorded response

2. Determining the contact area between the cantilever tip and sample from this fit

3. Extracting the reduced Young's modulus and damping parameters from the fit

[0013] Below, we detail the logic to determine material mechanical properties of the probed samples using the phenomenological parameters. To obtain the model of dynamic contact mechanics and extract a Green function to be employed in HAS-AFM, we must consider a novel model of dynamic contact mechanics that accounts, explicitly, for sample dynamics. We wish to solve the equation for wave propagation in a material of finite thickness and infinite lateral extent that is subjected to a local pressure distribution at its upper boundary, z = a, and held pinned at its lower boundary, z = 0:

where u ^m is the mth component of the displacement field at a point in the material, and where z is the coordinate running along the finite axis and parameterizes the position in the

lateral plane. The unit vector is normal to the surface, and thus dependent on

[0014] Define the operator

For an isotropic material, we have and

with all other components vanishing.

[0015] Under these symmetries, there are three equations of interest:

[0016] The pressure distribution is assumed to be cylindrically symmetric and monotonically decreasing with radial distance from the cylindrical axis of symmetry.

This forces u to have the form

where r is the horizontal distance from the center of the tip to the point of interest in the sample and r is the unit vector associated with the radial direction in cylindrical coordinates. An illustration of such a pressure distribution and the sample is given in Fig.4.

[0017] Casting Eqs.(4) in cylindrical coordinates reduces the problem from three coupled equations to two:

[0018] The sample is subject to the boundary conditions

[0019] It is useful to express the displacement as a 2-dimensional Fourier transform and Hankel transform:

where ] _a is the -th Bessel function of the first kind. For ί = r, a = 1; for ί = z, a = 0.

[0020] While the Fourier transforms diagonalize the t- and z- derivatives, they do not enforce the boundary conditions. Explicitly diagonalizing the matrix operators subject to the boundary conditions is an intractable problem; however, the method of images may be employed to enforce the boundary conditions while using the Fourier transforms.

[0021] It is easier to implement the method and visualize the image source distribution by leaving the location of the real pressure distribution at z = a and setting the free boundary condition at z = a + b. Image sources are inserted above the z = a + b plane to enforce the free condition on the upper boundary; image sources are inserted below the z = 0 plane to enforce the pinned condition on the lower boundary. Initially inserting the image sources to act in concert with the real source to enforce the boundary conditions generates a successive cascade of sources that must be inserted to counter the effects of each previous insertion. A visual depiction of the beginning of this cascade is given in Fig. (5). A source below the lower boundary must take the opposite sign of the source it counters; a source above the upper boundary must take the same sign as the source it counters. After the calculation is finished, b can be taken to 0.

[0022] Neglecting friction (P _r = 0), setting solving for performing the integrals

over m and k, and setting z = a yields

where Li _n is the n-th polylogarithm, and the has been written in terms of complex Lame parameters. Since the z-pressure was restricted to a Dirac delta at the center of the coordinate system, the resulting solution, Eq.(8), furnishes the continuous Green function for the sample in Fourier space, and thus we give the solution the variable name g .

[0023] This distributed surface Green function has a nontrivial action on the pressure distribution under integration, demonstrating that the spatial distribution of the force applied to the surface matters, not solely the total force. In effective approaches to AFM dynamics, such as that used in HAS-AFM, we care only about the total force acting at the sample surface, so the effective sample Green function must be some functional of the distributed surface Green function. To connect them, we must average the distributed surface Green function over an appropriate area with some particular weighting function. To determine the weighting function, we seek inspiration from known, solved problems in contact mechanics.

[0024] If an infinitely hard cylindrical probe impinges upon a sample surface, the displacement of the surface from equilibrium in the region of contact is uniform. That is,

[0025] Of course, beyond the contact region, the displacement is not uniform.

[0026] Since the Green function must be symmetric under averaging

over r with the weighting function in a function that is independent of r'

and thus has trivial action on the pressure distribution under integration for r' < R. If our pressure distribution has a finite, effective contact radius, it is natural, and indeed necessary, to choose R such that the integral subtends exactly this contact region. Since the resulting relationship between the average surface displacement and the pressure distribution is simply the convolution of an averaged Green function with the total force, the resulting averaged Green function thus corresponds to the sample Green function in the effective models employed in HAS-AFM.

[0027] Performing the averaging yields

Here, R is the (effective) contact radius between the tip and the sample surface. The first term is the usual contact mechanics term; the second term arises as a finite-thickness modification; the remaining terms arise from sample dynamics.

[0028] Eq.(10) is to be used in conjunction with the following equation for the cantilever's response while in contact with the sample

where is the cantilever's amplitude at a given driving frequency in free space, A _t is the

cantilever's amplitude at the given driving frequency while engaged with the sample, σ is simply a phenomenological parameter to fit to data is the cantilever Green function that is calibrated

before scanning, and is the "renormalized" sample Green function, in which K is another

parameter to be fit to data. The real and imaginary parts of the parameters λ and μ, which describe the elasticity and damping of a sample, can be fit to data.

[0029] The sample thickness, a, can usually be determined by alternative means (e.g. optical measurements, or specifying the sample dimensions during fabrication); the average density, p, can be determined by measuring the resonance frequency if the dimensions are known or by other means. Determining these parameters are easier than estimating the tip geometry and can be much more accurately measured. Once these parameters are known, the reduced Young's modulus and damping parameters can easily be determined point-by -point while a HAS-AFM scan is being conducted by fitting the functional form to the response of the cantilever at a plurality of frequencies. [0030] Traditional techniques that attempt to extract bulk properties of the sample consider only the first term in Eq.(10). This both incorrectly predicts the frequency dependence of that sample's response, and forces users of the technique to employ an additional technique to determine the contact radius between the tip and sample before the bulk parameters can be ascertained. This usually involves characterizing the tip geometry, which introduces a number of compounding uncertainties resulting in inaccurate bulk measurements and technical problems for the user. Fig. 6 depicts the predicted fractional change in the magnitude of the sample Green function as a function of frequency for a typical rat tail fibroblast. Interpreting the inverse of the Green function as the effective local stiffness of the cell is consistent with measurements using traditional methods.