純物質(zhì)分子受到的分子間引力及其應(yīng)用研究
——利用平均場(chǎng)理論計(jì)算

＊趙仁興

（河北科技大學(xué) 河北 050018）

分子間作用力是指兩個(gè)或多個(gè)分子之間除了離子鍵、共價(jià)鍵和金屬鍵外的其他弱相互作用力，常表現(xiàn)為范德華力、氫鍵等[1]。分子間作用力包括引力和斥力，但在絕大多數(shù)情況下，維系物質(zhì)存在的是和斥力平衡后的引力。物質(zhì)分子受到的引力大小決定了物質(zhì)的熔點(diǎn)、沸點(diǎn)、臨界溫度、汽化熱、熔解熱、密度、表面張力、黏度、擴(kuò)散系數(shù)等物理性質(zhì)，同時(shí)和不同物質(zhì)之間的溶解度、結(jié)晶、吸附等有關(guān)。

從1873年范德瓦爾斯建議引入引力、斥力將理想氣體狀態(tài)方程修正為VDW方程[2]以來(lái)，關(guān)于分子間作用力的研究一直持續(xù)不斷。1912年葛生（Keeson）研究了極性分子（永久偶極矩）之間的相互作用，1921年德拜（Debye）分析了極性分子與非極性分子之間（永久偶極矩和誘導(dǎo)偶極矩）的相互作用，1930年倫敦（London）研究了非極性（瞬間偶極矩與誘導(dǎo)偶極矩）分子之間的相互作用,以后又不斷推出氫鍵、靜電作用、π-π堆積及荷電基團(tuán)、疏水基團(tuán)等的相互作用[1]。我國(guó)學(xué)者金松壽等[3]20世紀(jì)80年代推出了選擇性分子間引力理論。

具有代表性的分子間作用力公式有兩個(gè)極性分子（永久偶極矩）之間吸引的平均作用能公式[4]和分子間總的相互作用勢(shì)能蘭納德—瓊斯公式[5]。

上述研究極大地推進(jìn)了物質(zhì)物性的估算、結(jié)構(gòu)的研究及新材料的制備等，但至今為止物質(zhì)的吸附、色譜分離、溶解度、沸點(diǎn)等仍不能用分子間引力及氫鍵能等直接進(jìn)行量化的計(jì)算，只能作些定性推斷。

量子化學(xué)也將分子間力的研究作為一個(gè)重要目標(biāo)。近二十年來(lái)量子化學(xué)的從頭計(jì)算法有了很大發(fā)展，可用來(lái)預(yù)測(cè)分子的結(jié)構(gòu)和能量，但由于分子之間作用復(fù)雜并形式多樣，限制了其計(jì)算的精度[6]。

我國(guó)利用分子間引力也開展了多方面的研究，其中包括制冷劑性能、汽化熱和沸點(diǎn)的計(jì)算、表面活性劑自組裝調(diào)控、金屬表面小分子吸附、共沸物的研究等[7-10]。

物質(zhì)是眾多分子的聚集，因此除需了解兩兩（例如偶極矩—偶極矩之間）分子間的作用力外，也需要得到分子在物質(zhì)中受到的總引力，并且只有在了解總引力的基礎(chǔ)上，才能更好的了解物質(zhì)本身。在充分了解不同基團(tuán)之間的分子間作用后就可有目標(biāo)的組合不同的基團(tuán)，生產(chǎn)更好的產(chǎn)品，供人類使用。

基于上述考慮，本研究采用和以前研究不同的方法計(jì)算分子間引力。物質(zhì)所表現(xiàn)的宏觀性能，應(yīng)該是分子在物質(zhì)內(nèi)部或表面所受到的綜合作用力決定的。為此該研究從分子熱運(yùn)動(dòng)的平均場(chǎng)出發(fā)，利用汽液平衡時(shí)的汽化熱，在不考慮不同引力種類的前提下，計(jì)算任意一個(gè)單分子在純物質(zhì)中受到其周圍分子的作用力之和，本文給出了常見(jiàn)單質(zhì)氣體純物質(zhì)分子受到的引力及相應(yīng)的位能，同時(shí)從引力出發(fā)給出了物質(zhì)臨界溫度、臨界壓力、表面張力的計(jì)算方法，并將計(jì)算結(jié)果和實(shí)測(cè)值進(jìn)行了比較，具有較高的一致性。本項(xiàng)目分子間引力的研究方法和傳統(tǒng)研究有很大不同，是一種全新的分子間相互作用的研究方法。

