奈米線尺寸效應

3.3.1 金塊材模擬

首先建立一個 15a x 15a x 15a 的的金塊材，如圖 3.3，在溫度為 1K 下執行 NVT 系踪，拉伸速率為 2x10⁷s^-1，施加 X 方向應變後所產 生之三方向應力_xx、_yy、_zz，藉由應力應變關係將這些值帶入(3.3.1) 式中可求得勁度矩陣(Stiffness matrix)中之 C11 、C21、與 C31 (在應變 為 0.01 時求得之彈性模數)。再分別依照此流程施加應變於 Y、Z 方

23

我們將初始奈米線結構在NPT系綜(constant pressure algorithm)，

溫度為1K，壓力為0.0 GPa，平衡時間為100ps的情況底下執行運算，

則我們可以得到X、Y、Z三方向的應力隨時間變化，此時三方向應力 收斂至零。同時我們可以得到一沿著軸向平衡應變量(Equilibrium strain)，此應變量與奈米線截面積大小有關。

24

在討論奈米線的尺寸效應時，我們採用四種不同大小的截面尺 寸，分別為6a x 6a、8a x 8a、12a x 12a、16a x 16a、20a x 20a (其中a 為晶格常數)，在拉伸前執行NPT系綜，使結構達到穩定鬆弛，6a及 8a平衡過程三軸應力及能量如圖3.4~3.7。及其平衡後結構圖3.8~3.11。

由結構圖可以發現在不同線寬下，其平衡穩定結果會略有所不 同，此結果為表面應力造成其平衡應變量(Equilibrium strain)[13][46]

的不同，我們可以得到一應變量對時間的關係圖，如圖3.12。可以明 顯發現，當線寬在六倍晶格常數時其結構在平衡時已經產生降伏，由 文獻[13]可知，奈米線的壓縮應力與其截面積有關，當截面積越小，

其壓縮應力會越大，而當線寬小於一臨界值時，奈米線承受不住此應 力大小而產生降伏。

3.3.3 奈米線拉伸與討論

拉伸結果的討論由3.3.2節平衡後得到的結構進一步做拉伸模 擬，排除了六倍晶格線寬，因其結構已產生降伏故我們不另做討論。

分析8a x 8a、12a x 12a、16a x 16a、20a x 20a 四種不同線寬在拉伸速 率為2E+07s^-1下的結果，其應力應變關係圖如3.13，數據如表 4。

此外，我們另外利用單軸NST的方法去做拉伸比較其結果，在8a x 8a的結構下可以得到降伏強度約為4.45GPa，降伏應變約為0.0927，

25

此結果與利用位移控制拉伸的方法誤差很小。

由應力應變圖可發現，當奈米線截面越大時，其結構楊氏係數會 越大，如圖3.14，因為奈米線側向自由表面部份佔整體結構的比例會 隨著截面積變大而變小，又因其特殊的環境，自由表面周遭的原子比 內部原子少而鍵結又較弱[7]，故其楊氏係數會上升並趨近於塊材的 值。另一方面，降伏應力則會下降，如圖3.15，其主要原因為截面越 大 ， 結 構 自 由 表 面 積 則 會 越 大 ， 所 能 提 供 的 差 排 成 核 機 會 越 高 [6][10]，造成其降伏應力及應變都會有下降的趨勢。

在拉伸的過程中，因兩側原子受到二端剛體的牽引，因此結構可 以保持較規則的排列，相對的，靠近中央的原子，因兩側的拉伸而必 需向內側填補，使材料形成頸縮(necking)的現象，結構的拉伸變形如 圖3.16。

其拉伸過程中，利用BDT local stress計算原子沿 X 軸方向局部 應力變化，觀察 Z=15Å (前後約三層原子厚度)的平面在降伏點附近 的應力變化，如圖3.17，在降伏前一瞬間(t=10ps)，結構達到一高應力 且 飽 和 的 狀 態 ， 隨 後 內 部 原 子 相 互 拉 扯 ， 在 降 伏 發 生 的 時 間 點 (t=20ps)，結構應力被釋放，直到t=50ps時，除了差排發生的附近會有 較高的應力，其餘部份相對較低。

由差排理論[47]可知，面心立方(FCC)金的滑動平面為{111}的最

26

緊密堆積平面，其主要原因為在此面上的單位面積原子密度最高，原 子間彼此間距小，鍵結力強，需要較大的力才能夠使其分開，故在此 方向上的原子只會有相對移動產生[48]。而每個滑移面分別搭配三個 滑移方向，因此總共會有12組滑移系統，如表13 。另一方面，利用 Schmid’s law 得到的 Schmid’s factor在[100]晶向最大值為0.408，如 表13，與模擬結果求得的值相同。應變率在 2x10⁷ s^-1、截面為 8a x 8a，如圖3.18。我們可以看到在 t=10 ps 到 t=20 ps 之間奈米線發生 了差排成核，進一步可以觀察在這段時間內的詳細情況，如圖3.19，

從圖(a)可以觀察在 t=18 ps 時，在奈米線自由表面附近開始有原子擾 動，這即是差排開始發生的位置，與文獻[13][44]敘述差排成核位置 由奈米線表面發生符合，而到圖(e) t=22ps 之後由表面慢慢往內部發 展形成滑移面。

在降伏區間內，應力成鋸齒狀分佈[4]，示意圖如圖3.20，彈性範 圍內，原子會些微偏離原平衡位置，但與其相鄰原子相對位置不會改 變，當超過彈性限度原子平面產生滑移時，結構形成差排(路徑1)，

結構從相對穩定的排列(a圖)至(b圖)，位勢能會上升同時應力上升達 到一高峰值，此為另一降伏應力，當滑移平面越過一個原子達位能最 低點時(路徑2)，拉伸造成的差排消失，系統處於另一穩定結構下，

同時應力值達到一低值，系統不斷重覆此步驟會形成鋸齒狀應力應變

27

關係圖，而當應變持續增加時，因頸縮的現象造成截面積會越來越 小，差排面的面積也越來越小，能形成的差排應力也相對變小，因此 能達到的高峰應力值會持續下降。此應力特性在此以應變率低於 10⁷~10⁸ s^-1 時較明顯，因當應變率太大時，因塑性變形滑移量大，無 法形成滑移和穩定的過程，此時應力曲線呈現平緩地下降。