A New Elementary Method for Determining the Tip Radius and Young’s Modulus in AFM Spherical Indentations

Atomic force microscopy (AFM) is a powerful tool for characterizing biological materials at the nanoscale utilizing the AFM nanoindentation method. When testing biological materials, spherical indenters are typically employed to reduce the possibility of damaging the sample. The accuracy of determining Young’s modulus depends, among other factors, on the calibration of the indenter, i.e., the determination of the tip radius. This paper demonstrates that the tip radius can be approximately calculated using a single force–indentation curve on an unknown, soft sample without performing any additional experimental calibration process. The proposed method is based on plotting a tangent line on the force indentation curve at the maximum indentation depth. Subsequently, using equations that relate the applied force, maximum indentation depth, and the tip radius, the calculation of the tip radius becomes trivial. It is significant to note that the method requires only a single force–indentation curve and does not necessitate knowledge of the sample’s Young’s modulus. Consequently, the determination of both the sample’s Young’s modulus and the tip radius can be performed simultaneously. Thus, the experimental effort is significantly reduced. The method was tested on 80 force–indentation curves obtained on an agarose gel, and the results were accurate.


Introduction
Atomic force microscopy (AFM) is a powerful tool that enables imaging and mechanical characterization of soft materials at the nanoscale [1,2].Using the AFM nanoindentation method, Young's modulus maps on biological materials can be created and used for the diagnosis of various diseases [2,3].In particular, groundbreaking research has shown that utilizing the AFM nanoindentation method can lead to the discrimination of cells as normal or cancerous [4]; to the characterization of human tissues as normal, benign or malignant [5,6]; to the early diagnosis of osteoarthritis [7]; to the mechanical characterization of proteins [8,9] and viruses [10]; and so on.The significant advantage of this approach lies in its user-independent disease diagnosis, which can be executed through mathematical criteria and automated computational processes.Nonetheless, several challenges remain to be addressed prior to the complete use of AFM technology in clinical activities [11].The principal goal regarding the AFM research on biological materials is to develop a reliable system used to characterize biological materials such as cells and tissues at the nanoscale and used for medical purposes.Towards this direction it is essential to develop intelligent systems in order to reduce the complexity and the experimental effort.An AFM nanoindentation experiment, requires the calibration of the AFM tip used for the experiments (i.e., the determination of the indenter's dimensions).In many cases, a spherical indenter is preferred, since it reduces the possibility of damaging the soft biological material.Thus, it is essential to determine the indenter's radius prior each experiment.When testing soft biological materials at the nanoscale using spherical indenters (and assuming that the indenter is orders of magnitude stiffer than the sample), the classical Hertz equation is commonly employed for data processing [12]: In Equation (1), F is the applied force on the sample, h is the indentation depth, R is the indenter's radius, and E, v are the Young's modulus and the Poisson's ratio of the material, respectively.However, Equation ( 1) is only valid for small indentation depths compared to the tip's radius (h R) [12].The accurate equation that relates the applied force to the indentation depth was firstly derived by Sneddon and is presented below [13]: In Equation ( 2), r c is the radius at contact depth (h c ) [14].In addition, the indentation depth is related to the contact radius with the following equation [13]: Equations ( 2) and ( 3) do not provide a direct relation between the applied force and the indentation depth.Thus, a new equation was recently derived [15].The idea was to write Equation (3) as follows: Subsequently, the r c R = f h R function was fitted to a simple equation of the form [15]: In Equation ( 5), c 1 , c 2 , . . ., c N are constants that depend on the h max /R ratio and can be found in [15].Equation (5) was substituted to the general differential equation that is valid for every axisymmetric indenter [16]: Subsequently, the solution of differential Equation ( 6) results in In Equation (7), Q is defined as follows [15]: It is important to note that Equations ( 2) and ( 7) yield identical results, as shown in Figure 1a,b.The functions 2) and ( 7) are presented for comparison (where, E * = E 1−v 2 is the sample's reduced modulus).For 0 ≤ h/R ≤ 1.32, N = 3 and c 1 = 1.022, c 2 = −0.1133,c 3 = −0.0742(Figure 1a) [15].For 0 ≤ h/R ≤ 4.9512, N = 6 and c 1 = 1.0100000, c 2 = −0.0730300,c 3 = −0.1357000,c 4 = 0.0359800, c 5 = −0.0040240,c 6 = 0.0001653 (Figure 1b) [15].It is important to further note that the latter case (N = 6) is applicable in any scenario [15].However, it is too complicated.Therefore, in conventional experiments, the values N = 3 and c 1 = 1.022, c 2 = −0.1133,c 3 = −0.0742are not only precise but also greatly reduce complexity.
Micromachines 2023, 14, x FOR PEER REVIEW 3 of 19 1.022,  2 = −0.1133, 3 = −0.0742are not only precise but also greatly reduce complexity.In addition, the F 2E * R 2 = f h R functions when using Equations ( 1) and ( 7) are shown comparatively in Figure 1c for the domain 0 ≤ h R ≤ 0.25 and in Figure 1d for the domain 0 ≤ h R ≤ 1.32.If Equation ( 1) is used instead of Equation ( 7) for h max /R = 1, the error in the Young's modulus calculation will be approximately 10%.Therefore, when conducting AFM nanoindentation tests on soft materials with spherical indenters, the most suitable equation for fitting the force-indentation data is Equation (7).However, it will be demonstrated that the preference for Equation (7) over Equation ( 1) is not solely based on avoiding errors in Young's modulus calculation.As already mentioned, the determination of the Young's modulus of the tested material necessitates knowledge of the tip radius.The tip radius is usually determined using scanning electron microscopy (SEM) imaging [17] or AFM gratings [18].This represents an extra experimental stage in AFM indentation experiments, considerably extending the time needed to determine the mechanical properties of the material.This is due to the necessity of using a new indenter for each experiment to prevent contamination or changes in the tip's shape or dimensions (hence, prior to any experiment, a new calibration process is essential).In addition, even one additional experimental step increases the possibilities of contaminating the AFM tip.
This paper will demonstrate that by employing a single force-indentation curve, it is possible to determine both the tip radius and the Young's modulus.Hence, the need for additional experimental processes for tip calibration can be readily circumvented.Therefore, the significant benefit of fitting the data to Equation (7), in addition to the offered accuracy in comparison to Equation (1), is the potential to conduct an AFM tip calibration without requiring any supplementary experimental procedures.