1.研究使用的公式和推導(dǎo)過(guò)程

(1)靜電力公式

有關(guān)研究認(rèn)為分子間力為長(zhǎng)程力[11]，作用范圍為0.3～0.5nm，但和水分子偶極矩尺寸6.20×10-23nm比較，在其作用范圍內(nèi)仍包含有大量的其他分子。因此本研究認(rèn)為純物質(zhì)中單個(gè)分子所受到的作用力，是其周圍無(wú)數(shù)分子的共同作用的結(jié)果。在汽、液相變過(guò)程研究范圍內(nèi)，由于單個(gè)分子占據(jù)的名義體積較大，分子間的斥力對(duì)物質(zhì)宏觀性質(zhì)尚不起決定性作用。考慮分子的熱運(yùn)動(dòng)引起的物質(zhì)同一相內(nèi)分子性能的均勻性，按平均場(chǎng)定義純物質(zhì)中單個(gè)分子受到周圍分子引力的總和為：

式中：

F—任一物質(zhì)分子在純物質(zhì)中受到周圍分子的引力標(biāo)量之和（不考慮方向）。

f—和純物質(zhì)分子特性有關(guān)的引力系數(shù)，f也可表示為f=a1×a2，a1為某分子外其他分子對(duì)其能產(chǎn)生引力的綜合因子；a2為某分子自身具有的能和周圍其他分子產(chǎn)生引力的綜合因子。本研究以f代替了諸如永久偶極矩、誘導(dǎo)偶極矩或瞬時(shí)偶極等因子，但應(yīng)該是某分子原子核的正電荷或核外電子與其他分子的原子核電荷、核外電子之間的相互作用系數(shù)，同時(shí)將4πε0也納入了f中，因此f值和以往研究的相應(yīng)數(shù)值無(wú)法進(jìn)行直接比較。

f′—某一對(duì)分子之間的引力系數(shù)。

n—某分子受到引力作用的分子對(duì)數(shù)。

r′—分子引力的作用半徑。以往分子相互作用的研究中，一般定義r′為兩個(gè)分子質(zhì)心間的距離。本文定義分子引力的作用半徑為某分子原子核和相鄰最近分子電子云之間的距離，計(jì)算公式如下：

式中：Vc—臨界點(diǎn)時(shí)的摩爾體積；

r—定義為分子的引力平衡半徑。任何物質(zhì)在壓力、溫度不變的情況下，其密度一般不會(huì)發(fā)生變化，因此可以認(rèn)為在一定壓力、溫度下，分子所占有的名義球體體積的半徑即是分子間引力平衡半徑。

在單一物質(zhì)均相無(wú)序狀態(tài)下，分子引力平衡半徑r的計(jì)算式如下：

式中：V—摩爾體積；N0—阿伏伽德羅常數(shù)。

必須說(shuō)明的是r值會(huì)隨物質(zhì)溫度、壓力而發(fā)生變化。

某物質(zhì)分子原子核與電子云之間的距離按該物質(zhì)臨界容積的0.25倍計(jì)算出的球體半徑。按羅于根《元素電離能與原子的實(shí)際半徑》中的計(jì)算，基態(tài)氫原子電子到核的距離為5.296×10－11m，氦原子為5.293×10－11m，兩者基本一致。按氫、氦臨界容積0.25倍計(jì)算出的氫分子、氦分子電子云半徑分別為4.81×10-10m、4.63×10-10m，約為按電離能算出原子半徑的10倍。

計(jì)算的中心分子和其外圍的第二層、第三層分子等的距離可以為等表示。

