羅強(qiáng) 楊恒 郭平? 趙建飛

1)(西南石油大學(xué)理學(xué)院,成都 610500)

2)(西南石油大學(xué),油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,成都 610500)

1 引 言

天然氣水合物作為一種清潔、高效的非常規(guī)能源,近年來(lái)得到了許多國(guó)家的高度重視和大力發(fā)展[1-3].目前,對(duì)水合物的研究主要集中在溫度[4]、壓強(qiáng)[5]和表面活性劑[6]及CO2對(duì)水合物的置換等方面的實(shí)驗(yàn)研究[7-10],并通過(guò)實(shí)驗(yàn)方法對(duì)3種類型的水合物結(jié)構(gòu)給出X-ray[11],Raman[12]譜.盡管實(shí)驗(yàn)研究取得了長(zhǎng)足進(jìn)展,但理論研究相對(duì)薄弱[13],在理論上,人們采用了分子動(dòng)力學(xué)(MD),蒙特卡羅(MC)等方法,模擬了溫度壓強(qiáng)對(duì)水合物結(jié)構(gòu)穩(wěn)定性、能量和熱力學(xué)性質(zhì)的影響[14-21],但難以給出其電子性質(zhì).為了深入研究水合物的微觀機(jī)理.人們采用密度泛函理論(DFT)對(duì)其做了研究,但主要針對(duì)結(jié)構(gòu)相對(duì)簡(jiǎn)單的I型水合物(512612)水合物的結(jié)構(gòu)和性質(zhì)進(jìn)行了研究,對(duì)II型及H型水籠子研究相對(duì)較少[22-30].

水合物是類冰狀的結(jié)晶物質(zhì),在一定條件下可以穩(wěn)定存在,結(jié)構(gòu)受客體分子和溫壓的影響,所以很多高壓低溫的星體,存在新結(jié)構(gòu)水合物的可能,探索水合物結(jié)構(gòu)可以有助于認(rèn)識(shí)外星球的組成.其次水合物有儲(chǔ)存氣體的能力,可以儲(chǔ)運(yùn)天然氣和埋存CO2,研究新結(jié)構(gòu)降低水合物形成能也是一個(gè)新的研究方向.Huang等[31]采用MC方法預(yù)測(cè)了III籠型冰相.曹瀟瀟等[32,33]研究了大客體烷烴分子與水籠子之間的穩(wěn)定性,也做出了新的預(yù)測(cè).

但我們發(fā)現(xiàn)給出的水合物結(jié)構(gòu)存在五重對(duì)稱軸,單個(gè)水籠子難以簡(jiǎn)單地在空間存在周期性排列,會(huì)導(dǎo)致運(yùn)算量大、計(jì)算成本高.所以理論上尋找結(jié)構(gòu)簡(jiǎn)單且具有高對(duì)稱的水合物新相具有重要意義.

本研究采用密度泛函理論[34,35]第一性原理方法,構(gòu)建出新型結(jié)構(gòu)的甲烷水合物,并對(duì)其結(jié)構(gòu)和電子性質(zhì)等進(jìn)行研究.通過(guò)得到結(jié)構(gòu)參數(shù)、形成能、XRD圖、態(tài)密度和能帶,為研究天然氣水合物的賦存、合成和儲(chǔ)運(yùn)等提供理論參數(shù)依據(jù).

2 計(jì)算模型與方法

2.1 模 型

模型是受到固體物理中面心立方晶體的第一布里淵區(qū)(具有周期性、高對(duì)稱性的截角八面體)啟發(fā).先做出一個(gè)截角八面體,將水分子中的氧原子置于其頂角上,甲烷分子中的碳原子放在中心,以此作為骨架,此時(shí)對(duì)稱群為IMM,如圖1(a)所示.采用自動(dòng)加氫方式填加氫原子,由8個(gè)在平行于{111}面上(位于過(guò)體對(duì)角線1/4處)的正六邊形和6個(gè)在{100}面上的正方形,構(gòu)成了具有截角八面體結(jié)構(gòu)的水籠子(4668),作為本研究的基本模型.通過(guò)調(diào)整氫原子的位置,再經(jīng)過(guò)計(jì)算機(jī)優(yōu)化后得到的結(jié)構(gòu)如圖1(b)所示,為了便于描述,稱其為N型甲烷水合物,晶格參數(shù)為7.700 ?,密度0.903 g/cm3,大于 I型的 0.813 g/cm3、II型的0.804 g/cm3和H型的0.768 g/cm3[36].

圖1 模型示意圖(a)N型甲烷水合物骨架模型;(b)優(yōu)化后的結(jié)構(gòu)Fig.1.Model diagram:(a)N-methane hydrate skeleton model;(b)optimized structure of N-methane hydrate.

