Forcite分子动力学模块计算温度对带隙和光谱影响.pptVIP

* 本课程学习内容 计算石英晶体晶胞体积与温度的关系，了解分子动力学的初步使用。 二. 计算硅晶体晶胞体积与温度的关系，温度对带隙和光谱的影响。 石英熔点～1800K 作业1：查水晶c轴的定义？出处？ 作业2：根据?的数据，估算1800K时，?－石英晶胞的体积是多少？与动力学计算值比较。 MM80：Mirwald and Masone (1980) 静压实验。 作业3：找Tridymite结构 一. 计算石英晶体晶胞体积与温度的关系，了解分子动力学的初步使用。 1. 在F盘中建立文件夹class0920 2. 打开Materials Studio，在class0920中新建名为sio2的Project。 3. 由内部数据库，输入?、?-quartz的结构 导入?、?两种石英晶体结构 先优化?-SiO2的结构。选择Quartz-beta为当前文件。 选择分子动力学计算模块Forcite。 从Forcite的下拉菜单中打开Calculation对话框，选择Setup选项卡，设置Task为结构优化，计算精度Quality为Ultra-fine。 单击More按钮，在打开的对话框中勾选Optimize cell，设置迭代次数为500。 4. 用分子动力学计算模块Forcite优化两种石英晶体的结构 打开Energy选项卡，设置力场Forcefield为COMPASS27，计算精度Quality为Ultra-fine。 打开Job Control选项卡，选用本地机My Computer进行计算，设置并行数＝2，单击Run进行结构优化。 同样对?-SiO2进行结构优化，Project中出现两个新文件夹。 由File/Save Project保存文件；再由Window/Close all关闭窗口文件。 双击四个结构文件 由Window/Tile Vertically排列四个结构图。1、3分别是实验确定的?-SiO2 和?-SiO2；2、4分别是结构优化的?-SiO2 和?-SiO2。 3 2 4 1 依次使这四个结构文件为当前文件，同时在右下角的Properties 中选 Lattice 3D，显示相应结构的体积和空间群。这些信息列于结构图的左右两侧。可以看出2、4的体积一样。 依次选中一个结构图，用Find Symmetry找其对称性，信息列于结构图的上下两侧。 比较四个图的信息可知，1、3未变，2、4的体积一样，空间群也一样。即结构优化后， ?-SiO2 变为?-SiO2。这样我们就在不同的温度下计算?-SiO2的性质。 结构优化后， ?-SiO2转变为 ?-SiO2，体积一样。 1 2 1 4 3 关闭其它文件，仅保留结构图4，即结构优化后的?-SiO2的结构。 按右图所示的步骤，Edit/Select All/Copy，选择结构优化的?-SiO2的结构。 鼠标移至Project的SiO2上，按右键，在下拉菜单中选择3D Atomistic Document，工作窗口中出现一个空的3D Atomistic .xsd。 5. 计算不同温度下?-SiO2晶胞的体积 由Edit / Paste，将结构优化后的?-SiO2的结构复制到新的3D Viewer上。 在3D Viewer的空白处单击左键，取消选择。 鼠标移至Project的3D Atomistic.xsd上，按右键，单击下拉菜单中的Rename，将文件名改为sio2_100.xsd。 同样，在Project下继续复制结构优化后的quartz_beta.xsd，分别将这些新的3D Atomistic.xsd重新命名为SiO2_200 ~ SiO2_2000.xsd。 课堂时间少，同学作4个： SiO2_300.xsd、 SiO2_800.xsd、 SiO2_1600.xsd 、SiO2_1800.xsd 收敛性测试～模拟计算时间与晶胞体积的关系 单击More按钮，选择Dynamics标签，按右图进行设置。注意温度固定在1200K，其余设置不变。 选择分子动力学计算模块Forcite 。从Forcite的下拉菜单中打开Calculation对话框，选择Setup标签，设置Task为动力学(Dynamics)，计算精度Quality为medium。 计算1200K下?-SiO2的体积，比较模拟时间对计算出的体积数值的影响，以设置合适的模拟时间。 结构文件为 sio2_1200.xsd，设为当前文件。 打开Job Control选项卡，选用本地机My Computer进行计算，设置并行数＝2，单击Run进行结构优化。 打开Energy选项卡，设置力场Forcefield为COMPASS27，计算精度Quality为Medium。 计算结束后，P