A New Model for Processing Force-Indentation Curves
When employing spherical indenters, the data adhere to Equations ( 7) and (8).However, in the majority of cases, the maximum indentation depth is h max ≤ R. It has been previously demonstrated that for h max R ≤ 1.32 [15]: As already mentioned in the introduction, the constants c 1 , c 2 , c 3 , for the domain h max /R ≤ 1.32, are c 1 = 1.022, c 2 = −0.1133,c 3 = −0.0742[15].The slope of the forceindentation curve at any given point is defined as the contact stiffness: The contact stiffness at the maximum indentation depth is presented below: Furthermore, the equation of the tangent line to the force-indentation curve at h = h max is given by: In Equation ( 12), b represents the point of intersection between the tangent line (12) and the force axis.By substituting h = h max into Equation (12), the value of b can be readily calculated: In addition, using also Equations ( 9) and (11), Equation ( 14) leads to an interesting conclusion; the Young's modulus can be easily calculated by plotting the tangent line of the force indentation curve at the point h = h max .By utilizing a linear fit to the tangent line, b, S can be determined as fitting coefficients.Subsequently, the Young's modulus can be easily calculated using Equation (14), assuming that the sample's Poisson's ratio is known.
The analysis can be extended by calculating the point of intersection between the tangent line (Equation ( 12)) and the indentation axis.At the point of intersection F = 0, thus By combining Equations ( 11), ( 14) and (15), it is concluded Equation ( 16) yields an intriguing outcome.If the point of intersection between the tangent line and the indentation axis is determined, the indenter's radius can be calculated using Equation (16).This outcome holds significance as the indenter's radius can be calculated without requiring knowledge of the sample's Young's modulus.Therefore, by employing Equation ( 16), the indenter's radius can be ascertained, followed by an easy calculation of the Young's modulus using Equation (14).The procedure is also presented in Figure 1a and is summarized as follows.If h max R ≤ 1.32, the force indentation data can be fitted to a simple equation of the form In Equation ( 17), a, b, c are fitting parameters.Subsequently, the tangent line (Equation ( 12)) is plotted at the point h = h max of the fitted curve.The factors b, S are determined as fitting parameters and the tip radius and the Young's modulus are calculated using the Equations ( 16) and ( 14), respectively.Equations ( 14) and ( 16) represent simplified forms of the general case for h max R ≤ 1.32 and are applicable to the majority of cases.In addition, it is straightforward to determine the suitability of the set of Equations ( 14) and ( 16), for the experiments by observing the nominal tip radius.For example, let us assume that the nominal tip radius provided by the manufacturer is 1 µm.If the maximum indentation depth is significantly smaller (e.g., ∼0.5 µm), then it is safe to use the simplified Equations ( 14) and ( 16).However, it is significant to note that the method has no restrictions regarding the maximum indentation depth and it can be also applied for any h max R ratio.The general case is presented below: For example, for h max R ≤ 4.9512, N = 6, c 1 = 1.0100000, c 2 = −0.0730300,c 3 = −0.1357000,c 4 = 0.0359800, c 5 = −0.0040240,and c 6 = 0.0001653 [15] (as also mentioned in the introduction).Thus, in this case, the contact stiffness at h = h max is given by the following equation: Furthermore, Equation ( 14) is also adjusted as shown below: Thus, the point of intersection between the tangent line and the indentation axis becomes Furthermore, Equation ( 21) can be expressed as follows: The graphical representation of Equation ( 22) is presented in Figure 2b.It is interesting to note that for big h max R ratios, h com.R tends to a limit value which its equal to 0.5.In addition, Equation ( 20) can be also written in the form Micromachines 2023, 14, x FOR PEER REVIEW 7 of 19 .The slope, S, and the point of intersection between the tangent line and the force axis, denoted as b, are determined as fitting parameters.Subsequently, the indenter's radius can be calculated using Equation (16).Lastly, the Young's modulus is determined using Equation ( 14).(b) The

