The ‘pop-in’ phenomenon

The nanoindentation (or instrumented or depth sensing indentation) is a variety of indentation hardness tests applied to small volumes. During nanoindentation, an indenter is brought into contact with a sample and mechanically loaded.

The following parts give a short overview of models existing in the literature used for the extraction of mechanical properties of homogeneous bulk materials from indentation experiments with conical indenters.

Parameters such as contact load \(F_\text{c}\) and depth of penetration \(h\) are recorded at a rapid rate (normally 10Hz) during loading and unloading steps of the indentation test. Usually, the depth resolution is around the fraction of \(\text{nm}\)-level and the resolution is around \(\text{nN}\)-level.

During the loading step of nanoindentation, a discontinuity in the measured depth is commonly referred to as a pop-in event.

The pop-in event

A pop-in (event) is a sudden (load or displacement) burst during the loading of an indenter on a sample. If the nanoindentation experiment is load-controlled, an horizontal plateau is also observed on the load-displacement curve, when a pop-in occurs at the critical load \(F_\text{crit}\) and critical displacement \(h_\text{crit}\) (see Figure 7). In the case of a displacement-controlled nanoindentation experiments, a vertical drop of the load is observed on the load-displacement curve.

_images/load-disp_curve_popin.png

Figure 7 Schematics of indentation load-displacement curve with a pop-in: a) load-controlled and b) displacement-controlled nanoindentation experiments.

\(F_\text{c,max}\) is the maximum applied load, \(h_\text{t}\) is the total indentation depth and \(h_\text{r}\) is the residual indentation depth.

Find here the Matlab function to quantify pop-in by peak detection: peakdet.m.

Warning

In this toolbox, only load-controlled nanoindentation experiments are analyzed.

Many authors observed pop-in events on metals or metallic thin films, ceramics, metallic glasses, semiconductors, hard brittle thin films deposited on a soft elastoplastic substrate…

The pop-in event is often explained by on of the following mechanisms as a function of the indented specimen (ceramic, metal, semiconductor, coated or multilayer specimen…) and the experimental conditions (time, temperature, geometry of the indenter…) :

It is really important to mention that pop-in events are function of the indenter shape [20], [9], the temperature [21] and the strain rate [20], [18].

In parallel, pop-in event is not always a perfect plateau on the load-displacement curve. Sometimes, sliding pop-in are occuring with no observable end point.

_images/load-disp_curve_popin_sliding.png

Figure 8 Schematics of indentation load-displacement curve with a sliding pop-in.

According to the litterature, such sliding pop-in are detected during experiments given the high-resolution of transducer [7].

Some authors proposed to describe the statistics of the pop-in event with a cumulative Weibull-type distribution [5] or with a cumulative fraction function based on a rate equation, when a time or a temperature dependence of the pop-in is demonstrated.

Weibull-type distribution

The cumulative Weibull distribution function with 2 parameters \((\lambda, m)\) is a continuous probability distribution and is often used in the description of particle size distribution and in survival or failure analysis [29].

(1)\[V\left(x, \lambda, m\right) = exp\left(\left(\frac{x}{\lambda}\right)^m\right)\]
(2)\[W\left(x, \lambda, m\right) = 1 - exp\left(\left(\frac{x}{\lambda}\right)^m\right)\]

With \(V\) the survival probability function and \(W\) the mortal probability function. \(m\) is a dimensionless material-dependent constant, often named the Weibull modulus (from 0 to usually 50) [1]. If \(m=1\), the rate of failure remains constant and there is random failure occurring.

In the case of several indentations performed on the same sample, \(x\) can be the distribution of the critical loads \(F_\text{crit}\) or the critical displacements \(h_\text{crit}\), at which the pop-in events appear on the corresponding load-displacement curves. In this case, higher is \(m\), more homogeneous is the distribution of the pop-in.

\(\lambda\) is the scale parameter of the cumulative Weibull distribution. In the case of indentations, \(\lambda\) is the mean critical load \(F^0_{crit}\) or the mean critical displacement \(h^0_{crit}\) at which the pop-in event appears for a given material.

Find here the Matlab function to calculate the cumulative mortal Weibull distribution: weibull_cdf.m.

Find here the Matlab function to calculate the complementary cumulative survival Weibull distribution: weibull_cdf_survival.m.

Chechenin et al. proposed to use a modified cumulative Weibull distribution function for the description of the statistics of the pop-in event [5]. This function is set to have a probability of 0.5, when \(F_{crit}\) (the critical load) is equal to \(F^0_{crit}\) (the mean critical load).

