孫凌濤，石東平

(重慶文理學(xué)院新材料技術(shù)研究院，重慶 402160)

Cu-Co雙金屬團(tuán)簇結(jié)構(gòu)演化及其性質(zhì)的分子動(dòng)力學(xué)模擬

孫凌濤，石東平

(重慶文理學(xué)院新材料技術(shù)研究院，重慶 402160)

采用分子動(dòng)力學(xué)結(jié)合嵌入原子方法對(duì)比研究了Co分布于Cu-Co團(tuán)簇不同層的結(jié)構(gòu)和性質(zhì).研究表明：Co原子分層摻雜可對(duì)團(tuán)簇的結(jié)構(gòu)轉(zhuǎn)變點(diǎn)和熔點(diǎn)進(jìn)行誘導(dǎo)控制；分層摻雜的Cu-Co團(tuán)簇第一相變是一種擴(kuò)散度較小的由立方八面體轉(zhuǎn)變?yōu)槎骟w的相變；Co原子易于向低能態(tài)團(tuán)簇的亞表層(111)面偏析，從而誘導(dǎo)團(tuán)簇結(jié)構(gòu)紊亂，造成其熔點(diǎn)差異.

團(tuán)簇; 結(jié)構(gòu); 偏析; 分子動(dòng)力學(xué)

1 引 言

介于原子、分子與塊體之間的一種物質(zhì)新層次，團(tuán)簇因較大的表體比引起的表面效應(yīng)，會(huì)表現(xiàn)出許多特殊的物理化學(xué)性質(zhì)[1].由兩種金屬原子聚集而成的雙金屬團(tuán)簇，可以通過異質(zhì)原子成分、尺寸和排列方式對(duì)團(tuán)簇的物理、化學(xué)和電磁性能進(jìn)行調(diào)節(jié)[2，3].在雙金屬納米團(tuán)簇展現(xiàn)的眾多特殊功能中，熱力學(xué)行為表現(xiàn)出許多奇特性質(zhì)，如升溫過程中團(tuán)簇表現(xiàn)出尺寸效應(yīng)[4]和分層預(yù)融[5]等.這使得雙金屬團(tuán)簇在傳感、納米催化和微電子等新材料領(lǐng)域中具有廣闊應(yīng)用前景，相關(guān)理論[6]、實(shí)驗(yàn)[7]研究成為近年熱點(diǎn).

對(duì)于雙金屬團(tuán)簇，為了釋放內(nèi)應(yīng)力，小原子比大原子有更強(qiáng)烈的向核分布趨勢(shì)，這種效應(yīng)在二十面體團(tuán)簇中較為強(qiáng)烈，而在FCC和十面體團(tuán)簇中較弱.另外，為了降低粒子的表面能，雙金屬團(tuán)簇中表面能較小的原子易于向表面偏析[8].因此任意兩種金屬原子由于原子尺寸和表面能的差異會(huì)使團(tuán)簇在熱力學(xué)過程中出現(xiàn)偏析現(xiàn)象，從而誘發(fā)團(tuán)簇結(jié)構(gòu)形貌、表面組分及性質(zhì)發(fā)生改變[9].最近的理論和實(shí)驗(yàn)研究發(fā)現(xiàn)，原子偏析會(huì)使一些雙金屬體系易于形成表面能小的原子包裹表能較大原子的核殼結(jié)構(gòu)[10，11].例如，通過溫控，Ag原子偏析會(huì)誘導(dǎo)Co-Ag團(tuán)簇形成具有近似hcp核的核殼結(jié)構(gòu)[12]，也可以形成具有Ag-Co-Ag的洋蔥型結(jié)構(gòu)[13].通過在Au團(tuán)簇表面包裹Ag原子，可以大幅降低Ag-Au團(tuán)簇馬氏體轉(zhuǎn)變溫度[14].以不同比例摻雜的Cu-Co雙金屬團(tuán)簇以一定冷卻速率降溫后，由于原子偏析會(huì)誘導(dǎo)其形成特定的可控結(jié)構(gòu)[15，16].以Co原子為基體的Co-Rh團(tuán)簇，部分Co原子偏析到表面，另外一部分仍占據(jù)于核心層，這將誘導(dǎo)團(tuán)簇形成洋蔥型結(jié)構(gòu)[17].目前原子偏析作為一個(gè)研究側(cè)入點(diǎn)，主要側(cè)重于對(duì)雙金屬團(tuán)簇降溫過程中團(tuán)簇結(jié)構(gòu)形貌的研究，但對(duì)升溫過程中原子偏析誘導(dǎo)的團(tuán)簇結(jié)構(gòu)、形貌、相變點(diǎn)變化的作用細(xì)節(jié)和機(jī)理還不夠深入.另外，目前的研究多采用表面能小的原子隨機(jī)摻雜或按層分布于表面能較大的基體團(tuán)簇，來考察表面能小的原子向雙金屬團(tuán)簇表面偏析對(duì)團(tuán)簇的結(jié)構(gòu)和性質(zhì)的影響，對(duì)表面能大的原子隨機(jī)或分層摻雜于表面能小的團(tuán)簇中研究還較少.

