范德瓦尔斯方程中有关分子吸引力引起的修正问题

范德瓦尔斯方程中有关分子吸引力引起的修正问题专业：电子信息与物理系姓名：郑凯学号：08405142指导教师：张栓柱摘要：本文用统计的方法导出范德瓦尔斯方程及修正值的具体表达式, 并讨论修正量与状态量的关系，对范德瓦尔斯方程的提出背景及其内涵作一系统介绍关键词：范德瓦尔斯方程统计解释历史分子力1873年以前, 人们一直沿用理想气体的状态方程PV=νRT, 此方程适用于低压、高温的条件下; 可是由实际气体的等温线可知, 在低温高压下, 实际气体与理想气体有很大不同。

因此, 需要对理想气体的状态方程作一定的修正, 1873年范德瓦尔斯对理想气体的两条假设作改进, 得出了能描述真实气体行为的范德瓦尔斯方程。

这一方程的建立, 对气体动理论与实际气体的状态更加的接近了, 范德瓦尔斯的贡献及历史沿革可以由以下几个方面论述。

1、气体动理论的发展范德瓦尔斯方程是1873年提出的，它是分气体动理论中最重要的一叶。

在此之前，气体动理论经历了大致三个阶段。

第一阶段是在17世纪。

意大利托里拆利研究大气压强，获得了托里拆利真空。

法国的帕斯卡研究了流体压强，建立了液体压强的有关定律。

在这基础上形成了朴素的热动力学概念。

1658年，波意耳通过实验，取得空气体积和压强成反比的定量关系。

波意耳认为：空气是由大量的小颗粒组成的，浮游在涡旋运动的宇宙物质中，所以这些粒子不停地杂乱无章地云顶。

这是最初的气体动理论，用它能解释压强是如何形成的。

第二阶段是187世纪到19世纪上半叶。

1738年，伯努利提出“弹子球”模型：把空气看成是“弹性流体”，一颗颗的“弹子球”与器壁碰撞, 由此可导出气体定律。

这个观点是气体动理论的基础, 也正是后来被称之为“理想气体”的模型。

第三阶段是在19世纪中叶后。

焦耳热功当量实验和空气自由膨胀实验, 推动了热动力学的发展, 使气体动理论的研究进入了一个新阶段。

作为热动力学奠基人之一的德国物理学家克劳休修, 第一个发表论文, 阐述了分子运动理论, 他还提出了气体分子平均自由程的概念, 创立了作为后来导出气体状态方程基础的“维里理论”。

这个时期, 英国的麦克思韦将统计方法引入到分子理论的研究中, 提出了分子速度的分布律。

之后, 玻尔兹曼又将重力场引入分子运动速度分布律中, 并予以熵以统计的意义。

至此, 气体动理论终于由一些定性的论据研究发展成为一门有严密体系的科学。

由于有了新的事实根据, 原来那种将气体看做是大小可以忽略不计的弹性模型已不能解释新的实验结果。

研究真实气体的性质, 建立相应的理论, 成为十九世纪后半叶物理学家的一个重要课题之一。

2、范德瓦尔斯的贡献范德瓦尔斯于于1837年出生于荷兰的一个木匠家庭, 由于家境贫寒, 上完小学后就没能再上中学；辍学在家的范德瓦尔斯凭着对科学的热爱, 靠自己的勤奋努力, 自学成材, 最后终于考入莱顿大学, 并在1873年年获博士学位。

他的博士论文就是《论气态和液态的连续性》，这篇论文考虑了气体的分子体积和分子间的吸引力的影响，推出了著名的物态方程: ( P+ a/V) (V-b) =RT。

后来被人们称为范德瓦尔斯方程。

该方程不仅能解释很多实验结果, 而且能以常数a、b值计算临界状态参数。

1910年, 范德瓦尔斯由于他的气态和液态方程获诺贝尔物理学奖。

对近代物理学中的气体动理论的发展做出了杰出的贡献。

3、推导范德瓦尔斯方程3.1范德瓦尔斯气体理想气体是一个近似模型, 它忽略了分子的体积和分子间的引力。

克劳修斯和范德瓦尔把气体看做有相互吸引力作用的刚性球即具有一定的体积且相互之间存在引力的作用。

这种模型与实际气体更接近，我们常把这种气体称作范德瓦尔斯气。

3.2刚性球分子引起的体积修正理想气体模型中把分子看成是没有体积的质点，所以每个分子可以自由活动的空间就是容器的容积v 。

当把它看做是有一定体积的刚性球分子时，每个分子能够活动的空间显然不再等于容器的容积v, 而应从v 中减去一个反映气体分子所占有体积的修正量b 。

设想在气体内除某一分子 在运动的过程中不断与它们相碰，若用d 表示分子的有效直径，当分子 与任一分子 相碰时，它们中心间的距离就是d ，如图1所示。

当分子 沿任意方向趋近其它任一分子时， 的中心将被排除于直径为2d0的球形区域外。

即分子背向球形区域的那一半不进入球形区域，如图2所示。

所以，修正量为（1）其中()d v 03430π=为一个分子的体积3.3分子间相互作用引起压强的修正在气体内部，由于每个分子受到各向的吸引，平均效果相互抵消。

但当分子靠近器壁时，一边受气体分子的吸引，而近壁一边则没有。

