Simulation and experimental analysis of tip response of tapping mode atomic force microscope
Anjie Peng^{1} , Yaxin Chen^{2} , Rui Guan^{3} , Zheng Wei^{4}
^{1, 2, 3, 4}College of Mechanical and Electrical Engineering, Beijing University of Chemical Technology, Beijing, China
^{1, 2, 3, 4}Beijing University of Chemical Technology, Beijing, China
^{4}Corresponding author
Vibroengineering PROCEDIA, Vol. 40, 2022, p. 5056.
https://doi.org/10.21595/vp.2022.22341
Received 14 December 2021; received in revised form 27 December 2021; accepted 3 January 2022; published 5 February 2022
In this study, the vertical deflection responses of tapping mode atomic force microscope (TMAFM) microcantilever tip are obtained by simulation and experiment. The results show that, under the blocking of the sample on one side, the steadystate response of the tip is still a sinusoidal form almost symmetrical about the equilibrium position. Furthermore, from the perspective of energy dissipation of the microcantilever system, the phases of two surfaces with different properties are simulated under different background dissipation. The result shows that eliminating partial background dissipation can increase the phase contrast between the two surfaces. These results are of significance for understanding the tip response and phase optimization in TMAFM.
 It is proved by experiment and simulation that in TMAFM, the displacement response of the tip is almost a sine wave symmetrical about the equilibrium position.
 The phase response of the tip on the surface with different properties is obtained by numerical simulation.
 The simulation results show that reducing the background dissipation can improve the phase contrast between the components of the sample, which has hardly been mentioned by previous studies.
Keywords: atomic force microscope, tapping mode, displacement response, phase response.
1. Introduction
Since the emergence of atomic force microscopes (AFM), it has become one of the most powerful tools in micro and nanotechnology because of its excellent resolution [1, 2]. The key component of AFM is a microcantilever with sharp a tip. The information of the sample is recorded by transforming the interaction forces between the tip and the sample into the displacement signal of the tip [2]. Tapping mode (TM) is one of the most commonly used operation modes in atomic force microscopes. It can obtain highresolution images of various information of samples, and can greatly reduce the wear of the tip and damages to the sample compared with contact mode [3]. Therefore, tapping mode atomic force microscope (TMAFM) is widely used in various fields, including biomolecules, polymers, and nanostructures [47].
During the operation of TMAFM, the microcantilever oscillates above the sample under the excitation of the piezoelectric actuator, and the tip contacts intermittently with the sample. The interactions between the tip and sample are very complex, including van der Waals interactions, shortrange repulsive, adhesion, air damping, and capillary forces [8], which may make the response of the tip very complex. Payam [9] obtained the time history responses of the microcantilever with different distances between the tip and sample through simulation. Cleveland et al. [10] showed through experiments that the time history response of the tip vertical deflection was nearly sinusoidal. However, few researchers have conducted detailed analysis on the displacement response of the microcantilever tip in simulation and experiment.
In addition, phaseimaging is an important part of TMAFM, which is very sensitive to the properties of samples. The different components and properties of the sample are distinguished by recording the difference between the tip response signal and the excitation signal. Phaseimaging has been widely used to study various heterogeneous materials and their viscoelasticity, wettability, and so on [1113]. The phase in TMAFM is considered to be related to the energy dissipation of the system [10, 14, 15]. At present, few studies have optimized the phase contrast in TMAFM from the perspective of reducing the background dissipation of the microcantilever system.
In this paper, the time history response of the tip displacement is obtained by numerical simulation and experiment. Then, from the perspective of energy dissipation, the phases of the tip working on two surfaces with different properties under different background dissipation is simulated, and the factors affecting the phase contrast between the two surfaces are analyzed.
2. Analysis of tip response in TMAFM
2.1. Model
The key component of TMAFM is a microcantilever with a sharp tip, which is sensitive to small forces. The piezoelectric actuator drives the substrate with sinusoidal displacement to make the microcantilever oscillate near its first natural frequency. The information of the sample can be obtained by detecting the change of the response at the tip. As shown in Fig. 1(a), the length of the microcantilever is $l$, the width is $b$, and the thickness is $d$, the displacement excitation ${z}_{f}=D\mathrm{s}\mathrm{i}\mathrm{n}\varpi t$, $D$ is the excitation amplitude, $\varpi $ is the excitation frequency. The differential equation of the microcantilever is as follows:
where $w$ is the absolute displacement of the microcantilever, $u$ is the relative displacement of the microcantilever relative to the substrate, $E$ and $I$ are the Young’s Modulus and section moment of inertia of the microcantilever respectively, $\rho $ is the density of the microcantilever, $m\text{'}$ is the mass of the tip, $c$ is the equivalent damping coefficient per unit length of the microcantilever, ${F}_{ts}$ represents the interaction forces between the tip and the sample, and $\delta \left(x\right)$ Dirac function.
Fig. 1. Microcantilever model of TMAFM
a) Continuous microcantilever model
b) Simplified model of microcantilever
In TMAFM, the surface information of the sample is imaged by monitoring the response at the tip of the microcantilever. At the same time, the microcantilever works near its firstorder natural frequency, so the tip position and firstorder mode can be used to simplify the microcantilever model. As shown in Fig. 1(b), the differential equation of motion of the needle tip can be written as:
where ${m}_{e}$ is the effective mass of the microcantilever, ${k}_{e}$ is the force constant of the microcantilever, ${c}_{0}$ is the equivalent damping coefficient representing the background dissipation, ${\omega}_{0}={\left({k}_{e}/{m}_{e}\right)}^{1/2}$ is the firstorder natural angular frequency, ${Q}_{b}={m}_{e}{\omega}_{0}/{c}_{0}$ is the effective quality factor representing the background dissipation of the microcantilever system. In tapping mode, the microcantilever has various energy dissipation paths when vibrates in air, such as air viscous dissipation, squeeze film damping, liquid bridge dissipation, and so on. This type of dissipation, which does not reflect the sample information is called background dissipation.
In AFM, the interaction between the tip and the sample is very complex, and the interaction forces are usually simplified in the response analysis. Assuming that the spherical tip interacts with the flat sample surface, van der Waals forces and DMT contact force are typically used to describe attraction and repulsion forces, respectively. The forces between the tip and sample can be expressed as follows [8, 15, 16]:
where $h$ is the distance between the tip and the sample at the equilibrium position, $H$ is Hamaker’s constant, $R$ is the radius of the tip, ${a}_{0}$ is the intermolecular distance, $\eta $ is the viscosity of the sample, and ${E}^{*}$ is the equivalent Young’s modulus, which can be expressed as follows:
where, ${E}_{t}$, ${E}_{s}$, ${\nu}_{t}$ and ${\nu}_{s}$ are the Young’s modulus and Poisson’s ratio of the tip and sample, respectively.
2.2. Analysis of time history response of the tip displacement
In order to simulate the response of the tip in TMAFM, the parameters of the real micro cantilever are used. The effective quality factor ${Q}_{b}$ used in this paper is obtained by processing the results of Tune experiment by halfpower method. The excitation frequency is 99.9 % of the system’s natural frequency. The specific parameters are shown in Table 1.
Table 1. Parameters used in simulation
Symbol