針對(duì)目前研究中存在的問題，本文以表面能相差較大的Cu(1592 mJ/m2)和Co(2197 mJ/m2)[18]原子，構(gòu)建了Co原子分布于基體Cu團(tuán)簇不同層的雙金屬團(tuán)簇.然后，采用分子動(dòng)力學(xué)模擬方法，并結(jié)合勢(shì)能-溫度曲線、團(tuán)簇快照?qǐng)D、對(duì)分布函數(shù)及Honeycutt-Andersen(HA)指數(shù)[19]等不同表征分析手段，對(duì)升溫過程中Co分層摻雜Cu-Co團(tuán)簇的結(jié)構(gòu)及性質(zhì)進(jìn)行了研究.

2 計(jì)算方法

團(tuán)簇模型采用五階Gear預(yù)報(bào)-校正算法積分牛頓方程，用標(biāo)度方法控制系統(tǒng)的溫度.模擬系統(tǒng)使用正則系綜(NVT)且無外力的周期性邊界條件.由于截?cái)喟霃竭x取的差異，模擬純金屬團(tuán)簇的勢(shì)函數(shù)不適合模擬異質(zhì)二元合金團(tuán)簇，本文使用由Zhou等[20，21]提出的一種合金原子間的相互作用—鑲嵌原子勢(shì)(EAM)，它采用了和嵌入原子法相同的分析形勢(shì)，可用來研究16種任意組合金的相互作用.團(tuán)簇?zé)崃W(xué)模擬過程將團(tuán)簇從300 K連續(xù)升溫至1400 K，升溫速率0.6 K/ps，模擬時(shí)間步長(zhǎng)0.5 fs.含有309個(gè)原子初始團(tuán)簇是從10a0×10a0×10a0(晶格常數(shù)a0=0.3615)的FCC塊體中截取，然后分別用Co原子替換立方八面體Cu團(tuán)簇亞表層及次表層的Cu原子所得.Co分布于Cu團(tuán)簇次亞表層(簡(jiǎn)稱：次表Co團(tuán)簇)和亞表層(簡(jiǎn)稱：亞表Co團(tuán)簇)的初始團(tuán)簇的結(jié)構(gòu)形貌及其Cu和Co原子的徑向密度函數(shù)圖1中(a)、(b)所示.

圖1 團(tuán)簇中Cu和Co原子的徑向分布密度及各層的表面形貌Fig.1 Atomic radius distributions of the clusters and their corresponding initial morphologies

3 結(jié)果與討論

首先給出了升溫過程中純Cu及其兩種摻雜團(tuán)簇溫度-能量曲線，如圖2所示.可以看出，隨著溫度的升高，三種組分的團(tuán)簇存在相同的變化趨勢(shì)，都會(huì)在低溫時(shí)出現(xiàn)勢(shì)能向下跳躍點(diǎn)，在高溫時(shí)出現(xiàn)勢(shì)能向上跳躍點(diǎn)，在跳躍點(diǎn)前后團(tuán)簇的能量與溫度呈近似線性關(guān)系.所不同的是純Cu團(tuán)簇和次表Co團(tuán)簇的第一個(gè)勢(shì)能拐點(diǎn)出現(xiàn)在300-400 K之間，而亞表Co團(tuán)簇的第一勢(shì)能拐點(diǎn)出現(xiàn)在400-500 K之間；純Cu團(tuán)簇、次表Co團(tuán)簇及亞表Co團(tuán)簇的第二個(gè)拐點(diǎn)分別出現(xiàn)在800-900 K、900-1000 K、1000-1100 K之間.從圖2縱向來看，純Cu團(tuán)簇的平均勢(shì)能整體高于亞表Co及次表Co團(tuán)簇，亞表Co團(tuán)簇能量相對(duì)較低.