(3)\[W\left(\frac{F_\text{crit}}{F^0_\text{crit}}\right) = 1 - exp\left(-ln2\left(\frac{F_\text{crit}}{F^0_\text{crit}}\right)^m\right)\]

The cumulative Weibull distribution [29] and the modified cumulative Weibull distribution [5] are implemented in the PopIn toolbox.

Find here the Matlab function to calculate the cumulative modified mortal Weibull distribution: weibull_modified_cdf.m.

Find here the Matlab function to calculate the complementary cumulative modified survival Weibull distribution: weibull_modified_cdf_survival.m.

Dislocations nucleation

According to many authors, pop-in corresponds to the nucleation of at least one dislocation, during indentation of crystalline materials (see Figure 9) [20], [21], [22], [14], [30], [17], [13], [19], [31] and [32]. Before the pop-in event, the indentation load-displacement obtained with a spherical indenter, can be fitted with the Hertz equation (see Figure 10) [10]:

(4)\[F_\text{c} = (4/3) E^{*} \sqrt{R^{*}} h_\text{crit}^{1.5}\]

With \(E^{*}\) the reduced Young’s modulus calculated from indentation test defined by (5) and \(R^{*}\) the reduced radius of the spherical indenter defined by (6).

(5)\[\frac{1}{E^{*}} = \frac{1}{E_\text{s}^{'}} + \frac{1}{E_\text{i}^{'}}\]

With \(E_\text{s}^{'}\) the reduced Young’s modulus of the sample and \(E_\text{i}^{'}\) the reduced Young’s modulus of the indenter.

(6)\[\frac{1}{R^{*}} = \frac{1}{R_\text{s}} + \frac{1}{R_\text{i}}\]

With \(R_\text{s}\) the radius of the sample (usually \(+\infty\)) and \(R_\text{i}\) the radius of the indenter.

Find here the Matlab function to calculate the load with the Hertz equation: elasticLoad.m.

Find here the Matlab function to calculate the displacement with the Hertz equation: elasticDisp.m.

Find here the Matlab function to calculate the contact radius with the Hertz equation: elasticRadius.m.

Find here the Matlab function to calculate the reduced value of a variable: reducedValue.m.

_images/popin_mechanisms_1.png

Figure 9 Schematics cross section of deformation profile of an elastic-plastic substrate under indentation : 1) elastic deformation, 2) elastoplastic deformation (nucleation of dislocation) and 3) transfer of dislocations across a grain boundary.

_images/load-disp_curve_two_popin_Hertzian_fit.png

Figure 10 Schematic of indentation load-displacement curve with two pop-in events (the 1st for the nucleation of dislocation and the 2nd for the strain transfer across a grain boundary).

The elastic-plastic transition is often studied statistically as a function of temperature \(T\) and indentation rate \(F(t)\). Schuh C.A. et al. observed a good consistency between experimental results and a thermally activated mechanism of incipient plasticity [20], [21], [22], [14]. They also proposed to use a statistical thermal activation model (see following equations) with a stress-biasing term. Based on this model, it is possible to extract the activation energy \(\epsilon-\sigma V\), activation volume \(V\), and attempt frequency for the rate-limiting event that controls yield \(W(F_\text{c})\). The activation energy \(\epsilon-\sigma V\) is separated into a term characterizing the activation enthalpy \(\epsilon\), and a term giving the stress bias \(\sigma V\), where \(\sigma\) is the biasing stress over the activation volume \(V\).