2.2 計(jì)算方法及參數(shù)設(shè)置

計(jì)算采用軟件是Materials Studio 2016的CASTEP[37,38]模塊,它是對(duì)固體材料進(jìn)行理論模擬研究的量子力學(xué)程序,由英國(guó)劍橋大學(xué)凝聚態(tài)理論研究組開(kāi)發(fā),其可靠性已經(jīng)為大量實(shí)際計(jì)算所驗(yàn)證[39,40].本文使用CASTEP計(jì)算時(shí),采用PBE形式的廣義梯度近似(generalized gradient approximation,GGA)[41]泛函處理交換關(guān)聯(lián)勢(shì)進(jìn)行計(jì)算,并考慮了Grimme色散修正.

經(jīng)收斂性測(cè)試及計(jì)算,平面波截?cái)嗄苋?00 eV,K點(diǎn)[42]網(wǎng)格尺寸采用3×3×3,選擇超軟贗勢(shì)處理電子離子間相互作用,單原子能量收斂精度為5.0×10—5eV/atom,作用于原子最大力收斂精度是0.01e V/?,最大壓力收斂精度是0.2 GPa,最大位移收斂精度是0.005 ?,能量收斂算法采用BFGS[43-46]優(yōu)化算法,SCF自洽精度設(shè)為單原子能量收斂至 5.0×10—7eV.

3 計(jì)算結(jié)果及分析

3.1 優(yōu)化后結(jié)構(gòu)參數(shù)

由于N型甲烷水合物是一個(gè)新結(jié)構(gòu),難以與前人的研究結(jié)果進(jìn)行對(duì)比,但是其結(jié)構(gòu)與I型甲烷水合物較為相似,所以我們和I型甲烷水合物結(jié)構(gòu)[47]進(jìn)行對(duì)比研究; 另一方面對(duì)于交換關(guān)聯(lián)泛函的處理,采用GGA進(jìn)行計(jì)算,并采用了Grimme色散修正,優(yōu)化前后晶格參數(shù)如表1所列.

表1 平均晶格參數(shù)Table 1. Average lattice parameters.

N型甲烷水合物水籠子是由6個(gè)四邊形和8個(gè)六邊形構(gòu)成,從表1中可以看出,優(yōu)化前的六邊形平均鍵長(zhǎng)為2.6804 ?,優(yōu)化后的六邊形平均邊長(zhǎng)分別為2.7231 ?,與I型甲烷水合物水籠子中的六邊形平均邊長(zhǎng)2.7152 ?較為接近; 優(yōu)化前四邊形平均邊長(zhǎng)為2.7656 ?,優(yōu)化后為2.7233 ?; 水分子的平均鍵長(zhǎng)在優(yōu)化前為0.994 ?,優(yōu)化后為1.0056 ?,優(yōu)化前后水分子鍵長(zhǎng)幾乎沒(méi)有發(fā)生變化,與已知I型甲烷水合物計(jì)算的0.993 ?水分子平均鍵長(zhǎng)亦很接近.水分子平均鍵角優(yōu)化前是109.406°,優(yōu)化后的鍵角變小,為 107.738°,與 I型甲烷水合物計(jì)算的106.629°接近.

甲烷分子平均鍵長(zhǎng)優(yōu)化前為1.1178 ?,優(yōu)化后為 1.0973 ?,比表1中的 I型甲烷水合物的1.0917 ?大,這是N型甲烷水合物水籠子比I型甲烷水合物水籠子大,水籠子與甲烷相互作用較弱導(dǎo)致的.甲烷分子的平均鍵角優(yōu)化前為109.454°,優(yōu)化后的結(jié)構(gòu)與I型甲烷水合物值相同,說(shuō)明優(yōu)化并不影響甲烷分子上氫原子的位置; 甲烷分子到水分子平均距離優(yōu)化前為4.278 ?,優(yōu)化后為4.2831 ?,比表1中的I型甲烷水合物的值3.7812 ?要大,說(shuō)明N型甲烷水合物水籠子空間更大,能容納更大體積的氣體分子.

氫鍵是水合物穩(wěn)定性的關(guān)鍵因素,優(yōu)化前的氫鍵平均鍵長(zhǎng)為1.7025 ?,優(yōu)化后為1.7293 ?,而表1中I型甲烷水合物參考值是1.7128 ?,I型甲烷水合物平均氫鍵鍵長(zhǎng)較長(zhǎng),這可能是由于N型甲烷水合物籠子更大,甲烷分子到水分子平均距離大,甲烷影響作用弱導(dǎo)致.