圖2 升溫過程中團(tuán)簇原子的溫度-能量曲線Fig.2 Variation of energy with temperature for clusters

圖2中勢(shì)能拐點(diǎn)位置及其能量縱向差異說明Co原子分層摻雜會(huì)直接影響到團(tuán)簇的結(jié)構(gòu)轉(zhuǎn)變、熔點(diǎn)及能量.用勢(shì)能-溫度曲線來確定團(tuán)簇的相變點(diǎn)非常直觀，但是由于勢(shì)能拐點(diǎn)的跳躍區(qū)間覆蓋了一個(gè)范圍，無法精準(zhǔn)鎖定某一溫度值.基于此，還需用溫度-熱熔曲線來精確化結(jié)構(gòu)相變點(diǎn)，如圖3所示.可以看出三種組分團(tuán)簇的相變點(diǎn)位置存在較大的差異，亞表Co團(tuán)簇、次亞Co團(tuán)簇及純Cu團(tuán)簇第一相變點(diǎn)分別為428 K、370 K和377 K，相應(yīng)的第二相變點(diǎn)分變?yōu)?008 K、934 K及853 K.

圖3 團(tuán)簇的熱容-溫度曲線Fig.3 Temperature dependence of the heating capacity for clusters

為了考察團(tuán)簇相變點(diǎn)的團(tuán)簇結(jié)構(gòu)變化，圖4給出了亞表Co團(tuán)簇分別在第一相變點(diǎn)前后的對(duì)分布函數(shù)(圖中a、b所示)和次表Co團(tuán)簇在第一相變點(diǎn)前后的對(duì)分布函數(shù)(圖中c、d所示).可以看出，亞表Co團(tuán)簇轉(zhuǎn)變前后對(duì)分布函數(shù)峰位、峰強(qiáng)度都有著明顯的差異，相對(duì)于轉(zhuǎn)變前來說，轉(zhuǎn)變后的第三峰明顯減小，從第四峰開始，轉(zhuǎn)變后的峰位開始錯(cuò)位并逐漸消失，說明轉(zhuǎn)變前后結(jié)構(gòu)發(fā)生了明顯變化，結(jié)合快照?qǐng)D考察，亞表Co團(tuán)簇由FCC結(jié)構(gòu)轉(zhuǎn)變?yōu)榫哂忻黠@五重對(duì)稱軸的二十面體結(jié)構(gòu).對(duì)于次表Co團(tuán)簇，其轉(zhuǎn)變前后的對(duì)分布函數(shù)分別與亞表Co團(tuán)簇轉(zhuǎn)變前后的對(duì)分布函數(shù)有著相同的趨勢(shì)，他們的峰位和峰強(qiáng)基本一致，這說明次表Co團(tuán)簇轉(zhuǎn)變前后也分別是FCC結(jié)構(gòu)和二十面體結(jié)構(gòu)，對(duì)應(yīng)的結(jié)構(gòu)快照?qǐng)D驗(yàn)證了這一點(diǎn).總的來看，亞表Co團(tuán)簇和次表Co團(tuán)簇，盡管都發(fā)生了從FCC到二十面體的結(jié)構(gòu)轉(zhuǎn)變，但其結(jié)構(gòu)中Co原子所處層位基本沒變化，說明這種結(jié)構(gòu)轉(zhuǎn)變是一種擴(kuò)散度較小的相變.另外，由于二十面體較立方八面體結(jié)構(gòu)更穩(wěn)定，所以團(tuán)簇在從立方八面體向二十面體結(jié)構(gòu)轉(zhuǎn)變的過程中勢(shì)能會(huì)降低，從而說明了圖2中一級(jí)相變點(diǎn)勢(shì)能向下跳躍的原因.

圖4 亞表Co團(tuán)簇和次表Co團(tuán)簇在第一相變點(diǎn)前后的對(duì)分布函數(shù)Fig.4 Pair correlation function of subsurface cluster and 3rd-shell cluster at the temperature before and after the first phase transition point