因而使得分子受到一个垂直于器壁指()νππ33403443421)1(N N N b d R=•≈••-=向气体内部的拉力作用。

在靠近器壁处取一厚度为引力有效作用距离s 的区域ABCD ，如图3所示。

假想分子在AB 处就和器壁相碰而被反射，则产生的压强为p bv RTp ∆--=, （2）但实际上，分子必须通过ABCD 区域才能碰到器壁，分子在这区域中受到的向内的拉力f 将使它在垂直于器壁方向上的动量减小。

因而器壁实际受到的压强要比(2)式值小p ∆，即（3） 压强的减小值为()()v p m 2∆⨯=∆相碰的分子数单位时间内与单位壁面 （4） 显然，()n mv ∞∆同时，单位时间内与单位面积相碰的分子数也与n 成正比，所以∞∆p vn221∞（5）p bv RTp ∆--=写成等式，得vap 2=∆ （6）其中a 为比例系数。

3.4范德瓦尔斯方程 把(6) 代入(3), 得()RT b v a p ab v RT p v v=-⎪⎪⎭⎫⎝⎛+--=22即 （7）上式就是1mol 范德瓦尔斯气体的状态方程。

在方程的推导过程中我们对体积和压强作了修正，得到反映压强和体积修正的两个数a 、b 。

这两个数在一般的教材中都认为是与其它因素无关的常数。

4、用统计方法导出范德瓦尔斯方程 4.1配分函数设1mol 的范德瓦尔斯气体由N 个单原子分子组成, 分子的质量为m 任一分子i 的动量为Pi,任一对分子i, j 间的相互作用能为()r ij u 则系统的能量为()∑∑+=r p ij iu miE 22（8）配分函数为()QN ddd d d dr zdydxdz d y d x dp iN dpdq N Z h m zyxz y x ep p p p p p ehh e NNNNu mNNEji ij NNN⎪⎭⎫⎝⎛⎰⎰⎰⎰⎰=∑∑=•=〈-2!1............!1!1223111233_1112βπββ（9）其中()zy xz y x eNNNd dd d d dr Q ij......111⎰⎰-=β（10）引入函数()1-=-efr ijuijβ (11)则(10)可变成()()ττNji j i ijij NNNijji d d d d d d d d Q f f fz y x z y x f......1......1 (1)111⎰∑∑∑⎰⎰∏⎰〈''〈+++=+= (12)其中z y x lll l d d d d =τ，上式作近似处理后可得，(13)把(13)代入( 9) ,得()⎥⎥⎦⎤⎢⎢⎣⎡+=⎰τl Nd r f v Q N v212()⎪⎪⎭⎫⎝⎛+=⎰⎪⎭⎫ ⎝⎛τβπl NNd r f v N Z N vh m 212!12232 (14)4.2范德瓦尔斯方程由热力学与统计物理学公式知,气体的压强可表示为vZp ∂∂=-ln 1β (15) 由(14)得()⎪⎪⎭⎫⎝⎛+++⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎰⎪⎭⎫ ⎝⎛τβπl Nd r f v v N N Z N h m 21ln ln 2!1ln ln 2232 上式最后一项对数展开后取首项,得 (16)(16)两边对v 求偏导, 得()()ττl l d r f v N d r f v v v v N v Z v N N ⎰⎰-=⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=∂∂22222ln ln (17)把(17)代入( 15), 得()⎪⎪⎭⎫⎝⎛-=⎰-τβl d r f v N p v N 2212 (18)用kT1=β 代入并整理,得 ()⎪⎭⎫⎝⎛-=⎰τl d r f v N NkT pv 21 (19)()τβπl Nd r f v v N N Z N h m ⎰⎪⎭⎫ ⎝⎛++⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=2ln 2!1ln ln 2232由于()()()d d r r r u n 00r ≤〉⎪⎩⎪⎨⎧∞-=常数，μμ所以()()()()()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-=⎪⎭⎫ ⎝⎛--=--=-⎰∑⎪⎪⎭⎫ ⎝⎛⎰⎰⎰⎰⎰∞=∞----d n N d dr r dr N dr N d N d r f N d kT d d r e r r e eonkT kT r u lr u l 00133002122033!03022412122γγβμττγγπππ (20) 把(20)代入( 19), 得(21) 即 (22)因为d v 06130π=，所以上式变为()⎪⎪⎭⎫ ⎝⎛∑∞=-⎪⎪⎭⎫ ⎝⎛+=d kT d N v n v kT v N NkT pv 033!02_411320μγγγγπ (23) 由于v N v 〈〈04，所以(24)()⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎪⎪⎪⎭⎫⎝⎛--•+=∑∞=13303!030211γμγγπd d d n kT n N v NkT pv ∑⎪⎪⎭⎫ ⎝⎛∞=-⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+=1230206141γγμππd kT N d n v kT v N NkT pv vN vN v v o411410-≈+把(24)代入( 23)并整理, 得()()NkT N v n kT p v vd kT d N =-⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-+∑⎪⎪⎭⎫ ⎝⎛∞=021324033!102γγμγγπ (25)上式就是用统计方法导出的范德瓦尔斯方程。

5、对范德瓦尔斯方程的讨论 5.1体积与压强的修正把(25) 与( 7) 比较可知,由于分子的固有体积引起的修正b 与由于分子间相互作用引起的修正a 为()⎪⎪⎭⎫ ⎝⎛∑∞=-=d kT d N n kT a 033!102132μγγγγπ （26）v o N b 4=5.2压强的修正与温度的关系5.2.1一般情况下从（26）式我们可以看到，a 是T 的函数，且随温度的升高而减小。

这说明，随着温度的升高，分子间的相互作用减弱,分子势能增大。

但通常认为,势能是由分子间的相对位置决定的, 即分子势能只与系统的体积有关,与温度无关. 实际上,这主要是由于我们没有注意到系统的统计特性造成的。