從結(jié)構(gòu)分析中可以得出,N型甲烷水合物水分子的鍵長(zhǎng)鍵角比初始結(jié)構(gòu)變化較小,這是由于水分子是極性分子,在形成水籠子的過(guò)程中出現(xiàn)極化,形成氫鍵所導(dǎo)致的; N型甲烷水合物中六邊形平均邊長(zhǎng),水分子平均鍵長(zhǎng)、鍵角,氫鍵平均鍵長(zhǎng)都比I型甲烷水合物的參數(shù)要大一些,而甲烷平均鍵角沒(méi)有發(fā)生變化,甲烷平均鍵長(zhǎng)增加,這是由于N型甲烷水合物籠子較I型甲烷水合物的大,甲烷分子與水籠子相互作用較弱,由范德瓦耳斯(van der Waals)力結(jié)合所導(dǎo)致的.

3.2 X射線衍射(XRD)

為了進(jìn)一步了解N型甲烷水合物結(jié)構(gòu),從理論上做出了其X射線衍射圖(XRD),如圖2所示.由于N型甲烷水合物是從理論上提出的全新甲烷水合物結(jié)構(gòu),實(shí)驗(yàn)數(shù)據(jù)相對(duì)缺乏,由于與I型甲烷水合物結(jié)構(gòu)比較接近,我們參考對(duì)比丁家祥等[48]I型甲烷水合物實(shí)驗(yàn)數(shù)據(jù).

從圖2可以看出,N型甲烷水合物在2θ值為16.5°,23.0°,28.5°,33.0°,37.0°,41.0°,44.0°和47.5°附近出現(xiàn)了 8個(gè)峰值,與 I型甲烷水合物XRD實(shí)驗(yàn)衍射圖峰值中的2θ值符合較好,但比實(shí)驗(yàn)中的峰值個(gè)數(shù)少.這是因?yàn)閷?shí)際水合物中可能存在雜質(zhì),丁家祥等[48]做的實(shí)驗(yàn)溫度是 —90 ℃,而我們采用的是理想晶格,所選的溫度趨于絕對(duì)零度,且結(jié)構(gòu)本身與I型甲烷水合物結(jié)構(gòu)又有所不同,這些均導(dǎo)致了峰值個(gè)數(shù)會(huì)出現(xiàn)差異.

圖2 N型甲烷水合物結(jié)構(gòu)X射線衍射圖Fig.2.X-ray diffraction of N-methane hydrate structure.

3.3 形成能

人們認(rèn)為水合物中客體分子與水籠子相互作用,是由范德瓦耳斯力作用導(dǎo)致的,為了考察N型甲烷水合物形成的難易程度,采用GGA進(jìn)行優(yōu)化計(jì)算,同時(shí)考慮Grimme色散修正,計(jì)算公式如下[45]:

式中 ΔEform為N型甲烷水合物的形成能,Etotal為N型甲烷水合物的總能量,Ecage為N型甲烷水合物水籠子的總能量,ECH4為甲烷分子的能量,n為晶格中甲烷分子個(gè)數(shù).根據(jù)(1)式計(jì)算出GGA近似情況下的形成能,列入表2中 ΔEform行.

表2 GGA近似下的形成能Table 2. Formation energy of GGA approximation.

表2中可看出,優(yōu)化后形成能為 —0.247 eV,較I型甲烷水合物的形成能 —0.581 eV絕對(duì)值更小,表明N型甲烷水合物易形成,但穩(wěn)定性比I型甲烷水合物差; 同時(shí)也反映了甲烷分子形成N型甲烷水合物的作用非常微弱,它們之間靠范德瓦耳斯力的結(jié)合形成.

3.4 態(tài)密度及分波態(tài)密度

為更好地理解N型甲烷水合物形成的機(jī)理,對(duì)其電子性質(zhì)作了進(jìn)一步研究,分析態(tài)密度(DOS)及分波態(tài)密度(PDOS)是研究電子性質(zhì)的主要手段之一.具體針對(duì)單個(gè)甲烷分子及水籠子中的甲烷分子、水籠子與N型甲烷水合物的水籠子的態(tài)密度及分波態(tài)密度進(jìn)行對(duì)比分析,其中能量為零處為費(fèi)米能級(jí),態(tài)密度的展寬寬度為0.2 eV.

3.4.1 甲烷分子態(tài)密度及分波態(tài)密度

