余藝, 程慶軍, 郭詩(shī)奇, 楊劍瑜

?

合金納米粒子生長(zhǎng)形態(tài)的分子動(dòng)力學(xué)模擬

余藝, 程慶軍, 郭詩(shī)奇, 楊劍瑜

(湖南工程學(xué)院理學(xué)院, 湖南湘潭, 411104)

選用Cu–Ag納米粒子作為研究對(duì)象, 采用分子動(dòng)力學(xué)方法和嵌入原子模型模擬了將Cu原子注入到由Ag原子構(gòu)成的納米粒子的過(guò)程, 研究了注入能量(10, 30, 50 eV)與納米粒子結(jié)構(gòu)的關(guān)系。采用共近鄰分析技術(shù)對(duì)合金納米的結(jié)構(gòu)演變進(jìn)行了研究。模擬結(jié)果表明: 當(dāng)入射原子動(dòng)能為50 eV時(shí), 合金納米粒子從截角八面體轉(zhuǎn)變?yōu)槎骟w結(jié)構(gòu)。這種轉(zhuǎn)變是因?yàn)槿肷湓优c基底原子碰撞, 促使系統(tǒng)溫度升高。

Cu–Ag合金納米粒子; 分子動(dòng)力學(xué); 近鄰分析技術(shù); 截角八面體; 二十面體

