《《CASTEP概述》PPT课件.ppt》由会员分享,可在线阅读,更多相关《《CASTEP概述》PPT课件.ppt(73页珍藏版)》请在金锄头文库上搜索。
1、CASTEP概述 CASTAP是特别为固体材料学而设计的一个现代的量子力学基本程序 其使用了密度泛函 DFT 平面波赝势方法 进行第一原理量子力学计算 以探索如半导体 陶瓷 金属 矿物和沸石等材料的晶体和表面性质 典型的应用包括表面化学 键结构 态密度和光学性质等研究 CASTAP也可用于研究体系的电贺密度和波函数的3D形式 此外 CASTAP可用于有效研究点缺陷 空位 间隙和置换杂质 和扩展缺陷 如晶界和位错 的性质 MaterialStudio使用组件对话框中的CASTAP选项允许准备 启动 分析和监测CASTAP服役工作 计算 允许选择计算选项 如基集 交换关联势和收敛判据 作业控制和文
2、档控制 分析 允许处理和演示CASTAP计算结果 这一工具提供加速整体直观化以及键结构图 态密度图形和光学性质图形 关于CASTAP CASTAP的任务 CASTAP计算是要进行的三个任务中的一个 即单个点的能量计算 几何优化或分子动力学 可提供这些计算中的每一个以便产生特定的物理性能 性质为一种附加的任务 允许重新开始已完成的计算以便产生最初没有提出的额外性能 在CASTAP计算中有很多运行步骤 可分为如下几组 结构定义 必须规定包含所感兴趣结构的周期性的3D模型文件 有大量方法规定一种结构 可使用构建晶体 BuildCrystal 或构建真空板 BuildVacuumStab 来构建 也可
3、从已经存在的的结构文档中引入 还可修正已存在的结构 注意 CASTAP仅能在3D周期模型文件基础上进行计算 必须构建超单胞 以便研究分子体系 提示 CASTAP计算所需时间随原子数平方的增加而增加 因此 建议是用最小的初晶胞来描述体系 可使用Build Symmetry PrimitiveCell菜单选项来转换成初晶胞 计算设置 合适的3D模型文件一旦确定 必须选择计算类型和相关参数 例如 对于动力学计算必须确定系综和参数 包括温度 时间步长和步数 选择运行计算的磁盘并开始CASTAP作业 结果分析 计算完成后 相关于CASTAP作业的文档返回用户 在项目面板适当位置显示 这些文档的一些进一步
4、处理要求获得可观察量如光学性质 CASTAP中选择一项任务 1从模块面板 ModuleExplorer 选择CASTAP Calculation 2选择设置表 3从任务列表中选择所要求的任务 CASTAP能量任务 CASTAP能量任务允许计算特定体系的总能量以及物理性质 除了总能量之外 在计算之后还可报告作用于原子上的力 也能创建电荷密度文件 利用材料观测仪 MaterialVisualizer 允许目测电荷密度的立体分布 还能报告计算中使用的Monkhorst Park的k点的电子能量 因此在CASTAP分析中可生成态密度图 对于能够得到可靠结构信息的体系的电子性质的研究 能量任务是有用的
5、只要给定应力性质 也可用于计算没有内部自由度的高对称性体系的状态方程 即压力 体积 能量 体积关系 注意 具有内部自由度的体系中 利用几何优化 GeometryOptimization 任务可获得状态方程 CASTAP中能量的缺损单位是电子伏特 eV 各种能量单位的换算关系见Mohr P J 2000 1eV 0 036749308Ha 23 0605kcal mole 96 4853kJ mole CASTAP几何优化任务 CASTAP几何优化任务允许改善结构的几何 获得稳定结构或多晶型物 通过一个迭代过程来完成这项任务 迭代过程中调整原子坐标和晶胞参数使结构的总能量最小化 CASTAP几何
6、优化是基于减小计算力和应力的数量级 直到小于规定的收敛误差 也可能给定外部应力张量来对拉应力 压应力和切应力等作用下的体系行为模型化 在这些情况下反复迭代内部 几何优化处理产生的模型结构与真实结构紧密相似 利用CASTAP计算的晶格参数精度列于右图 应力张量直到与所施加的外部应力相等 状态方程计算 在所施加静压力下几何优化可用于确定材料的体模量B和对压力的导数B dB dP 过程包括计算理论状态方程 EOS 该方程描述单胞体积于外部静压力的关系 工艺非常类似于真实实验 使用几何优化对话框中的应力列表将外部压力固定 通过进行几何优化可以找到在此压力下的单胞体积 随后的P V数据分析与实验研究精确
7、一致 描述EOS选择分析表达式 其参数适于计算数据点 最流行的EOS形式是三阶Birch Murnaghan方程 式中V0为平衡体积 Cohen等进行了EOS各种解析式的的详细比较研究 注意 从相应实验中获得的B和B 值依赖于计算使用的压力值范围 利用金刚石压砧获得的实验值通常在0 30GPa范围内 因此推荐理论研究也在这个范围内 在研究中避免使用负压力值也很重要 此外 用于生成P V数据序列的压力值可能是不均匀的 在低压力范围要求更精确采样以便获得体模量精确值 P V 几何优化方法 在缺损条件下 CASTAP使用BFGS几何优化方法 该方法通常提供了寻找最低能量结构的最快途径 这是支持CAS
8、TAP单胞优化的唯一模式 衰减分子动力学 Dampedmoleculardynamics 方法是另一种可以选择的方法 该方法对具有平滑势能表面的体系如分子晶体或表面分子与BFGS同样有效 CASTAP动力学任务 CASTAP动力学任务允许模拟结构中原子在计算力的影响下将如何移动 在进行CASTAP动力学计算以前 可以选择热力学系综和相应参数 定义模拟时间和模拟温度 选择热力学系综 对牛顿运动定律积分允许探索体系恒值能量表面 NVE动力学 然而 在体系与环境进行热交换条件下发生最本质的现象 使用NVT系综 或者是确定性的Nos 系综或者是随机性的Langevin系综 可模拟该条件 定义时间步长
9、timestep 在积分算法中重要参数是时间步长 为更好利用计算时间 应使用大的时间步长 然而 如果时间步长过大 则可导致积分过程的不稳定和不精确 典型地 这表示为运动常数的系统偏差 注意 量子力学分子动力学计算要求比力场动力学使用更小的时间步长 动力学过程的约束 CASTAP支持LangevinNVT或NVE动力学过程的线性约束 然而 借助MaterialStudio界面可以近似使用以下两种更基本的约束 质心固定 单个原子固定 使用seedname cell文档可以利用更复杂的约束 CASTAP性质任务 CASTAP性质任务允许在完成能量 几何优化或动力学运行之后求出电子和结构性质 可以产生
10、的性质如下 态密度 DOS 利用原始模拟中产生的电荷密度和势能 非自恰计算价带和导带的精细Monkhorst Pack网格上的电子本征值 带结构 利用原始模拟中产生的电荷密度和势能 非自恰计算价带和导带的布里渊区高对称性方向电子本征值 光学性质 计算电子能带间转变的矩阵元素 CASTAP分析对话可用于生成包含可以测得的光学性质的网格和图形文件 布局数分析 进行Mulliken分析 计算决定原子电荷的键总数和角动量 以及自旋极化计算所需的磁矩 任旋地 可产生态密度微分计算所要求的分量 应力 计算应力张量 并写入seedname castep文档 如果要进行单胞参数固定时进行几何优化运行和要检查点
11、阵偏离平衡的程度 这些信息是有用的 例如 可进行符合于给定体系理论基态的固定单胞的点缺陷的超晶胞研究 几何优化后的应力值显示了与超单胞近似相关联的弹性效应 注意 为计算某种性质 从适当模拟得到的结果文档必须以当前的文件夹形式出现 用第一原理预测AIAs的晶格参数 本指南主要是阐明在MaterialsStudio当中如何运用量子力学来测定物质的晶体结构 你将从中学到如何构建晶体结构以及如何设置CASTEP几何优化运行和分析结果 本指南的内容如下 1构建AlAs的晶体结构2设置和运行CASTEP中的计算3分析结果4比较实验数据和结构注意如果你的服务器没有足够快的CPU 本指南限制使用CASTEP进
12、行几何优化计算 因为它会占用相当长的时间 1构建AlAs的晶体结构 为了构建晶体结构 我们需要知道你想要构建的晶体的空间群信息 晶格参数以及它的内部坐标 以AlAs为例 它的空间群是F 43m或空间群数字是216 它有两种基本元素Al和As 其分数坐标分别为 000 和 0 250 250 25 它的晶格参数为5 6622埃 第一步是构建晶格 在Projectexplorer的跟目录上右键单击 选中New 3DAtomisticDocument 接着在3DAtomisticDocument右键单击 把它更名为AlAs 从菜单中选择Build Crystals BuildCrystal 然后显示
13、出BuildCrystal对话框 如下 在Entergroup中选择F 43m或在Entergroup中单击 然后键入216 再按下TAB键 空间群信息框中的信息也随着F 43m空间群的信息而发生变化 选择LatticeParameters标签 把a的数值从10 00改为5 662 单击Build按钮 一个没有原子的晶格就在3Dmodeldocument中显示出来 现在我们就可以添加原子了 从菜单栏中选择Build AddAtoms 通过它 我们可以把原子添加到指定的位置 其对话框如下 在AddAtoms对话框中选择Options标签 确定Coordinatesystem为Fractional
14、 如上所示 选择Atoms标签 在Element文 本框中键入Al 然后按下Add按钮 铝原子就添加到结构中了 在Element文本框中键入As 在a b c文本框中键入0 25 按Add按钮 关闭对话框 原子添加完毕 我们再使用对称操作工具来构建晶体结构当中剩余的原子 这些原子也显示在邻近的单胞中 当然 我们也可以通过重新建造晶体结构来移去这些原子 从菜单栏中选择Build Crystals RebuildCrystal 按下Rebuild按钮 在显示出的晶体结构中那些原子就被移走了 我们可以把显示方式变为BallandStick 在模型文档中右键单击 选择DisplayStyles 按下B
15、allandstick按钮 关闭对话框 在3D视窗中的晶体结构是传统的单胞 它显示的是格子的立方对称 如果存在的话 CASTEP使用的则是格子的全部对称 既包含有两个原子的原胞和包含有8个原子的单胞是相对应的 不论单胞如何定义 电荷密度 键长 每一类原子的总体能量都是一样的 并且由于使用了较少的原子 使计算时间得以减少 从菜单栏中选择Build Symmetry PrimitiveCell 在模型文档中显示如下 AlAs的原胞 2设置和运行CASTEP中的计算 从工具栏中选择CASTEP工具 再选择Calculation或从菜单栏中选择Modules CASTEP Calculation CA
16、STEPCalculation对话框如下 下面我们将要优化它的几何结构 把Task改为GeometryOptimization 把Quality改为Fine 优化当中的默认设置是优化原子坐标 尽管如此 在本例中我们不仅要优化原子坐标也要优化晶格 按下Task右侧的More 按钮 选中OptimizeCell 关闭对话框 当我们改变Quality时 其他的参数也会有所改变来反映Quality的改变 选择Properties标签 可从中指定我们想要计算的属性 选中Bandstructure和Densityofstates 另外 我们也可以具体指明jobcontrol选项 例如实时更新等 选择JobControl标签 选中More 按钮 在CASTEPJobControlOptions对话框中 把Update的时间间隔改为30秒 关闭对话框 按下Run按钮 关闭对话框 几秒钟之后 在ProjectExplorer中出现一个新的文件 它包含所有的运行结果 一个工作日志窗口也会出现 它包含工作的运行状 态 我们也可以从JobExplorer中得到工作运行状况的信息 从菜单栏中选择View Exp