鍵對(duì)分析方法已被廣泛用于監(jiān)測(cè)液體、玻璃態(tài)及非晶結(jié)構(gòu)，這種技術(shù)主要采用原子與近鄰原子的關(guān)系(HA指數(shù))來考察團(tuán)簇升溫過程中整體的結(jié)構(gòu)變化，圖5、6分別給出了亞表Co團(tuán)簇和次表Co團(tuán)簇的HA指數(shù)隨溫度變化的關(guān)系曲線，其中1421鍵表征面心立方結(jié)構(gòu)，1422鍵表征體心立方結(jié)構(gòu)，1551鍵表征二十面體結(jié)構(gòu)，1541健表征缺陷及無序結(jié)構(gòu).圖5中在300 K時(shí)1551鍵值為0，而1421鍵值為0.73.在370 K左右，亞表Co團(tuán)簇的HA指數(shù)有明顯的變化，其中1421鍵值變?yōu)?.29，團(tuán)簇已不再具有FCC結(jié)構(gòu)，相應(yīng)1551鍵值由0變?yōu)?.03，此時(shí)團(tuán)簇已由FCC結(jié)構(gòu)轉(zhuǎn)變?yōu)槎骟w結(jié)構(gòu)，這和圖4結(jié)論一致.從300 K到熔點(diǎn)934 K之后，團(tuán)簇的1551、1541、1421及1422鍵值相差幅度由大變小，最后基本接近，這是由于從低溫到高溫團(tuán)簇由固體變?yōu)橐后w熔化所造成的.圖6中亞表Co團(tuán)簇HA指數(shù)在升溫區(qū)間內(nèi)也存在明顯的兩次跳躍，但其變化區(qū)間在400-500 K及1000-1100 K之間.在一級(jí)相變點(diǎn)，亞表Co團(tuán)簇在轉(zhuǎn)變前1421鍵值為0.67，轉(zhuǎn)變后1421明顯降低，此時(shí)1551鍵值為0.03團(tuán)簇此時(shí)已轉(zhuǎn)變?yōu)槎骟w.在熔點(diǎn)后，團(tuán)簇各HA指數(shù)大小接近，此時(shí)團(tuán)簇已熔化.

圖5 亞表Co團(tuán)簇HA指數(shù)隨溫度的變化Fig.5 HA index as a function of temperature for the subsurface cluster

圖6 次表Co團(tuán)簇HA指數(shù)隨溫度的變化Fig.6 HA index as a function of temperature for the 3rd-shell cluster

為了研究團(tuán)簇分層摻雜產(chǎn)生的結(jié)構(gòu)轉(zhuǎn)變點(diǎn)及熔點(diǎn)的差異，對(duì)Co分布于團(tuán)簇中不同位置點(diǎn)的能量進(jìn)行了分析.主要的技術(shù)手段是用一個(gè)Co原子替換純Cu二十面體團(tuán)簇中每一層不同位置(棱邊、頂點(diǎn)、(111)面)的Cu原子，然后將替換后的團(tuán)簇在0K下弛豫1ns得到團(tuán)簇每個(gè)原子的平均能量(見圖7).從圖中可以發(fā)現(xiàn)，Co原子分布于亞表層(111)面時(shí)能量最低，并且其處于表層時(shí)的原子能量遠(yuǎn)高于其它各層.Co原子位于亞表層時(shí)，在棱邊的能量最高，其次是頂點(diǎn)，最低的是(111)面；Co原子位于次表層時(shí)，Co在頂點(diǎn)和棱邊的能量相同且略低于(111)面.Co原子這種不同分布，會(huì)造成三個(gè)結(jié)果，一是Co原子分布在Cu團(tuán)簇的亞表層(111)面的能量低于Co分布在次表層(111)面，這會(huì)導(dǎo)致要使亞表Co團(tuán)簇熔化需要更高的溫度，這與圖3所得結(jié)果取得一致.二是隨著溫度的升高，亞表及次表Co團(tuán)簇中的Co原子位置在各層頂點(diǎn)、棱邊和(111)面變化會(huì)導(dǎo)致團(tuán)簇能量不斷變化，這與圖2曲線取得一致.三是隨著溫度的升高，次亞表層Co團(tuán)簇中的Co原子易于偏析到團(tuán)簇的亞表層，而亞表層Co團(tuán)簇中的Co原子位置易于從頂點(diǎn)和棱邊向(111)面變化.

圖7 團(tuán)簇在0 K弛豫1 ns后原子的平均能量及此表Co團(tuán)簇升溫過程中的原子能量和Cu原子與團(tuán)簇質(zhì)心的間距Fig.7 The average atom energy of clusters at 0 K and relaxation after 1 ns，the atomic energy in the process of Co clusters on the second floor heating and the distance from mass center of Cu atom to cluster