納米合金是由2種以上納米粒子構(gòu)成的金屬[1–2], 因其獨(dú)特的物理、化學(xué)性質(zhì), 在材料科學(xué)中有著重要的地位及廣泛的應(yīng)用。納米合金的性質(zhì)與其結(jié)構(gòu)(包括構(gòu)型、尺寸、金屬組元及分布、成分比例等)有密切關(guān)系。根據(jù)應(yīng)用需要調(diào)控其結(jié)構(gòu)是納米合金制備技術(shù)發(fā)展的方向。實(shí)驗(yàn)制備納米合金, 一般是讓其在一定基底上生長(zhǎng)成形, 其結(jié)構(gòu)隨生長(zhǎng)條件(溫度、壓強(qiáng)、組元成分比例、基底構(gòu)型與尺寸等)的變化而演變[3]。如何實(shí)現(xiàn)生長(zhǎng)條件可控、協(xié)調(diào)變化, 進(jìn)而對(duì)納米合金結(jié)構(gòu)進(jìn)行調(diào)控, 是有重要理論價(jià)值與實(shí)踐意義的課題。合金納米粒子的結(jié)構(gòu)源于形成合金納米粒子的組元表面能、原子半徑的相對(duì)大小、合金形成能力的強(qiáng)弱、環(huán)境溫度、粒子尺寸與合金成分比例等多種因素。目前國(guó)內(nèi)外的研究幾乎都是從靜態(tài)角度, 利用能量?jī)?yōu)化的方法研究合金納米粒子的結(jié)構(gòu)。Ferrando課題組應(yīng)用這一方法研究了同為面心立方結(jié)構(gòu)的2種金屬組成的合金納米粒子結(jié)構(gòu)[4–6], 發(fā)現(xiàn)由于表面聚集, 這些合金納米粒子都有形成核殼結(jié)構(gòu)的趨勢(shì), 但其具體構(gòu)型與粒子尺寸、成分比例有密切關(guān)系; Aguado課題組研究了2種體心立方結(jié)構(gòu)的金屬元素組成的合金納米粒子的結(jié)構(gòu), 發(fā)現(xiàn)尺寸較小時(shí), 合金納米粒子不是核殼結(jié)構(gòu), 而是與純?cè)丶{米粒子相似的多面體結(jié)構(gòu)[7]。如上所述, 實(shí)驗(yàn)與理論都表明, 合金納米粒子制備過(guò)程中, 其結(jié)構(gòu)與吸附原子在納米粒子表面的擴(kuò)散與生長(zhǎng)行為有密切聯(lián)系。在原子尺度上對(duì)其結(jié)構(gòu)演變、生長(zhǎng)過(guò)程進(jìn)行研究有助于理解其性質(zhì), 并達(dá)到人為控制其生長(zhǎng)形態(tài)的目的。這些研究結(jié)果說(shuō)明, 合金納米粒子結(jié)構(gòu)與合金成分比例、組元與溫度等有密切關(guān)系。雖然能量?jī)?yōu)化的方法被廣泛應(yīng)用, 但這種方法有2個(gè)缺陷: 第一, 能量局部最小的數(shù)量隨納米粒子包含原子數(shù)的增加呈指數(shù)遞增, 因此要找出全局最小需要極大的計(jì)算量, 這就使得此方法只適用于小尺度系統(tǒng)[6]; 第二, 合金納米粒子制備是一個(gè)動(dòng)力學(xué)過(guò)程, 制備條件對(duì)納米粒子結(jié)構(gòu)影響非常大, 能量最小對(duì)應(yīng)的合金納米粒子結(jié)構(gòu)僅是理想組態(tài), 因此用能量最優(yōu)化方法得到的合金納米粒子結(jié)構(gòu)與實(shí)驗(yàn)結(jié)果有較大差別。由于實(shí)驗(yàn)只是調(diào)節(jié)某一個(gè)生長(zhǎng)條件, 也沒(méi)有考慮到制備過(guò)程中生長(zhǎng)條件是動(dòng)態(tài)變化的, 這些實(shí)驗(yàn)也就不能實(shí)現(xiàn)生長(zhǎng)條件的可控、協(xié)調(diào)變化, 而且生長(zhǎng)條件影響納米合金結(jié)構(gòu)的物理機(jī)制也還有待深入研究。因此, 已有實(shí)驗(yàn)還只是一些零星的、非系統(tǒng)性的研究結(jié)果[3, 8–9]。理論研究方面, 分子動(dòng)力學(xué)能跟蹤原子運(yùn)動(dòng)軌跡[10–12], 已經(jīng)被廣泛應(yīng)用于研究合金納米粒子的生長(zhǎng)形態(tài)。為研究Cu–Ag合金納米粒子生長(zhǎng)形態(tài)的調(diào)控和合金納米粒子結(jié)構(gòu)隨溫度、尺寸、合金成份比例與組元結(jié)構(gòu)等生長(zhǎng)條件的變化規(guī)律, 將Cu沉積原子注入到Ag原子納米尺度的基底中, 基底的動(dòng)能設(shè)置為0, 用高能粒子(Cu原子)的碰撞使基底轉(zhuǎn)變生長(zhǎng)方式。為實(shí)現(xiàn)這一目標(biāo), 本文從理論上闡明影響合金納米粒子結(jié)構(gòu)的因素及演變規(guī)律,著重探討注入能量這一因素對(duì)截角八面體結(jié)構(gòu)的201個(gè)Ag原子生長(zhǎng)形態(tài)的影響。

1 模型與方法

利用分子動(dòng)力學(xué)進(jìn)行計(jì)算機(jī)模擬時(shí), 準(zhǔn)確描述原子間作用力勢(shì)函數(shù)是模擬成功與否的關(guān)鍵[12–14]。嵌入原子方法(EAM)能夠很好地描述金屬原子間多體相互作用。表面能和晶格常數(shù)2個(gè)物理量是納米粒子生長(zhǎng)模擬的重要因素, 其中晶格常數(shù)的實(shí)驗(yàn)值被作為輸入?yún)?shù)來(lái)確定EAM模型參數(shù)。因此, 應(yīng)用EAM模型計(jì)算了2個(gè)元素的表面能()可以驗(yàn)證模型的正確性。在一般情況下, 表面能隨著表面原子的配位數(shù)減少而相應(yīng)減少, 因此對(duì)于面心立方結(jié)構(gòu)的金屬, 其表面能的順序應(yīng)該是110>100>111(表1)。EAM模型準(zhǔn)確地預(yù)測(cè)了表面能的這一趨勢(shì), 計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果和第一原理結(jié)果也比較接近。因此, 本模型能準(zhǔn)確描述Cu–Ag納米粒子。