本研究為簡(jiǎn)化計(jì)，不考慮周圍密接分子數(shù)及不同距離的分子數(shù)，所有其他因素綜合為f值。必須注意的是，因?yàn)槲覀儾捎玫氖瞧療嵊?jì)算出來(lái)的f值，因此它是周圍分子對(duì)某一分子引力的綜合系數(shù)，是由周圍具有相同距離的分子及多個(gè)倍數(shù)距離分子共同作用的結(jié)果，這種假定的不足之處來(lái)源于：由于不同元素或化合物外層電子的多少，0.25倍Vc發(fā)生的變化；不同分子層引力的距離與單一距離的不同。

(2)汽、液相變過(guò)程中的蒸發(fā)熱和單分子引力系數(shù)的計(jì)算

①汽、液相變過(guò)程中的能量

相變是物質(zhì)在自然界存在的一種最重要的過(guò)程，其中包括汽—液、汽—固、液—固相之間的相變及固體結(jié)晶內(nèi)部不同相的變化。

汽、液相變過(guò)程中摩爾分子吸收的熱量[12]為；

②由汽化熱計(jì)算單分子引力系數(shù)的公式

“汽化熱包括兩個(gè)組成部分，P（Vg-Vm）是液體汽化時(shí)擴(kuò)大體積所需對(duì)外作的等壓功，而Ug,m—UL,m是液體分子變?yōu)橄嗤瑴囟认碌臍怏w分子所需克服周圍分子吸引力作的功”[13]。在汽液平衡時(shí)，氣相和液相分子均已有相同的溫度和壓力，因此在分子由氣相分子凝結(jié)成液相分子或由液相分子汽化為氣相分子時(shí)（按分子運(yùn)動(dòng)論，分子的動(dòng)能和壓力、溫度有關(guān)），不會(huì)有壓力、溫度的變化，分子的動(dòng)能不會(huì)再發(fā)生變化，汽體、液體內(nèi)分子動(dòng)能是相等。因此在汽、液相變化中需要克服分子間引力的作用、需吸收外環(huán)境的熱量或放出熱量。

根據(jù)引力公式可得：

式中：ΔHm、ΔUm、Hgm、HLm、Ugm、ULm、Vgm、VLm—分別是汽化時(shí)摩爾焓、內(nèi)能的變化及液體、氣體的摩爾焓、摩爾內(nèi)能、摩爾體積。P為汽化時(shí)的壓力。公式（7）進(jìn)行單位換算后，可用于實(shí)際計(jì)算引力系數(shù)的公式為：

式中：f—純物質(zhì)中單個(gè)分子的引力系數(shù)（kJ·m）；

M—分子量；

N0—阿伏伽德羅常數(shù)（6.023×1023）；

△HM—汽化熱（kJ/kg）；

P—汽化時(shí)的壓力（MPa）；

Vg—飽和氣體比容（m3/kg）；

(3)單分子引力計(jì)算

根據(jù)引力系數(shù)計(jì)算公式，單個(gè)分子受到的引力計(jì)算公式如下：

式中：Vg′—飽和氣體比容（m3/moL）；VL′—飽和液體比容（m3/moL）；FL、Fg—飽和液體或氣體中單個(gè)分子受到的引力（N）。

(4)氣相、液相摩爾分子間的引力位能計(jì)算公式

根據(jù)位能的加和性能，摩爾分子間的引力位能Eg、EL的計(jì)算公式：

式中：Eg、EL—摩爾分子引力位能（kJ/mol）；

其余符號(hào)同前。

2.常見(jiàn)氣態(tài)元素引力系數(shù)的計(jì)算結(jié)果和討論

(1)常見(jiàn)氣態(tài)元素的物理性質(zhì)

本文計(jì)算用的10種元素?zé)崃W(xué)數(shù)據(jù)除氯外均來(lái)源于參考文獻(xiàn)[13-14]。

(2)常見(jiàn)氣態(tài)元素分子間引力系數(shù)和溫度之間的關(guān)系