為了探究團(tuán)簇熔點(diǎn)差異的原因，圖7附圖還給出了Co原子與團(tuán)簇質(zhì)心的平均距離隨溫度變化函數(shù).由于原子偏析會(huì)使團(tuán)簇局部晶格錯(cuò)配，且處在亞表層的Co原子發(fā)生偏析時(shí)其位置很難精確定位，于是附圖中只給出了次表Co團(tuán)簇中Co原子溫度-離質(zhì)心平均距離關(guān)系，其中Co原子與質(zhì)心的平均距離越大說明偏析過程越明顯.可以看出，團(tuán)簇在300-400 K之間有一個(gè)向下的跳躍，但這種轉(zhuǎn)變量較小，說明此時(shí)Co原子位置有一個(gè)較小的變化，這和圖4結(jié)論“Co原子所處層位基本沒變化，說明這種結(jié)構(gòu)轉(zhuǎn)變是一種擴(kuò)散度較小的相變”取得一致.在400-900 K之間團(tuán)簇中Co的原子的相對(duì)位置基本沒變，團(tuán)簇的形貌基本保持二十面體結(jié)構(gòu).在900-1000 K之間，隨著溫度升高，Co原子離質(zhì)心的距離有一個(gè)突然的激增，表明大量Co原子開始從次亞表層向亞表層偏析，且在934 K團(tuán)簇由二十體結(jié)構(gòu)轉(zhuǎn)變?yōu)闊o序結(jié)構(gòu)，表面團(tuán)簇已經(jīng)熔化，這說明Co原子的偏析伴隨于團(tuán)簇熔化過程中，是造成團(tuán)簇熔點(diǎn)差異的主要誘因.

4 結(jié) 論

本文采用分子動(dòng)力學(xué)結(jié)合嵌入原子方法，對(duì)比研究了升溫過程中亞表及次表Co團(tuán)簇的結(jié)構(gòu)和性質(zhì).研究結(jié)果表明：分層摻雜會(huì)使純Cu團(tuán)簇、次表Co團(tuán)簇及亞表Co團(tuán)簇的結(jié)構(gòu)轉(zhuǎn)變點(diǎn)和熔點(diǎn)有差別，說明分層摻雜可以對(duì)團(tuán)簇的結(jié)構(gòu)轉(zhuǎn)變點(diǎn)和熔點(diǎn)進(jìn)行誘導(dǎo)，這為制備可控熔點(diǎn)的雙金屬團(tuán)簇提供一條有效途徑；分層摻雜的Cu-Co團(tuán)簇第一相變是一種擴(kuò)散度較小的由立方八面體轉(zhuǎn)變?yōu)槎骟w結(jié)構(gòu)的相變；由于分布在團(tuán)簇各層不同頂點(diǎn)、棱邊、(111)面上的Co原子能量的不同，會(huì)使Co原子易于向低能態(tài)的團(tuán)簇亞表層(111)面偏析，從而誘導(dǎo)團(tuán)簇結(jié)構(gòu)紊亂，造成其熔點(diǎn)差異.

[1] Wang G H.Clusterphysics[M].Shanghai: Shanghai Scientific & Technical Publisher, 2003(in Chinese)[王廣厚. 團(tuán)簇物理學(xué)[M]. 上海： 上?？茖W(xué)技術(shù)出版社，2003]

[2] Ferrando R，Jellinek J，Johnston R L. Nanoalloys： from theory to applications of alloy clusters and nanoparticles[J].Chem.Rev., 2008，108(3)： 845.

[3] Mariscal M M，Dassie S A，Leiva E P M. Collision as a way of forming bimetallic nanoclusters of various structures and chemical compositions[J].J.Chem.Phys., 2005，123(3)： 184505.

[4] Kart H H，Yildirim H，Ozdemir K S，etal.Physical properties of Cu nanoparticles：A molecular dynamics study[J].Mater.Chem.Phys., 2014，147(1-2)： 204.

[5] Lewis L J，Jensen P，Barrat J L. Melting，freezing，and coalescence of gold nanoclusters[J].Phys.Rev. B，1997，56(4)： 2248.

[6] Baletto F，F(xiàn)errando R. Structural properties of nanoclusters： Energetic，thermodynamic，and kinetic effects[J].Rev.Mod.Phys., 2005，77(1)： 371.

