Nanoindentation, also called instrumented indentation testing,[1] is a variety of indentation hardness tests applied to small volumes. Indentation is perhaps the most commonly applied means of testing the mechanical properties of materials. The nanoindentation technique was developed in the mid-1970s to measure the hardness of small volumes of material.[2]
In a traditional indentation test (macro or micro indentation), a hard tip whose mechanical properties are known (frequently made of a very hard material like diamond) is pressed into a sample whose properties are unknown. The load placed on the indenter tip is increased as the tip penetrates further into the specimen and soon reaches a user-defined value. At this point, the load may be held constant for a period or removed. The area of the residual indentation in the sample is measured and the hardness,
H
Pmax
Ar
H= | Pmax |
Ar |
.
For most techniques, the projected area may be measured directly using light microscopy. As can be seen from this equation, a given load will make a smaller indent in a "hard" material than a "soft" one.
This technique is limited due to large and varied tip shapes, with indenter rigs which do not have very good spatial resolution (the location of the area to be indented is very hard to specify accurately). Comparison across experiments, typically done in different laboratories, is difficult and often meaningless. Nanoindentation improves on these macro- and micro-indentation tests by indenting on the nanoscale with a very precise tip shape, high spatial resolutions to place the indents, and by providing real-time load-displacement (into the surface) data while the indentation is in progress.
In nanoindentation small loads and tip sizes are used, so the indentation area may only be a few square micrometres or even nanometres. This presents problems in determining the hardness, as the contact area is not easily found. Atomic force microscopy or scanning electron microscopy techniques may be utilized to image the indentation, but can be quite cumbersome. Instead, an indenter with a geometry known to high precision (usually a Berkovich tip, which has a three-sided pyramid geometry) is employed. During the course of the instrumented indentation process, a record of the depth of penetration is made, and then the area of the indent is determined using the known geometry of the indentation tip. While indenting, various parameters such as load and depth of penetration can be measured. A record of these values can be plotted on a graph to create a load-displacement curve (such as the one shown in Figure 1). These curves can be used to extract mechanical properties of the material.[3]
The slope of the curve,
dP/dh
S
Er
E | ||||
|
\sqrt{\pi | |
Ap(hc)
hc
\beta
Ap(hc)
Ap(hc)=C0h
2+C | |
1h |
1+C | |
2h |
1/2 | |
c |
+C3h
1/4 | |
c |
+\ldots+C8h
1/128 | |
c |
C0
Er
Es
1/Er=(1-\nu
2)/E | |
i+(1-\nu |
2)/E | |
s. |
Here, the subscript
i
\nu
Ei
\nui
\nus
There are two different types of hardness that can be obtained from a nano indenter: one is as in traditional macroindentation tests where one attains a single hardness value per experiment; the other is based on the hardness as the material is being indented resulting in hardness as a function of depth.
H= | Pmax |
Ar |
.
The hardness is given by the equation above, relating the maximum load to the indentation area. The area can be measured after the indentation by in-situ atomic force microscopy, or by 'after-the event' optical (or electron) microscopy. An example indentation image, from which the area may be determined, is shown at right.
Some nanoindenters use an area function based on the geometry of the tip, compensating for elastic load during the test. Use of this area function provides a method of gaining real-time nanohardness values from a load-displacement graph. However, there is some controversy over the use of area functions to estimate the residual areas versus direct measurement. An area function
Ap(hc)
h
The strain-rate sensitivity of the flow stress
m
m= | \partialln{\sigma |
\sigma=\sigma(\varepsilon |
)
\varepsilon |
m
dln{H}=mdln{
\varepsilonp |
The subscripts
p
V*
V*=9kBT
| |||
T
m
V*\propto(Hm)-1
The construction of a depth-sensing indentation system is made possible by the inclusion of very sensitive displacement and load sensing systems. Load transducers must be capable of measuring forces in the micronewton range and displacement sensors are very frequently capable of sub-nanometre resolution. Environmental isolation is crucial to the operation of the instrument. Vibrations transmitted to the device, fluctuations in atmospheric temperature and pressure, and thermal fluctuations of the components during the course of an experiment can cause significant errors.
Dynamic nanoindentation or continuous stiffness measurement (CSM, also offered commercially as CMX, dynamics...), introduced in 1989,[4] is a significant improvement over the quasi-static mode described above. It consists into overlapping a very small, fast (> 40 Hz) oscillation onto the main loading signal and evaluate the magnitude of the resulting partial unloadings by a lock-in amplifier, so as to quasi-continuously determine the contact stiffness. This allows for the continuous evaluation of the hardness and Young's modulus of the material over the depth of the indentation, which is of great advantage with coatings and graded materials. The CSM method is also pivotal for the experimental determination of the local creep and strain-rate dependent mechanical properties of materials, as well as the local damping of visco-elastic materials. The harmonic amplitude of the oscillations is usually chosen around 2 nm (RMS), which is a trade-off value avoiding an underestimation of the stiffness due to the "dynamic unloading error"[5] or the "plasticity error"[6] during measurements on materials with unusually high elastic-to-plastic ratio (E/H > 150), such as soft metals.
The ability to conduct nanoindentation studies with nanometre depth, and sub-nanonewton force resolution is also possible using a standard AFM setup. The AFM allows for nanomechanical studies to be conducted alongside topographic analyses, without the use of dedicated instruments. Load-displacement curves can be collected similarly for a variety of materials – provided that they are softer than the AFM tip – and mechanical properties can be directly calculated from these curves.[7] Conversely, some commercial nanoindentation systems offer the possibility to use a piezo-driven stage to image the topography of residual indents with the nanoindenter tip.
Utilizing optical fiber Fabry-Perôt interferometry, nanoindentation studies can be performed with unparalleled precision, achieving micro-mechanical characterization of soft biomaterials.[8] Optical interferometry, allows for nanomechanical studies of biomaterials alongside topographic analyses without the need for dedicated instruments. This technology is particularly adept at micro-mechanical analysis of soft and living materials. Load-displacement curves can be gathered for a diverse range of materials, and their mechanical properties can be directly inferred from these curves. Additionally, some advanced systems offer the capability to integrate optical imaging with micro-mechanical characterization, enabling a comprehensive understanding of the relationship between structure and stiffness in biomaterials.
The indentation curves have often at least thousands of data points. The hardness and elastic modulus can quickly be calculated by using a programming language or a spreadsheet. Instrumented indentation testing machines come with the software specifically designed to analyze the indentation data from their own machine. The Indentation Grapher (Dureza) software is able to import text data from several commercial machines or custom made equipment.[9] Spreadsheet programs such as MS-Excel or OpenOffice Calculate do not have the ability to fit to the non-linear power law equation from indentation data. A linear fit can be done by offset
(-h(P0))
XY
The Martens hardness,
HM
hmax
Pmax
HM= | Pmax |
As |
.
The displacement is used to calculate the contact surface area,
As
As=24.5
2 | |
h | |
max |
The indentation hardness,
HIT
H | ||||
|
.
Here, the hardness is related to the projected contact area
Ap
As the indent size decreases the error caused by tip rounding increases. The tip wear can be accounted for within the software by using a simple polynomial function. As the indenter tip wears the
C1
C0
C1
Ap=C0
2 | |
h | |
max |
+C1hmax.
Calculating the elastic modulus with software involves using software filtering techniques to separate the critical unloading data from the rest of the load-displacement data. The start and end points are usually found by using user defined percentages. This user input increases the variability because of possible human error. It would be best if the entire calculation process was automatically done for more consistent results. A good nanoindentation machine prints out the load unload curve data with labels to each of the segments such as loading, top hold, unload, bottom hold, and reloading. If multiple cycles are used then each one should be labeled. However mores nanoindenters only give the raw data for the load-unload curves. An automatic software technique finds the sharp change from the top hold time to the beginning of the unloading. This can be found by doing a linear fit to the top hold time data. The unload data starts when the load is 1.5 times standard deviation less than the hold time load. The minimum data point is the end of the unloading data. The computer calculates the elastic modulus with this data according to the Oliver—Pharr (nonlinear). The Doerner-Nix method is less complicated to program because it is a linear curve fit of the selected minimum to maximum data. However, it is limited because the calculated elastic modulus will decrease as more data points are used along the unloading curve. The Oliver–Pharr nonlinear curve fit method to the unloading curve data where
h
hf
k
m
k
hf
m
dP/dh
P=k\left(h-hf\right)m.
An image of the indent can also be measured using software. The atomic force microscope (AFM) scans the indent. First the lowest point of the indentation is found. Make an array of lines around the using linear lines from indent center along the indent surface. Where the section line is more than several standard deviations (> 3
\sigma
For nanoindentation experiments performed with a conical indenter on a thin film deposited on a substrate or on a multilayer sample, the NIMS Matlab toolbox[10] is useful for load-displacement curves analysis and calculations of Young's modulus and hardness of the coating.[10] In the case of pop-in, the PopIn Matlab toolbox[11] is a solution to analyze statistically pop-in distribution and to extract critical load or critical indentation depth, just before pop-in.[11] Finally, for indentation maps obtained following the grid indentation technique, the TriDiMap Matlab toolbox[12] offers the possibility to plot 2D or 3D maps and to analyze statistically mechanical properties distribution of each constituent, in case of a heterogeneous material by doing deconvolution of probability density function.[12]
Molecular dynamics (MD) has been a very powerful technique to investigate the nanoindentation at atomic scale. For instance, Alexey et al [13] employed MD to simulate the nanoindentation process of a titanium crystal, dependence of deformation of the crystalline structure on the type of the indenter is observed, which is very hard to harvest in experiment. Tao et al [14] performed MD simulations of nanoindentation on Cu/Ni nanotwinned multilayers films using a spherical indenter and investigated the effects of hetero-twin interface and twin thickness on hardness. Recently, a review paper by Carlos et al [15] is published upon the atomistic studies of nanoindentation. This review covers different nanoindentation mechanisms and effects of surface orientation, crystallography (fcc, bcc, hcp, etc), surface and bulk damage on plasticity. All of the MD-obtained results are very difficult to be achieved in experiment due to the resolution limitation of structural characterization techniques. Among various MD simulation software, such as GROMACS, Xenoview, Amber, etc., LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator), which is developed by Sandia National Laboratories, is the most widely used for simulation. An interaction potential and an input file including information of atom ID, coordinates, charges, ensemble, time step, etc are fed to the simulator, and then running could be executed. After specified running timesteps, information such as energy, atomic trajectories, and structural information (such as coordination number) could be output for further analysis, which makes it possible to investigate the nanoindentation mechanism at atomic-scale. Another interesting Matlab toolbox called STABiX has been developed to quantify slip transmission at grain boundaries by analyzing indentation experiments in bicrystal.[16]
Nanoindentation is a robust technique for determination of mechanical properties. By combining the application of low loads, measuring the resulting displacement, and determining the contact area between the tip of the indenter and the sample a wide range of mechanical properties are able to be measured.[17] [18] The application that drove the innovation of the technique is testing thin film properties for which conventional testing are not feasible. Conventional mechanical testing such as tensile testing or dynamic mechanical analysis (DMA) can only return the average property without any indication of variability across the sample. However, nanoindentation can be used for determination of local properties of homogeneous as well as heterogeneous materials.[19] [20] The reduction in sample size requirements has allowed the technique to become broadly applied to products where the manufactured state does not present enough material for microhardness testing. Applications in this area include medical implants, consumer goods, and packaging.[21] Alternative uses of the technique are used to test MEMs devices by utilizing the low-loads and small scale displacements the nanoindenter is capable of.[22]
Conventional nanoindentation methods for calculation of Modulus of elasticity (based on the unloading curve) are limited to linear, isotropic materials.
Problems associated with the "pile-up" or "sink-in" of the material on the edges of the indent during the indentation process remain a problem that is still under investigation. It is possible to measure the pile-up contact area using computerized image analysis of atomic force microscope (AFM) images of the indentations.[23] This process also depends on the linear isotropic elastic recovery for the indent reconstruction.
Nanoindentation of soft material has intrinsic challenges due to adhesion, surface detection and tip dependency of results. There is an ongoing research to overcome such problems.[24]
Two critical issues need to be considered when attempting nanoindentation measurements on soft materials: stiffness and viscoelasticity.
The first is the requirement that in any force-displacement measurement platform the stiffness of the machine (
kmachine
ksample
kmachine
kmachine
The stiffness of a sample is given by
ksample
a
Esample
a
E
kmachine
kmachine
ksample
Esample
To measure
Esample
kmachine
a
Esample
kmachine
a
The second issue concerning soft samples is their viscoelasticity. Methods to handle viscoelasticity include the following.
In the classical treatment of viscoelasticity, the load-displacement (P-h) response measured from the sample is fitted to predictions from an assumed constitutive model (e.g. the Maxwell model) of the material comprising spring and dashpot elements.[26] Such an approach can be very time consuming, and cannot in general prove the assumed constitutive law in an unambiguous manner.
Dynamic indentation with an oscillatory load can be performed, and the viscoelastic behavior of the sample is presented in terms of the resultant storage and loss moduli, often as variations over the load frequency.[27] However, the storage and loss moduli obtained this way are not intrinsic material constants, but depend on the oscillation frequency and the indenter probe geometry.
A rate-jump method can be used to return an intrinsic elastic modulus of the sample that is independent of the test conditions.[28] In this method, a constitutive law comprising any network of (in general) non-linear dashpots and linear elastic springs is assumed to hold within a very short time window about the time instant tc at which a sudden step change in the loading rate is applied on the sample. Since the dashpots are described by relations of the form
\epsilon |
\epsilon |
\sigma
\sigma
\sigma |
\sigma |
\epsilon |
\epsilon
\sigma
\epsilon |
∆
\epsilon |
\sigma |
The last equation indicates that the fields ∆
\sigma
\epsilon |
Specific equations from this rate-jump method have been developed for specific test platforms. For example, in depth-sensing nanoindentation, the elastic modulus and hardness are evaluated at the onset of an unloading stage following a load-hold stage. Such an onset point for unloading is a rate-jump point, and solving the equation
\epsilon
\sigma
1 | |
Se |
=
1 | |
2Era |
=
| |||
|
=
1 | |
S |
-
| |||||
|
h |
h
P |
u
Se
Er
a
f(hc)=\pia2
hc
S
Se
hc=hmax-\epsilon
Pmax | |
Se |
=hmax-\epsilonPmax\left(
1 | - | |
S |
| |||||
|
\right)
\epsilon
While nanoindentation testing can be relatively simple, the interpretation of results is challenging. One of the main challenges is the use of proper tip depending on the application and proper interpretation of the results. For instance, it has been shown that the elastic modulus can be tip dependent.[24]
Indent depths during nanoindentation can vary from a few nm up to around a micron. Over this range, there are strong "size effects" – ie inferred mechanical properties exhibit a dependence on depth. Several reviews[30] [31] [32] cover these effects. They usually take the form of the material apparently becoming harder for shallower depths. For example, the hardness of pure gold has been found[31] to vary from about 2 GPa for a depth of 5 nm to 0.5 GPa for a depth of 100 nm, while the "correct" value for large scale indentation of such gold is about 0.1 GPa. There have been many investigations of the causes of this effect. Postulated explanations include the need to create very high gradients of plastic strain with small indents, requiring "geometrically necessary dislocations".[33] [34] [35] Another suggestion[36] [37] [38] is that there may be no dislocations in the region being deformed, with the need for their nucleation creating a requirement for higher stresses to allow plastic deformation to start (leading to a "pop-in" feature on a load-displacement plot). However, no systematic, universal correction can be made for such "size effects" and it is not normally possible with nanoindenters to deform a volume that is large enough to be representative of the bulk material. For typical polycrystalline samples, such volumes must contain a relatively large number of grains, so as to capture the effects of grain size, texture, grain boundary structure etc. In practice, this usually requires the dimensions of the deformed region to be of the order of hundreds of microns. It may also be noted that the fine scale of nanoindentation can make the outcome sensitive to surface roughness[39] and to the presence of oxide layers and other surface contamination.