依據(jù)式（2）定義的引力計(jì)算半徑和式（8）計(jì)算出的常見(jiàn)氣態(tài)元素單分子引力系數(shù)和溫度之間的關(guān)系見(jiàn)圖1～圖5。

圖1 氦3、氦4引力系數(shù)和溫度關(guān)系曲線

圖2 氫、仲氫、氘引力系數(shù)和溫度關(guān)系

圖3 惰性元素引力系數(shù)和溫度關(guān)系

圖4 氮、氧、氟引力系數(shù)和溫度關(guān)系

圖5 氯元素引力系數(shù)和溫度關(guān)系

上述10種元素的偶極距均為零，說(shuō)明上述元素的分子之間引力主要來(lái)源于色散力。隨溫度上升，不同元素的引力系數(shù)變化各有不同。第一周期氦、氫元素，無(wú)論是氦3或氦4、標(biāo)準(zhǔn)氫或仲氫、氘分子間的引力系數(shù)均隨溫度的上升而直線上升，說(shuō)明第一周期元素色散力會(huì)隨溫度的上升而增加，相對(duì)偏差較大。

其他周期元素的變化相對(duì)平穩(wěn)，熔點(diǎn)到臨界點(diǎn)溫度范圍內(nèi)引力系數(shù)與平均值相對(duì)偏差的平均值為0.34%～1.34%。不同元素的平均引力系數(shù)及相對(duì)偏差如下：氦3（2.23E-35，13.65%）、氦4（4.41E-35，9.14%）、仲氫（4.01E-34，7.81%）、標(biāo)準(zhǔn)氫（4.13E-34，6.92%）、氘（5.475E-34，4.70%），氖（6.33E-34，1.34%）、氟（2.66E-33，0.44%）、氮（2.52E-33，0.61%）、氬（2.71E-33，0.68%）、氧（2.97E-33，0.70%）、氪（4.01E-33，0.70%）、氙（6.07E-33，0.58%）、氯（9.66E-33，0.73%）。引力系數(shù)較大的元素，如氧、氮、氯、氪、氙等在臨界點(diǎn)附近，引力系數(shù)有明顯下降。分子引力系數(shù)越大的元素，引力系數(shù)下降越明顯，可能是靠近臨界點(diǎn)時(shí)氣相中分子的聚集[15]行為所造成的。本文認(rèn)為通過(guò)一個(gè)公式能將元素不同氣、液相飽和溫度下4種數(shù)值相差很大的數(shù)據(jù)，計(jì)算出一個(gè)表征元素分子間引力并相對(duì)一致的數(shù)值，即使計(jì)算出的引力系數(shù)不是絕對(duì)正確，但這一相對(duì)穩(wěn)定的數(shù)值體現(xiàn)了該元素的一種本質(zhì)屬性，應(yīng)具有較大的物理意義。

(3)不同元素引力系數(shù)的比較

不同元素引力系數(shù)變化的規(guī)律為：

①在周期表同一周期中外層電子數(shù)飽和的惰性元素引力系數(shù)最小，其他元素的引力系數(shù)相近，但隨外層電子數(shù)增加依次有所增加；

②在周期表同一列中，從上到下外層電子數(shù)越多，引力系數(shù)越大，也說(shuō)明了分子間引力是原子核正電荷與核外電子的相互作用所引起的。

(4)單個(gè)分子的引力、克分子引力位能的計(jì)算結(jié)果