表1 利用第一原理(FP)與嵌入原子方法(EAM)計(jì)算的Cu和Ag表面能及其實(shí)驗(yàn)值(Exp) /( mJ·m-2)

研究發(fā)現(xiàn)[1–2], 對(duì)于面心立方結(jié)構(gòu)(FCC)金屬, 截角八面體(TO)與二十面體(ICO)結(jié)構(gòu)的納米粒子是最穩(wěn)定的結(jié)構(gòu)。因此, 選擇201個(gè)原子的截角八面體(TO201)和147個(gè)原子的二十面體(ICO147)作為研究對(duì)象, 這2種納米粒子的構(gòu)型如圖1所示。截角八面體由8個(gè){111}晶面和6個(gè){100}晶面構(gòu)成, 需要指出的是{111}與{100}面是面心立方結(jié)構(gòu)的密排面, 也是最穩(wěn)定的表面, 截角八面體構(gòu)型與球形比較接近, 這也使得它的結(jié)構(gòu)更為穩(wěn)定。二十面體則是以2個(gè)五角形金字塔為基礎(chǔ), 它的表面只有排列最為緊密的{111}晶面, 從這點(diǎn)來(lái)說(shuō), 它也是一種穩(wěn)定結(jié)構(gòu)。另外, 為便于說(shuō)明這2種納米粒子的結(jié)構(gòu)特點(diǎn), 分別在TO201與ICO147納米粒子表面上描繪出一個(gè)正六邊形與正五邊形, 如圖1中的白色多邊形所示。

圖1 201個(gè)原子的截角八面體和147個(gè)原子的二十面體結(jié)構(gòu)

在生長(zhǎng)模擬中, 利用分子動(dòng)力學(xué)方法模擬將具有一定動(dòng)能的Cu原子沉積到截角八面體的Ag納米粒子上, 沉積原子的能量范圍是10～50 eV。本文以注入能量為10, 30, 50 eV為例來(lái)說(shuō)明模擬結(jié)構(gòu)。生長(zhǎng)模擬在等容等能(NVE)與等容等溫(NVT)的正則系統(tǒng)中進(jìn)行。首先在NVE系統(tǒng)中, 模擬具有一定動(dòng)能的Cu原子的注入沉積過(guò)程, 經(jīng)過(guò)0.3 ns的模擬, 系統(tǒng)基本達(dá)到平衡狀態(tài)。然后經(jīng)過(guò)2.7 ns的NVT系統(tǒng)模擬, 將系統(tǒng)平衡在溫度為100 K的狀態(tài)。最后記錄所有原子的位置和速度, 以便于做結(jié)構(gòu)和動(dòng)力學(xué)特性分析。在分子動(dòng)力學(xué)模擬中, 利用預(yù)測(cè)–校正算法求解原子的牛頓運(yùn)動(dòng)方程, 時(shí)間步長(zhǎng)取2.0 fs。為保證系統(tǒng)的動(dòng)量為0, 3個(gè)Cu原子被均等地沉積到Ag基底上, 其中每個(gè)Cu原子的運(yùn)動(dòng)方向都指向Ag納米粒子中心。

如圖2所示, 當(dāng)入射的Cu原子動(dòng)能為50 eV時(shí), 在開(kāi)始階段(NVE系統(tǒng)), 系統(tǒng)溫度(約800 K)遠(yuǎn)高于設(shè)定的溫度(100 K), 這是因?yàn)锳g基底原子獲得了入射Cu原子的動(dòng)能。然后在NVT系統(tǒng)下, Cu–Ag系統(tǒng)溫度逐漸降低, 最后平衡在設(shè)定的100 K左右。所有物理量都在100 K下統(tǒng)計(jì)獲得。