Contact Point Determination
In AFM indentation experiments on biological samples, a critical aspect affecting result accuracy is the determination of the contact point between the tip and the sample.For precise determination of the contact point, the AtomicJ software (https://sourceforge.net/projects/jrobust/) was utilized [19].The procedure is straightforward: each point of the curve is taken as a trial contact point, a polynomial is fitted to the precontact section, and the suitable contact model is applied to the force-indentation data [19].The tested point that resulted in the lowest total sum of squares is considered as the contact point [19].

Measurements
Spherical indenters (borosilicate glass spheres with Young's modulus 64 GPa) were employed for the AFM indentation experiments (CP-PNPL-BSG-A, sQube, obtained from NanoAndMore GMBH, Wetzlar, Germany).The nominal tip radius of such an indenter, .The slope, S, and the point of intersection between the tangent line and the force axis, denoted as b, are determined as fitting parameters.Subsequently, the indenter's radius can be calculated using Equation (16).Lastly, the Young's modulus is determined using Equation ( 14).
function is presented in Figure 2c.For big h max R ratios the factor b 2E * R 2 tends to −0.5.However, it is crucial to note that attaining the limit values is challenging in real experiments.For example, the probability of a plastic deformation is high in such cases.However, they present significant mathematical interest and can be used for a deeper understanding of the underlying theory.For example, assume that h max R = 5.
In this case, h com.R = 0.5.This result can be used as a paradigm to realize that the relation between h com.and h max does not depend on the Young's modulus value.This outcome is of paramount significance, as the method's reliability is not contingent on the sample type when it comes to calculating R. It is also noteworthy that the independence of the relation between h com.and h max from the Young's modulus applies to any indentation depth, as indicated by Equation (21).

AFM Indentation Experiments on an Agarose Gel 2.2.1. Contact Point Determination
In AFM indentation experiments on biological samples, a critical aspect affecting result accuracy is the determination of the contact point between the tip and the sample.For precise determination of the contact point, the AtomicJ software (https://sourceforge.net/ projects/jrobust/) was utilized [19].The procedure is straightforward: each point of the curve is taken as a trial contact point, a polynomial is fitted to the precontact section, and the suitable contact model is applied to the force-indentation data [19].The tested point that resulted in the lowest total sum of squares is considered as the contact point [19].

Measurements
Spherical indenters (borosilicate glass spheres with Young's modulus 64 GPa) were employed for the AFM indentation experiments (CP-PNPL-BSG-A, sQube, obtained from NanoAndMore GMBH, Wetzlar, Germany).The nominal tip radius of such an indenter, as specified by the manufacturer, is 1 µm with a deviation of ±10% (i.e., 0.9 µm ≤ R ≤ 1.1 µm).The indenters were calibrated prior to the experiments using the AFM test grating TGT1 from NT-MDT Instruments.For precise quantitative measurements, it is necessary to calibrate the probe parameters.To calibrate the laser detection system's sensitivity in terms of nanometer deflection per volt signal, a force vs. distance curve on mica was firstly obtained [20].By positioning two cursors on the contact section of the force vs. distance curve, the deflection sensitivity was determined [20].The spring's constant determination was performed using the thermal noise method.The experiments were performed on agarose gels with concentration 2.5% in a 35 mm petri dish.The Poisson's ratio of the agarose gel can be considered equal to v = 0.5 due to the high-water content.The Young's modulus was calculated through conventional fitting procedures using Equation (9) (and using also the measured R-value), as well as employing the method proposed in this paper (i.e., fitting the tangent line to Equation (12) and subsequently employing Equation ( 16) to determine the tip radius and ( 14) to determine the Young's modulus).

Results
An example illustrating the application of the proposed method is provided below.The experiment was performed using a spherical tip with a nominal radius equal to 1 µm.The tip was calibrated using an AFM grating as described in Section 2.2.2, and the result was R meas.= 0.92 µm (this result is within the ±10% range provided by the manufacturer, 0.9 µm ≤ R meas.≤ 1.1 µm).The calibration of the spherical indenter using the AFM grating is presented in Figure 3.The force-indentation data are shown in Figure 4a.The data were fitted to Equation ( 17): The R-squared coefficient resulted in R 2 s.c.= 0.9808.The force-indentation data and the fitted curve (Equation ( 24)) are presented comparatively in Figure 4a.Subsequently, the tangent line at h = h max = 432 nm was plotted.
The tangent line is described by the following Equation (Figure 4b): Thus, S = 0.2505 N/m, b = −3.363•10−8 N and h com.= − b S = 1.342•10 −7 m.In addition, Equation ( 16) is written as follows: Equation ( 26) can be easily solved using any basic software (e.g., Matlab).The tip radius resulted in R = 0.921•10 −6 m = 0.921 µm.This result is nearly identical to the value that was measured using the AFM grating.
A graphical solution of Equation ( 26) is also presented in Figure 4c.In particular, the functions y 1 (R) = 1.342•10 −7 m and were plotted within the domain 0.8•10 −6 m ≤ R ≤ 1.2•10 −6 m.The common point of the two functions is R = 0.921•10 −6 m.In addition, it is also easy to calculate the Young's modulus using Equation ( 14): By solving Equation ( 27), the Young's modulus resulted in E = 132.5 kPa.Three additional paradigms are also presented in Figure 5.To validate the accuracy of the method, an additional 76 force-indentation curves were also processed.The results are shown in Figure 6. Figure 6a displays the outcomes concerning the tip radius.The mean ± standard deviation value resulted in 0.9184 µm ± 0.0135 µm, which is in agreement with the result obtained using the AFM grating.A histogram of the R-values is also presented in Figure 6b.In Figure 6c, the Young's modulus values as calculated using the proposed by this paper approach (Equations ( 14) and ( 16)) and using a classic fitting procedure are presented for comparison.The results are nearly identical.
Micromachines 2023, 14, x FOR PEER REVIEW 8 of 19 1.1 μm).The indenters were calibrated prior to the experiments using the AFM test grating TGT1 from NT-MDT Instruments.For precise quantitative measurements, it is necessary to calibrate the probe parameters.To calibrate the laser detection system's sensitivity in terms of nanometer deflection per volt signal, a force vs. distance curve on mica was firstly obtained [20].By positioning two cursors on the contact section of the force vs. distance curve, the deflection sensitivity was determined [20].The spring's constant determination was performed using the thermal noise method.The experiments were performed on agarose gels with concentration 2.5% in a 35 mm petri dish.The Poisson's ratio of the agarose gel can be considered equal to  = 0.5 due to the high-water content.The Young's modulus was calculated through conventional fitting procedures using Equation ( 9) (and using also the measured R-value), as well as employing the method proposed in this paper (i.e., fitting the tangent line to Equation ( 12) and subsequently employing Equation ( 16) to determine the tip radius and ( 14) to determine the Young's modulus).

Results
An example illustrating the application of the proposed method is provided below.The experiment was performed using a spherical tip with a nominal radius equal to 1 μm.The tip was calibrated using an AFM grating as described in Section 2.2.2, and the result was  .= 0.92 μm (this result is within the ±10% range provided by the manufacturer, 0.9 μm ≤  .≤ 1.1 μm).The calibration of the spherical indenter using the AFM grating is presented in Figure 3.The force-indentation data are shown in Figure 4a.The data were fitted to Equation ( 17):  = 281ℎ 3/2 − 23340ℎ 2 − 1.019 • 10 10 ℎ 3 (.. ) The R-squared coefficient resulted in  .. 2 = 0.9808.The force-indentation data and the fitted curve (Equation ( 24)) are presented comparatively in Figure 4a.Subsequently, the tangent line at ℎ = ℎ  = 432 nm was plotted.16) is written as follows:

Discussion
In this paper, a new method is presented that allows for the determination of the Young's modulus of soft biological materials without the need for an experimental calibration procedure for the AFM tip.When using spherical indenters, the radius of the indenter can be readily calculated by employing the tangent line of the fitted curve of the force-indentation data at h = h max .More specifically, the point of intersection between the tangent line and the indentation axis can unveil the value of R (as defined in Equations ( 16) and ( 21)).Subsequently, by utilizing the point of intersection of the tangent line with the force axis, the Young's modulus can be readily determined (as described in Equations ( 14) and ( 20)).The proposed approach is reliable, as the results were nearly identical to those obtained using conventional methods, such as tip calibration through an AFM grating and Young's modulus determination using traditional fitting procedures.The application of the new method requires solving Equations ( 16) and ( 14) (or ( 21) and ( 20) in the general case).Furthermore, utilizing Equation ( 22), an equation that establishes a numerical relationship between the ratio h com.
h max and the ratio h max R can be derived.The h max data are presented in Figure 7a.This is a noteworthy finding, as the data can be fitted to a polynomial curve to derive a straightforward equation that establishes a relationship between h com.
h max with h max R within the specified domain of interest.For example, the data h max R = f h com.
h max were fitted to the function (for the domain 0 ≤ h/R ≤ 1.32): The fit was perfect (R 2 s.c.= 1.0000).The fitting coefficients are as follows: p 4 = −8149, p 3 = 8816, p 2 = −3592 and p 1 = 637.3and p 0 = −39.27.This is a significant result, as Equation ( 28) can be used to determine the tip radius after graphically estimating the point h com.For example, let us consider the case introduced at the beginning of the results section (depicted in Figure 4).In this example, h com.= 134.2nm and h max = 432 nm.Thus, h com.
h max = 0.3106.By employing Equation (28), it is straightforward to calculate that h max R = 0.4691.Thus, R = 0.921 µm.Furthermore, to enhance the ease of applying the proposed method, two tables presenting the data for h max R = f h com.h max are also provided.In Table 1, the values of h max are presented within the domain of 0.01 ≤ h/R ≤ 1.32, and in Table 2, within the domain 0.05 ≤ h/R ≤ 5.00.It is also noteworthy to emphasize that a 'rational approach' would involve concurrently determining the tip radius and the sample's Young's modulus by employing a simple fit to Equation (9), under the assumption of h max /R ≤ 1.32.The reason is that Equation ( 9) can be written as follows: By combining Equations ( 17) and (29), it is concluded Thus, given that the coefficients c 1 , c 2 , c 3 are known, it may be assumed that the Young's modulus can be calculated using Equation (31), and subsequently, employing Equation (30) or Equation (32), the tip radius can also be determined.However, this is not accurate; there exist various combinations of the fitting coefficients a, b, and c that can result in the same curve.This fact can be also proved using Equation (24).In this case, b = 23, 340 N m 2 .Therefore, utilizing Equation ( 31), E = 2.398•10 5 Pa.In addition, a = 281 N m 3/2 .Therefore, by applying Equation (30), R = 0.225•10 −6 m, which is approximately four times smaller than the actual value (0.921 µm).Hence, it is imperative to adhere to the procedure outlined in this paper for the calculation of R and E.
The precise equations concerning deep spherical indentations are typically circumvented in experiments involving soft biological materials, as fitting the data to Equation ( 1) is simpler and the errors in Young's modulus calculations for h max /R < 1 are not substantial.However, this paper demonstrates that Equations ( 7) and ( 8) can provide significantly more options compared to Equation (1).Indeed, these equations can lead to the determination of the indenter's radius using a simple force-indentation curve.If we use Equation ( 1) instead of Equation ( 7), the contact stiffness becomes Subsequently, the point of intersection between the tangent line and the indentation axis can be calculated: Hence, when utilizing Equation ( 1), it becomes impossible to calculate R, as the values for h com.consistently equate to h max 3 .Therefore, in cases for which h max R and Equation (1) accurately describe the data, the method cannot be applied.
It is also crucial to highlight the significant reliability of the proposed method.The experiments were conducted on an agarose gel at arbitrarily selected points.The calculated Young's moduli were within the range of 102 kPa ≤ E ≤ 174 kPa (see also Figure 6c).However, despite the significant variation in the Young's modulus, there is only a slight variation in the calculation of the indenter's radius (as clearly depicted in Figure 6a).If we were to test a hypothetically perfect elastic half-space, the force-indentation data and the fitted curve would be identical, and the calculation of R would be consistent across all force-indentation curves.However, for real soft samples, the data do not perfectly follow Equation (7).Consequently, errors in h com.and b would emerge, leading to variations in R calculations.Nevertheless, the main outcome is that even though there exists a disparity between the fitted curve and the data in all instances, the inaccuracies pertaining to the tip radius calculation were exceedingly minor.Specifically, the computed values in each measurement, as depicted in Figure 6a, closely resembled the measured values obtained through the AFM grating (as shown in Figure 3).This is a logical outcome when we consider Equation (22).This equation suggests that there is a specific relationship between h com.and h max that can lead to the calculation of the ratio h max /R and, as a result, to the calculation of the indenter's radius (see also Tables 1 and 2).Importantly, this relationship remains uninfluenced by the Young's modulus.When collecting force-indentation data, h max is a known parameter and determining h com.becomes straightforward by utilizing the tangent line to the force-indentation curve at the maximum indentation depth.Thus, the ratio h com./h max unveils the ratio h max R , consequently determining the radius of the indenter.Furthermore, it is worth noting that by multiplying h com./h max with h max /R in Tables 1 and 2, a new table can be generated that establishes a relationship between h com.and R.
In addition, it is also important to note that most of the cells and biological tissues present a viscoelastic behavior [21,22].However, for small indentation rates, the elastic models derived from Hertzian mechanics can be employed for data fitting [22].From this perspective, our method remains valid for viscoelastic materials.By utilizing an extremely small indentation rate, we can determine the indenter's radius (R).Subsequently, we can employ the R-value for dynamic loading at higher indentation rates, allowing us to extract the viscoelastic properties of the material.However, a very interesting question is whether it is possible to calculate R using force-indentation curves for different loading conditions [23,24].This constitutes a fascinating task for future research.
It is also noteworthy to emphasize that the development of algorithms for intelligent micro-and nano-systems is crucial for technological progress and its applications in medicine and biology.AFM processes represent cutting-edge research in today's context, offering numerous possibilities for potential clinical applications in the future.However, in order to attain this objective, it is crucial to develop straightforward and automated procedures for data processing.Simplifying experimental procedures through algorithms based on mathematical criteria is crucial for the utilization of AFM technology in medical applications, such as disease diagnosis.

Conclusions
This paper introduces a novel method for calibrating the indenter in AFM nanoindentation experiments involving soft materials using spherical indenters.The calibration of the indenter is founded on processing the force-indentation curve employing rigorous mathematical criteria.Hence, it becomes possible to calculate the Young's modulus and the AFM tip radius using the force-indentation data without requiring any additional experimental procedures.The fundamental steps of the method are outlined as follows:

•
Plot the tangent line of the fitted curve at the maximum indentation depth.

•
Determine the point of intersection between the tangent line and the indentation axis and solve either (16) or (21).Alternatively, for simplification, employ Equation (28) or refer to Tables 1 and 2.

•
To calculate the Young's modulus at the tested point, identify the intersection point between the tangent line and the force axis, and solve either Equations (14) or (20).
The proposed approach can be integrated into typical AFM equipment to automate and streamline the experimental procedures.The development of intelligent AFM systems increases the potential for utilizing AFM processes in practical clinical applications, such as disease diagnosis.

Figure 2 .
Figure 2. (a) The blue curve represents the force-indentation curve.The red dotted line represents the tangent line at the point ℎ = ℎ  .The slope of the tangent line equals to the contact stiffness at ℎ = ℎ  (i.e.,  = tan ()).The slope, S, and the point of intersection between the tangent line and the force axis, denoted as b, are determined as fitting parameters.Subsequently, the indenter's radius can be calculated using Equation(16).Lastly, the Young's modulus is determined using Equation(14).(b) The

Figure 2 .
Figure 2. (a)The blue curve represents the force-indentation curve.The red dotted line represents the tangent line at the point h = h max .The slope of the tangent line equals to the contact stiffness at h = h max (i.e., S = tan(θ)).The slope, S, and the point of intersection between the tangent line and the force axis, denoted as b, are determined as fitting parameters.Subsequently, the indenter's radius can be calculated using Equation(16).Lastly, the Young's modulus is determined using Equation (14).

Figure 3 .
Figure 3. Calibration of the spherical tip using the AFM grating.The tip radius resulted in R = 0.92 µm.

Figure 4 .
Figure 4.A paradigm of the proposed method.(a) The force-indentation data were fitted to Equation (17) (F = 281h 3/2 + 23340h 2 + 1.019•10 10 h 3 , R 2 s.c.= 0.9808).(b) The tangent line at the maximum indentation depth h = h max = 432 nm was plotted.The point of intersection between the line and the F-axis is b = −3.363•10−8 N and the point of intersection between the line and the h-axis is h com.= 1.342•10 −7 m.(c) A graphical solution of Equation (26).

Figure 5 .
Figure 5. Three additional paradigms of the proposed method.Figure 5. Three additional paradigms of the proposed method.

Figure 5 .
Figure 5. Three additional paradigms of the proposed method.Figure 5. Three additional paradigms of the proposed method.

Figure 6 .
Figure 6.Evaluating the reliability of the proposed method.(a) The determination of the indenter's radius using 80 measurements.The mean ± standard deviation value obtained was 0.9184 μm ± 0.0135 μm, which agrees with the tip radius measurement obtained using the AFM grating.(b) A histogram constructed using the values presented in (a).(c) The Young's modulus was calculated using the method proposed in this paper and a conventional fitting procedure.The results from the two methods were nearly identical.

Figure 6 .
Figure 6.Evaluating the reliability of the proposed method.(a) The determination of the indenter's radius using 80 measurements.The mean ± standard deviation value obtained was 0.9184 µm ± 0.0135 µm, which agrees with the tip radius measurement obtained using the AFM grating.(b) A histogram constructed using the values presented in (a).(c) The Young's modulus was calculated using the method proposed in this paper and a conventional fitting procedure.The results from the two methods were nearly identical.

Figure 7 .
Figure 7. (a) The data for h max R = f h com.h max in the domain 0.05 ≤ h/R ≤ 5.00.(b) The same data within the domain 0 ≤ h/R ≤ 1.32.The data were fitted to a fourth-degree polynomial function (Equation (28)).

Table 1 .
The h com./h max ratio for the domain 0.01 ≤ h max /R ≤ 1.00.

Table 2 .
The h com./h max ratio for the domain 0.05 ≤ h max /R ≤ 5.00.