依據(jù)公式（9～12），不同元素計(jì)算溫度范圍內(nèi)，飽和液相、汽相單個(gè)分子受到的引力和的克分子引力位能（焦耳）如下：氦3（32-44.7，3.63-35.3）、氦4（85-103.4，11.1-61.8）、仲氫（809.9-919.4，60.8-511）、標(biāo)準(zhǔn)氫（703.62-854.52，62.8-561）、氘（1099-1232，107.3-652），氖（1241-1743.7，161-911.4）、氟（4819-7209，91-2761）、氮（3537-5792，345-2947）、氬（4223-6407.8，547-3232）、氧（547-3232，85-3575）、氪（5929-8930，743-4423）、氙（8741-12423，1046-5657）、氯（20861.7-12974，865-8774）。在氣液平衡狀態(tài)下不同溫度時(shí)氣體或液體內(nèi)的單分子受到的引力是不同的，在所計(jì)算的物質(zhì)中氦3分子受到的引力最小，氯分子受到的引力最大；依據(jù)臨界點(diǎn)附近的數(shù)據(jù)外推得到臨界位能和臨界溫度的關(guān)系是正相關(guān)，相關(guān)系數(shù)達(dá)到0.998699。由于溫度是分子動(dòng)能的量度，在分子的動(dòng)能和位能達(dá)到某一平衡時(shí)會(huì)出現(xiàn)臨界狀態(tài)。

圖6、圖7給出了氧元素汽、液相中克分子引力位能及單個(gè)分子受到引力隨溫度的變化。隨溫度上升，液相中的單個(gè)分子受到的引力的下降速率大于引力位能的下降速率；汽相中單個(gè)分子受到的引力的上升速率小于引力位能的上升速率。

圖6 汽液相中克分子氧的引力位能隨溫度變化

圖7 汽液相中氧分子受到的引力隨溫度變化

將同一溫度下汽、液兩相克分子位能進(jìn)行平均，其平均值和溫度線性相關(guān)系數(shù)基本接近于1，說(shuō)明汽、液平衡時(shí)分子動(dòng)能和分子間的位能應(yīng)滿足一定的關(guān)系，圖8給出了氦4的示例。

圖8 氦4液汽相克分子位能的平均值和溫度的關(guān)系

(5)和其他方法計(jì)算出的分子間力的比較

《建立分子簡(jiǎn)化模型用經(jīng)典電磁理論探析分子間作用力》說(shuō)明了分子間引力、斥力、分子電流安培力之間的關(guān)系[16]，同時(shí)計(jì)算了氦分子的分子間力，本計(jì)算和其計(jì)算出的引力比較見(jiàn)表4，差別未超過(guò)數(shù)量級(jí)，計(jì)算結(jié)果有可比性，同時(shí)說(shuō)明了本文計(jì)算出的分子間引力具有一定的可靠性。

表1 氦分子間引力計(jì)算結(jié)果和其他研究結(jié)果的比較

3.利用分子間引力計(jì)算物質(zhì)的臨界溫度和壓強(qiáng)

(1)分子間引力和臨界溫度、壓強(qiáng)關(guān)系式推導(dǎo)

對(duì)式（8）進(jìn)行變換，增加一個(gè)不可壓縮的量則可得到下式：

經(jīng)推導(dǎo)到臨界點(diǎn)：

式中：E—平均公斤引力系數(shù)（kJ·m·kg-4/3），E=f/

Vc—臨界容積（m3/kg）；

β—臨界態(tài)修正系數(shù)，無(wú)量綱；

Pc—臨界壓強(qiáng)（MPa）。

該式說(shuō)明臨界點(diǎn)的氣壓和物質(zhì)引力系數(shù)與（Vc）^（3/4）的倒數(shù)乘積成正比。

式中：R—8.3145（kJ·mol-1·K-1）；Z—壓縮因子。

(2)計(jì)算結(jié)果和討論

①相關(guān)方程計(jì)算結(jié)果

圖9、圖10分別給出了第一周期五種元素、同位素和其余元素分子間引力系數(shù)（沸程內(nèi)的平均值）和臨界體積4/3方倒數(shù)的乘積與臨界壓強(qiáng)的關(guān)系曲線，R2分別達(dá)到0.9679和0.9805，利用相關(guān)系數(shù)反算出的臨界壓力見(jiàn)表2，計(jì)算結(jié)果表明誤差在可接受范圍內(nèi)。圖11、圖12、圖13、圖14給出了分子間引力系數(shù)（沸程內(nèi)的平均值）和臨界體積^（1/3）方倒數(shù)的乘積與臨界溫度、沸點(diǎn)關(guān)系，同樣具有良好的正相關(guān)性。結(jié)果證實(shí)了本文提出引力關(guān)系式的推論，間接證實(shí)了利用汽化熱計(jì)算分子間力的可行性，同時(shí)可以看到分子間引力系數(shù)越大，物質(zhì)的臨界點(diǎn)、沸點(diǎn)溫度越高。

