用第一性原理计算C(碳)的结构与物理参数.docxVIP

用第一性原理计算C(碳)的结构与物理参数1. 实验目的掌握第一性原理和密度泛涵的计算方法；学会使用Visualizer 的各种建模和可视化工具；熟悉CASTEP模块的功能。2. 实验原理CASTEP是基于密度泛涵理论平面波赝势基础上的量子力学计算。密度泛涵理论的基本思想是原子、分子和固体的基本物理性质可以用粒子密度函数进行描述。可以归纳为两个基本定理：定理1：粒子数密度函数是一个决定系统基态物理性质的基本参量。定理2：在粒子数不变的条件下能量对密度函数变分得到系统基态的能量。不计自旋的全同费米子的哈密顿量为：其中动能项为：库仑作用项为：V为对所有粒子均相同的局域势u(r)表示的外场影响：粒子数密度函数为：对于给定的，能量泛函定义为：；系统基态的能量：3. 实验内容A, 碳的弹性模量计算B,计算材料的晶格参数，并画出所计算材料的能带图 C,计算碳的电子结构计算； 4. 实验设备和仪器硬件：多台PC机和一台高性能计算服务器。软件：主要利用Materials studio软件包里的Materials Visualizer和CASTEP模块。5. 实验步骤一,建立所研究材料的结构模型①按照所研究材料的晶胞参数建立晶体结构。方法是可以自己建立模型，也可以在软件自带的结构库中直接读入结构文件。如：在菜单栏中选择File | Import，进入structures/文件夹，选择c.msi文件。②按照计算需要对晶胞内的原子进行替换，并用Build | Symmetry | Primitive Cell将模型设置为原胞形式。二,设置并运行量子力学计算在工具栏中选择CASTEP模块。其对话框如下图。A, 碳的弹性模量计算1：几何优化①选择Setup项，将Task项设置为 Geometry Optimization，标准设为 Fine,点击此对话框中的 More按钮，在Optimize Cell项前打钩，关闭对话框。 ②选择 Job Control按钮，点击此对话框中的 More按钮，设置Gateway，选择在本机计算还是在服务器上计算。③选择 Properties 按钮，确定没有选择任何项。按下Run按钮，开始运行计算。关闭对话框。优化成功：优化曲线：2：判断计算结果是否正确①在 Project Explorer工作栏内,双击.castep文件，将其激活。 ②在菜单栏中选择Edit | Find 键入 completed successfully?查找到此文字，找到后向上几行，有一个两行的表格，如果两行显示都为“yes”，说明计算正常结束。3：计算弹性常数.在CASTEP对话框中选择 Setup 项，在Task 的下拉单中选择 Elastic Constants，点击对话框中的 More按钮，设置每次拉伸的步数， 按下Run 按钮。正在计算：计算完成： 3计算结果分析在计算结果文件夹中激活.xsd文件。在工具栏中选择 CASTEP工具，然后选择Analysis，选择Electron density项，按下Import按钮，电子的等值面就显示在结构中。 打开CASTEP的Analysis对话框，选择Elastic Constants项，将CASTEP Cij文件夹中的.xsd文件激活，按下Calculate按钮。Cij和Sij参数：弹性模量结果：B, 计算材料的晶格参数，并画出所计算材料的能带图 1：几何优化①选择Setup项，将Task项设置为 Geometry Optimization，标准设为 Fine,点击此对话框中的 More按钮，在Optimize Cell项前打钩，关闭对话框。 ②选择 Job Control按钮，点击此对话框中的 More按钮，设置Gateway，选择在本机计算还是在服务器上计算。③选择 Properties 按钮在Band structure选项前打钩。按下Run按钮，开始运行计算。关闭对话框。计算完成：优化结果图一：优化结果图二：2：判断计算结果是否正确①在 Project Explorer工作栏内,双击.castep文件，将其激活。 ②在菜单栏中选择Edit | Find 键入 completed successfully?查找到此文字，找到后向上几行，有一个两行的表格，如果两行显示都为“yes”，说明计算正常结束。 3计算结果分析在计算结果文件夹中激活.xsd文件。在工具栏中选择 CASTEP工具，然后选择Analysis，选择Electron density项，按下Import按钮，电子的等值面就显示在结构中。 打开CASTEP的Analysis对话框，选择Band structure项，按下View按钮；晶格参数：能带结果图：C,计算碳的电子结构计算；1：几何优化①选择Setup项，将Task项设置为 G