圖2 Cu原子入射能量為50 eV時(shí)基底納米粒子溫度隨時(shí)間變化

2 結(jié)果與討論

當(dāng)把Cu原子注入到Ag納米粒子TO201, 納米合金首先在NVE中模擬沉積過(guò)程(見(jiàn)圖2), 入射能量等于50 eV時(shí), Cu原子的入射能量傳遞到基底, 基底原子之間劇烈碰撞, 使系統(tǒng)溫度變?yōu)榧s850 K。NVE中模擬的850 K高溫會(huì)導(dǎo)致Ag納米粒子基底開(kāi)始快速扭曲變形, 然后經(jīng)過(guò)2.7 ns的NVT模擬, 系統(tǒng)平衡溫度冷卻至100 K, 納米粒子由高溫(850 K)時(shí)的液態(tài)凝固成一定結(jié)構(gòu)的固態(tài)。通常的分析技術(shù), 如徑向分布函數(shù)和靜態(tài)結(jié)構(gòu)因子, 只能從整體上對(duì)結(jié)構(gòu)作有序與無(wú)序的判斷。為更細(xì)致地研究納米合金的結(jié)構(gòu)變化過(guò)程, 采用共近鄰分析技術(shù)(CNA)確定凝固的Cu–Ag合金納米粒子結(jié)構(gòu)。共近鄰分析技術(shù)通過(guò)3位數(shù)(,,)描述原子所屬的狀態(tài),表示成鍵的2原子的共有最近鄰原子數(shù),表示共有最近鄰原子之間的成鍵數(shù),表示共有最近鄰原子所成鍵鏈中最長(zhǎng)鏈上的鍵數(shù)。當(dāng)將Cu原子入射到Ag基底中, 入射能量為10, 30, 50 eV時(shí), 共近鄰分析的結(jié)果列于表2。

表2 Cu注入到Ag的TO201中(3, 1, 1), (4, 2, 1),(4, 2, 2)和(5, 5, 5)等鍵對(duì)占總鍵對(duì)的百分比 %

根據(jù)共近鄰分析方法, 每一種鍵對(duì)都對(duì)應(yīng)著特定的原子空間幾何構(gòu)型。其中FCC型晶體中以(4, 2, 1)鍵為特征鍵對(duì); HCP型晶體則同時(shí)以(4, 2, 1)和(4, 2, 2)為特征鍵對(duì)。(3, 1, 1)是納米粒子的表面原子, (5, 5, 5)是二十面體5重對(duì)稱性的晶序特性。因此僅需要(3, 1, 1), (4, 2, 1), (4, 2, 2), (5, 5, 5)這幾種晶序就能區(qū)分TO201, ICO147結(jié)構(gòu)。

如表2所示, 當(dāng)Cu原子以10 eV和30 eV的能量入射TO201時(shí), 新的Cu–Ag合金納米粒子中沒(méi)有(4, 2, 2)和(5, 5, 5)構(gòu)型, 僅有(3, 1, 1)與(4, 2, 1)鍵, 且這2種鍵型所占百分比與Ag的TO201類似, TO201的Ag基底具有{111}和{100}面, 經(jīng)CNA分析, (4, 2, 2)和(5, 5, 5)鍵對(duì)數(shù)為0。這表明當(dāng)Cu原子以較低能量入射時(shí), Ag基底仍保持著原來(lái)的結(jié)構(gòu)而沒(méi)有發(fā)生變化。當(dāng)入射的Cu原子動(dòng)能增加大到50 eV時(shí), (4, 2, 2)鍵所占百分比提高到了18%, (4, 2, 1)鍵的比例減少到41%。需要指出的是, 在此情形下, 有少量的(5, 5, 5)鍵型。這種鍵型特征與ICO147類似。共近鄰分析方法說(shuō)明, 由于高能量的Cu原子注入, Ag納米粒子由TO轉(zhuǎn)變?yōu)镮CO結(jié)構(gòu)。用Cu原子注入后, 合金納米粒子結(jié)構(gòu)演變過(guò)程圖更直接地說(shuō)明這種轉(zhuǎn)變。

