Cu薄膜三維生長(cháng)的Monte Carlo模擬
利用蒙特卡羅(Monte Carlo)方法模擬了Cu薄膜在四方基底上的三維生長(cháng)過(guò)程。模型中考慮了三個(gè)主要的原子熱運動(dòng)過(guò)程:原子沉積、原子擴散、原子脫附,各過(guò)程發(fā)生的概率是由各運動(dòng)的速率來(lái)決定的。討論了基底溫度、沉積速率及原子覆蓋度對Cu薄膜的表面形貌及表面粗糙度的影響。模擬結果表明:隨基底溫度升高或沉積速率下降,島的平均尺寸增大,數目減少,薄膜以層狀生長(cháng)方式生長(cháng);Cu薄膜表面粗糙度隨溫度的升高而減小,當基底溫度處于某一臨界溫度之內時(shí),表面粗糙度隨沉積速率的變化很大,但當基底溫度超過(guò)臨界溫度時(shí),表面粗糙度隨沉積速率的變化很小;薄膜的粗糙度與薄膜亞單層的形核密切相關(guān)。
關(guān)鍵詞：粗糙度;Monte Carlo模型;薄膜生長(cháng);形貌
隨著(zhù)微電子行業(yè)的高速發(fā)展, 作為主要互連線(xiàn)的Cu 薄膜應用越來(lái)越廣泛。而Cu 薄膜互連線(xiàn)的性能主要依賴(lài)于Cu 薄膜的微觀(guān)結構, 因此從原子尺度上研究Cu 薄膜生長(cháng)的微觀(guān)機制, 從而控制薄膜的性能, 提高薄膜制備質(zhì)量具有重要的研究意義。薄膜生長(cháng)過(guò)程是一個(gè)隨機動(dòng)力學(xué)過(guò)程, 因此,Monte Carlo 方法很自然應用于研究模擬這一過(guò)程。國內外研究者用Monte Carlo 方法從不同的方面研究了薄膜生長(cháng)的微觀(guān)過(guò)程。劉祖黎等[1] 用動(dòng)力學(xué)晶格Monte Carlo 方法模擬了Cu 薄膜的三維生長(cháng), 研究了濺射沉積條件下粒子的沉積角度、沉積速率以及入射能量對Cu 薄膜生長(cháng)的影響, 但沒(méi)有研究溫度的變化對薄膜生長(cháng)的影響。在他們的模型中通過(guò)控制沉積及擴散的時(shí)間來(lái)實(shí)現薄膜的三維生長(cháng)模擬。吳子若等[2] 利用動(dòng)力學(xué)晶格Monte Carlo 方法模擬了Cu 薄膜的三維生長(cháng)過(guò)程, 討論了基底溫度、沉積速率及原子覆蓋率對Cu 原子遷移、成核和表面島生長(cháng)等微觀(guān)生長(cháng)機制的影響, 獲得了Cu 薄膜的表面形貌圖并計算了表面粗糙度。在他們的模型中, 關(guān)于原子的擴散只考慮了單原子遷移、雙原子遷移及臺階原子的遷移, 更為重要的是他們采用的三種遷移勢壘都是常數, 也就意味著(zhù)在研究擴散的時(shí)候沒(méi)有考慮擴散原子周?chē)拥姆植记闆r。鄭小平等[3]建立了一個(gè)比較合理的三維模型, 并通過(guò)模擬成像和定量計算研究了薄膜生長(cháng)過(guò)程中的兩個(gè)重要問(wèn)題, 早期成核與表面粗糙度。他們運用勢函數計算出遷移速率和脫附速率, 再將每一個(gè)事件組成一個(gè)列表, 然后從表中隨機地抽取一個(gè)事件執行, 抽取的依據是某一事件的發(fā)生幾率與事件列表中所有事件幾率之和的比值最大。由于在計算擴散時(shí)考慮了擴散路徑的影響, 因此將會(huì )導致較大的計算量。
本文基于作者的Monte Carlo 模型模擬了在四方形基底上Cu 薄膜的三維生長(cháng)過(guò)程[4] , 研究了沉積速率、基底溫度及原子覆蓋度對Cu 薄膜表面形貌和表面粗糙度的影響。在計算過(guò)程中, 采用Born-Mayer 勢來(lái)計算Cu 原子的擴散速率, 在擴散過(guò)程中只考慮了原子的向下遷移, 并忽略了脫附過(guò)程對薄膜生長(cháng)的影響。
本文采用Monte Carlo 方法模擬Cu 薄膜的生長(cháng)過(guò)程并通過(guò)改變基底溫度、沉積速率及覆蓋度等來(lái)觀(guān)察各種因素對薄膜粗糙度的影響。得出如下結論:
(1) 在一定的沉積速率和覆蓋度的條件下, 薄膜的表面粗糙度會(huì )隨著(zhù)基底溫度的升高而減小, 薄膜的生長(cháng)方式由島狀生長(cháng)轉為層狀生長(cháng), 直到達到臨界溫度后, 薄膜的表面粗糙度隨著(zhù)基底溫度的升高而保持不變。
(2) 當覆蓋度給定, 基底溫度在一定范圍內時(shí),薄膜的表面粗糙度會(huì )隨著(zhù)沉積速率的減小而減小,薄膜的生長(cháng)方式向層狀生長(cháng)轉化, 但超過(guò)一定溫度后, 薄膜的表面粗糙度隨沉積速率的變化較小。
(3) 當基底溫度及沉積速率一定時(shí), 薄膜表面的原子覆蓋度越大粗糙度越大, 基底溫度較低時(shí), 覆蓋度的變化對粗糙度的影響較大, 而基底溫度較高時(shí),覆蓋度對粗糙度的影響變化不大; 薄膜表面的粗糙度與亞單層生長(cháng)時(shí)薄膜原子的形核有直接關(guān)系。
Abstract: The Cu film growth on surfaces of square-lattice substrate was modeled and simulated in 3-dimensional Monte Carlo method.The model considers three major thermodynamic processes:the deposition,diffusion and desorption of individual Cu atom,with the probabilities depending on the deposition rate,diffusion coefficient,and desorption rate,respectively.The impacts of the film growth conditions,such as the substrate temperature,deposition rate,and Cu coverage,on the microstructures of the film were calculated.The results show that the substrate temperature and deposition rate strongly affects the microstructures of the film.For instance,as the substrate temperature rises up,or the deposition rate decreases,the layer-by-layer growth mode dominates,with small islands coalescing into bigger ones,resulting in smoother surfaces.Below the critical temperature,an increase of the deposition rate roughens the surfaces;whereas above the critical temperature,the deposition rate weakly influences the surface roughness.Possible mechanisms were also tentatively discussed.
Keywords: Roughness,Monte Carlo model,Thin film growth,Morphologies
基金項目: 國家重點(diǎn)基礎研究發(fā)展計劃資助項目(2005CB321704);; 國家杰出青年科學(xué)基金項目(10925209)
參考文獻：
[1]劉祖黎,張雪峰,姚凱倫,等.濺射沉積Cu膜生長(cháng)的Monte Carlo模擬[J].真空科學(xué)與技術(shù)學(xué)報,2005,25(2):83-87
[2]吳子若,程鑫彬,王占山.動(dòng)力學(xué)晶格蒙特卡洛方法模擬Cu薄膜生長(cháng)[J].光子學(xué)報,2010,39(1):62-66
[3]鄭小平,張佩峰,范多旺,等.薄膜生長(cháng)的計算機模擬[J].材料研究學(xué)報,2005,19(2):170-178
[4]朱國.薄膜生長(cháng)初期的蒙特卡羅模型[J].大連理工大學(xué)學(xué)報,2009,49(6):781-785
[5]劉祖黎,魏合林,王漢文,等.薄膜生長(cháng)的隨機模型[J].物理學(xué)報,1999,48(7):1302-1308
[6]Bruschi P,Cagnoni P,Nannini A.Temperature DependentMonte Carlo Simulations of Thin Metal Film Growth and Per-colation[J].Phys Rev,1997,B55(12):7955-7963
[7]朱國.薄膜生長(cháng)的三維蒙特卡羅模型[J].材料研究學(xué)報,2009,23(6):640-645
[8]Wu Feng-Min,Lu Hang-Jun,Wu Zi-Qin.Simulation of Multi-layer Homoepitaxial Growth on Cu(100)Surface[J].ChinesePhysics,2006,15(4):807-812
[9]Elsholz F,Scholl E,Rosenfeld A.Control of Surface Rough-ness in Amorphous Thin-Film Growth[J].Appl Phys Lett,2004,84(21):4167
[10]王志永,朱志立,谷錦華.不同沉積速率微晶硅薄膜生長(cháng)模式的蒙特卡洛模擬研究[J].真空科學(xué)與技術(shù)學(xué)報,2010,30(6):632-635
[11]雒向東.磁控濺射Cu膜的表面形貌演化研究[J].工藝技術(shù)與材料,2007,32(2):138-141
