余藝, 程慶軍, 郭詩奇, 楊劍瑜

?

合金納米粒子生長形態(tài)的分子動力學模擬

余藝, 程慶軍, 郭詩奇, 楊劍瑜

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

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

Cu–Ag合金納米粒子; 分子動力學; 近鄰分析技術; 截角八面體; 二十面體

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

1 模型與方法

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

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

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

圖1 201個原子的截角八面體和147個原子的二十面體結構

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

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

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

2 結果與討論

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

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

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

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

如圖3所示, 動能為10 eV的Cu原子入射到Ag的TO基底上, Ag基底幾乎保持了原來的TO結構, 入射的Cu只是嵌入到Ag基底中。最終平衡后, 合金納米粒子保持了原來的六邊對稱性。

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

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

圖3 注入能量為10 eV的Cu原子沉積到結構為TO的Ag基底中納米粒子構型演變

圖4 注入能量為50 eV的Cu原子沉積到結構為TO的Ag基底中納米粒子構型演變

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

為了解釋這一模擬結果, 進行了一系列不同尺寸的TO和ICO納米粒子的MD模擬實驗, 證明每個原子結合能的大小與原子尺寸相關。

3 結論

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

參考文獻：

[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] 文玉華, 朱如曾, 周富信. 分子動力學模擬的主要技術[J]. 力學進展, 2003, 33(1): 67–71.

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

[12] 周健, 陸小華, 王延儒, 等. 分子動力學模擬[J].南京化工大學學報, 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.

(責任編校:劉剛毅)

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

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