范德华力的量化

1 范德华力的计算模型

1.1 Lennard-Jones理论

范德华力由三个部分组成：取向力、诱导力和色散力。取向力，即极性分子晶体的结合力，是由于极性分子在永久点偶极矩而产生的。诱导力是诱导偶极矩于极性分子偶极矩之间的作用力。非极性分子之间存在着瞬间、周期变化的偶极矩，这种瞬间偶极矩间的相互作用产生了非极性分子晶体的结合力，即色散力。

一般情况下，对于非极性分子，或惰性元素的单原子分子来说，只存在色散力。而对于极性分子来说，则上述三种力都存在。但是不管哪一种力，引起的吸引能均与分子间距r的6次方成反比，即：

(1) 当两分子靠得很近时，同样也会参数排斥力，排斥能与分子间距r的12次方成反比，即：

(2) 因此，两分子间的吸引力和排斥力产生的能力可表示为：

(3) 而常用的表示方式为：

(5) 式中，

，

式(5)就是Lennard-Jones势能。

忽略高阶项，分子间的相互作用势能可表示为：

(6) 其中，Λ为常数，z为分子间的距离。

如图1，范德华力的总势能表示为：

(7)

图1 两个粒子间的相互作用示意图

其中，n表示单位体积中包含的分子数。上式可进一步表示为：

(8)

式中，是Hamaker。

上式中，假设，如图2，就能得到粒子和平面之间相互作用的势能和范德华力：

(9)

(10) 式中，d=2r。

图2 粒子和平面相互作用示意图。

如果z相对于r足够小，上式可近似为：

(11)

1.2 针尖与样品接触模型的选择

探针针尖与样品接触过程中，会产生变形，弹性变形理论中涉及到表面力，有三种理论模型：JKSR模型、DMT模型、MYD/BHW模型。一般只选择JKSR模型和DMT 模型。JKSR模型适合于大半径、高粘附能和低模量体系。DMT模型适用于曲率半径小、粘附能低、弹性模量高的体系。

DMT模型中，弹性球体和刚性平面的作用力为

(12)

H是Hamaker长沙，d为弹性球体的直径，z0是弹性球体与刚性平面表面刚接触式的分子间距，h是球体和平面接触后的分子间距，K是等效刚度，E和σ是杨氏模量和泊松比。

我们知道探针针尖也是由很多原子组成的。分析探针于样品原子的作用力时，样品是由很多原子组成，那么这些原子肯定都会对探针产生力的作用，如图3。所以，我们可以把针尖和样品的作用简化成探针最尖端的原子和平面的作用。

图3 探针和样品作用示意图

(1) 当针尖远离样品时，

(13)

(2) 当针尖接触到样品时，，根据DMT接触模型

(14)

是针尖和样品的等效刚度，

(15)

，和，分别表示针尖和样品的杨氏模量和泊松比。是针尖和样品接触时分子间的距离。

2 悬臂的模型

悬臂的模型为

m——悬臂质量，ω0——悬臂固有频率，Q——品质因数，K——悬臂刚度系数，z——是针尖在z方向位移，zpie——压电陶瓷的z向输出位移