水系統(tǒng)結(jié)構(gòu)和動(dòng)力學(xué)特性的計(jì)算機(jī)模擬研究.pdf

1、本文結(jié)合計(jì)算機(jī)應(yīng)用技術(shù)和物理學(xué)原理,對(duì)于水系統(tǒng)的結(jié)構(gòu)和動(dòng)力學(xué)特性進(jìn)行了計(jì)算機(jī)模擬仿真研究,實(shí)現(xiàn)了計(jì)算機(jī)方法在物理研究中的實(shí)踐應(yīng)用,并得到了一些令人滿意的結(jié)果。
　　 本論文主體可分為五大部分。前兩章為基礎(chǔ)知識(shí)和原理的介紹,后三章針對(duì)研究課題,闡述開展的主要工作。
　　 第一章先介紹了水系統(tǒng)的微觀結(jié)構(gòu)和性質(zhì),然后闡述了水的固體相——冰的結(jié)構(gòu)和物理特征。在此基礎(chǔ)上,本章節(jié)對(duì)于水系統(tǒng)在兩種特殊條件下形成的奇異結(jié)構(gòu)——?dú)怏w水化合

2、物和受限空間水進(jìn)行了詳細(xì)介紹。
　　 第二章主要闡述計(jì)算機(jī)模擬方法在物理學(xué)中的應(yīng)用發(fā)展以及一些具體的計(jì)算機(jī)模擬理論。其中,對(duì)晶格動(dòng)力學(xué)模擬方法和分子動(dòng)力學(xué)模擬方法進(jìn)行了詳細(xì)的介紹。
　　 第三章主要內(nèi)容是在前兩章的基礎(chǔ)上,根據(jù)水系統(tǒng)的物理特征,采用計(jì)算機(jī)模擬的方法,對(duì)于描述一型、二型氣體水化合物和受限于碳納米管水的物理模型進(jìn)行優(yōu)化選擇,然后對(duì)需要分析的物理參量進(jìn)行程序編碼,實(shí)現(xiàn)從物理語言描述到計(jì)算機(jī)語言描述的轉(zhuǎn)換。通過編

3、制程序在計(jì)算機(jī)中的運(yùn)行,把實(shí)際物理問題模擬抽象成二進(jìn)制數(shù)據(jù)。這一章中根據(jù)水系統(tǒng)中的“冰規(guī)則”,我們對(duì)水合物晶格上氫原子坐標(biāo)的無序放置提出了自己有效的算法思想。通過對(duì)幾種不同的水分子模型的計(jì)算,根據(jù)輸出優(yōu)劣性的比較;選擇了TIP4P水分子模型作為我們所研究的水系統(tǒng)中的基本單元模型。在此基礎(chǔ)上,我們定義了氣體水合物體系和碳納米管水體系的模擬細(xì)節(jié),對(duì)研究體系進(jìn)行完整全面的建模操作。要從量化角度界定物理體系的一些特征,通常需要考量一些物理過程和

4、物理參量。這一章節(jié)的后半部分則是根據(jù)一些物理過程和物理參量的定義,將其用計(jì)算機(jī)語言表述出來。這些表述既要符合物理規(guī)律的真實(shí)含義,再現(xiàn)物理過程,同時(shí)又要求符合計(jì)算仿真的算法要求,能夠在設(shè)備上有效運(yùn)行,得到精確的輸出結(jié)果。這為后面章節(jié)中對(duì)兩種特殊水系統(tǒng)的模擬仿真研究,提供了計(jì)算機(jī)技術(shù)支持。
　　 第四章主要內(nèi)容是使用晶格動(dòng)力學(xué)計(jì)算方法對(duì)氣體水化合物結(jié)構(gòu)和動(dòng)力學(xué)特性的模擬結(jié)果的討論,并將計(jì)算結(jié)果與已有的實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比分析。在上一章中

5、已經(jīng)對(duì)于氣體水化合物的模擬體系進(jìn)行了建模,并對(duì)一些描述氣體水化合物物理特性的參量進(jìn)行了計(jì)算機(jī)語言的闡述。在此基礎(chǔ)上,這一章將對(duì)于不同晶格結(jié)構(gòu)的一型和二型氣體水化合物進(jìn)行比較性的研究,首先對(duì)于結(jié)構(gòu)比較簡單的單原子分子稀有氣體水合物進(jìn)行研究。討論了同一種氣體分子在分別填入一型和二型氣體水化合物之后其動(dòng)力學(xué)特性出現(xiàn)的差異,指出十二面體小籠結(jié)構(gòu)在整個(gè)氣體水化合物體系中扮演了重要角色,它與客體分子的耦合作用最為強(qiáng)烈,對(duì)于結(jié)構(gòu)不同的一型和二型氣體水

