畢業(yè)論文-Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)研究.doc

平頂山學(xué)院2014屆本科生畢業(yè)論文 Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)研究 劉倩倩Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)研究摘要本文是從銅鐵原子總數(shù)目和銅鐵原子個(gè)數(shù)比兩個(gè)維度來(lái)研究Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)。用Materials Studio(簡(jiǎn)稱MS)軟件建立初始模型，選取三種不同銅鐵原子的尺寸(m+n)和六種不同銅鐵原子組成(m/n)的Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)后，從MS里導(dǎo)出的數(shù)據(jù)，用LAMMPS編程進(jìn)行計(jì)算研究，并在VMD中觀察結(jié)果。結(jié)果表明，相同原子尺寸中銅鐵原子組成越小，即鐵原子含量越多，其原子勢(shì)能越低，合金團(tuán)簇結(jié)構(gòu)越穩(wěn)定；相同原子組成中原子尺寸越大，即原子總數(shù)目越大，其原子勢(shì)能越小，合金團(tuán)簇越穩(wěn)定。關(guān)鍵詞：納米團(tuán)簇，合金，分子動(dòng)力學(xué)模擬，Materials Studio Molecular Dynamics Research of Cu - Fe Binary Alloy Nano Clusters Structure Abstract This article is about to study the molecular-dynamics simulation of Cu-Fe binary alloy nanocluster structure from the aspects of the total number of the Copper and iron atoms and the proportion of the copper and iron atoms. With Materials Studio software (MS), the initial model is established. After the total number of the three different copper and iron atoms (m + n) and six different proportion of the copper and iron atoms (m/n)are selected, data is derived through the MS.Then using LAMMPS programming calculated is conducted and research is observed in VMD. The results show that in the same total number of the atoms copper and iron atoms, the more content of iron atoms is, the lower atomic potential energy is ,that is to say , the more stable alloy clusters structure is; Meanwhile in the same proportion of the copper and iron atoms, the larger total number of the Copper and iron atoms is, the lower atomic potential energy is, namely,the greater the stability of the alloy clusters is.Key Words:nanocluster,alloy,Molecular dynamics simulation ,Materials Studio 目 錄1 緒論12 納米團(tuán)簇結(jié)構(gòu)22.1 合金概述22.2 納米團(tuán)簇介紹23 分子動(dòng)力學(xué)基本理論93.1 理論基礎(chǔ)93.2 分子動(dòng)力學(xué)的實(shí)現(xiàn)及常用的技術(shù)113.3原子間相互作用勢(shì)124 Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)模擬方法144.1 相關(guān)軟件介紹144.2 模擬計(jì)算過(guò)程155 Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)研究16結(jié)果及分析16參考文獻(xiàn)19致謝201 緒論在過(guò)去的二十幾年中，納米科學(xué)的研究發(fā)展非常迅猛，在這個(gè)研究領(lǐng)域中，人們發(fā)現(xiàn)納米材料在數(shù)據(jù)的儲(chǔ)存、生物醫(yī)學(xué)、電磁波吸收與屏蔽以及提高催化活性等領(lǐng)域有潛在的應(yīng)用前景。在納米材料中，團(tuán)簇研究可追溯到50年代后期Pecker等用超聲噴注法獲得團(tuán)簇。之后，法國(guó)科學(xué)家Leleyter和Joyes在研究濺射過(guò)程中發(fā)現(xiàn)各種帶電和中性團(tuán)簇。但直到70年代末仍處于零星分散的狀態(tài)。8O年代國(guó)際上團(tuán)簇研究有了迅速的發(fā)展，取得了令人觸目的進(jìn)展，其中最為突出的是，1984年美國(guó)加州大學(xué)伯克利分校Knight等發(fā)現(xiàn)超聲膨脹產(chǎn)生Na團(tuán)簇具有幻數(shù)結(jié)構(gòu)，與其價(jià)電子結(jié)構(gòu)呈殼層分布相對(duì)應(yīng)。接著發(fā)現(xiàn)C60籠形團(tuán)簇及其大量制備的簡(jiǎn)單方法，引起科學(xué)界的轟動(dòng)。一些著名大學(xué)和研究機(jī)構(gòu)都積極開(kāi)展團(tuán)簇研究，如瑞士無(wú)機(jī)分析和物化研究所，美國(guó)能源部、阿貢國(guó)立實(shí)驗(yàn)室、海軍研究所，芝加哥大學(xué)、德國(guó)馬一普研究所、薩爾蘭大學(xué)，日本分子科學(xué)研究所和東京大學(xué)等。召開(kāi)了一系列以團(tuán)簇為中心議題的國(guó)際會(huì)議。例如，從1976年開(kāi)始已連續(xù)召開(kāi)八屆的小顆粒和無(wú)機(jī)團(tuán)簇國(guó)際會(huì)議(ISSPIC)，小團(tuán)簇物理和化學(xué)國(guó)際會(huì)議(NATO，1986)，元素和分子團(tuán)簇國(guó)際會(huì)(Italy，1987)，微團(tuán)簇的第一屆NEC會(huì)議(日本，1987)等。l989年美國(guó)能源部材料科學(xué)委員會(huì)、基礎(chǔ)能源局和材料科學(xué)部聯(lián)合主持召開(kāi)來(lái)自學(xué)術(shù)界和工業(yè)界的物理、化學(xué)和材料等各方面專家對(duì)“團(tuán)簇和團(tuán)簇合成材料的研究機(jī)遇”進(jìn)行戰(zhàn)略討論，分析了當(dāng)前研究狀況提出了進(jìn)一步發(fā)展該領(lǐng)域的必要性、規(guī)劃和藍(lán)圖。標(biāo)志著團(tuán)簇科學(xué)研究由初創(chuàng)時(shí)期的分散孤立狀態(tài)向有目的地組織跨學(xué)科協(xié)作以便建立新型學(xué)科體系的方向發(fā)展，由簡(jiǎn)單體系和單一特性的純基礎(chǔ)研究向復(fù)雜系統(tǒng)和綜合性質(zhì)以及基礎(chǔ)研究和應(yīng)用開(kāi)發(fā)密切結(jié)合的方面開(kāi)拓。本文用Materials Studio(簡(jiǎn)稱MS)軟件建立初始模型，選取108個(gè)、365個(gè)和864個(gè)三種銅鐵原子總數(shù)目和0.0、0.2、0.4、0.6、0.8、1.0六種銅鐵原子個(gè)數(shù)比例的Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)后，從MS里導(dǎo)出的數(shù)據(jù)，用LAMMPS編程進(jìn)行計(jì)算研究，并在VMD中觀察結(jié)果。本文是從銅鐵原子總數(shù)目和銅鐵原子個(gè)數(shù)比兩個(gè)維度來(lái)研究Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)，得到直觀且有意義的結(jié)論。2 納米團(tuán)簇結(jié)構(gòu)2.1 合金概述合金是由兩種或兩種以上的元素組成，其中至少有一種為金屬，經(jīng)過(guò)熔煉燒結(jié)或其他方法組合形成具有金屬性的材料。合金存在的狀態(tài)通常由合金的成分、溫度和壓力三個(gè)因素確定，合金的化學(xué)成分變化時(shí)，則合金中所存在的相及相的相對(duì)含量也隨之發(fā)生變化，同樣，當(dāng)溫度和壓力發(fā)生變化時(shí)，合金所存在的狀態(tài)也要發(fā)生變化，因此合金的狀態(tài)可由合金的成分和溫度兩個(gè)因素確定。與純金屬相比，固溶體合金結(jié)晶有如下特點(diǎn)：第一，固溶體合金結(jié)晶同純金屬一樣也遵循著形核與核長(zhǎng)大規(guī)律，不同的是形核時(shí)除結(jié)構(gòu)起伏，能量起伏外還要有成份起伏。第二，純金屬凝固是在恒溫下，合金凝固則是在變溫條件下進(jìn)行（或在一個(gè)溫度區(qū)間）。第三，純金屬無(wú)成份變化，固溶體結(jié)晶過(guò)程自始至終伴隨著異類原子的互擴(kuò)散，發(fā)生著溶質(zhì)在液固兩相中的重新分配。二元合金的過(guò)冷度基本機(jī)理在理論和實(shí)驗(yàn)上經(jīng)過(guò)了嚴(yán)格的驗(yàn)證；在凝固模擬過(guò)程中，微觀模型可以預(yù)測(cè)微觀結(jié)構(gòu)和微觀偏析隨著工藝過(guò)程和合金成分的變化，宏觀模型可以預(yù)測(cè)鑄件充型和凝固過(guò)程中溫度場(chǎng)、流場(chǎng)和濃度場(chǎng)的變化，相圖計(jì)算可以提供多元相平衡信息，為了實(shí)現(xiàn)全面地預(yù)測(cè)熔體充型過(guò)程、鑄件微觀結(jié)構(gòu)和微觀偏析一這一凝固模擬最終目標(biāo)，必須將相圖計(jì)算、微觀模型計(jì)算和宏觀模型計(jì)算整合在一起，目前在這方面還需要開(kāi)發(fā)用戶友好的相圖計(jì)算引擎和高效省時(shí)的數(shù)值計(jì)算方法。2.2 納米團(tuán)簇介紹原子簇研究是化學(xué)中引人關(guān)注的領(lǐng)域。它可以幫助人們認(rèn)識(shí)原子團(tuán)簇是怎樣作為基本單元，構(gòu)造起微觀原子與宏觀材料間的橋梁。盡管等離子體質(zhì)譜實(shí)驗(yàn)?zāi)芴峁└鞣N原子簇生成的信息，但仍然需要理論計(jì)算去研究其結(jié)構(gòu)和能量。早期的理論計(jì)算假設(shè)原子簇是某種結(jié)構(gòu)（例如fcc和bcc）大塊固體的碎片，進(jìn)一步研究發(fā)現(xiàn)原子簇以多面體方式生長(zhǎng)，非晶態(tài)的結(jié)構(gòu)占較為重要地位（即含有5次對(duì)稱軸）。我們應(yīng)用分子動(dòng)力學(xué)(MD)方法研究了金銅合金體系。比較原子簇增大過(guò)程中非晶態(tài)與晶態(tài)結(jié)構(gòu)穩(wěn)定性差異。在分子動(dòng)力學(xué)模擬中，描述原子間相互作用關(guān)系的勢(shì)函數(shù)是整個(gè)模擬的關(guān)鍵。傳統(tǒng)的兩體勢(shì)在研究材料的彈性系數(shù)性質(zhì)和預(yù)言金屬的結(jié)合能及空位形成能時(shí)難以獲得準(zhǔn)確的結(jié)果。80年代以來(lái)，各種經(jīng)驗(yàn)或半經(jīng)驗(yàn)的多體模型得到迅速發(fā)展，特別是嵌入原子法(EAM)，即克服了兩體勢(shì)的缺陷，又不會(huì)使計(jì)算量太大。由幾個(gè)至上千個(gè)原子、分子結(jié)合成相對(duì)穩(wěn)定的微觀和亞微觀聚集體叫做團(tuán)簇。它是介于原子、分子和宏觀固體之間的一個(gè)新層次，它的許多物理、化學(xué)性能既不同于單個(gè)原子、分子，又不同于常規(guī)固體，是凝聚態(tài)物質(zhì)中的一種新結(jié)構(gòu)。而納米團(tuán)簇的組裝則是指納米團(tuán)簇或納米團(tuán)簇化合物在一維、二維或三維空間中以規(guī)則的順序排列成有特定性質(zhì)和功能的點(diǎn)陣的過(guò)程。組裝成的結(jié)構(gòu)有的保持了納米團(tuán)簇或納米團(tuán)簇化合物的原有特性，有的則產(chǎn)生了新的特性和功能。通過(guò)對(duì)納米粒子的修飾、納米粒子間分子繩索的結(jié)構(gòu)和大小以及組裝載體的性質(zhì)的調(diào)整可以獲得各種性能的結(jié)構(gòu)。2.2.1 納米團(tuán)簇的特性現(xiàn)在普遍認(rèn)為直徑在1100nm尺寸的顆粒屬納米粒子的范疇。這段尺寸的粒子的物理和化學(xué)性質(zhì)與大于100nm以上的粒子有著明顯的區(qū)別，但對(duì)其性質(zhì)遠(yuǎn)沒(méi)有深入研究。迄今人工合成的最新枝狀化合物的最大尺寸還只能達(dá)到10nm，而光刻的最小尺寸也只能接近100nm(Intel公司PentiumIII微處理器使用的光刻技術(shù)達(dá)到180nm)，膠體粒子和納米團(tuán)簇的尺寸大體位于這一間隙。因此納米團(tuán)簇的發(fā)現(xiàn)正為填補(bǔ)這段間隙的研究架起了橋梁。納米團(tuán)簇粒子的幾個(gè)主要的特性表現(xiàn)為: 表面效應(yīng)。當(dāng)固體粒子直徑小于100nm時(shí)，固體表面的特殊性質(zhì)開(kāi)始表現(xiàn)出來(lái)，這主要是由于表面原子數(shù)目占主體原子數(shù)目的比例開(kāi)始明顯升高，粒子越小，表面原子數(shù)占主體原子數(shù)的比例就越高。納米粒子表面的許多原子處于多個(gè)方位無(wú)原子接近的狀態(tài)，所以活性很高，易發(fā)生位置的移動(dòng)或與周圍的其他物質(zhì)發(fā)生作用。納米粒子的這一性質(zhì)已被廣泛用于催化反應(yīng)和表面修飾的研究中。久保(Kubo)效應(yīng)。Kubo效應(yīng)是指當(dāng)粒子尺寸下降到一定值時(shí)，費(fèi)米能級(jí)附近的電子能級(jí)由準(zhǔn)連續(xù)變?yōu)殡x散能級(jí)的現(xiàn)象。Kubo 最早提出了著名的公式:=(4/3)EF/N，式中為能級(jí)間距，EF為費(fèi)米能級(jí)，N為總電子數(shù)。對(duì)宏觀物體而言，N趨向于無(wú)限大，由上式可得0，即對(duì)宏觀物體能級(jí)間距趨向于零；而對(duì)納米粒子，所包含的原子數(shù)有限，N值很小，這就導(dǎo)致有一定的值，即能級(jí)間距發(fā)生分裂而呈現(xiàn)出Kube效應(yīng)。在金屬中，費(fèi)米能級(jí)附近的電子能級(jí)是連續(xù)的，隨著顆粒尺寸的減小，能級(jí)發(fā)生分裂，金屬納米粒子的性質(zhì)趨向于半導(dǎo)體。對(duì)能量的吸收和輻射表現(xiàn)出量子效應(yīng)。如對(duì)光的吸收，隨著粒子尺寸的減小，能級(jí)差增大，吸收光的波長(zhǎng)向短波方向移動(dòng)即所謂的藍(lán)移。穩(wěn)定性與幻數(shù)。許多原子團(tuán)簇中原子的個(gè)數(shù)不是任意的。對(duì)某一特定原子團(tuán)簇而言，只有原子數(shù)為某些固定值時(shí)，團(tuán)簇才有最高的穩(wěn)定性，原子團(tuán)簇以這些數(shù)目出現(xiàn)的概率最高。人們把相對(duì)穩(wěn)定的團(tuán)簇中所包括的原子個(gè)數(shù)(n)稱為幻數(shù)。不同元素的團(tuán)簇其幻數(shù)是不同的，即使同一元素，若制備方法不同其幻數(shù)也可能不一樣。2.2.2 團(tuán)簇組裝(1) 溶劑揮發(fā)自組裝法溶劑揮發(fā)自組裝法是利用納米團(tuán)簇或團(tuán)簇化合物間的弱的分子間力、靜電引力或它們與基體的作用力通過(guò)溶劑的揮發(fā)使其自發(fā)地組裝成規(guī)則、有序的結(jié)構(gòu)或超晶格。膠體在一定的條件下，有自組織或自組裝成有序結(jié)構(gòu)的能力。生命體中的許多結(jié)構(gòu)和實(shí)體均是由納米級(jí)的生物大分子組裝而成，如生物膜及生物膜中蛋白質(zhì)的自組裝、細(xì)胞中微管蛋白質(zhì)組裝成纖維狀的微管結(jié)構(gòu)等等。Bawendi等將修飾上三辛基膦的CdSe納米粒子在一定壓力和溫度下溶于辛烷和辛醇的混合液中，然后，通過(guò)降壓逐漸地?fù)]發(fā)辛烷，由于CdSe納米粒子在辛醇中的溶解度較小使得它從溶液中逐步析出，由此組裝出結(jié)構(gòu)非常勻稱、有規(guī)則的納米超晶格。Korgel等將表面修飾有C12巰基化合物的58nm銀膠溶液滴加在碳和云母基體上， 溶劑緩慢揮發(fā)后制得規(guī)則的單層膜。Motte等研究了Ag2S納米膠體在碳膜上進(jìn)行二維和三維自組裝，通過(guò)控制滴加和蒸發(fā)Ag2S膠體的速度，可以獲得非常規(guī)整的超晶格。Kiely等制備出表面組裝癸烷巰基化合物的金膠，將金膠溶液滴加在涂碳的銅網(wǎng)上，溶劑揮發(fā)后，通過(guò)掃描電鏡觀察到金膠化合物自組裝的3種排列方式:一是大小不同的顆?；旌暇鶆虻嘏帕谐呻p峰陣列；二是大顆粒和小顆粒各自自動(dòng)分開(kāi)形成緊密的六邊形結(jié)構(gòu)；三是大小不同的顆粒隨機(jī)地占據(jù)六邊形晶格的位點(diǎn)。Sun等在油酸和油烯基胺中合成直徑為310nm的鐵磁性納米粒子PtFe，調(diào)節(jié)Pt和Fe的比例可控制納米粒子的大小。將分散的納米粒子溶液滴加在涂有銅膜的SiO網(wǎng)膜上，溶劑揮發(fā)后，在TEM上觀察到三維的超晶格。這種材料可以支撐高密度的磁性反轉(zhuǎn)。(2)模板法模板法是利用納米粒子與模板間特定的作用力，如化學(xué)鍵、分子間作用力和靜電引力等對(duì)納米粒子進(jìn)行有序組裝。模板可以是固體(基片和多孔材料)、高分子膜、單分子膜、離子、分子、超分子和生物分子等。下面介紹其中的幾種模板法。分子模板法在基質(zhì)上進(jìn)行單分子層自組裝的研究在各方面都有比較深入的研究和成熟的方法。由于分子自組裝有較高的有序度和密度(特別是長(zhǎng)鏈分子的自組裝)，而且表層可以形成各種均勻的活性基團(tuán)，在此基礎(chǔ)上進(jìn)行納米粒子的自組裝自然是簡(jiǎn)便而又有效。自組裝的分子一般為雙官能團(tuán)分子即分子的兩端各有一特定功能的活性基團(tuán)，一端鍵合在基質(zhì)上，另一端與納米粒子鍵合。早在1992年，Colvin等通過(guò)雙官能團(tuán)分子將具有半導(dǎo)體性質(zhì)的CdS膠體組裝在涂有金或鋁的玻璃片上。在鋁基體上可以自組裝含羧基的分子，分子的另一端為巰基可與CdS鍵合。他們測(cè)定了表面覆蓋度以及組裝層的厚度，并研究了二維組裝膜的性質(zhì)。Andres等制備出平均粒徑為3.7nm的十二烷基硫醇組裝的金納米粒子，將金溶膠在MoS2基底上組裝，能自發(fā)形成(111)面的單層平行取向的膜。將此膜放入芳基雙硫醇的溶液中，芳基雙硫醇能將鄰近的納米粒子聯(lián)接起來(lái)形成二維超晶格結(jié)構(gòu)。通過(guò)改變金納米粒子的尺寸、組成、聯(lián)接分子的長(zhǎng)度和化學(xué)結(jié)構(gòu)以及基底的性質(zhì)，可以獲得很大范圍的電導(dǎo)變化。半導(dǎo)體CdS是一種光敏材料，在光催化、光電池等方面有廣泛的應(yīng)用。Hickey等用硅氧烷巰基化合物把5nm的CdS組裝在氧化錫電極上，并研究了光電流的性質(zhì)，測(cè)定了半導(dǎo)體CdS中空穴和電子的傳遞速率以及量子限和電位對(duì)電荷傳遞速度常數(shù)的影響。用烷基硫醇和巰丙酰甘氨酸修飾的金納米團(tuán)簇通過(guò)Cu2+的羧基化合物可以進(jìn)行可逆的組裝和聚集，這為制備新功能材料提供了新的方法。胡效亞等使用半胱胺分子(一端為巰基，另一端為氨基)在金電極表面組裝單分子層，再在此基礎(chǔ)上組裝單層金膠，研究了金膠表面的陰離子吸附行為和電化學(xué)性質(zhì)，在宏觀上觀察到納米金膠電化學(xué)性質(zhì)與金的區(qū)別，金膠能將金的許多吸附性質(zhì)“放大”。由于金膠與金電極接觸緊貼，所以在電化學(xué)性質(zhì)研究中有非常好的電子傳遞速度，從而獲得較理想的多方面的電化學(xué)信息。固態(tài)模板法Zhou等在室溫避光的高真空條件下(1.33Pa)，將孔徑為3nm的MCM241硅石在羰基釕原子簇的二氯甲烷溶液中浸泡72h，MCM241的顏色變成紅色或棕色，形成排列非常有序的一維陣列，對(duì)己烯和環(huán)辛烯的加氫反應(yīng)有極高的催化活性，并有良好的電學(xué)性質(zhì)。Blaaderen等用電子束在高分子薄膜上打出規(guī)則排布的孔洞，這些洞的深度和直徑與被組裝粒子相匹配，將分散于溶液中的微米粒子組裝到這些高分子薄膜的微孔中，通過(guò)選擇適當(dāng)混合溶劑和調(diào)節(jié)離子強(qiáng)度而使粒子一層層沉積在模板上，形成三維有序結(jié)構(gòu)。離子模板法離子模板法又稱離子自組裝，它是指帶正(負(fù))電荷的納米粒子或納米粒子化合物與帶負(fù)(正)電荷的離子、超分子、生物分子、納米粒子或納米粒子化合物通過(guò)靜電引力交替組裝成的超晶格。Liu等在合成金膠過(guò)程中加入陽(yáng)離子聚合物，制備出覆蓋陽(yáng)離子聚合物的納米金粒。在基片上先組裝一層陽(yáng)離子納米金粒，然后通過(guò)靜電引力再組裝一層陰離子聚合物，如此反復(fù)共組裝了15層，厚度為5nm，該膜的導(dǎo)電率與金的導(dǎo)電率在一個(gè)數(shù)量級(jí)上，而另外一些方法組裝的膜其導(dǎo)電率遠(yuǎn)低于它，這充分表現(xiàn)出超晶格的特性。這樣的超晶格在微電子器件制作方面有應(yīng)用背景。Liu等還將帶正電荷的TiO2納米粒子與帶負(fù)電荷的聚合物交替地、一層一層地組裝在玻璃和硅片上，制成的材料有望在光催化和非線性光學(xué)材料方面得到應(yīng)用。生物分子模板法將生物分子修飾到納米粒子表面，利用生物分子之間的識(shí)別功能，實(shí)現(xiàn)納米粒子的自組裝。Mirkin等52將一種低聚核苷酸修飾到一組金膠表面上，再將另一種低聚核苷酸修飾到另一組金膠表面上。當(dāng)兩組溶液混合時(shí)，兩種低聚核苷酸能相互識(shí)別和配對(duì)，形成二維或三維的組裝體系。當(dāng)溫度升高時(shí)，這種體系還能分解到最初狀態(tài)，也就是說(shuō)這是一種可逆的自組裝。排列點(diǎn)陣法使用掃描隧道顯微鏡(STM)和掃描探頭顯微鏡(SPMs)的探頭，在基片上堆積和排列納米粒子。這樣排列的納米粒子點(diǎn)陣具有排列規(guī)整，顆粒大小可以控制等特點(diǎn)，可以制成由幾個(gè)原子組成的原子簇。Kolb等將銅原子先還原在STM的探針上，當(dāng)探針接觸到金基片后，就在該處留下一粒銅的原子簇，反復(fù)操作就可以在金基片上有序地排列銅原子簇的點(diǎn)陣。每個(gè)原子簇高0.8nm，相當(dāng)于4層銅原子，用此方法還可以制作2層原子的點(diǎn)陣。Sugimura等在硅片上用掃描探頭顯微鏡(SPMs)制作有機(jī)硅烷單層的點(diǎn)陣，在此點(diǎn)陣上還可以自組裝其他分子，制成分子級(jí)器件。2.2.3 團(tuán)簇組裝體系的表征(1) X-射線衍射法在介孔材料的表征手段中，最常用的是X-射線衍射(XRD)技術(shù)。通過(guò)對(duì)X-射線衍射圖的解析，可獲得有關(guān)晶體的物質(zhì)組成、結(jié)構(gòu)(原子的三維立體坐標(biāo)、化學(xué)鍵、分子立體構(gòu)型和構(gòu)象、價(jià)電子云密度等)及分子間相互作用的信息。X-射線衍射也是測(cè)量納米晶尺寸的常用手段。它不僅可確定試樣物相及其相含量，還可判斷晶粒尺寸大小。在小角度散射區(qū)域內(nèi)(210)出現(xiàn)的衍射峰是確認(rèn)介孔結(jié)構(gòu)存在的有力判據(jù)之一。(2) 電鏡分析掃描電子顯微鏡(SEM)是觀察微小顆粒固體表面形態(tài)、結(jié)構(gòu)、物相及尺寸分布的有力工具。利用它可以檢查各步反應(yīng)前后材料的顯微鏡變化及微晶顆粒的破壞、腐蝕等情況。透射電子顯微鏡(TEM)是研究物質(zhì)結(jié)構(gòu)最直觀的方法，通過(guò)透射電子顯微鏡人們可以直接觀察物質(zhì)的微觀結(jié)構(gòu)。在這里利用透射電子顯微鏡可以觀察介孔材料中生成的納米團(tuán)簇的尺寸、形態(tài)及分布。對(duì)于介孔材料納米團(tuán)簇，高分辨電子顯微鏡(HREM)是最有力的觀察工具，它的分辨率為1.9，對(duì)孔徑只有幾個(gè)埃的沸石分子篩都能進(jìn)行有效的觀察，結(jié)合計(jì)算機(jī)模擬，形成原子像，可以得到介孔材料中納米團(tuán)簇的精細(xì)結(jié)構(gòu)。(3) 光譜分析紅外光譜(IR)是在電磁波紅外區(qū)(15000-10 cm-1)觀察物質(zhì)吸收和發(fā)射，以研究分子的振動(dòng)和轉(zhuǎn)動(dòng)光譜的譜學(xué)方法。它根據(jù)譜帶的特征頻率研究未知物成分(定性)，根據(jù)譜帶強(qiáng)度確定樣品中某個(gè)組分含量(定量)，它還可研究分子結(jié)構(gòu)(如官能團(tuán)、化學(xué)鍵)、鑒定異構(gòu)體判斷化合物結(jié)構(gòu)，利用譜帶變化還可研究分子間的相互作用。紅外光譜已廣泛地應(yīng)用于研究介孔材料及介孔材料中組裝的化合物作為紅外光譜的互補(bǔ)，拉曼光譜也是一種常用測(cè)試手段。由于拉曼散射信號(hào)弱，又常受到熒光干擾，因此提高拉曼光譜的檢測(cè)信噪比是拉曼光譜技術(shù)發(fā)展的重要途徑。表面增強(qiáng)拉受光譜法(SERS)解決了這個(gè)問(wèn)題。拉曼的表面增強(qiáng)，指的是由樣品分子處于某種特殊制備的金屬表面或表面附近而引起拉曼散射增加的特殊效應(yīng)。同傳統(tǒng)的拉曼散射譜比較SERS具有高達(dá)104，106倍的拉曼散射強(qiáng)度。除了普遍使用IR/Raman光譜外，還有熒光光譜、紫光一可見(jiàn)光譜等常見(jiàn)光譜學(xué)方法，普遍地應(yīng)用于納米材料、納米摻雜和納米表面修飾材料的組分鑒定，以及納米晶/介孔組裝體光學(xué)性質(zhì)的研究。(4) 其他表征方法擴(kuò)展X-射線吸收精細(xì)結(jié)構(gòu)光譜(EXAFS)近年來(lái)已成為分析缺少長(zhǎng)程有序體系的先進(jìn)技術(shù)。它能提供X-射線吸收邊之外所發(fā)射的精細(xì)光譜，對(duì)于難以得到單晶的材料尤其顯示獨(dú)到的優(yōu)越。用這種方法可得到有關(guān)配位原子種類、配位數(shù)、鍵長(zhǎng)或原子間距等X-射線吸收的原子化學(xué)環(huán)境方面的數(shù)據(jù)。EXAFS的最大優(yōu)點(diǎn)是對(duì)樣品有廣泛的適用性:晶態(tài)、非晶態(tài)、氣-液-固態(tài)物質(zhì)、玻璃態(tài)、表面分散態(tài)的配位結(jié)構(gòu)及表面吸附態(tài)等均可應(yīng)用，因此它是結(jié)構(gòu)表征的重要手段之一。在一些復(fù)合型介孔氧化物中，不同制備條件(如無(wú)機(jī)物種、合成溫度、原料配比等)會(huì)造成原子的組合方式發(fā)生一定變化，導(dǎo)致晶型錯(cuò)位、結(jié)構(gòu)重組。為了研究介孔結(jié)構(gòu)金屬原子的配位情況及組合方式，常借助于核磁共振(NRM)對(duì)材料進(jìn)行分析。2.2.4 應(yīng)用組裝納米團(tuán)簇的主要應(yīng)用有: 制作微小的電子器件。在集成線路以及其他微小電子器件中，目前使用的是光刻技術(shù)來(lái)制作元件和導(dǎo)線，對(duì)于尺寸遠(yuǎn)小于光波長(zhǎng)的器件，光刻技術(shù)已無(wú)法發(fā)揮作用。用幾個(gè)納米大小的顆粒排列組裝成的器件和導(dǎo)線(量子點(diǎn)和量子線)以及用納米碳管制作的分子器件，將是這一領(lǐng)域里最有希望和前景的方法和技術(shù)。這不僅可以提高運(yùn)算速度，而且可大大降低功耗。有序納米碳管的陣列組裝具有良好的場(chǎng)發(fā)射性質(zhì)，用它做成的顯示器與液晶顯示器相比具有顯示速度快、對(duì)比度高、不同角度觀察效果一致等優(yōu)點(diǎn)，被認(rèn)為是下一代的顯示器。制備超晶格材料。超晶格(Su2perlattices)材料是指在原有晶體晶格的周期結(jié)構(gòu)上疊加一個(gè)周期，這種周期結(jié)構(gòu)的勢(shì)阱(Potential well)區(qū)厚度小于電子平均自由程，勢(shì)壘(Potential hill)區(qū)足夠窄，以致相鄰勢(shì)阱中的電子波函數(shù)能夠互相耦合。晶體與超晶格的關(guān)系就如同蛋白質(zhì)的三級(jí)結(jié)構(gòu)與四級(jí)結(jié)構(gòu)的關(guān)系，超晶格與組成它的基本單元晶體相比具有許多新的特性和功能。提供新型的光學(xué)材料和非線性光學(xué)材料。如微波吸收性能的改變、增強(qiáng)拉曼光譜的散射等。制備新型的陶瓷材料、磁性材料和其他特異性材料。構(gòu)建新型催化材料。如組裝納米金電催化氧化CO。納米粒子表面的高活性是制備新型催化材料的基礎(chǔ)。設(shè)計(jì)和制備新型化學(xué)傳感器。能源的貯存。如電池的電極材料，貯存氫的納米碳管等。2.2.5 目前団簇科學(xué)的研究主要方向（1）研究團(tuán)簇的組成及電子構(gòu)型的規(guī)律、幻數(shù)和幾何結(jié)構(gòu)、穩(wěn)定性的規(guī)律；（2）研究團(tuán)簇的成核和形成過(guò)程及機(jī)制，研究團(tuán)簇的制備方法、尤其是獲取尺寸均一與可控的團(tuán)簇束流；（3）研究金屬、半導(dǎo)體及非金屬和各種化合物團(tuán)簇的光、電、磁、力學(xué)、化學(xué)等性質(zhì)，它們與結(jié)構(gòu)和尺寸的關(guān)系，及向大塊物質(zhì)轉(zhuǎn)變的關(guān)節(jié)點(diǎn)；（4）研究團(tuán)簇材料的合成和性質(zhì)；（5）探索新的理論，不僅能解釋現(xiàn)有團(tuán)簇的效應(yīng)和現(xiàn)象，而且能解釋和預(yù)知團(tuán)簇的結(jié)構(gòu)，模擬團(tuán)簇動(dòng)力學(xué)性質(zhì)，指導(dǎo)實(shí)驗(yàn)；（6）發(fā)展新的方法對(duì)團(tuán)簇表面進(jìn)行修飾和控制。3 分子動(dòng)力學(xué)基本理論3.1 理論基礎(chǔ)分子動(dòng)力學(xué)方法(Molecular Dynamics Method，簡(jiǎn)稱MD)是最常見(jiàn)的一種計(jì)算模擬方法。它主要是根據(jù)牛頓力學(xué)來(lái)模擬原子或者分子體系的運(yùn)動(dòng)，假定原子的運(yùn)動(dòng)遵循經(jīng)典的力學(xué)定律，即：在N個(gè)原子所構(gòu)成的體系里，每個(gè)原子遵循牛頓第二定律， (2.1)是原子質(zhì)量，是其加速度，是由其他原子作用在該原子的力。利用原子或分子的自然運(yùn)動(dòng)在相空間中抽取樣本，從而進(jìn)行統(tǒng)計(jì)計(jì)算。計(jì)算體系的相積分，并進(jìn)一步計(jì)算體系的熱力學(xué)量和其他宏觀性質(zhì)19。它的結(jié)果主要取決于原子間相互作用勢(shì)的選擇20。所以，MD方法是一種確定性算法，給定原子的初始位置和速度后，就能完全確定體系隨時(shí)間的演變。MD在不同的科學(xué)領(lǐng)域里特別是在材料科學(xué)的研究中均有廣泛的應(yīng)用，扮演了非常重要的角色。例如在研究金屬的固液界面、固固界面熱力學(xué)和動(dòng)力學(xué)方面，都是有力的工具21,22。宏觀條件下，具有完全相同的性質(zhì)和結(jié)構(gòu)、處于各種運(yùn)動(dòng)狀態(tài)卻又各自獨(dú)立的系統(tǒng)的集合，稱為系綜。系綜是應(yīng)用統(tǒng)計(jì)方法描述熱力學(xué)系統(tǒng)的統(tǒng)計(jì)規(guī)律而引入的概念。按照約束條件的不同，系綜可以分為微正則系綜(NVE)，正則系綜(NVT)，等溫等壓系綜(NPT)，等壓等焓系綜(NPH)和巨正則系綜(VT )幾類。下面對(duì)分子動(dòng)力學(xué)模擬中一些常用的系綜進(jìn)行簡(jiǎn)要的介紹。1. 微正則系綜在 NVE中，整個(gè)體系與外界無(wú)能量交換，也無(wú)粒子交換。體系的粒子數(shù)N，體積V和能量E守恒。因此，微正則系綜所描述的體系就是一個(gè)孤立體系。體系達(dá)到平衡后，其總能量也恒定，溫度和壓強(qiáng)值在一定范圍內(nèi)波動(dòng)23。平衡狀態(tài)下的孤立體系中，一切可能的粒子微觀運(yùn)動(dòng)狀態(tài)出現(xiàn)的幾率都相同，且不隨時(shí)間而變。這個(gè)假設(shè)被稱為等概率原理。在微正則系綜中，微觀狀態(tài)出現(xiàn)在相空間各個(gè)區(qū)域的概率和為l。由于無(wú)法用某一個(gè)微觀狀態(tài)量來(lái)表示體系的微觀狀態(tài)量，因此只能用微觀量在一切可能的微觀狀態(tài)上的統(tǒng)計(jì)平均值來(lái)表示。所以，若我們假設(shè)微觀狀態(tài)處在相空間某個(gè)區(qū)域的概率為 (p，g，f)，對(duì)于一個(gè)我們想求的微觀量X，它在一切可能的微觀狀態(tài)上的平均值就等于： (2.2)其中，它表示相空間中的單位體積元，表示微觀狀態(tài)分布函數(shù)，表示對(duì)應(yīng)于相空間中每一點(diǎn)的微觀量。利用等概率原理，只要知道體系的微觀狀態(tài)數(shù)，就可以求出體系的熵。原則上我們可以計(jì)算出包括體系能量、壓強(qiáng)、溫度在內(nèi)的全部熱力學(xué)函數(shù)。2. 正則系綜具有確定粒子數(shù)N，體積V以及溫度T的體系被稱為正則系綜(NVT)。為了保持系統(tǒng)的溫度恒定，若假設(shè)系統(tǒng)與一巨大熱源相接觸，與之發(fā)生能量交換，直至達(dá)到熱平衡。由于熱源足夠大，其溫度始終保持恒定，并不因和體系間發(fā)生了熱交換而發(fā)生改變。由此可知，正則系綜不是一個(gè)孤立體系，其中可能的粒子微觀運(yùn)動(dòng)狀態(tài)出現(xiàn)的幾率是不相同的，因此等概率原理不適用。但是，由于體系和熱源所形成的系統(tǒng)是一個(gè)能量確定的孤立系統(tǒng)，等概率原理對(duì)這個(gè)復(fù)合系統(tǒng)是適用的。所以，我們可以利用類似于微正則系綜中處理孤立體系的做法，求出復(fù)合系統(tǒng)的微觀狀態(tài)分布和熵，進(jìn)一步求得所研究體系的微觀狀態(tài)分布。因此，正則系綜中體系的熱力學(xué)函數(shù)也可以求得。3. 等溫等壓系綜確定粒子數(shù)N，壓強(qiáng)P以及溫度T的體系被稱為等溫等壓系綜(NPT)。所研究的體系與外界之間沒(méi)有粒子交換。為了保持系統(tǒng)溫度恒定，我們采用類似正則系綜中的處理方法，假設(shè)體系和一個(gè)巨大的熱源相接觸并與之發(fā)生能量轉(zhuǎn)換，直至達(dá)到熱平衡；同時(shí)，等溫等壓系綜中還要求體系的壓強(qiáng)恒定，所以還需要假設(shè)體系和一個(gè)壓浴相耦合，并與之發(fā)生壓強(qiáng)傳遞直至平衡。等溫等壓體系的研究對(duì)實(shí)驗(yàn)來(lái)說(shuō)相當(dāng)重要，因?yàn)楹芏鄬?shí)驗(yàn)都是在等溫等壓的情況下進(jìn)行的。本論文就是采用NPT處理數(shù)據(jù)，進(jìn)而研究銅原子填充碳納米管復(fù)合結(jié)構(gòu)的穩(wěn)定性。3.2 分子動(dòng)力學(xué)的實(shí)現(xiàn)及常用的技術(shù)分子動(dòng)力學(xué)模擬步驟大致有以下幾步：1.讀取模擬參數(shù)，如溫度，壓強(qiáng)，時(shí)間步長(zhǎng)，模擬時(shí)間總長(zhǎng)度。值得說(shuō)明的是，時(shí)間步長(zhǎng)，就是積分的間隔，因而時(shí)間步長(zhǎng)的選取對(duì)動(dòng)力學(xué)模擬非常重要。太長(zhǎng)的時(shí)間步長(zhǎng)會(huì)造成原子或分子間的激烈碰撞，導(dǎo)致體系數(shù)據(jù)溢出；太短的時(shí)間步長(zhǎng)會(huì)降低模擬過(guò)程搜索相空間的能力。因此，一般選取的時(shí)間步長(zhǎng)為體系各個(gè)自由度中最短運(yùn)動(dòng)周期的十分之一。2.初始化體系。其中包括周期性邊界條件，體系原子數(shù)，原子初始位置和速度。3.結(jié)合原子勢(shì)計(jì)算每個(gè)原子受到的作用力。4.使用數(shù)值積分方法，求解牛頓運(yùn)動(dòng)方程。5.統(tǒng)計(jì)體系的熱力學(xué)量。MD的核心部分是對(duì)運(yùn)動(dòng)方程的時(shí)間積分。其基本思想是利用有限差分法對(duì)二階常微分方程進(jìn)行積分。通常把模擬時(shí)間離散為有限格點(diǎn)，相鄰的格點(diǎn)距離稱為模擬的時(shí)間步長(zhǎng)。初始時(shí)刻t的原子坐標(biāo)以及速度，加速度已知，就可以計(jì)算出下一時(shí)刻t+的坐標(biāo)以及速度等量值。MD方法中常見(jiàn)的數(shù)值積分方法有Verlet算法24，蛙跳算法，以及預(yù)測(cè)-校正算法。MD模擬中常用的物理量包括體系的動(dòng)能，勢(shì)能，總能，溫度和壓強(qiáng)等。1.周期性邊界條件在MD計(jì)算中，人們往往利用周期性邊界條件來(lái)解決模擬的體系尺寸遠(yuǎn)小于實(shí)際物質(zhì)的尺寸的問(wèn)題。在使用周期性邊界條件時(shí)，粒子被限制在一個(gè)模擬盒子中(稱為模擬元胞)，假設(shè)這個(gè)模擬元胞沿三個(gè)坐標(biāo)方向無(wú)限復(fù)制，直至完全包括整個(gè)空間。模擬元胞中的 粒子不僅可以與其他粒子 發(fā)生作用，也和近鄰元胞中的粒子發(fā)生作用，即相互作用可以穿透元胞的邊界。顯然，應(yīng)用了周期性邊界條件之后，消除了表面效應(yīng)，并且此時(shí)模擬元胞的邊界位置沒(méi)有作用了。2.優(yōu)化算法在MD計(jì)算中，經(jīng)常需要利用體系的勢(shì)能極小值。因此，人們常用優(yōu)化算法來(lái)找尋體系的最優(yōu)結(jié)構(gòu)，搜索體系的勢(shì)能極小值。目前，兩種優(yōu)化算法應(yīng)用較為廣泛，其一是最速下降法，另一則為共軛梯度算法。本文中使用的是最速下降法，即沿著力的方向向更低能量態(tài)進(jìn)行搜索，直到體系的能量最終收斂到極小值。在最速下降法中，每次選擇的搜索方向，僅是由當(dāng)前函數(shù)的狀態(tài)決定的。3.3原子間相互作用勢(shì)在一定的物理模型的基礎(chǔ)上發(fā)展相應(yīng)的原子間相互作用勢(shì)，進(jìn)而研究材料的性質(zhì)和不同狀態(tài)下的行為，成為材料研究中一種必要的研究手段。在分子動(dòng)力學(xué)模擬中，粒子之間的相互作用力是由粒子之間的勢(shì)函數(shù)對(duì)位置求偏導(dǎo)得到的。因此，作用勢(shì)的選擇與MD計(jì)算的關(guān)系極為密切。選擇不同的作用勢(shì)，體系的勢(shì)能面會(huì)有不同的形狀，動(dòng)力學(xué)計(jì)算所得的分子運(yùn)動(dòng)和分子內(nèi)部運(yùn)動(dòng)的軌跡也會(huì)不同，進(jìn)而影響到抽樣的結(jié)果和依據(jù)抽樣結(jié)果的勢(shì)能計(jì)算。在MD模擬發(fā)展的初期，人們普遍使用的是對(duì)勢(shì)。它只描述兩個(gè)原子之間的相互作用，主要應(yīng)用在液體、氣體的模擬計(jì)算中，另外也經(jīng)常用于對(duì)分子動(dòng)力學(xué)擬技術(shù)本身的研究中。近期的工作大多選取多體勢(shì)。與對(duì)勢(shì)僅考慮兩兩原子間的相互作用不同，多體勢(shì)認(rèn)為：在多原子系統(tǒng)中，一個(gè)原子的位置改變，必將影響到整體系統(tǒng)中一定范圍的電子云分布，從而影響到其他原子間的相互作用。因此，多原子體系的勢(shì)函數(shù)用多體勢(shì)來(lái)表述更為準(zhǔn)確。多體勢(shì)函數(shù)在描述金屬及共價(jià)鍵方面，有十分明顯的優(yōu)勢(shì)。通過(guò)前面介紹的MD模擬方法可以看出，勢(shì)能函數(shù)的選取是材料的計(jì)算機(jī)模擬的核心問(wèn)題。勢(shì)能函數(shù)選取的好壞直接關(guān)系到模擬結(jié)果的正確與否。勢(shì)能函數(shù)總是表示成粒子坐標(biāo)的數(shù)學(xué)函數(shù)形式，它給出體系總能量與粒子坐標(biāo)的關(guān)系。如果這個(gè)函數(shù)對(duì)體系的描述越準(zhǔn)確，那么用它來(lái)計(jì)算所得到的物理體系的性質(zhì)就與真實(shí)情況卻接近。一般在描述原子間短程作用勢(shì)時(shí)要采用多體勢(shì)，因?yàn)楝F(xiàn)代的半經(jīng)驗(yàn)多體勢(shì)己經(jīng)需要考慮在短程的范圍內(nèi)環(huán)境對(duì)原子的總體作用，比如Airebo勢(shì)為CC成鍵相互作用，截至距離為2；而對(duì)于需要考慮長(zhǎng)程作用的原子體系來(lái)說(shuō)，還要在短程的基礎(chǔ)上結(jié)合兩體勢(shì)，比如經(jīng)典的LJ勢(shì)25,26和后來(lái)的Exp-6勢(shì)。C-C原子間的相互作用以Airebo勢(shì)來(lái)描述，Cu-Cu原子間的相互作用以EAM多體勢(shì)來(lái)描述，C-Cu原子的相互作用以LJ勢(shì)來(lái)描述。下面重點(diǎn)介紹LJ勢(shì)和EAM勢(shì)：1.LJ勢(shì) 碳原子和銅原子之間的作用為范德瓦耳斯力(Van der Waals Force)，Lellnard-Jones(LJ)勢(shì)是計(jì)算化學(xué)中用來(lái)模擬兩分子間勢(shì)能的一個(gè)函數(shù)，最早是由John Edward Lennard-Jones在1931年提出的。它是對(duì)勢(shì)模型中的典型，在理論和實(shí)驗(yàn)研究領(lǐng)域都具有重要意義。由于其解析形式簡(jiǎn)單而被廣泛使用，特別是用來(lái)描述惰性氣體分子間的相互作用尤為精確。LJ勢(shì)以兩體距離為唯一變量，包含兩個(gè)參數(shù)。其形式為： (2.3) 其中， 等于勢(shì)能井的深度， 是相互作用的勢(shì)能正好為零時(shí)的兩體距離。在實(shí)際應(yīng)用中， 、參數(shù)往往通過(guò)擬合已知實(shí)驗(yàn)數(shù)據(jù)或精確量子計(jì)算結(jié)果而確定。2.EAM勢(shì)多數(shù)針對(duì)金屬的多體勢(shì)是以有效介質(zhì)為基礎(chǔ)建立的。常用的包括： 有效介質(zhì)理論(EMT)、嵌入原子方法(EAM)、緊束縛方法(TB)、二階動(dòng)量近似(SMA)和膠體模型(Glue Model)等。這其中以EAM的應(yīng)用最為廣泛。本論文采用的即是EAM勢(shì)。下面我們將對(duì)EAM勢(shì)進(jìn)行簡(jiǎn)單介紹。EAM是一類多體相互作用的統(tǒng)稱。它的物理思想起源于Friedel提出的原子嵌入能概念。原子嵌入能的基本思想是：原子的凝聚能主要取決于該原子所占據(jù)位置的局域電子密度。EAM勢(shì)是Daw和BaSkes在結(jié)合電子的密度泛函理論和有效介質(zhì)理論發(fā)展起來(lái)的一種構(gòu)造原子間相互作用勢(shì)的方法。在EAM下，可以將體系的總能分為兩個(gè)部分：一部分是所研究的原子的嵌入能；另一部分是其余的所有構(gòu)成固體的原子之間的斥能。 (2.4) (2.5) 公式(2-4)中，第一項(xiàng)是原子的嵌入能(多體相互作用)，第二項(xiàng)是原子 和之間兩體相互作用部分。得到EAM能量表達(dá)式的過(guò)程中最主要用到了兩個(gè)近似：1.體系的總能量可以表示為嵌入能和二體斥能之和；2.體系的電荷分布近似為各個(gè)原子電荷的簡(jiǎn)單線性疊加。此近似忽略了原子間成鍵時(shí)電荷的轉(zhuǎn)移和重新分布，但是對(duì)于大多數(shù)金屬而言，此關(guān)系近似成立。Foiles等人對(duì)金屬EAM勢(shì)的構(gòu)造做了重要改進(jìn)，并給出了Au、Ag、Cu、Pt、Pd、Ni等金屬的EAM勢(shì)，使之不僅適用于純金屬材料的研究，而且可以對(duì)復(fù)合材料給出比較理想的描述。在EAM模型中，各參數(shù)是用晶格常數(shù)、內(nèi)聚能、彈性常數(shù)和空位形成能等實(shí)驗(yàn)可測(cè)的量擬合得到的。 本文中的研究體系所選取的原子勢(shì)是EAM勢(shì)。采用EAM勢(shì)的優(yōu)點(diǎn)是其更適合對(duì)金屬及合金的模擬，并且，由于EAM勢(shì)易于分析，為得到理想的熱力學(xué)性質(zhì)進(jìn)行調(diào)整也十分方便。4 Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)模擬方法4.1 相關(guān)軟件介紹1. MS簡(jiǎn)介Materials Studio是專門為材料科學(xué)領(lǐng)域研究者開(kāi)發(fā)的一款可運(yùn)行在PC上的模擬軟件。它可以幫助研究者們能更方便地建立三維結(jié)構(gòu)模型，并對(duì)各種晶體、無(wú)定型以及高分子材料的性質(zhì)及相關(guān)過(guò)程進(jìn)行深入的研究。Materials Studio生成的結(jié)構(gòu)、圖表及視頻片斷等數(shù)據(jù)可以及時(shí)地與其它PC軟件共享，使任何一位研究者都能充分享用Materials Studio軟件所帶來(lái)的先進(jìn)技術(shù)。Materials Studio軟件能使任何研究者達(dá)到與世界一流研究部門相一致的材料模擬的能力。模擬的內(nèi)容包括了催化劑、聚合物、固體及表面、晶體與衍射、化學(xué)反應(yīng)等材料和化學(xué)研究領(lǐng)域的主要課題。2. LAMMPS計(jì)算軟件LAMMPS(Large-scale Atomic/Molecular Massively Parallel Simulation)一般來(lái)講，只要是分子動(dòng)力學(xué)所涉及到的領(lǐng)域，LAMMPS代碼也都涉及到了。LAMMPS可以支持包括氣態(tài)、液態(tài)或者固態(tài)相形態(tài)和各種系綜下百萬(wàn)級(jí)的原子分子體系，并提供多種勢(shì)函數(shù)支持。LAMMPS的另一個(gè)優(yōu)點(diǎn)是具有良好的并行擴(kuò)展性。模擬in文件共包括四部分：初始化(initialization)，建模(lattice configuration setup)，模擬控制(simulation control)以及輸出控制(output control)。3. VMD軟件分子動(dòng)力學(xué)可視化VMD(Visual MD)是分子動(dòng)力學(xué)可視化程序。它使用Open GL提供高質(zhì)量的3D分子圖形，顯示、動(dòng)畫(huà)演示和分析大的原子、分子體系。在內(nèi)存允許的范圍內(nèi)，它可以無(wú)限制的顯示原子、分子數(shù)量，以及動(dòng)畫(huà)畫(huà)面的數(shù)量。VMD提供多種分子顯示方式和著色方法而且具有立體顯示功能。此外，在我們的模擬過(guò)程中，VMD能幫助我們更好的分析銅原子鏈填充碳納米管復(fù)合結(jié)構(gòu)的穩(wěn)定性。4.2 模擬計(jì)算過(guò)程本研究先在MS里建立銅原子鏈填充碳納米管模型，用LAMMPS軟件計(jì)算模擬，分析結(jié)果，并在VMD中改變?cè)O(shè)置，觀察圖像。本文為研究管直徑、管長(zhǎng)度、銅原子數(shù)目對(duì)復(fù)合結(jié)構(gòu)穩(wěn)定性的影響，建立了(8,0)、(9,0)、(10,0)和(11,0)四種碳管，106.50、200.22、298.20碳管長(zhǎng)度及不同銅原子數(shù)目的模型進(jìn)行研究。本文采用MD研究銅填充碳管的過(guò)程及填充后碳管中的銅原子排布情況。其中，C-C原子間的相互作用以Airebo勢(shì)來(lái)描述，Cu-Cu原子間的相互作用以EAM多體勢(shì)來(lái)描述，C-Cu原子的相互作用以LJ勢(shì)來(lái)描述。本文在等溫等壓系綜(NPT)下進(jìn)行，溫度先是在2000K平衡10000步，再?gòu)?000K到300K以5K/500ps速度緩慢降溫步，時(shí)間步長(zhǎng)為1fs，每1000步記錄一次構(gòu)型；再由300K到0K以5K/500ps速度最速下降，每一溫度平衡500ps，將數(shù)據(jù)整理分析對(duì)比，得出結(jié)果。 5 Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)的分子動(dòng)力學(xué)研究結(jié)果及分析5.1 Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)圖1、圖2和圖3是在VMD中觀察到的不同尺寸的銅鐵原子中不同的原子比例的團(tuán)簇結(jié)構(gòu)的側(cè)視圖，斜視圖，如圖所示，圖1是尺寸為108的銅鐵原子中銅鐵原子的組成分別為0.0、0.2、0.4、0.6、0.8和1.0時(shí)的Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)，圖2是尺寸為365的銅鐵原子中銅鐵原子的組成分別為0.0、0.2、0.4、0.6、0.8和1.0時(shí)的Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)，圖1是尺寸為864的銅鐵原子中銅鐵原子的組成分別為0.0、0.2、0.4、0.6、0.8和1.0時(shí)的Cu-Fe二元合金納米團(tuán)簇結(jié)構(gòu)。 圖一（1） 108-0.0 圖一（2） 365-0.0 圖一（3） 864-0.0 圖二（1） 108-0.2 圖二（2） 365-0.2 圖二（3） 864-0.2 圖三（1） 108-0.4 圖三（2） 365-0.4 圖三（3） 864-0.4 圖四（1） 108-0.6 圖四（2） 365-0.6 圖四（3） 864-0.6 圖五（1） 108-0.8 圖五（2） 365-0.8 圖五（3） 864-0.8 圖六（1） 108-1.0 圖六（2） 365-1.0 圖六（3） 864-1.0 不同的銅鐵原子尺寸中各個(gè)原子組成的銅鐵二元納米團(tuán)簇結(jié)構(gòu)的正視圖m-n中，m表示銅鐵原子尺寸，n表示銅鐵原子組成由圖一、圖二和圖三可知，銅鐵原子二元合金團(tuán)簇結(jié)構(gòu)的穩(wěn)定結(jié)構(gòu)呈立方體形狀，銅鐵原子分布局部集中，交替排列；隨著原子尺寸的增大，原子排列切面與合金團(tuán)簇結(jié)構(gòu)的表面夾角逐漸增大。0.00.20.40.60.81.0108-379.63339-388.37112-396.93775-413.46008-425.22492-441.18676365-1286.4441-1305.7704-1346.9768-1370.2737-1380.2392-1501.4639864-3035.5968-3095.3011-3198.3809-3310.6201-3408.8707-3549.4957表1 不同銅鐵原子尺寸（m+n）中不同原子組成(m/n)與原子勢(shì)能的關(guān)系其函數(shù)曲線圖一如下：由表1和曲線圖一可知，相同原子尺寸中銅鐵原子組成越小，及鐵原子含量越多，其原子勢(shì)能越小，合金團(tuán)簇結(jié)構(gòu)越穩(wěn)定；相同原子組成中原子尺