圖9 第一周期元素引力系數(shù)/臨界體積^（4/3）的關(guān)系

圖10 其余元素引力系數(shù)/臨界體積^（4/3）的關(guān)系

表2 利用回歸曲線計(jì)算出的臨界壓力和臨界溫度

圖11 引力系數(shù)/臨界體積^（1/3）與臨界溫度關(guān)系

圖12 引力系數(shù)/臨界體積^（1/3）與沸點(diǎn)溫度關(guān)系

(3)通過(guò)引力公式直接計(jì)算臨界壓力、臨界溫度

由公式（14）（15）可知，在由分子間力計(jì)算臨界壓力、臨界溫度時(shí)需要確定β值。在不作β修整時(shí)，采用的引力系數(shù)是元素沸程范圍內(nèi)的引力系數(shù)的平均值，每個(gè)元素代入壓縮因子Z為實(shí)際壓縮因子條件下計(jì)算出的臨界溫度、臨界壓力和實(shí)際值相差較大，不同元素臨界壓力相對(duì)誤差絕對(duì)值為16.81%～53.78%；臨界溫度為8.24%～53.86%。但考慮到是由引力直接計(jì)算出來(lái)的，臨界壓力、溫度的誤差在50%以下，應(yīng)當(dāng)說(shuō)引力的假設(shè)是基本成立的。每個(gè)元素求得的臨界壓力、臨界溫度修正值β，除氦4外，溫度修正的β和壓力修正的β值基本一致，氬、氪、氙、氮、氧、氟、氯為-0.1351～-0.2158。第一周期元素的β變化較大，標(biāo)準(zhǔn)氫、仲氫、氘、氦4、氦3臨界壓強(qiáng)的β值分別為0.08896、0.0791、-0.0407、0.2109、0.319，臨界溫度的β值分別為0.0672、0.0791、-0.0407、-0.1747、0.3193。

4.由分子間引力計(jì)算飽和液體的表面張力

(1)公式推導(dǎo)

至今尚未有直接利用分子間力計(jì)算表面張力方法，現(xiàn)擬采用計(jì)算出的引力計(jì)算不同溫度下的表面張力。在液相和氣相的邊界層被認(rèn)為是具有和氣相、液相完全不同性質(zhì)的第三相。由于邊界層（表面相）兩側(cè)承受到的吸引力不同，表面處于張緊狀態(tài)，并傾向于保持最小的表面積。表面張力一般用γ表示，定義為表面上單位長(zhǎng)度的力，達(dá)因/厘米。純質(zhì)液體表面張力的估算公式有：Quayle的等張比容法、Sastri的對(duì)比關(guān)聯(lián)式等，產(chǎn)生的誤差一般在5%以內(nèi)[17]。

汽液兩相接觸時(shí)，由于兩相內(nèi)部引力不同，液體表面自由能ΔG=γdS。在等溫等壓下，液體形成新表面時(shí)吉布斯函數(shù)的變化應(yīng)等于除體積功外的其他功。依據(jù)前面的公式，該功應(yīng)是分子間引力所作的功。根據(jù)推理，形成新的汽液表面時(shí)，分子間引力作的功為：

公式的前一部分是指的液體分子由內(nèi)部進(jìn)入表面時(shí)克服液體分子間引力所作的功；后一部分為氣體分子和液體分子接觸形成表面時(shí)，氣體分子引力所作的功。