[7] Schmidt M，Kusche R，Issendorff B V，etal. Irregular variations in the melting point of size-selected atomic clusters[J].Nature，1998，393(6682)： 238.

[8] Yang J Y，Hu W Y，Wu Y R，etal. Substrate dependence of growth configurations for Co-Cu bimetallic clusters[J].Cryst.GrowthDes., 2012，12(6)： 2978.

[9] Chushak Y G，Bartell L S. Freezing of Ni-Al bimetallic nanoclusters in computer simulations[J].J.Phys.Chem. B，2003，107(16)： 3747.

[10] Song P X，Wen D S. Molecular dynamics simulation of a Core-Shell structured metallic nanoparticle[J].J.Phys.Chem. C，2010，114(19)： 8688.

[11] Ferrer D，Torres-Castro A，Gao X，etal. Three-layer core/shell structure in Au-Pd bimetallic nanoparticles[J].Nano.Lett., 2007，7(6)： 1701.

[12] Dorfbauer F，Schrefl T，Kirschner M，etal. Nanostructure calculation of CoAg core-shell clusters[J].J.Appl.Phys., 2006，99(8)： 08G706.

[13] Wang Q，Li G J，Li D G，etal. Evolution of three-shell onion-like and core-shell structures in(AgCo)201 bimetallic clusters[J].Chin.Phys. B，2009，18(5)： 1843.

[14] Chen F Y，Johnston R L.Martensitic transformation in Ag-Au bimetallic core-shell nanoalloys[J].Appl.Phys.Lett., 2008，92(2)： 023112.

[15] Zhang Y J，Li Y Q，Xiao X Y，etal. Effect of Cu atomic segregation on the frozen structures of Co-Cu bimetallic clusters[J].NANO：BriefReportsandReviews，2012，7(6)： 1250047.

[16] Li G J，Wang K，Wang Q，etal. Formation of icosahedral and hcp structures in bimetallic Co-Cu clusters during the freezing processes[J].Mater.Lett., 2012，88： 126.

[17] Fromen M C, Morillo J，Casanove M J，etal. Structure and chemical order in Co-Rh nanoparticles[J].Epl-Europhys.Lett., 2006，73(1)： 885.

[18] Nanda K K，Sahu S N，Behera S N. Liquid-drop model for the size-dependent melting of low-dimensional systems[J].Phys.Rev. A，2002，66(1)： 013208.

[19] Honeycutt J D，Andersen H C. Molecular dynamics study of melting and freezing of small Lennard-Jones clusters[J].J.Phys.Chem., 1987，91(19)： 4950.

[20] Zhou X W，Wadley H N G，Johnson R A，etal. Atomic scale structure of sputtered metal multilayers[J].Acta.Mater., 2001，49(19)： 4005.

[21] Wadley H N G，Zhou X，Johnson R A，etal. Mechanisms，models and methods of vapor deposition[J].Prog.Mat.Sci., 2001，46(5)： 329.

Molecular dynamics simulation on the structure evolution of bimetallic Cu-Co clusters and their properties

SUN Ling-Tao，SHI Dong-Ping

(Research Institute for New Materials Technology，Chongqing University of Arts and Sciences，Chongqing 402160，China)

A comparison study on the structures and properties of the Co atom distribution in different layers of Cu-Co clusters was performed by the molecular dynamics combining with atom-embedded method.Research results showed that the layer-doping of Co atom can induce and control the structural transformation point and melting point of clusters.The first phase transition of layer-doped Cu-Co clusters is a smaller diffusion degree of phase transformation from cubic octahedron into icosahedron.Additionally，Co atom is easy to have the surface segregation tendency to the subsurface(111)of the low-energy state of clusters，further resulting in the structural disorder of the cluster and the difference of the melting point.

Clustes; Structure; Segregation; Molecular dynamics

2014-12-09

重慶文理學(xué)院校級(jí)科研項(xiàng)目(Y2013CJ26); 重慶市高校微納米材料工程與技術(shù)重點(diǎn)實(shí)驗(yàn)室度開放課題(KFJJ1404)

孫凌濤(1987—), 男, 寧夏中衛(wèi)人, 碩士, 助教, 主要研究方向?yàn)椴牧衔锢砼c化學(xué).E-mail: cquptslt@163.com

103969/j.issn.1000-0364.2015.08.010

O561

A

1000-0364(2015)08-0586-05