图4-2 表面原子在2600k下随时间变化原子均方位移改变的MSD曲线
1表面1 2表面2 3表面3 4表面4
5表面5 6表面6 7表面7 8表面8
图4-1-2本实验模型，我们一共建立八层面，考虑到前面失败的原因和各种干扰因素，我们对模型进行了如下改善为，其中后二层固定，前六层为自由面，升高温度，考察其表面的熔化情况，在2600K时，我们对模型进行优化，释放该模型中的因为快速加热带来的应力，然后用MSD对运算结果进行统计分析为。由图可知，在2600K时，原子的第一层至第六层MSD随时间的几乎没变化。这表明原子层没有扩散运动，只是一种热振动，还保持相对完整的晶体状态。
31
12108MDS(A) 1 2 3 4 5 6 7 826420020406080100120t(ps)
图4-3 表面原子在2700k下随时间变化原子均方位移改变的MSD曲线
1表面1 2表面2 3表面3 4表面4 5表面5 6表面6 7表面7 8表面8
图4-1-3是在2700K时，原子第一层的MSD随时间的变化发生显著的增加。第二层及以下几层随时间增加变化不明显。这表明第一层原子具有一定的扩散运动，其余各层保持的只是热振动状态，这说明HfO2虽然试验熔点是3030K,但在2700K左右就开始熔化。
32
图4-4 表面原子在2800k下随时间变化原子均方位移改变的MSD曲线
1表面1 2表面2 3表面3 4表面4
5表面5 6表面6 7表面7 8表面8
图4-1-3是在2800K时，原子前五层的MSD随时间的变化发生显著的增加。第六层及以下几层随时间增加变化不明显。这表明前五层原子具有一定的扩散运动，其余各层保持的只是热振动状态，说明HfO2试验熔点是3030K,在2800K左右较2700K熔化程度加大。
虽然本实验的目的是为了研究在金属氧化物HfO2的熔化过程中是否像水一样在其表面有“类液”情况的出现，但因时间的限制，在本实验阶段我们还无法能把本实验做完整，还有待于模拟运算一下2900K,3000K的情况。基于此原因，我们只能算出以上温度点的熔化情况。
33
参考文献
[1] G E Moore. Electronics 38:114, 1965.
[2] The International Technology Roadmap for Semiconductors. 2008. [3] P Darmawan et al., Appl. Phys. Lett. 93:062901, 2008.
[4] 方泽波等. Al2O3高k栅介质的可靠性[J] .无机材料学报 23:357, 2008.
[5] 卢红亮等. 原子层淀积 Al2O3薄膜的热稳定性研究[J].无机材料学报 21:1217, 2006. [6] C J Forst et al., Nature 427: 53, 2004.
[7] K Y Tse and J Robertson, Phys. Rev. Lett. 99:086805, 2007.
[8] 卢红亮等. Mn/GaAs(100)界面电子结构的同步辐射光电子能谱研究.物理学报[J].55:1374, 2006.
[9] E. Botcharova, M. Heilmaier, J. Freudenberger, G. Drew, D. Kudashow, U. Martin and
L. Schultz.Supersaturated solid solution of niobium in copper by mechanical alloying Journal of Alloys and Compounds Volume 351, Issues 1-2, (10 March 2003) Pages 119-125 [10] 李新宇.纳米金属多层膜中离子束混合相变机制的研究.[博士学位论文]. 北京：清华大学
材料科学与工程系，2006 年
[11] Tsaur B Y, Mayer J W, Tu K N. Ion-beam-induced metastable Pt2Si3 phase: Iformation,
structure and properties. J. Appl. Phys., 1980, 51: 5326~5333
[12] Zhang Z J, Liu B X, Thermodynamics and crystallographic mechanism of metastable
phase formation in the Nb-Fe, Nb-Co and Nb-Ni systems by ion mixing. J. Appl. Phys., 1994, 75: 4948~4952
[13] 邵阳．超导和庞磁阻材料的计算机模拟研究：[硕士学位论文]．北京：清华大学，2004
34
[14] 罗熙淳,梁迎春,董 申. 分子动力学在纳米机械加工技术中的应用[J] . 中国机械工
程,1999 ,10 (6) :692.
[15] 胡壮麒,王鲁红,刘 轶. 电子和原子层次材料行为的计算机模拟[J] . 材料研究学
报,1998 ,12 (1) :1.
[16] 卢百平,张丰收,傅恒志. 材料科学中的分子动力学模拟[A] . 编者西北工业大学2002 年
材料论坛论文集[ C] . 西安: ,2002. 124.
[17] Daw N S ,Baskes M I. Embedded atom method derivationand application to impurities ,
surfaces , and defect s in metals[J ] . Physical Review B , 1984 ,29 (12) :848628495. [18] 郭俊梅,邓德国,潘建生,等. 计算材料与材料设计[J ] . 贵金属,1999 ,20 (4) :62. [19] 董科军,刘让苏,郑采星,等. 液态金属凝固过程分子动力学模拟的并行算法研究[J] . 计算
机科学与工程,2000 ,22(3) :45.
[20] 王鲁红,刘浩哲,胡壮麒. 金属的液态结构及快速凝固模拟[J ] . 金属学报,1998 ,34
(9) :933.
[21] 吴兴惠，项金钟．现代材料计算与设计教程[M]．北京：电子工业出版社，2002．298~300 [22] Stanek C R, Atomic scale disorder in fluorite and fluorite related oxides. Ph.D. thesis, Imperial College, 2003
[23] Pauling L and Wilson E B, Introduction of quantum mechanics, with applications to chemistry, McGraw Hill, New York, 1985
[24] 陈舜麟．计算材料学[M]．北京：化学工业出版社，2005．135~143
[25] 黄昆（原著），韩汝奇（改编）．固体物理学[M]．北京：高等教育出版社，1988．49~55 [26] Ahn K H and Millis A J, Effects of uniaxial strain in LaMnO3. Phys[J]. Rev. B 2001, 64:115103(12p)
35