圖3為單個(gè)甲烷分子及水籠子中的甲烷分子的態(tài)密度及分波態(tài)密度,從圖3可以看出,水籠子中甲烷分子的態(tài)密度較處于自由狀態(tài)的單個(gè)甲烷分子,其能量分布整體向左平移; 從s態(tài)和p態(tài)電子分波態(tài)密度來(lái)看,在費(fèi)米能級(jí)附近及以下有類似的變化趨勢(shì); 在能量約為5—17 eV區(qū)間,出現(xiàn)了明顯的sp軌道雜化,且總體變化趨勢(shì)朝能量更低的方向移動(dòng),這也表明了水籠子中甲烷分子能量更低,在接近絕對(duì)零度時(shí),易形成甲烷水合物.

3.4.2 水籠子態(tài)密度及分波態(tài)密度

為了進(jìn)一步理解甲烷與N型水籠子的相互作用,采用GGA近似下的態(tài)密度及分波態(tài)密度對(duì)水籠子做了分析,如圖4所示.N型水籠子與空籠子的態(tài)密度圖形大部分重合,主要在能量5—17 eV區(qū)間附近態(tài)密度發(fā)生了較微弱的變化,結(jié)合分波態(tài)密度來(lái)看,s態(tài)電子變化可以忽略,主要是氧原子的p態(tài)電子引起的微弱變化,這也反映出甲烷分子和水籠子的作用非常微弱,它們之間靠范德瓦耳斯力結(jié)合.

圖3 甲烷分子態(tài)密度和分波態(tài)密度Fig.3.Density of states and partial density of states in methane molecules.

圖4 水籠子的態(tài)密度及分波態(tài)密度Fig.4.Density of state and partial density of state in water cage.

圖5 能帶結(jié)構(gòu)Fig.5.Energy band structures.

3.5 能 帶

能帶結(jié)構(gòu)是反映物質(zhì)電子性質(zhì)的另一重要概念,通過(guò)理論計(jì)算,給出了GGA近似下的能帶結(jié)構(gòu),以及與之對(duì)應(yīng)的N型甲烷水合物和水籠子中甲烷態(tài)密度,如圖5所示.

費(fèi)米能級(jí)(能量為零)以上的是導(dǎo)帶,以下為價(jià)帶,能隙(禁帶寬度)為導(dǎo)帶低與價(jià)帶頂?shù)哪芰恐?N型甲烷水合物的能隙為5.397 eV,價(jià)帶頂與導(dǎo)帶底位于G點(diǎn)處,可看為絕緣體; 結(jié)合N型甲烷水合物和水籠子中的甲烷態(tài)密度看,價(jià)帶頂與導(dǎo)帶底附近的能態(tài)密度由水籠子貢獻(xiàn),再結(jié)合圖4中水籠子的分波態(tài)密度來(lái)看,組成水籠子中的水分子,其氧原子的p態(tài)電子對(duì)禁帶寬度起了決定作用,而甲烷分子的態(tài)密度偏離帶低及帶頂附近,混合在水籠子的能帶中,并未進(jìn)入禁帶,導(dǎo)致能隙發(fā)生改變,從這個(gè)角度來(lái)看,也體現(xiàn)了甲烷與水籠子的微弱作用,在宏觀層面表現(xiàn)出是范德瓦耳斯作用.

4 結(jié) 論

基于密度泛函理論,對(duì)首次提出的N型甲烷水合物,運(yùn)用MS軟件中CASTEP模塊研究了其微觀機(jī)理.考慮Grimme色散修正,采用GGA進(jìn)行對(duì)比計(jì)算,分析其結(jié)構(gòu)、XRD、形成能、電子態(tài)密度和能帶結(jié)構(gòu)等,得到了如下結(jié)論.

1)N型甲烷水合物結(jié)構(gòu)為截角八面體型的水籠子(4668),對(duì)稱群為IMM,為目前已知的對(duì)稱性最高的水合物水籠子,晶格參數(shù)為7.700 ?,密度為0.903 g/cm3;

2)N型甲烷水合物中六邊形平均邊長(zhǎng),水分子平均鍵長(zhǎng)、鍵角,氫鍵平均鍵長(zhǎng)都比I型甲烷水合物的參數(shù)要大,這是由N型甲烷水合物水籠子比I型甲烷水合物水籠子大所導(dǎo)致; 而甲烷平均鍵角沒(méi)有發(fā)生變化,甲烷平均鍵長(zhǎng)增加,甲烷分子與水籠子相互作用力為范德瓦耳斯力; 結(jié)合XRD來(lái)看,N型甲烷水合物與I型甲烷水合物結(jié)構(gòu)較為接近;

3)N型甲烷水合物形成能、電子態(tài)密度及分波態(tài)密度均表明了甲烷與N型水籠子相互作用微弱,靠分子力作用;

4)N型甲烷水合物為絕緣體材料,能隙大于5 eV.