薄膜的形成-14-2012.ppt

薄膜的形成——薄膜形成过程的计算机模拟 ★ 薄膜形成过程的计算机模拟 Monte Carlo方法计算机模拟 蒙特卡罗原本是摩纳哥的一个赌城的名字，赌博中的“掷骰子”当然具有随机性，而由此抽象出来的用于科学研究的蒙特卡罗方法便是一种统计模拟随机抽样的方法，已广泛用于模拟各种输运现象，解决各种物理问题。 它的基本思想是，为了求解数学、物理、工程技术以及受产管理等方面的问题，首先建立一个概率模型或随机过程，位它的参数等了问题的解；然后通过对模型或过程的观察或抽样试验来计算所求参数的统计特征，最后给出所求解的近似值。 实例：薄膜二维生长的 Monte Carlo模拟 薄膜的形成——薄膜形成过程的计算机模拟 建模及其模拟 材料性质分析 模拟结果分析 Monte Carlo 模拟的过程 薄膜的形成——薄膜形成过程的计算机模拟 本模拟采用Monte Carlo方法模拟二维平面上粒子的单层膜的生长成膜情况。所考虑的成膜过程是大量微观粒子在给定宏观约束条件下的集体行为,但就每一个粒子的行为而言,则是随机的。 结合粒子的随机运动,根据过程的物理特性设定一些概率规则,理论模拟结果与实际情况进行比较。 用计算机模拟研究了薄膜生长机理,给出了微观生长过程的详细图像。 薄膜的形成——薄膜形成过程的计算机模拟 1. 利用周期性边界条件建立一个100×100的二维方格点阵。 2. 每次有一个粒子随机落到衬底上。若粒子落到第二层上则根据下落的位置,按照一定规律扩散到衬底上。 落下的粒子只能位于某个阵点所在的坐标处 薄膜的形成——薄膜形成过程的计算机模拟 3. 考虑到衬底点阵原子势能的极小值应当位于各阵点的正上方,所以落下的粒子只能位于四个最近邻位置的某个阵点所在的坐标处,而不会落在阵点之间的位置。当新位置的能量高于原来位置的能量时,新位置就有一定的概率被接受。 沉积粒子及其4个最近邻位置 薄膜的形成——薄膜形成过程的计算机模拟 薄膜的形成——薄膜形成过程的计算机模拟 薄膜的形成——薄膜形成过程的计算机模拟 薄膜的形成——薄膜形成过程的计算机模拟 分子动力学计算机模拟 Moleculer dycamics的基本思想： 首先须建立一组分子的运动方程，通过直接求解体系中的一个个分子的运动方程得到各个分子的位置和动量，即得到在相空间的轨迹，从而得到体系的静态和动态性质，进而得到体系的宏观性质。 实例：纳米磁性薄膜的分子动力学模拟 薄膜的形成——薄膜形成过程的计算机模拟 基本假设： 1. 在二维情况下，即薄膜的厚度相比面积而言可忽略。 2. 所有的分子团簇的形状假定为圆形。 3. 在L-J势的前提下，忽略排斥作用，只有吸引作用。 4. 在某次模拟薄膜生长的过程中，没有外界的分子团簇供给。 5. 只考虑最近邻相互作用。 薄膜的形成——薄膜形成过程的计算机模拟 1. 本模型模拟的是纳米磁性薄膜的生长过程，主要是在二维的平面上粒子数一定，温度恒定，压强不变的条件下系综的模拟。在模拟过程中，将粒子理想化地看成是理想的刚性小球，在二维平面上即为圆形。 本模型的建模过程: 2. 模拟的初始条件：在模拟的初始态，认为在200*200像素的坐标平面内随几分布了N（100） 个 点 。刚性小球的半径也是在5nm范围内随机产生的分布。（用一个象素表示1nm） 3. 小球的相互作用势为L-J势，为了简化，不考虑排斥作用，只考虑了吸引作用，且只考虑最近邻原子间的相互吸引来决定系统的能量，所以两个团簇是相向运动的。 薄膜的形成——薄膜形成过程的计算机模拟 模拟过程中,为了问题的简化及编程的实现，两个大小任意的分子团簇在体积（二维下模拟时为面积）相等的情况下相遇时将合二为一，而且运用了Monte Carlo方法。 5. 在合并的过程中，有Monte Carlo 过程来决定是否合并为单个分子团簇，这通过计算转换概率w来实现，在Monte Carlo 过程中，由于现在模拟能力的限制，无法计算出能量变化量，所以先假定能量变化量为已知的常数，根据能量最低原理计算！ 6. 分子及分子团簇在合并过程的运动有分子动力学中的L-J对势来决定。同时将合并的过程时序上离散为一定间隔的步长.在每一个离散态上观察分子及分子团簇的状态。 薄膜的形成——薄膜形成过程的计算机模拟 分子作用势的选取：引用的是L-J势，主要考虑到在本次的模拟中，由于计算和程序能力的限制，我们模拟的n100个点的情况。所以在模拟的二维平面上可以看成是很稀薄的分子分布，这正好符合L-J势的适用范围；可以得到相对简单的势函数形式，简化了程序和计算。 薄膜的形成——薄膜形成过程的计算机模拟 具体的模拟步骤： 确定要模拟的分子团簇数目，将分子团簇随机分布在限定区域内。