黃美嬌臺灣大學：機械工程學研究所洪榮泰Hong, Rong-TaiRong-TaiHong2007-11-282018-06-282007-11-282018-06-282006http://ntur.lib.ntu.edu.tw//handle/246246/61446本論文以平行分子動力學模擬，模擬在物理氣相沉積下， 銅金屬原子以濺鍍法方式沉積在導線凹槽內之銅大馬士格製程。 分別研究大馬士格製程中銅金屬沉積在凹槽的形貌及微結構之變化。 所探討之凹槽內擴散阻絕層結構分別為: 鈦單晶六面密堆疊表面指數(0001)， 以及鉭單晶體心面之表面指數(100)、 (110)及(111)。在形貌方面 ， 主要探討基板溫度、 入射動能及凹槽之展弦比參數下， 對不同阻絕層之材料下之覆蓋百分比。 在微結構方面， 對鉭擴散 阻絕層不同的晶格表面下， 分析銅對鉭擴散阻絕層之沉積紋路及銅原子與擴散阻絕層混合情形。 在分子動力學模擬中 ， 所使用之原子間勢能函數為Johnson的解析鄰近嵌入式原子方法。 基於勢能函數轉換不變性的理論下， 來處理本文銅原子與擴散阻絕層材料Cu/Ta及Cu/Ti之合金勢能之計算。 由結果可知 ， 增加基材溫度及入射銅原子之動能 ，可提昇覆蓋程度。由於鉭相對於鈦具有較佳熱穩定性及較深的勢能井，銅原子之沉積覆蓋相對較高於以鈦為主之擴散阻絕層。此外 ， 鈦之合金的形成相較於鉭為主來得嚴重。 對鉭不同表面指標下，具有表面開口結構之(100)鉭晶格下， 其覆蓋程度相對較差， 這對在高展弦比下所形成之效果更為明顯。對於沉積紋路所形成方向 ， 在較高表面能量之擴散阻絕層下， 所形成之紋路較傾向水平，反之則以垂直纖維形成。 此外 ， 基材溫度也會影響沉積紋路之方向。 本文使用叢集電腦環境來發展凹槽沉積平行分子動力學模擬:空間分割及原子分割法之平行程式碼， 比較上述平行算則，在不同的叢集電腦系統之適用性及效能。由本文結果， 由於空間分割法在資料交換的次數較低， 其效能大於原子分割法。The trench filling and microstructures of depositing copper atoms on the titanium and tantalum diffusion barrier layers in a damascene process have been studied using parallel molecular dynamics simulation with the embedded atom method (EAM) as interaction potential for the present alloy metal system which is based on the invariance transformation of alloy system. The effects of different process parameters including incident energies of depositing atoms, substrate temperature, the different surface indices of tantalum and the aspect ratio of trench were investigated. The surface indices of tantalum layer considered were (100), (110) and (111) and the aspect ratios were 1, 1.5 and 2, respectively. The coverage of trench filling is improved as increasing the substrate temperature and incident copper atoms energy. Comparing with titanium and tantalum diffusion barrier layers under the same process parameters, it is found that due to a better thermal stability significant improvement in coverage percentage can be obtained using the tantalum barrier layer, especially at higher substrate temperatures. The enhancement of trench filling is studied in case of (110) and (111) crystal structure of tantalum because there are fewer open structures in sidewall of diffusion barrier layers. In microstructure of trench filling, the alloy formation is serious in the case of titanium than that of tantalum diffusion barrier layers. The orientations of textures are changed with different surface energies and the deposited copper within trench has fiber structure with lower surface energy. To treat a larger system, parallel molecular dynamics simulation is also implemented and we compare the efficiency of different PC cluster computing environments including homogeneous and heterogeneous computer systems. The clarification of the suitability is studied between such systems. In this work, the different parallel schemes are developed including atomic and spatial decomposition methods. Comparing the efficiencies between the two schemes, the spatial decomposition performs better in this work.1 緒論 1.1 引言 . . . . . . . . . . . . . . . . . . . . . . . . 1 1.2 大馬士革製程. . . . . . . . . . . . . . . . . . . 2 1.3 研究動機. . . . . . . . . . . . . . . . . . . . . . 4 1.4 文獻回顧 . . . . . . . . . . . . . . . . . . . . 5 1.5 …本文研究範圍. . . . . . . . . . . . . . . . 8 2 分子動力學模擬 2.1 歷史回顧. . . . . . . . . . . . . . . . . . . . . . 10 2.2 理論背景. . . . . . . . . . . . . . . . . . . . . . 11 2.3 作用勢能函數 . . . . . . . . . . . . . . . . . . . 13 2.4 數值方法. . . . . . . . . . . . . . . . . . . . . . 13 2.5 模擬方法. . . . . . . . . . . . . . . . . . . . . . 15 3 嵌入原子法 3.1 原子勢能理論 . . . . . . . . . . . . . . . . . . . 20 3.2 嵌入原子法 . . . . . . . . . . . . . . . . . . . . 21 3.3 Johnson嵌入原子法 . . . . . . . . . . . . . . . 23 4 平行分子動力學模擬 4.1 平行演算法. . . . . . . . . . . . . . . . . . . . 33 4.2 系統配置 . . . . . . . . . . . . . . . . . . . . . 36 4.3 效能分析. . . . . . . . . . . . . . . . . . . . . . 38 5.模擬模式 5.1 作用勢能模式. . . . . . . . . . . . . . . . . . . 47 5.2 表面能量之計算. . . . . . . . . . . . . . . . . . 50 5.3 徑向密度分佈函數. . . . . . . . . . . . . . . . . 52 5.4 導線凹槽沉積模式. . . . . . . . . . . . . . . . . 54 6 導線沉積模擬結果與討論 6.1 鈦擴散阻絕層沉積形貌. . . . . . . . . . . . . . . 59 6.2 鉭擴散阻絕層. . . . . . . . . . . . . . . . . . . 67 6.3 鉭擴散阻絕層(110)及(111)表面 . . . . . . . . 76 6.4 鉭擴散阻絕層高展弦比銅之沉積. . . . . . . . . . . 85 6.5 導線沉積微結構. . . . . . . . . . . . . . . . . . 95 7 結論與展望 7.1 結論 . . . . . . . . . . . . . . . . . . . . . . . 104 7.2 展望 . . . . . . . . . . . . . . . . . . . . . . . 105 附錄A ...................................116 附錄B ...................................119 附錄 C ..................................12122076645 bytesapplication/pdfen-US大馬士革製程擴散阻絕層平行分子動力學模擬Damascene processdiffusion barrier layerparallel molecular dynamics simulationthesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/61446/1/ntu-95-D87522002-1.pdf