如圖3所示, 動(dòng)能為10 eV的Cu原子入射到Ag的TO基底上, Ag基底幾乎保持了原來(lái)的TO結(jié)構(gòu), 入射的Cu只是嵌入到Ag基底中。最終平衡后, 合金納米粒子保持了原來(lái)的六邊對(duì)稱性。

當(dāng)入射的Cu原子動(dòng)能為50 eV時(shí)(圖4), 與10 eV相比, 50 eV的Cu原子與基底的Ag原子強(qiáng)烈作用, 粒子的運(yùn)動(dòng)過(guò)程要更加復(fù)雜, 溫度急劇升高促使結(jié)構(gòu)發(fā)生變化, 由6重對(duì)稱演變成5重對(duì)稱。Ag納米粒子出現(xiàn)了5重對(duì)稱結(jié)構(gòu), 并出現(xiàn)了幾個(gè){100}面, 這樣的結(jié)構(gòu)更類似于ICO構(gòu)型。另外, 因Ag原子表面能較低、原子半徑較大, 最終導(dǎo)致50 eV入射的Cu沉積原子在進(jìn)入TO后偏析于Ag納米粒子表面。

如圖5所示, 對(duì)尺寸較小的Cu納米粒子, ICO比TO穩(wěn)定, 因此TO比ICO更容易演變成其它結(jié)構(gòu)。對(duì)所有尺寸的Ag納米粒子, TO比ICO更穩(wěn)定, 這一結(jié)果與其它理論的計(jì)算結(jié)果一致。當(dāng)Cu注入到Ag的TO201基底中, 因TO很穩(wěn)定, 其很難演變成ICO。因此, 只有當(dāng)以較大動(dòng)能的Cu原子注入到Ag基底中, 使基底原子發(fā)生較大變化, 才能使Ag基底從TO轉(zhuǎn)化為ICO。

圖3 注入能量為10 eV的Cu原子沉積到結(jié)構(gòu)為T(mén)O的Ag基底中納米粒子構(gòu)型演變

圖4 注入能量為50 eV的Cu原子沉積到結(jié)構(gòu)為T(mén)O的Ag基底中納米粒子構(gòu)型演變

圖5 (Cu, Ag)TO和ICO結(jié)構(gòu)結(jié)合能隨原子尺寸的變化

為了解釋這一模擬結(jié)果, 進(jìn)行了一系列不同尺寸的TO和ICO納米粒子的MD模擬實(shí)驗(yàn), 證明每個(gè)原子結(jié)合能的大小與原子尺寸相關(guān)。

3 結(jié)論

用分子動(dòng)力學(xué)和原子嵌入模型研究了能量在10～50 eV的Cu原子沉積到TO201構(gòu)型的Ag納米粒子基底的過(guò)程, 分析了生長(zhǎng)出的Cu–Ag納米合金結(jié)構(gòu)的早期演變過(guò)程。應(yīng)用共近鄰分析技術(shù), 發(fā)現(xiàn)當(dāng)入射Cu原子的能量為10 eV和30 eV時(shí), 生長(zhǎng)出的納米合金粒子基本保持了原有TO形態(tài); 在50 eV時(shí), 納米粒子由六重對(duì)稱性變?yōu)槲逯貙?duì)稱性。因沉積原子的能量被均勻地分配到基底每個(gè)原子, 基底原子不斷擠壓碰撞導(dǎo)致了系統(tǒng)溫度的升高, 溫度的升高進(jìn)而導(dǎo)致了納米粒子結(jié)構(gòu)的轉(zhuǎn)變。

參考文獻(xiàn)：