Value

Symbol

Value


Effective mass (kg)

${m}_{e}$

9.13×10^{12}

Tip radius (nm)

$R$

6

Force constant (N/m)

${k}_{e}$

42

Equivalent Young’s modulus (GPa)

${E}^{*}$

10.08

Effective quality factor

${Q}_{b}$

500

Hamaker’s constant (J)

$H$

2×10^{19}

Natural frequency (kHz)

${f}_{0}$

341.36

Intermolecular distance (nm)

${a}_{0}$

0.38

Drive amplitude (nm)

$D$

0.1

Sample viscosity (Pa·s)

$\eta $

20

Drive frequency (kHz)

$f$

341.02

By substituting the parameters in Table 1 into Eq. (13), and solving Eq. (2) can obtain the tip response under the sample limitation. It is noted that the interaction forces between the tip and the sample in Eq. (3) are piecewise nonlinear. Therefore, Eq. (2) is a piecewise nonlinear secondorder differential equation. Here, the fourthorder RangeKutta algorithm is used to solve Eq. (2).
When the microcantilever is far away from the sample, the interaction forces between the tip and sample can be ignored, and Eq. (2) can be considered as a secondorder linear differential equation. The time history response of the vertical deflection of the tip is shown in Fig. 2(a). The tip oscillates from the equilibrium position, and finally, the vertical deflection of the tip tends toward a steady state. The fast Fourier transform (FFT) of the steadystate response after 5 ms is shown in Fig. 2(b). It can be seen from the FFT diagram that there is only one frequency component which is the same as the excitation frequency. These results show that the steadystate response of the tip displacement when the microcantilever is far away from the sample is an absolute sine wave symmetrical about the equilibrium position. At this time, the amplitude of the steadystate response of the tip without sample constraint is ${A}_{0}$.
Fig. 2. Time history response and FFT of the tip vertical deflection when the microcantilever is far away from the sample
a) Timehistory response of the tip vertical deflection
b) FFT of the steadystate response
Fig. 3. Time history response and FFT of the tip vertical deflection at different initial distances between tip and sample
a) Timehistory response of the tip vertical deflection when $h=$ 0.9A0
b) FFT of the steadystate response when $h=$ 0.9A0
c) Timehistory response of the tip vertical deflection when $h=$ 0.1${A}_{0}$
d) FFT of the steadystate response when $h=$ 0.1${A}_{0}$
In tapping mode, the amplitude of the tip is generally 85 %90 % of that when oscillating in air. It is assumed that the sample surface is smooth and uniform. When the initial distance between the needle tip and the sample $h=$ 0.9${A}_{0}$, the time history response and FFT of the tip vertical deflection are shown in Fig. 3(a) and 3(b), respectively. The magenta dotted line represents the position of the sample, and the orange dotted line represents the symmetrical position of the sample with respect to the equilibrium position. The FFT of steadystate response is shown in Fig. 3(b), in which there are many very small frequency doubling components and a very small direct current (DC) component in addition to the fundamental frequency. This is due to the tip intermittently contacts with the sample, and the system is no longer linear, but piecewise nonlinear. Because the frequency doubling components and DC component of the response are very small compared with the fundamental frequency, the steadystate time history response of the tip displacement can be considered as a standard sinusoidal signal symmetrical about the equilibrium position.
In a more extreme case, when the initial distance between the needle tip and the sample $h=$ 0.1${A}_{0}$, the time history response and FFT of the tip vertical deflection are shown in Fig. 3(c) and 3(d), respectively. As can be seen from the FFT diagram, the DC component and frequency doubling components are very obvious at this time. This is because the depth of the tip pressing into the sample is large, and the interaction time between the tip and the sample becomes longer in one cycle. In this case, although the DC component and frequency doubling components are obvious, they are very small compared with the fundamental frequency, so that the steadystate response can also be regarded as a sine wave.
The experiment is carried out on a Bruker’s Dimension^{}Icon AFM. The sample is a highly oriented pyrolytic graphene (HOPG). We select TappingMode to scan the sample. During operation, the vertical deflection signal of the tip is recorded by HighSpeed Data Capture (HSDC) function. Fig. 4(a) shows the time history response of the tip vertical deflection when the amplitude setpoint ${A}_{sp}=$0.889${A}_{0}$. The magenta dotted line represents the amplitude setpoint, and the orange dotted line represents the symmetrical position of the amplitude setpoint with respect to the equilibrium position. FFT is performed on the time history response data of a certain point of the sample, as shown in Fig 4(b). The results show that the vertical deflection of the tip is a sine wave almost symmetrical about the equilibrium position. The experimental and simulation results show that the steadystate response of the needle tip in intermittent contact with the sample is a sinusoidal form symmetrical about the equilibrium position.
Fig. 4. Experimental results of time history response and FFT of the tip vertical deflection when the amplitude setpoint ${A}_{sp}=$0.889${A}_{0}$
a) Timehistory response of the tip vertical deflection
b) FFT of the steadystate response
2.3. Influence of background dissipation on phase contrast in TMAFM
In TMAFM, the phase image is performed by monitoring the phase shift between the tip response and excitation. The phase reflects the energy dissipation of the microcantilever system. During the operation of TMAFM, there are various energy dissipation, such as air viscous dissipation, squeeze film dissipation, liquid bridge dissipation, and so on [18].
Fig. 5. Simulation results of phase and phase contrast of the twocomponent sample under different background dissipation
Such dissipation that cannot reflect the sample information is called background dissipation, and the background dissipation of the system can be expressed by the quality factor ${Q}_{b}$. When the tip works in two areas with different surface properties, the phase shifts ${\phi}_{1}$ and ${\phi}_{2}$ of the two areas are different, and the difference between the twophase shifts is the phase contrast.
As shown in Fig. 5, the sample is composed of two components with different properties. The viscosity ${\eta}_{1}$ of surface І is 800 Pa·s, and the viscosity ${\eta}_{2}$ of surface Ⅱ is 200 Pa·s. The phase shifts of the two surfaces and the phase contrast results between the two surfaces under different background dissipation are obtained by simulation. It can be seen from the figure that the phase contrast $\mathrm{\Delta}\phi $ between the two surfaces increases with the increase of the quality factor representing background dissipation. The results show that reducing the background dissipation of the system can improve the phase contrast between the components of the sample.
3. Conclusions
In this work, the displacement response of the tip in tapping mode is obtained by simulation and experiment. The results show that when one side is constrained by the sample, the steadystate time history response of the tip vertical displacement is still a sinusoidal form symmetrical about the equilibrium position, which is inconsistent with our perceptual understanding. In addition, the phase response of the tip on two surfaces with different properties under different background dissipation is simulated. The results show that reducing the background dissipation can improve the phase contrast between the components of the sample, which has hardly been mentioned by previous studies. The results of this study are of significance to understanding the response of intermittent contact in micronanoscale and provide a reference for the imaging of TMAFM.
References
 G. Binnig, C. F. Quate, and C. Gerber, “Atomic Force Microscope,” Physical Review Letters, Vol. 56, No. 9, pp. 930–933, Mar. 1986, https://doi.org/10.1103/physrevlett.56.930 [Publisher]
 R. Garcia, “Nanomechanical mapping of soft materials with the atomic force microscope: methods, theory and applications,” Chemical Society Reviews, Vol. 49, No. 16, pp. 5850–5884, Aug. 2020, https://doi.org/10.1039/d0cs00318b [Publisher]
 W. Xiang, Y. Tian, and X. Liu, “Dynamic analysis of tapping mode atomic force microscope (AFM) for critical dimension measurement,” Precision Engineering, Vol. 64, pp. 269–279, Jul. 2020, https://doi.org/10.1016/j.precisioneng.2020.03.023 [Publisher]
 Y. Gan, “Atomic and subnanometer resolution in ambient conditions by atomic force microscopy,” Elsevier BV, Surface Science Reports, Mar. 2009. [Publisher]
 Y. F. Dufrêne et al., “Imaging modes of atomic force microscopy for application in molecular and cell biology,” Nature Nanotechnology, Vol. 12, No. 4, pp. 295–307, Apr. 2017, https://doi.org/10.1038/nnano.2017.45 [Publisher]
 J. Melcher et al., “Origins of phase contrast in the atomic force microscope in liquids,” Proceedings of the National Academy of Sciences, Vol. 106, No. 33, pp. 13655–13660, Aug. 2009, https://doi.org/10.1073/pnas.0902240106 [Publisher]
 D. Wang and T. P. Russell, “Advances in Atomic Force Microscopy for Probing Polymer Structure and Properties,” Macromolecules, Vol. 51, No. 1, pp. 3–24, Jan. 2018, https://doi.org/10.1021/acs.macromol.7b01459 [Publisher]
 R. García, “Dynamic atomic force microscopy methods,” Elsevier BV, Surface Science Reports, Sep. 2002. [Publisher]
 A. F. Payam, “Dynamic modeling and sensitivity analysis of dAFM in the transient and steady state motions,” Ultramicroscopy, Vol. 169, pp. 55–61, Oct. 2016, https://doi.org/10.1016/j.ultramic.2016.05.011 [Publisher]
 J. P. Cleveland, B. Anczykowski, A. E. Schmid, and V. B. Elings, “Energy dissipation in tappingmode atomic force microscopy,” Applied Physics Letters, Vol. 72, No. 20, pp. 2613–2615, May 1998, https://doi.org/10.1063/1.121434 [Publisher]
 M. Stark, C. Möller, D. J. Müller, and R. Guckenberger, “From Images to Interactions: HighResolution Phase Imaging in TappingMode Atomic Force Microscopy,” Biophysical Journal, Vol. 80, No. 6, pp. 3009–3018, Jun. 2001, https://doi.org/10.1016/s00063495(01)762662 [Publisher]
 A. Gil, J. Colchero, M. Luna, J. GómezHerrero, and A. M. Baró, “Adsorption of Water on Solid Surfaces Studied by Scanning Force Microscopy,” Langmuir, Vol. 16, No. 11, pp. 5086–5092, May 2000, https://doi.org/10.1021/la9913304 [Publisher]
 J. Tamayo and R. Garcı́a, “Effects of elastic and inelastic interactions on phase contrast images in tappingmode scanning force microscopy,” Applied Physics Letters, Vol. 71, No. 16, pp. 2394–2396, Oct. 1997, https://doi.org/10.1063/1.120039 [Publisher]
 J. Tamayo, “Energy dissipation in tappingmode scanning force microscopy with low quality factors,” Applied Physics Letters, Vol. 75, No. 22, pp. 3569–3571, Nov. 1999, https://doi.org/10.1063/1.125391 [Publisher]
 B. Vasić, A. Matković, and R. Gajić, “Phase imaging and nanoscale energy dissipation of supported graphene using amplitude modulation atomic force microscopy,” Nanotechnology, Vol. 28, No. 46, p. 465708, Nov. 2017, https://doi.org/10.1088/13616528/aa8e3b [Publisher]
 B. V. Derjaguin, V. M. Muller, and Y. P. Toporov, “Effect of contact deformations on the adhesion of particles,” Journal of Colloid and Interface Science, Vol. 53, No. 2, pp. 314–326, Nov. 1975, https://doi.org/10.1016/00219797(75)900181 [Publisher]
 M. Abbasi, “A simulation of atomic force microscope microcantilever in the tapping mode utilizing couple stress theory,” Micron, Vol. 107, pp. 20–27, Apr. 2018, https://doi.org/10.1016/j.micron.2018.01.008 [Publisher]
 M. Imboden and P. Mohanty, “Dissipation in nanoelectromechanical systems,” Elsevier BV, Physics Reports, Jan. 2014. [Publisher]