式中：f—分子間的引力系數(shù)。n—某進(jìn)入液體表面的中心分子最鄰近的分子數(shù)，為滿足表面自由能最小的條件，每個(gè)表面分子一般只失去一個(gè)鄰近分子，因此為1/n，實(shí)際引起表面張力不僅是第一層分子，也有其他層的分子，因此n一修正量，dS/π。r2—單位表面積含有的分子數(shù)。r—液體分子名義半徑；積分第一項(xiàng)為分子引力由液體內(nèi)半徑為r進(jìn)入表面到半徑為過(guò)程中作的功；后一項(xiàng)為dS面積上氣體相對(duì)于液體分子數(shù)的比例及氣體分子由氣體內(nèi)進(jìn)入到汽液表面氣體分子內(nèi)引力作的功，實(shí)際氣體也有多層分子做功為此加入第二修正量m。汽液界面上氣體分子和液體分子受到的引力應(yīng)該相等，因此：

積分后公式為：

對(duì)公式進(jìn)行單位換算，則得汽液相平衡時(shí)液相表面張力如下：

式中：γ—表面張力（μN(yùn)/m）；VL、Vg、VC—液體、氣體、臨界比容（m3/kg）；M—分子量；n、m—修正系數(shù)。

(2)計(jì)算結(jié)果和討論

由不同物質(zhì)表面張力的計(jì)算結(jié)果可知，查表結(jié)果和計(jì)算結(jié)果有一定誤差，但未超過(guò)數(shù)量級(jí)。

不同元素的表面張力誤差的計(jì)算結(jié)果見(jiàn)表3。結(jié)果表明計(jì)算不同元素表面張力的誤差不同，其中氖和氬元素的誤差最大，其他一般在5%～15%之間。不同元素的修正值不同，但多種元素的N值4.5～5，M值在0.5～0.61。不同元素M、N值的不同或計(jì)算和查表結(jié)果不能完全符合的原因可能和表面分子之間的結(jié)構(gòu)有關(guān)。

表3 其他元素表面張力計(jì)算的誤差

5.結(jié)論

（1）利用汽化熱計(jì)算和利用經(jīng)典電磁理論計(jì)算得出的氦分子間引力基本一致，說(shuō)明提出的方法適用于飽和態(tài)汽、液分子間引力的計(jì)算。（2）第一周期元素的分子間引力系數(shù)隨溫度上升而增加，其他周期元素的引力系數(shù)隨溫度上升有不同的變化但和平均值的偏差多數(shù)在1%以內(nèi)，引力系數(shù)越大的元素，到臨界點(diǎn)附近引力系數(shù)有明顯下降。在同一周期中外層電子數(shù)飽和的惰性元素引力系數(shù)最小，其他元素的引力系數(shù)相近，但隨外層電子數(shù)增加依次有所增加；同一列中，從上到下外層電子數(shù)越多，引力系數(shù)越大。原子核內(nèi)的中子數(shù)量對(duì)引力系數(shù)有影響，電子的自旋方向?qū)σο禂?shù)的影響較小。（3）計(jì)算元素臨界點(diǎn)的單分子受到引力為10-13～10-11N，克分子引力位能為39～10000J；元素克分子引力位能和元素的臨界溫度線性相關(guān)，相關(guān)系數(shù)達(dá)到0.9987。元素臨界壓力與平均引力系數(shù)和臨界體積4/3次方倒數(shù)的乘積正相關(guān)，元素臨界溫度與平均引力系數(shù)和臨界體積的1/3次方倒數(shù)乘積正相關(guān)。利用推導(dǎo)出的公式可以計(jì)算出元素的臨界壓力和臨界溫度。（4）利用計(jì)算出的分子間力可通過(guò)推導(dǎo)出的公式計(jì)算出飽和液體的表面張力，部分元素的計(jì)算結(jié)果可滿足工程要求，部分元素計(jì)算出的結(jié)果和查表值誤差較大，原因可能在于不同元素汽、液的分子結(jié)構(gòu)不同有關(guān)。（5）上文是分子間力定量化計(jì)算的一種有意義的方法，將會(huì)促進(jìn)分子間引力研究的發(fā)展，為改善物質(zhì)性能或生成新物質(zhì)提供科學(xué)依據(jù)。