(7)\[\dot{n} = \eta exp\left(-\frac{\epsilon-\sigma V}{kT}\right)\]
(8)\[\dot{N} = \eta exp\left(-\frac{\epsilon}{kT}\right)\cdot\int\int\limits_\Omega\int exp \left(\frac{\sigma V}{kT}\right)d\Omega\]
(9)\[\dot{F}(t) = \left[1-F(t)\right]{N}(t)\]
(10)\[F(t) = 1- exp\left(-\int_{0}^{t}{\dot{N}(t')dt'}\right)\]
(11)\[\tau_\text{max} = 0.31p_0 = 0.47p_\text{max} = \left(\frac{0.47}{\pi}\right)\left(\frac{4E^{*}}{3R^{*}}\right)F^{1/3}_{crit}\]
(12)\[\Omega \approx Ka^3 = K\left(\frac{3F_{crit}R^{*}}{4E^{*}}\right)\]
(13)\[F_\text{c} = \dot{F_\text{c}}\cdot t\]

With \(\dot{n}\) is the local rate at which the critical event occurs per unit volume of material, \(\eta\) the pre-exponential frequency factor, \(k\) the Boltzmann constant, \(\tau_\text{max}\) the maximun shear stress obtained for the maximum pressure \(p_\text{max}\), at a single point beneath the indenter given for an elastic Hertzian contact [10]. \(K\) is a proportionality constant of order \(\pi\).

Find here the Matlab function to calculate the maximum shear stress: maxShearStress.m.

Find here the Matlab function to calculate the maximum pressure: maxPressure.m.

Find here the Matlab function to calculate the mean pressure: meanPressure.m.

Finally, a statistical expression of the onset of plasticity can be formulated by combining previous equations, giving the cumulative fraction function \(W(F_\text{c})\) :

(14)\[W(F_\text{c}) = 1-exp\left(-\frac{9KR^{*}\eta}{4E^{*}\dot{F_\text{c}}\alpha^6}exp\left(-\frac{\epsilon}{kT}\right)\left(\beta(\alpha,F_\text{c})\right)\right)\]
(15)\[\alpha = \left(\frac{0.47}{\pi}\right)\left(\frac{4E^{*}}{3R^{*}}\right)^{2/3}\frac{V}{kT}\]
(16)\[\beta = 120exp(-F_\text{c}^{1/3}\alpha) + F_\text{c}^{5/3}\alpha^5 - 5F_\text{c}^{4/3}\alpha^4 + 20F_\text{c}\alpha^3-60F_\text{c}^{2/3}\alpha^2+120F_\text{c}^{1/3}\alpha-120\]

Find here the Matlab function to calculate the cumulative survival distribution as a function of the load rate and the temperature of nanoindentation tests: mason_cdf.m.

Find here the Matlab function to calculate the \(\alpha\) function: alphaMason.m.

Find here the Matlab function to calculate the \(\beta\) function: betaMason.m.

Note

The surface mechanical state (presence of dislocations after polishing steps or surface free of dislocation) can modified the statistics of pop-in behavior [27].

Note

See this Github repository for the plot of stress distributions at the surface and along the axis of symmetry, caused by Hertz pressure acting on a circular area radius.

In the case of prismatic dislocation loops beneath the indenter generated during pop-in, an energy balance model was also proposed by Wang D. et al. in 2021 [28].

Finally, after a short survey of the literature, it appears that pop-in or strain burst events are observed as well, on the load-displacement or stress-strain curves obtained during compression tests of metallic nano- or micro-pillar [3], [6] and [11]. This strain burst is most of the time, attributed to the nucleation of dislocations into the pillar, initially dislocation-free.

Strain transfer across grain boundaries in metals

Some authors observed sometimes two pop-in on the load-displacement curve during indentation performed close to a grain boundary (see Figure 9 and Figure 10). The first pop-in is usually attributed to the nucleation of dislocation in a metallic material (see previous section of this documentation), and the second pop-in is related to the presence nearby of the grain boundary.

The occurrence of such a strain burst is found to be related to the slip activity (function of the phase material and the grain orientation), to the grain boundary resistance (function of the grain boundary misorientation [26], to the local chemistry [4] (impurities, embrittled hydrogen…), to the prior plastic deformation [4], to the distance between the indenter and the grain boundary [26], or other experimental parameters like the shape of the indenter and the grain boundary inclination…

Recently, the STABiX Matlab toolbox was developed to analyse in simple way slip transmission in a bicrystal [15], [24] and [23].

Rupture of a hard brittle film on an elastic-plastic substrate

In the case of indentation made into a hard brittle film (e.g.: native or thermally/anodically grown oxide, ALD coatings …) on an elastic-plastic (ductile) substrate, pop-in were observed experimentally and linked to the fracture of the brittle film [12], [2], [16] and [25].

_images/load-disp_curve_popin_load_sub.png

Figure 11 Schematic of indentation load-displacement curve with a pop-in event and the plot of the load carried by the elastic-plastic substrate.

Some authors explained that a circumferential crack appears at the location of the elastic–plastic boundary in the substrate [2] and [16]. The radius \(c\) of this plastic zone in the ductile substrate is defined by the following equation:

(17)\[c = \sqrt{\frac{3F_\text{crit,s}}{2\pi\sigma_\text{e}}}\]

With \(F_\text{crit,s}\) the load carried by the substrate at the critical indentation depth at which pop-in occurs, obtained usually with the power law relationship (18). \(\sigma_\text{e}\) is the yield stress of the ductile substrate [8].

(18)\[F_\text{crit,s} = Kh_\text{crit}^n\]

With \(h_\text{crit}\) the critical displacement at which the pop-in appears. The two constants \(K\) and \(n\) are obtained from a fitting procedure of the load-displacement curve obtained from nanoindentation tests performed on the substrate without the brittle film.

_images/popin_mechanisms_2.png

Figure 12 Schematic cross section of deformation profile of a hard brittle film on an elastic-plastic substrate under indentation.

In the Figure 12, \(R\) is the radius of the spherical indenter, \(t\) is the thickness of the thin film, \(F_\text{c}\) is the applied load, \(h\) is the indentation displacement, \(c\) the radius of the plastic zone in the substrate and \(a_\text{c}\) the contact radius between the indenter and the thin film.

Find here the Matlab function to calculate the plastic radius \(c\): plasticRadius.m.

Find here the Matlab function to fit a load-displacement curve: load_displacement_fit.m.

References

[1]Afferante L. et al., “Is Weibull’s modulus really a material constant? Example case with interacting collinear cracks” (2006).
[2](1, 2) Bahr D.F. et al., “Indentation induced film fracture in hard film – soft substrate systems” (2003).
[3]Bei H. et al., “Effects of pre-strain on the compressive stress–strain response of Mo-alloy single-crystal micropillars” (2008).
[4](1, 2) Britton T.B. et al., “Nanoindentation study of slip transfer phenomenon at grain boundaries” (2011).
[5](1, 2, 3) Chechenin N.G. et al., “Nanoindentation of amorphous aluminum oxide films II. Critical parameters for the breakthrough and a membrane effect in thin hard films on soft substrates.” (1995).
[6]Crosby T. et al., “The origin of strain avalanches in sub-micron plasticity of fcc metals” (2015).
[7]Fu K. et al., “Toughness Assessment and Fracture Mechanism of Brittle Thin Films Under Nano-Indentation” (2016).
[8]Hainsworth S.V. et al., “Analysis of nanoindentation load-displacement loading curves.” (1996).
[9]Jang J.-I. et al., “Rate-dependent inhomogeneous-to-homogeneous transition of plastic flows during nanoindentation of bulk metallic glasses: Fact or artifact?” (2007).
[10](1, 2) Johnson K.L., “Contact Mechanics” (1987).
[11]Kraft O. et al., “Plasticity in Confined Dimensions” (2010).
[12]Kramer D.E. et al., “Surface constrained plasticity: Oxide rupture and the yield point process” (2001).
[13]Lu J.-Y. et al. “Thermally activated pop-in and indentation size effects in GaN films” (2012).
[14](1, 2) Mason J. et al., “Determining the activation energy and volume for the onset of plasticity during nanoindentation” (2006).
[15]Mercier D. et al. “A Matlab toolbox to analyze slip transfer through grain boundaries” (2015).
[16](1, 2) Morash K.R. and Bahr D.F., “An energy method to analyze through thickness thin film fracture during indentation” (2007).
[17]Morris J.R. et al., “Size Effects and Stochastic Behavior of Nanoindentation Pop In” (2011).
[18]Nieh T.G. et al., “Strain rate-dependent deformation in bulk metallic glasses” (2002).
[19]Ramalingam S. et al. “Determining Activation Volume for the Pressure-Induced Phase Transformation in β-Eucryptite Through Nanoindentation” (2012).
[20](1, 2, 3, 4) Schuh C.A. and Lund A.C., “Application of nucleation theory to the rate dependence of incipient plasticity during nanoindentation” (2004).
[21](1, 2, 3) Schuh C.A. et al., “Quantitative insight into dislocation nucleation from high-temperature nanoindentation experiments” (2005).
[22](1, 2) Schuh C.A., “Nanoindentation studies of materials” (2006).
[23]STABiX toolbox repository
[24]STABiX toolbox documentation
[25]Stauffer D.D. et al., “Plastic response of the native oxide on Cr and Al thin films from in situ conductive nanoindentation” (2012).
[26](1, 2) Wang M.G. and Ngan A.H.W., “Indentation strain burst phenomenon induced by grain boundaries in niobium” (2004).
[27]Wang L., “Influences of sample preparation on the indentation size effect and nanoindentation pop-in on nickel” PhD Thesis (2012).
[28]Wang D. et al., “Understanding the hydrogen effect on pop-in behavior of an equiatomic high-entropy alloy during in-situ nanoindentation”, Journal of Materials Science & Technology (2021).
[29](1, 2) Weibull W., “A statistical distribution function of wide applicability”, J. Appl. Mech.-Trans. ASME (1951), 18(3).
[30]Wo P.C. et al., “Time-dependent incipient plasticity in Ni3Al as observed in nanoindentation” (2006).
[31]Wu D. et al., “Effect of tip radius on the incipient plasticity of chromium studied by nanoindentation” (2014).
[32]Wu D. and Nieh T.G., “Incipient plasticity and dislocation nucleation in body-centered cubic chromium” (2014).