[1] Baletto F, Ferrando R. Structural properties of nanoclusters: Energetic, thermodynamic, and kinetic effects [J]. Rev Mod Phys, 2005, 77(1): 371–423.

[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–910.

[3] Chaudhuri R G, Paria S. Core/shell nanoparticles: classes, properties, synthesis mechanisms, characterization, and applications [J]. Chem Rev, 2012, 112(4): 2 373–2 433.

[4] Rapallo A, Rossi G, Ferrando R, et al. Global optimization of bimetallic cluster structures. I. Size-mismatched Ag-Cu, Ag-Ni, and Au-Cu systems [J]. J Chem Phys, 2005, 122(19): 194308(1–13).

[5] Rapallo A, Rossi G, Ferrando R, et al. Global optimization of bimetallic cluster structures. II. Size-matched Ag-Pd, Ag-Au, and Pd-Pt systems [J]. J Chem Phys, 2005, 122(19): 194309(1–9).

[6] Negreiros F R, Barcaro G, Kuntova Z. et al. Structures of gas-phase Ag-Pd nanoclusters: A computational study [J]. J Chem Phys, 2010, 132(23): 234703(1–8).

[7] Aguado A, Lopez J M. Structure determination in 55-atom Li-Na and Na-K nanoalloys [J]. J Chem Phys, 2010, 133(9): 094302(1–10).

[8] Serpell C J, Cookson J, Ozkaya D, et al. Core@shell bimetal-lic nanoparticle synthesis via anion coordination [J]. Nat Chem, 2011, 3(6): 478–483.

[9] Radilovic V, Ophus C, Marquis E A, et al. Highly monodisperse core-shell particles created by solid-state reactions [J]. Nat Mater, 2011, 10(9): 710–715.

[10] 文玉華, 朱如曾, 周富信. 分子動(dòng)力學(xué)模擬的主要技術(shù)[J]. 力學(xué)進(jìn)展, 2003, 33(1): 67–71.

[11] 蔡錫年. 分子動(dòng)力學(xué)和物理力學(xué)[M]. 北京: 科學(xué)出版社, 1986: 100–110.

[12] 周健, 陸小華, 王延儒, 等. 分子動(dòng)力學(xué)模擬[J].南京化工大學(xué)學(xué)報(bào), 1997, 19(2): 20–25.

[13] Vitos L, Ruban A V, Skriver H L. et al. The surface energy of metals [J]. Surf Sci, 1998, 411(1–2): 186–202.

[14] Tyson W R, Miller W A. A Comprehensive Search for Stable Pt-Pd Nanoalloy Configurations and Their Use as Tunable Catalysts [J]. Surf Sci, 1977, 62(1): 267–276.

(責(zé)任編校:劉剛毅)

The growth configuration of nanoalloy: a molecular dynamics simulation

Yu Yi, Cheng Qingjun, Guo Shiqi, Yang Jianyu

(Science College, Hunan Institute of Engineering, Xiangtan 411104, China)

The impact depositions of Cu atoms onto the Ag truncated octahedron are simulated by molecular dynamics along with the embedded atom method. The dependence of incident energies (10, 30 and 50 eV) and the nanoalloy configuration is found. The common neighbor analysis is used to identify the structural evolutions. A structural transformation from truncated octahedral to icosahedral arrangement is observed at incident energy of 50 eV. The structural transformation can be induced by the increase of system temperature, because the incident energy of the deposited atom is delivered to the substrate atoms as a result of collisions.

Cu-Ag nanoalloy; molecular dynamics; common neighbor analysis; truncated octahedron; icosahedron

10.3969/j.issn.1672–6146.2015.01.005

O 781

1672–6146(2015)01–0017–04

楊劍瑜, hnieyjy@aliyun.com; 余藝, 1847997096@qq.com。

2014–11–15

國(guó)家自然科學(xué)基金(50971058), 2014年湖南省大學(xué)生科技創(chuàng)新項(xiàng)目。