6、化合物的穩(wěn)定機(jī)制作用明顯。在同一種水化合物的晶格結(jié)構(gòu)中分別填入直徑不同的單原子氣體分子,考察客體分子大小在整個(gè)體系穩(wěn)定中所起到的作用,總結(jié)出當(dāng)客體分子與優(yōu)勢籠的有效直徑比值在0.745以上時(shí),整個(gè)水化合物的水分子晶格結(jié)構(gòu)就可以穩(wěn)定存在。當(dāng)甲烷分子填入二型氣體水化合物時(shí),恰恰符合以上總結(jié)的規(guī)律。在此基礎(chǔ)上初步預(yù)測了二型甲烷氣體水化合物的存在。然后進(jìn)一步對(duì)雙原子分子氮?dú)庑纬傻臍怏w水合物的動(dòng)力學(xué)特性以及能量穩(wěn)定性進(jìn)行模擬,發(fā)現(xiàn)了雙原子分子內(nèi)部

7、振動(dòng)的激發(fā)峰值,這是由于非單原子氣體的復(fù)雜內(nèi)部結(jié)構(gòu)而產(chǎn)生的,這為下一步研究甲烷分子在水化合物晶格內(nèi)部振動(dòng)的動(dòng)力學(xué)特性提供了依據(jù)。同時(shí)對(duì)氮?dú)夥肿釉诨\型結(jié)構(gòu)中處于不同取向的晶格能量進(jìn)行計(jì)算,發(fā)現(xiàn)氮?dú)夥肿釉诨\型結(jié)構(gòu)內(nèi)部沒有最低勢能取向,并且得到了氮?dú)馑衔锏淖畲筇畛渎?證明了氮?dú)夥肿与p填充結(jié)構(gòu)的穩(wěn)定存在。在單原子分子氣體水合物和雙原子分子氣體水合物的研究基礎(chǔ)上,最后討論了結(jié)構(gòu)比較復(fù)雜的甲烷氣體水合物,對(duì)甲烷氣體水合物的氣體填充率和晶格能量關(guān)

8、系進(jìn)行討論,預(yù)測二型甲烷氣體水合物在自然界中存在的可能以及在大籠結(jié)構(gòu)中雙填充的可能。
　　 第五章討論了使用分子動(dòng)力學(xué)計(jì)算方法對(duì)受限于單壁碳納米管中的水分子性質(zhì)的模擬結(jié)果。在第三章中對(duì)于碳納米管及其周圍水分子體系的建模基礎(chǔ)之上,從管徑和手性兩個(gè)因素出發(fā)來討論它們對(duì)于碳納米管中水分子靜態(tài)性質(zhì)和動(dòng)態(tài)性質(zhì)的影響。對(duì)于碳納米管水的相變這一靜態(tài)性質(zhì),我們主要從其結(jié)構(gòu)角度入手,考量水分子在常壓和高壓兩種情況下在不同直徑和手性的碳納米管中的相

9、變結(jié)構(gòu)。發(fā)現(xiàn)碳管內(nèi)水分子禁區(qū)的存在,并且隨著壓強(qiáng)的增大這一禁區(qū)從2.8(A)減小至2.4(A)左右。在相同壓強(qiáng)下,隨著碳納米管直徑的增大,其形成的有序冰納米管結(jié)構(gòu)呈現(xiàn)更加復(fù)雜的結(jié)構(gòu),并且相變溫度隨管子直徑的增大而降低,出現(xiàn)了與理論預(yù)測完全相反的結(jié)果。對(duì)于相同的碳納米管,壓強(qiáng)較大的情況下,其內(nèi)部水分子呈現(xiàn)更加密集的分布結(jié)構(gòu),對(duì)于(9,9)和(10,10)兩種較粗的管子,高壓下出現(xiàn)了冰納米管與冰鏈的雙層分布結(jié)構(gòu)。通過對(duì)水分子偶極矩分布的分析

10、,發(fā)現(xiàn)對(duì)于碳納米管水的靜態(tài)結(jié)構(gòu)分布,管壁與水分子之間的相互作用比水分子之間的氫鍵相互作用對(duì)其產(chǎn)生的影響更大。對(duì)于不同手性的碳納米管,只要直徑尺度相似,則其內(nèi)部水分子的相變性質(zhì)是一樣的。此外,還對(duì)各種不同手性和直徑碳納米管中水分子的平均氫鍵數(shù)進(jìn)行了計(jì)算,發(fā)現(xiàn)其平均氫鍵數(shù)均小于大塊水的狀態(tài),而且管徑越小,氫鍵結(jié)構(gòu)的缺陷越多。也就是說只有碳納米管的直徑對(duì)水的相變性質(zhì)產(chǎn)生影響,碳納米管的手性不對(duì)其產(chǎn)生影響。對(duì)不同直徑的armchair型和zig

11、zag型碳納米管內(nèi)水分子的遷移率進(jìn)行了計(jì)算,比較后發(fā)現(xiàn)對(duì)于碳納米管水的動(dòng)態(tài)輸運(yùn)性質(zhì),水分子之間的氫鍵相互作用比管壁與水分子之間的相互作用對(duì)其影響更大。隨著管徑的增大,相同手性碳管內(nèi)水分子的遷移通量增大,但是遷移率卻減小。在相同直徑的碳納米管內(nèi),水分子在armchair型管內(nèi)部的遷移率要比zigzag型的遷移率大,即碳納米管水的動(dòng)力學(xué)特性不僅僅與碳管的直徑有關(guān),還與碳管的手性有關(guān)。這些研究對(duì)于碳納米管制備中不同直徑和手性碳管分離的必要性和