多模态原子力显微镜的空气压膜阻尼效应

多模态原子力显微镜的空气压膜阻尼效应 Effect of air squeeze film damping in multi-mode atomic force microscopy 1 corresp first-author E-mail： zhaoyang@ahjzu.edu.cn E-mail： zhaoyang@ahjzu.edu.cn 赵 阳（1981-），女，安徽霍邱人，副教授，2003、2006和2018年于合肥工业大学分别获得学士、硕士和博士学位，主要研究方向为纳米扫描探针技术及应用、微纳米测量技术等。E-mail： zhaoyang@ahjzu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037040&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037047&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037045&type= 22.01332855 32.00399780 赵 阳 （1981-），女，安徽霍邱人，副教授，2003、2006和2018年于合肥工业大学分别获得学士、硕士和博士学位，主要研究方向为纳米扫描探针技术及应用、微纳米测量技术等。E-mail： zhaoyang@ahjzu.edu.cn zhaoyang@ahjzu.edu.cn 2 黄强先（1968-），男，教授，博士生导师，1990、1993和1998年于合肥工业大学分别获得学士、硕士和博士学位，主要研究方向为微纳米三维测量技术、纳米扫描探针技术研究及应用、仪器精度理论研究及应用等。E-mail： huangqx@hfut.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037051&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037057&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037053&type= 22.01332855 32.00400162 黄强先 （1968-），男，教授，博士生导师，1990、1993和1998年于合肥工业大学分别获得学士、硕士和博士学位，主要研究方向为微纳米三维测量技术、纳米扫描探针技术研究及应用、仪器精度理论研究及应用等。E-mail： huangqx@hfut.edu.cn huangqx@hfut.edu.cn 安徽建筑大学 电子与信息工程学院，安徽 合肥 230601 安徽建筑大学 电子与信息工程学院，安徽 合肥 230601 School of Electronic and Information Engineering， Anhui Jianzhu University， Hefei 230601， China School of Electronic and Information Engineering， Anhui Jianzhu University， Hefei 230601， China 合肥工业大学 仪器科学与光电工程学院，安徽 合肥 230009 合肥工业大学 仪器科学与光电工程学院，安徽 合肥 230009 School of Instrument Science and Opto-electronics Engineering， Hefei University of Technology， Hefei 230009， China School of Instrument Science and Opto-electronics Engineering， Hefei University of Technology， Hefei 230009， China 原子力显微镜处于轻敲模式工作于大气环境时，由于悬臂和被测试样间的距离很小，悬臂的振动会产生相应的压膜阻尼，影响悬臂的动态特性。基于欧拉-伯努利梁方程对悬臂梁进行建模，结合雷诺方程建立悬臂阻尼作用的分析模型，研究悬臂压膜阻尼效应的主要影响因素，分析不同谐振模态下压膜阻尼效应对悬臂阻尼系数、动态原子力显微镜的测量特性的影响，并基于多模态原子力显微镜进行了实验和测试。理论及实验结果表明，空气压膜阻尼对于基础谐振模态原子力显微镜的测量特性有一定的影响，在悬臂试样20~2 μm逼近过程中，悬臂品质因数明显降低，两种悬臂振幅分别减少了7.8%和20.6%，悬臂宽度增大则受到的压膜阻尼影响也增加；相同实验条件下，二阶谐振模态下的悬臂品质因数和振幅均没有明显变化。在基础谐振模态下压膜阻尼效应会引起微悬臂品质因数的降低，从而导致原子力显微镜测量分辨率的减小和测量速度的降低；多模态原子力显微镜的高阶谐振模态可显著降低空气压膜阻尼效应对悬臂品质因数及测量特性的影响。 When tapping mode atomic force microscopy （AFM） is operated in air， the distance between the cantilever and the sample is so small that cantilever oscillation induces squeeze film air damping in the gap. This study aims to model and analyze the effects of squeeze film damping by using the Euler-Bernoulli and Reynolds equations. Theoretical analysis shows that squeeze film damping can decrease the damping coefficient depending on the distance between the cantilever and the sample and the cantilever width. Then， resonant amplitude and quality factor （ Q ） of the cantilever in multi-mode AFM are tested. Experimental results show that the Q of the fundamental resonant cantilever decreases obviously. The amplitudes of two cantilevers with widths of 20-2 μm decrease by 7.8% and 20.6%. Squeeze film damping can affect the cantilever with a large width to a greater extent than the cantilever with a small width. This phenomenon can also influence AFM dynamic characteristics； in particular， it can decrease flexural sensitivity and scanning speed. However， this phenomenon can only slightly affect the Q and amplitude of the cantilever. The impact of squeeze film damping of the cantilever is significantly reduced when AFM is operated in a higher-order resonant mode. 原子力显微镜 压膜阻尼 品质因数 高阶谐振 atomic force microscopy squeeze film damping quality factor higher-order resonant 1 引 言 1 原子力显微镜（Atomic Force Microscopy， AFM）具有原子尺寸量级的高分辨率，是纳米技术研究的重要工具，被广泛用于纳米尺度成像、测量和操纵等领域 ［ 1 1 - 7 7 1-7 ］ 。AFM微悬臂工作在轻敲模式时，其探针针尖保持轻敲在试样表面，此时微悬臂与试样表面间距离非常小。在微悬臂的周期性振动过程中，它与试样间隙的气体相应地流入或流出，产生一定的压膜阻尼效应。由于微悬臂对力的敏感传递以及阻尼效应对于悬臂动态特性的影响 ［ 8 8 ］ ，空气压膜阻尼效应不能被忽略。空气阻尼对微悬臂的动态特性参数有较大影响，从而会进一步影响悬臂的测量特性。 Hoummady观察到空气压膜阻尼效应的实验现象 ［ 9 9 ］ ：当悬臂下方放置和不放置试样时微悬臂的振幅会有一定的偏移；改变微悬臂试样间的距离，其振幅和谐振频率皆会相应地偏移。Gunther 针对扫描探针显微镜中的音叉臂振动中空气阻尼进行了研究 ［ 10 10 ］ ，将探针末端和被测试样间的阻尼简化为球体和平板间的阻尼模型，分析了振动过程中音叉臂末端的探针和试样间产生的空气阻尼，对于振动过程中音叉臂和试样间的空气阻尼则未做探讨。Leveque，Girard等结合基于连续梁系统的欧拉-伯努利方程和流体力学理论分析了微悬臂探针低频振动时的空气阻尼效应 ［ 11 11 ］ ，进而分析了微悬臂受到的大气阻尼效应对振动幅值的影响，而振动过程中阻尼对于阻尼系数和品质因数的影响则未加分析。Zheng Wei团队进一步分析了不同形状的微悬臂在压膜阻尼影响下品质因数的变化 ［ 12 12 ］ 。 以上研究皆为基础谐振模态下压膜阻尼的影响，未考虑压膜阻尼对于微悬臂高阶谐振模态的影响。Ashok Kumar Pandey仿真了不同阶次谐振模态下微悬臂的空气压膜阻尼系数 ［ 13 13 ］ ，对于同一悬臂的一阶、二阶、三阶谐振模态下的空气压膜阻尼系数进行实验，结果表明，随着谐振阶次的提高，其空气压膜阻尼对于微悬臂的影响相应减小，品质因数相应增加。但是，其测量的品质因数其实是包含微悬臂振动与大气摩擦所引起的阻尼和压膜阻尼之和，且未对压膜阻尼对于微悬臂振动状态的影响做进一步的研究和分析。 本文基于欧拉-伯努利梁方程和雷诺方程对微悬臂的阻尼效应进行分析，研究基础谐振模态和高阶谐振模态微悬臂的压膜阻尼效应对微悬臂动态特性参数、动态AFM测量特性的影响，并围绕系统的阻尼效应进行了实验测试。 2 压膜阻尼效应 1 AFM微悬臂示意图如 图1 图1 所示。矩形悬臂为均质等截面，其长度、宽度和厚度分别为 L ， b 和 W ，密度为 <math id="M1"><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037064&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037061&type= 2.11666679 3.64066648 ，微悬臂的弹性模量为 E ，微悬臂的截面面积为 A （ A = Wb ），截面惯性矩为 I （ I = W 3 b /12）。 由于探针尺寸、质量和微悬臂相比数值微小，故本文忽略探针的作用，将运动的矩形微悬臂简化为矩形平板。忽略惯性效应且不考虑微悬臂在横向的位移，微悬臂和试样间空气压膜效应可用雷诺方程表示为： <math display="block" id="M2"><mtable><mtr><mtd><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><mfenced separators="|"><mrow><mfrac><mrow><mi>p</mi><msup><mrow><mi>h</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mi>μ</mi></mrow></mfrac><mfrac><mrow><mo>∂</mo><mi>p</mi></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mrow></mfenced><mo>+</mo><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac><mfenced separators="|"><mrow><mfrac><mrow><mi>p</mi><msup><mrow><mi>h</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mi>μ</mi></mrow></mfrac><mfrac><mrow><mo>∂</mo><mi>p</mi></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac></mrow></mfenced><mo>=</mo></mtd></mtr><mtr><mtd><mn mathvariant="normal">12</mn><mfenced separators="|"><mrow><mi>h</mi><mfrac><mrow><mo>∂</mo><mfenced separators="|"><mrow><mi>p</mi></mrow></mfenced></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mi>p</mi><mfrac><mrow><mo>∂</mo><mfenced separators="|"><mrow><mi>h</mi></mrow></mfenced></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac></mrow></mfenced><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037078&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037076&type= 46.73600006 23.11400032 其中： x 为悬臂长度方向， y 为微悬臂宽度方向； p 为板间气体压膜压强， μ 为气体黏度系数， h 为臂与试样间变化的距离，即流体膜的厚度。由于微悬臂振动频率远小于压膜阻尼的中介频率 ［ 14 14 ］ ，两板间的气体视为不可压缩气体，同时忽略温度变化，可将 式（1） 式（1） 简化为： <math display="block" id="M3"><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><mfenced separators="|"><mrow><msup><mrow><mi>h</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mfrac><mrow><mo>∂</mo><mi>p</mi></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mrow></mfenced><mo>+</mo><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac><mfenced separators="|"><mrow><msup><mrow><mi>h</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mfrac><mrow><mo>∂</mo><mi>p</mi></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac></mrow></mfenced><mo>=</mo><mn mathvariant="normal">12</mn><mi>μ</mi><mfenced separators="|"><mrow><mfrac><mrow><mo>∂</mo><mi>h</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037085&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037081&type= 55.37200165 9.90600014 . 进一步推导出压膜阻尼系数为： <math display="block" id="M4"><mi>c</mi><mo>=</mo><mfrac><mrow><mi>μ</mi><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><msup><mrow><mi>h</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037091&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037088&type= 11.93799973 8.38199997 . 考虑动态AFM微悬臂在大气环境下振动时，存在着微悬臂振动和大气摩擦产生的能量耗散以及压膜阻尼效应产生的能量耗散。设空气对微悬臂振动导致附加阻尼系数为 γ air ，则分布的外阻尼力为 <math id="M5"><msub><mrow><mi>γ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">r</mi></mrow></msub><mfrac><mrow><mo>∂</mo><mi>z</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037095&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037093&type= 8.63599968 7.95866632 ， f （ x ， t ）为微悬臂受到的 z 向力。微悬臂振动方程为： <math display="block" id="M6"><mtable columnalign="left"><mtr><mtd><mi>E</mi><mi>I</mi><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mi>z</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>t</mi></mrow></mfenced></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow></mfrac><mo>+</mo><mi>ρ</mi><mi>A</mi><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>z</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>t</mi></mrow></mfenced></mrow><mrow><mo>∂</mo><msup><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>γ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">r</mi></mrow></msub><mfrac><mrow><mo>∂</mo><mi>z</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>t</mi></mrow></mfenced></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037099&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037096&type= 48.34466553 18.62666702 其中： <math display="block" id="M7"><mi>z</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>ϕ</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mi>Z</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037103&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037101&type= 31.41133499 3.97933316 ， <math display="block" id="M8"><mtable columnalign="left"><mtr><mtd><mi>ϕ</mi><mfenced separators="|"><mrow><mi>x</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">h</mi><mtext> </mtext><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>x</mi><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>x</mi><mo>-</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>D</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">h</mi><mtext> </mtext><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>x</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>x</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037107&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037105&type= 48.93733215 10.24466610 <math display="block" id="M9"><msub><mrow><mi>D</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mfrac><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>L</mi><mo>+</mo><mi mathvariant="normal">c</mi><mi>h</mi><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>L</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>L</mi><mo>+</mo><mi mathvariant="normal">s</mi><mi>h</mi><msub><mrow><mi>α</mi></mrow><mrow><mi>n</mi></mrow></msub><mi>L</mi></mrow></mfrac><mo>,</mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mi>n</mi><mo>=</mo><mn mathvariant="normal">1</mn><mo>,</mo><mo> </mo><mo> </mo><mn mathvariant="normal">2</mn><mo>,</mo><mo> </mo><mo> </mo><mn mathvariant="normal">3</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037111&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037109&type= 63.92333221 8.46666718 ， 式中： A 0 为微悬臂振动的振幅值，各阶谐振模态的特征参数 α n 可由特征方程求出。 设 Z i （ t ）为AFM微悬臂的第 i 阶振型， ƞ i （ t ）为相应的广义坐标，则利用振型正交性质，并对所得方程进行解耦 ［ 15 15 ］ 可得AFM悬臂第 n 阶振动在广义坐标下的微分方程： <math display="block" id="M10"><msubsup><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><msub><mrow><mover><mrow><mi>η</mi></mrow><mrow><mo>⋅</mo><mo>⋅</mo></mrow></mover></mrow><mrow><mi>n</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>+</mo><mi>a</mi><msubsup><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><msub><mrow><mover><mrow><mi>η</mi></mrow><mrow><mo>⋅</mo></mrow></mover></mrow><mrow><mi>n</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>+</mo><msubsup><mrow><mi>ω</mi></mrow><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msubsup><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><msub><mrow><mi>η</mi></mrow><mrow><mi>n</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>=</mo><msubsup><mrow><mi>F</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037115&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037113&type= 66.80200195 6.60400009 ， <math display="block" id="M11"><msubsup><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><mo>=</mo><msubsup><mo>∫</mo><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi>L</mi></mrow></msubsup><mrow><mi>ρ</mi><mi>A</mi><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mfenced separators="|"><mrow><mi>x</mi></mrow></mfenced></mrow><mi mathvariant="normal">d</mi><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037119&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037117&type= 31.07266808 7.87400007 ， <math display="block" id="M12"><msubsup><mrow><mi>F</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>=</mo><msubsup><mo>∫</mo><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi>L</mi></mrow></msubsup><mrow><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>n</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi></mrow></mfenced><mi>f</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>t</mi></mrow></mfenced></mrow><mi mathvariant="normal">d</mi><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037125&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037123&type= 43.18000031 7.87400007 ， 其中： <math id="M13"><msubsup><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037131&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037127&type= 4.14866638 3.97933316 为第 n 阶广义质量， <math id="M14"><msubsup><mrow><mi>F</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037137&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037134&type= 8.63599968 3.97933316 为第 n 阶广义荷载， ω n 为第 n 阶振动时的角频率，且对方程解耦连时设定 γ air =aρA ， a 为待定常数，可根据试验确定。由 式（9） 式（9） 可根据阻尼比定义求取第 n 阶振动的阻尼比为： <math display="block" id="M15"><msub><mrow><mi>γ</mi></mrow><mrow><mi>n</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>c</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">r</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><msubsup><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><msub><mrow><mi>ω</mi></mrow><mrow><mi>n</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037141&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037139&type= 18.62666702 8.63599968 ， 式中： c air 为微悬臂所受的总阻尼系数，是微悬臂压膜阻尼系数和微悬臂内部能量耗散引起的阻尼系数之和。将一阶谐振模态的 α 1 =1.875，二阶谐振模态的 α 2 =4.694代入 式（9） 式（9） 计算，可得二阶谐振模态和一阶谐振模态微悬臂的阻尼比的关系为： <math display="block" id="M16"><mfrac><mrow><msub><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>=</mo><mn mathvariant="normal">3</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037145&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037143&type= 10.58333302 8.80533314 . 由此可见，和基础谐振模态相比，处于二阶谐振模态的微悬臂探针，其阻尼比相应地下降，从而减小了阻尼对于微悬臂振动状态的影响。 由于微悬臂工作在近谐振点处，其最大动力位移（即振幅）与静力位移之比为动力放大系数 δ ，即： <math display="block" id="M17"><mi>δ</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi>γ</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037147&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=36037146&type= 11.34533310 8.72066593 ， 其中 γ 为微悬臂谐振系统的阻尼比。可见微悬臂的振动幅值由阻尼比决定。在微悬臂和试样接近的过程中，二者距离的改变会导致其阻尼比的变化，进而导致微悬臂振幅的变化。因此，可以通过观察悬臂振幅的变化来观察压膜阻尼的变化情况。 3 实验结果及分析 1 3.1　实验系统 2 基于AFM工作原理，自制了多模态原子力显微镜 ［ 16 16 ］ 。该系统可以工作于接触模式、非接触模式，以及轻敲模式。轻敲模式下，系统可工作于基础谐振模态或者高阶谐振模态下。系统示意图如 图2 图2 所示。采用具有较高谐振频率的压电陶瓷片励振微悬臂，微悬臂变形检测系统采用光杠杆法检测其 Z 向变形信息。 探针-样品逼近系统由远距离逼近系统和近距离逼近系统两部分组成，位于多模态原子力显微镜系统的基座部分。测头及基座部分如 图3 图3 所示。远距离逼近系统通过精密微型直流电机驱动螺旋测微头旋转实现探针和试样的大范围、远距离逼近。逼近系统采用千分尺头，可保证逼近过程的平稳性；逼近选用微型直流电机、减速箱减速方式。直流电机通过联轴器与千分尺头相连接，带动千分尺头旋转，沿着导杆的方向上下运动，实现探针针尖与样品台 Z 方向距离的精密调节。 近距离逼近系统由封装在基座内部工作台下方的三维压电陶瓷扫描器组成，采用两段单管型压电扫描器控制三维位移量。控制其中一个未四分的压电陶瓷管，分别在其内、外壁施加电压，产生轴向伸缩，从而提供 Z 向位移，实现探针-试样间距离的小范围、近距离逼近。 3.2　压膜阻尼对悬臂品质因数的影响 2 选择Budget Sensor的ContAl探针， 使用公称弹性常数为3 N/m的微悬臂，微悬臂的几何尺寸（长，宽，高）为225 μm×28 μm×3 μm。激励悬臂工作于一阶谐振模态及二阶谐振模态，一阶谐振频率为82.30 kHz，二阶谐振频率为522.40kHz。保持其他条件不变，仅改变微悬臂-试样间距离，每次距离改变后进行扫频测试，通过幅频特性曲线中的谐振中心频率及半功率带宽计算品质因数。 悬臂与试样间距离控制的具体实现过程为：先给 Z 向压电陶瓷外壁供120 V电压，其伸长量约为2 μm；然后，驱动电机转动带动试样工作台逼近悬臂，直至探针近似轻敲至试样表面，此时停止电机，撤去 Z 向压电陶瓷电压，则探针-试样间距离为2 μm。以此为起始位置，通过给样品台中的 Z 向压电陶瓷的内外壁施加不同电压使其产生一定的变形量，从而改变微悬臂和试样表面间的距离。 品质因数测量数据如 表1 表1 所示。整个测试过程，探针-试样逐渐远离，使用电感测微仪测量其间距离约从2 μm逐渐增加至10 μm。当探针试样间距离减小时，压膜阻尼效应会导致一阶谐振微悬臂品质因数的明显降低，且随着距离的减小，其压膜阻尼相应增加，品质因数总减小量达到未受压膜阻尼时品质因数的18%。而二阶谐振微悬臂在探针试样间距离减小的过程中，品质因数没有明显的规律性变化。 Serial number Voltage of theinner wall/V Voltage of theouter wall/V Q （the first-order resonant mode） Q （the second-order resonant mode） Distance between the tip and the sample/μm 1 0 120 169.62 537.71 2.00 2 0 80 170.02 537.87 2.14 3 0 40 170.36 537.79 2.48 4 0 0 176.34 537.81 3.03 5 0 -40 179.92 537.79 3.44 6 0 -80 181.13 537.62 3.86 7 0 -120 184.02 537.80 4.42 8 0 -160 185.75 537.54 5.20 9 0 -198 189.23 537.12 5.93 10 40 -198 190.17 536.58 6.74 11 80 -198 192.43 536.87 7.57 12 120 -198 194.77 537.56 8.40 13 160 -198 195.27 537.27 9.05 14 198 -198 197.17 536.92 10.0 15 2 下方无试样 206.82 537.88 ≫20 3.3　压膜阻尼效应对AFM测量特性的影响 2 3 . 3 . 1　基础谐振模态下压膜阻尼效应 3 在微悬臂与样品间距离为20~2 μm区间进行了逼近并连续测试悬臂振幅。采用公称弹性常数为0.2 N/m的探针，几何尺寸（长，宽，高）为450 μm×50 μm×2 μm。使用相同的光强，分别驱动微悬臂工作于接触模式和一阶谐振模态的轻敲模式，测试微悬臂在此范围内向试样表面逼近过程中的力曲线。 图4 图4 （a）与 图4 图4 （b）分别为接触模式和轻敲模式下得到的测试曲线。 1,2; 在接触模式下，微悬臂无振动，向试样表面靠近，在逼近过程中其偏转量并未发生变化，说明在逼近过程中没有受到外力的作用。而动态轻敲模式下，微悬臂处于近谐振状态，振动频率为15.43 kHz。当微悬臂和试样表面距离较远（大于15 μm）时，逼近过程中微悬臂振动幅值基本不变；当微悬臂和试样表面间距离进一步减小时，微悬臂振动幅值随着微悬臂和试样表面间距离的减小而显著减小。整个逼近过程中振幅减小量约为180 mV，微悬臂自由振幅为870 mV，减小量达到微悬臂自由振幅的20.6%。 由于范德华力的作用区域在1 μm以内，实验中的微悬臂与试样间距离是超出范德华力作用范围的，接触模式的测试反映了逼近过程中没有其他外力的作用，轻敲模式下逼近过程中微悬臂受到的力应为由于微悬臂振动引起的压膜阻尼力。 实验还使用公称弹性常数为3 N/m的探针，在相同的实验条件下进行了压膜阻尼特性测试，测试结果如 图5 图5 （a）所示。 1;2; 从 图5 图5 （a）可以看到，其逼近力曲线的变化规律和0.2 N/m微悬臂相同，微悬臂自由振幅为3 250 mV，振幅减小量约为250 mV，约占微悬臂自由振幅的7.8%。进一步进行近距离的力曲线测试，结果如 图5 图5 （b）所示，可以看到，在微悬臂受到范德华力作用而振幅快速减小之前，随着微悬臂和试样间距离的减小，微悬臂振幅持续减小，减小量约为200 mV。 由实验结果可见，微悬臂处于基础谐振模态下，当针尖距离试样表面较近时（15 μm以内），微悬臂的振动会产生压膜阻尼效应，随着微悬臂试样间距离的减小，微悬臂的总阻尼系数和阻尼比增加，使微悬臂的动力放大系数减小，进而导致微悬臂振幅相应减小。同时，宽度较大的微悬臂受到的压膜阻尼影响更大，当微悬臂-试样间距离从20 μm减小到2 μm，宽度为28 μm的3 N/m微悬臂在接近过程中振幅减小了约7.8%，而宽度为50 μm的0.2 N/m微悬臂在接近过程中振幅减小了约20.6%。微悬臂与试样间距离和微悬臂宽度是影响压膜阻尼的主要因素。 3 . 3 . 2　高阶谐振模态下压膜阻尼效应 3 仍使用公称弹性常数为3 N/m的探针，驱动微悬臂工作于二阶谐振模态，其他实验条件保持和一阶谐振模态实验相同，测试微悬臂逼近力曲线，测试结果如 图6 图6 所示。 由 图6 图6 可知，二阶谐振模态下的微悬臂在20 μm范围内逼近过程中其幅值并无发生明显的变化，直至它受到范德华力作用而快速减小，说明在逼近过程中几乎没有受到外力的作用。高阶振动模式时，可以克服空气阻尼作用和流体力阻尼作用的交迭，由微悬臂振动过程中引起的空气阻尼作用占主导，而流体力的阻尼作用可以忽略。由此可见，相较于基础谐振模态，高阶谐振模态更有利于轻敲式AFM减小阻尼对于测试的影响。 AFM在工作于基础谐振模态下扫描试样形貌测试的过程中，微悬臂和试样间的距离非常小，附加产生的空气压膜阻尼效应则会引起微悬臂所受总阻尼的增加，进而引起微悬臂品质因数的减小。对于AFM而言，其微悬臂品质因数的减小一方面会导致悬臂最小可探测力梯度的增加，AFM测量分辨率的降低；另一方面会引起悬臂响应时间的增加，AFM扫描速度的减小。由此可见，压膜阻尼效应的产生会直接影响悬臂动态特性及AFM系统的测量特性。 4 结 论 1 本文对动态AFM微悬臂的空气压膜阻尼特性进行了研究，分析了影响空气压膜阻尼的因素以及空气压膜阻尼对于AFM系统测量特性的影响。基于多模态原子力显微镜系统，进行了幅值反馈模式下的基础和高阶谐振模态微悬臂及AFM系统测量特性的实验研究。测试结果表明，当微悬臂下方有试样且距离较近（10 μm以内）时，基础谐振模态微悬臂品质因数会由于压膜阻尼的影响而明显降低，且随着微悬臂和试样间距离的减小，其压膜阻尼相应增加，品质因数总减小量达到未受压膜阻尼时品质因数的18%。在微悬臂振幅测试中，所采用的两种微悬臂振幅均有明显减小。当微悬臂-试样间距离从20 μm减小到2 μm，宽度为28 μm的3 N/m微悬臂在接近过程中振幅减小了约7.8%，而宽度为50 μm的0.2 N/m微悬臂在接近过程中振幅减小了约20.6%。可见，微悬臂试样间距离和悬臂宽度是影响压膜阻尼的主要因素。压膜阻尼引起了微悬臂品质因数和振幅的明显下降，并进一步影响悬臂动态特性及AFM系统测量特性。 在相同的实验条件下，高阶谐振模态微悬臂品质因数及振幅没有发生明显的变化，说明采用高阶谐振模态可以有效地减少压膜阻尼对系统测量特性的影响。本文的研究为进一步研究动态AFM的阻尼效应提供了理论和实验基础。 参考文献 HEATH G R ， KOTS E ， ROBERTSON J L ， et al . Localization atomic force microscopy ［J］. Nature ， 2021 ， 594 （ 7863 ）： 385 - 390 . doi: 10.1038/s41586-021-03551-x http://dx.doi.org/10.1038/s41586-021-03551-x Localization atomic force microscopy SADEGH HASSANI S ， DARAEE M ， SOBAT Z . Application of atomic force microscopy in adhesion force measurements ［J］. Journal of Adhesion Science and Technology ， 2021 ， 35 （ 3 ）： 221 - 241 . doi: 10.1080/01694243.2020.1798647 http://dx.doi.org/10.1080/01694243.2020.1798647 Application of atomic force microscopy in adhesion force measurements YAN Y D ， HE Y ， GENG Y Q ， et al . Characterization study on machining PMMA thin-film using AFM tip-based dynamic plowing lithography ［J］. Scanning ， 2016 ， 38 （ 6 ）： 612 - 618 . doi: 10.1002/sca.21308 http://dx.doi.org/10.1002/sca.21308 Characterization study on machining PMMA thin-film using AFM tip-based dynamic plowing lithography 陈代谢 ， 殷伯华 ， 林云生 ， 等 . 大范围高速原子力显微镜的前馈反馈混合控制方法 ［J］. 光学 精密工程 ， 2011 ， 19 （ 4 ）： 836 - 843 . doi: 10.3788/ope.20111904.0836 http://dx.doi.org/10.3788/ope.20111904.0836 大范围高速原子力显微镜的前馈反馈混合控制方法 CHEN D X ， YIN B H ， LIN Y SH ， et al . Feed-forward and feed-back controller for large-range and high-speed AFM ［J］. Opt. Precision Eng. ， 2011 ， 19 （ 4 ）： 836 - 843 . （in Chinese） . doi: 10.3788/ope.20111904.0836 http://dx.doi.org/10.3788/ope.20111904.0836 Feed-forward and feed-back controller for large-range and high-speed AFM WASTL D S ， WEYMOUTH A J ， GIESSIBL F J . Atomically resolved graphitic surfaces in air by atomic force microscopy ［J］. ACS Nano ， 2014 ， 8 （ 5 ）： 5233 - 5239 . doi: 10.1021/nn501696q http://dx.doi.org/10.1021/nn501696q Atomically resolved graphitic surfaces in air by atomic force microscopy 吴俊杰 ， 刘俭 ， 魏佳斯 ， 等 . 双测头复合型微纳米测量仪的研制 ［J］. 光学 精密工程 ， 2020 ， 28 （ 2 ）： 415 - 423 . 双测头复合型微纳米测量仪的研制 WU J J ， LIU J ， WEI J S ， et al . Development of double probe composite micro-and nano measuring instrument ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 2 ）： 415 - 423 . （in Chinese） Development of double probe composite micro-and nano measuring instrument BALANTEKIN M . High-speed dynamic atomic force microscopy by using a Q -controlled cantilever eigenmode as an actuator ［J］. Ultramicroscopy ， 2015 ， 149 ： 45 - 50 . doi: 10.1016/j.ultramic.2014.11.016 http://dx.doi.org/10.1016/j.ultramic.2014.11.016 High-speed dynamic atomic force microscopy by using a Q-controlled cantilever eigenmode as an actuator CASTRO GARCÍA R . Amplitude Modulation atomic Force Microscopy ［M］. Weinheim ： Wiley-VCH ， 2010 . doi: 10.1002/9783527632183 http://dx.doi.org/10.1002/9783527632183 Amplitude Modulation atomic Force Microscopy HOUMMADY M ， FARNAULT E . Enhanced sensitivity to force gradients by using higher flexural modes of the atomic force microscope cantilever ［J］. Applied Physics A ， 1998 ， 66 （ 1 ）： S361 - S364 . doi: 10.1007/s003390051164 http://dx.doi.org/10.1007/s003390051164 Enhanced sensitivity to force gradients by using higher flexural modes of the atomic force microscope cantilever GÜNTHER P ， FISCHER U C ， DRANSFELD K . Scanning near-field acoustic microscopy ［J］. Applied Physics B ， 1989 ， 48 （ 1 ）： 89 - 92 . doi: 10.1007/bf00694423 http://dx.doi.org/10.1007/bf00694423 Scanning near-field acoustic microscopy LÉVÊQUE G ， GIRARD P ， BELAIDI S ， et al . Effects of air damping in noncontact resonant force microscopy ［J］. Review of Scientific Instruments ， 1997 ， 68 （ 11 ）： 4137 - 4144 . doi: 10.1063/1.1148358 http://dx.doi.org/10.1063/1.1148358 Effects of air damping in noncontact resonant force microscopy WEI Z ， LIU J ， ZHENG X T ， et al . Influence of squeeze film damping on quality factor in tapping mode atomic force microscope ［J］. Journal of Sound and Vibration ， 2021 ， 491 ： 115720 . doi: 10.1016/j.jsv.2020.115720 http://dx.doi.org/10.1016/j.jsv.2020.115720 Influence of squeeze film damping on quality factor in tapping mode atomic force microscope PANDEY A K ， PRATAP R . Effect of flexural modes on squeeze film damping in MEMS cantilever resonators ［J］. Journal of Micromechanics and Microengineering ， 2007 ， 17 （ 12 ）： 2475 - 2484 . doi: 10.1088/0960-1317/17/12/013 http://dx.doi.org/10.1088/0960-1317/17/12/013 Effect of flexural modes on squeeze film damping in MEMS cantilever resonators BAO M H ， YANG H . Squeeze film air damping in MEMS ［J］. Sensors and Actuators A： Physical ， 2007 ， 136 （ 1 ）： 3 - 27 . doi: 10.1016/j.sna.2007.01.008 http://dx.doi.org/10.1016/j.sna.2007.01.008 Squeeze film air damping in MEMS R . 克拉夫， J . 彭津 . 结构振动力学［M］. 第二版 . 北京 ： 高教出版社 ， 2006 . doi: 10.15385/yb.miracle.2006 http://dx.doi.org/10.15385/yb.miracle.2006 彭津 CLOUGH R ， PENZIEN J . Dynamics of Structures ［M］. 2nd Edition . Beijing ： Higher Education Press ， 2006 . （in Chinese） . doi: 10.15385/yb.miracle.2006 http://dx.doi.org/10.15385/yb.miracle.2006 Dynamics of Structures 黄强先 ， 张蕤 ， 刘凯 ， 等 . 多模态动态原子力显微镜系统 ［J］. 光学 精密工程 ， 2017 ， 25 （ 2 ）： 401 - 407 . doi: 10.3788/ope.20172502.0401 http://dx.doi.org/10.3788/ope.20172502.0401 多模态动态原子力显微镜系统 HUANG Q X ， ZHANG R ， LIU K ， et al . Multi-mode dynamic atomic force microscope system ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 2 ）： 401 - 407 . （in Chinese） . doi: 10.3788/ope.20172502.0401 http://dx.doi.org/10.3788/ope.20172502.0401 Multi-mode dynamic atomic force microscope system 10.37188/OPE.20223019.2362 TH89 A 1004-924X（2022）19-2362-08 赵阳,黄强先.多模态原子力显微镜的空气压膜阻尼效应[J].光学精密工程,2022,30(19):2362-2369. 赵阳,黄强先.多模态原子力显微镜的空气压膜阻尼效应[J].光学精密工程,2022,30(19):2362-2369. ZHAO Yang,HUANG Qiangxian.Effect of air squeeze film damping in multi-mode atomic force microscopy[J].Optics and Precision Engineering,2022,30(19):2362-2369. ZHAO Yang,HUANG Qiangxian.Effect of air squeeze film damping in multi-mode atomic force microscopy[J].Optics and Precision Engineering,2022,30(19):2362-2369. 安徽省高校自然科学基金资助项目（No.KJ2019A0765）；国家自然科学基金资助项目（No.50975075） 2022-03-16 2022-05-30 2022-10-10 2022-10-24 光学精密工程 19 30