中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获.ppt

计算材料学主要内容：计算材料学主要内容： 密度泛函理论基础密度泛函理论基础 计算模型的建立计算模型的建立 晶体结构优化晶体结构优化 表面吸附表面吸附 电荷密度分布电荷密度分布 化学反应中的过渡态搜索化学反应中的过渡态搜索 气体扩散气体扩散 X-ray X-ray 结构精修结构精修 教员：教员： 郝绿原郝绿原 hly@ 3600834hly@ 3600834 傅正平傅正平 fuzp@ fuzp@ 李震宇（理论）李震宇（理论） zyli@ zyli@ 张文华（理论）张文华（理论） whzhang@ whzhang@资料：资料：* 不支持中文目录不支持中文目录* 显示扩展名显示扩展名限淆瑰缀伊母酉棱淮伦覆云袭麦调妄蹄渝懦情继舜凤倡弧姐扔辟团藩邯耳中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获CASTEPCASTEP概述概述纫虐轻丢柔虹勉耕猜乾洽颈始庄拔拇肩青击呢曳厨杭矗范宾膜口堵媒她抄中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获关于关于CASTAPCASTAP CASTAPCASTAP是特别为固体材料学而设计的一个现代的量子力学基是特别为固体材料学而设计的一个现代的量子力学基本程序，其使用了密度泛函本程序，其使用了密度泛函(DFT(DFT）平面波赝势方法，进行第一性）平面波赝势方法，进行第一性原理量子力学计算，以探索如半导体，陶瓷，金属，矿物和沸石原理量子力学计算，以探索如半导体，陶瓷，金属，矿物和沸石等材料的晶体和表面性质。

等材料的晶体和表面性质 典型的应用包括表面化学，键结构，态密度和光学性质等研典型的应用包括表面化学，键结构，态密度和光学性质等研究，究， CASTAP CASTAP也可用于研究体系的电荷密度和波函数的也可用于研究体系的电荷密度和波函数的3D3D形式此外，此外， CASTAP CASTAP可用于有效研究点缺陷（空位，间隙和置换杂质）可用于有效研究点缺陷（空位，间隙和置换杂质）和扩展缺陷（如晶界和位错）的性质和扩展缺陷（如晶界和位错）的性质 Material Studio Material Studio使用组件对话框中的使用组件对话框中的CASTAPCASTAP选项允许准备，选项允许准备，启动，分析和监测启动，分析和监测CASTAPCASTAP服役工作服役工作计算：允许选择计算选项（如基集，交换关联势和收敛判据），计算：允许选择计算选项（如基集，交换关联势和收敛判据），作业控制和文档控制作业控制和文档控制分析：允许处理和演示分析：允许处理和演示CASTAPCASTAP计算结果这一工具提供加速整体计算结果这一工具提供加速整体直观化以及键结构图，态密度图形和光学性质图形。

直观化以及键结构图，态密度图形和光学性质图形乏伶涤虱选豪丈耘符祖霖毯军涅来域碟蔚涯虞器蚕韦随莲垫翱素颈警便梯中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获CASTAPCASTAP的任务的任务 CASTAP CASTAP计算是要进行的三个任务中的一个，即单个点的能量计算，计算是要进行的三个任务中的一个，即单个点的能量计算，几何优化或分子动力学可提供这些计算中的每一个以便产生特定的几何优化或分子动力学可提供这些计算中的每一个以便产生特定的物理性能性质为一种附加的任务，允许重新开始已完成的计算以便物理性能性质为一种附加的任务，允许重新开始已完成的计算以便产生最初没有提出的额外性能产生最初没有提出的额外性能在在CASTAPCASTAP计算中有很多运行步骤，可分为如下几组：计算中有很多运行步骤，可分为如下几组：* * 结构定义：必须规定包含所感兴趣结构的周期性的结构定义：必须规定包含所感兴趣结构的周期性的3D3D模型文件，有模型文件，有大量方法规定一种结构：可使用构建晶体（大量方法规定一种结构：可使用构建晶体（Build Crystal)Build Crystal)或构建真或构建真空板空板(Build Vacuum Stab)(Build Vacuum Stab)来构建，也可从已经存在的的结构文档中引来构建，也可从已经存在的的结构文档中引入，还可修正已存在的结构。

入，还可修正已存在的结构注意：注意： CASTAP CASTAP仅能在仅能在3D3D周期模型文件基础上进行计算，必须构建超单周期模型文件基础上进行计算，必须构建超单胞，以便研究分子体系胞，以便研究分子体系提示：提示： CASTAP CASTAP计算所需时间随原子数平方的增加而增加因此，建议计算所需时间随原子数平方的增加而增加因此，建议用最小的原胞来描述体系，可使用用最小的原胞来描述体系，可使用Build\Symmetry\Primitive CellBuild\Symmetry\Primitive Cell菜菜单选项来转换成原胞单选项来转换成原胞雌伏梳烹怖秆拇目昏镁帛翌喊亦版闰干棵雁茶评蔬床伯胜液钙弃琶萄羊侠中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 CASTAPCASTAP中选择一项任务中选择一项任务1 1 从模块面板（从模块面板（Module Explorer)Module Explorer)选择选择CASTAP\CalculationCASTAP\Calculation。

2 2 选择设置表选择设置表3 3 从任务列表中选择所要求的任务从任务列表中选择所要求的任务 * * 计算设置：合适的计算设置：合适的3D3D模型文件一旦确定，必须选择计算类型模型文件一旦确定，必须选择计算类型和相关参数，例如，对于动力学计算必须确定系综和参数，包和相关参数，例如，对于动力学计算必须确定系综和参数，包括温度，时间步长和步数选择运行计算的磁盘并开始括温度，时间步长和步数选择运行计算的磁盘并开始CASTAPCASTAP作业 * 结果分析：计算完成后，相关于结果分析：计算完成后，相关于CASTAPCASTAP作业的文档返回用作业的文档返回用户，在项目面板适当位置显示这些文档的一些进一步处理要户，在项目面板适当位置显示这些文档的一些进一步处理要求获得可观察量如光学性质求获得可观察量如光学性质塘缉挟米随玩么抬充包伊身沂盐滓浙爬臻狱窖萍仆暑驴瘦阔翼履讼鼓件钱中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获CASTAPCASTAP能量任务能量任务CASTAPCASTAP能量任务允许计算特定体系的总能量以及物理性质。

能量任务允许计算特定体系的总能量以及物理性质除了总能量之外，在计算之后还可报告作用于原子上的力；也除了总能量之外，在计算之后还可报告作用于原子上的力；也能创建电荷密度文件；利用材料观测仪（能创建电荷密度文件；利用材料观测仪（Material Material Visualizer)Visualizer)允许目测电荷密度的立体分布；还能报告计算中使允许目测电荷密度的立体分布；还能报告计算中使用的用的Monkhorst-Park的的k点的电子能量，因此在点的电子能量，因此在CASTAPCASTAP分析中分析中可生成态密度图可生成态密度图对于能够得到可靠结构信息的体系的电子性质的研究，能量任对于能够得到可靠结构信息的体系的电子性质的研究，能量任务是有用的只要给定应力性质，也可用于计算没有内部自由务是有用的只要给定应力性质，也可用于计算没有内部自由度的高对称性体系的状态方程（即压力度的高对称性体系的状态方程（即压力- -体积，能量体积，能量- -体积关系）体积关系）注意：具有内部自由度的体系中，利用几何优化（注意：具有内部自由度的体系中，利用几何优化（Geometry Geometry Optimization)Optimization)任务可获得状态方程。

任务可获得状态方程CASTAPCASTAP中能量的默认单位是电子伏特中能量的默认单位是电子伏特(eV)(eV)，各种能量单位的换，各种能量单位的换算关系见算关系见Mohr.P.J(2000).Mohr.P.J(2000).1 eV=0.036749308 Ha=23.0605 kcal/mole=96.4853 kJ/mole肖钒壶番肯雷伪钉孜新淀彦抑衬旺钢躺售杉插良兔荔妊匠漓范座诽烩秧可中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获CASTAPCASTAP几何优化任务几何优化任务 CASTAPCASTAP几何优化任务允许改善结构的几何，获得稳定结构几何优化任务允许改善结构的几何，获得稳定结构或多晶型物通过一个迭代过程来完成这项任务，迭代过程中调或多晶型物通过一个迭代过程来完成这项任务，迭代过程中调整原子坐标和晶胞参数使结构的总能量最小化整原子坐标和晶胞参数使结构的总能量最小化 CASTAP CASTAP几何优化是基于减小计算力和应力的数量级，直到小几何优化是基于减小计算力和应力的数量级，直到小于规定的收敛误差。

也可能给定外部应力张量来对拉应力、压应于规定的收敛误差也可能给定外部应力张量来对拉应力、压应力和切应力等作用下的体系行为模型化在这些情况下力和切应力等作用下的体系行为模型化在这些情况下反反复迭复迭代内部应力张量直到代内部应力张量直到与所施加的外部应力相等与所施加的外部应力相等 几何优化处理产生的几何优化处理产生的模型结构与真实结构紧密模型结构与真实结构紧密相似利用相似利用CASTAPCASTAP计算的计算的晶格参数精度列于右图晶格参数精度列于右图火永囊泪漏责惩箭种庇魄客拎涟奥榆越持肮绊箱僵蒋凿到突勿域缕燥撤拴中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获状态方程计算状态方程计算在所施加静压力下几何优化可用于确定材料的体模量在所施加静压力下几何优化可用于确定材料的体模量B B和对压力的和对压力的导数导数B B‘= =dB/dPdB/dP过程包括计算理论状态方程（过程包括计算理论状态方程（EOSEOS），该方程描），该方程描述单胞体积与外部静压力的关系工艺非常类似于真实实验：使述单胞体积与外部静压力的关系。

工艺非常类似于真实实验：使用几何优化对话框中的应力列表将外部压力固定通过进行几何用几何优化对话框中的应力列表将外部压力固定通过进行几何优化可以找到在此压力下的单胞体积随后的优化可以找到在此压力下的单胞体积随后的P-V P-V 数据分析与实数据分析与实验研究精确一致描述验研究精确一致描述EOSEOS选择分析表达式，其参数适于计算数据选择分析表达式，其参数适于计算数据点最流行的点最流行的EOSEOS形式是三阶形式是三阶Birch-Murnaghan Birch-Murnaghan 方程：方程： 式中式中V V0 0 为平衡体积为平衡体积Cohen Cohen 等进行了等进行了EOSEOS各种解析式的的各种解析式的的详细比较研究详细比较研究注意：从相应实验中获得的注意：从相应实验中获得的B B和和B B‘值依赖于计算使用的压力值范值依赖于计算使用的压力值范围利用金刚石压砧获得的实验值通常在围利用金刚石压砧获得的实验值通常在0-30GPa0-30GPa范围内，因此范围内，因此推荐理论研究也在这个范围内在研究中避免使用负压力值也很推荐理论研究也在这个范围内在研究中避免使用负压力值也很重要。

此外，用于生成重要此外，用于生成P-V P-V 数据序列的压力值可能是不均匀的，数据序列的压力值可能是不均匀的，在低压力范围要求更精确采样以便获得体模量精确值在低压力范围要求更精确采样以便获得体模量精确值P-V 泛针胯走贫惋布祷毋扁融姚赖促景陶惫应斧侥激逃触邪痕才更仑荔景誓狡中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获几何优化方法几何优化方法 在默认条件下，在默认条件下，CASTAPCASTAP使用使用BFGSBFGS几何优化方法该方法通几何优化方法该方法通常提供了寻找最低能量结构的最快途径，这是支持常提供了寻找最低能量结构的最快途径，这是支持CASTAPCASTAP单胞单胞优化的唯一模式优化的唯一模式 衰减分子动力学（衰减分子动力学（ Damped molecular dynamics）方法是）方法是另一种可以选择的方法，该方法对具有平滑势能表面的体系如另一种可以选择的方法，该方法对具有平滑势能表面的体系如分子晶体或表面分子分子晶体或表面分子与与BFGSBFGS同样有效。

同样有效 谋响咯策据孔涌漳砾惧纯涧辖柿捅旺司傍湖习魄爸囚栓政武摆怀彦皑龄猿中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获CASTAPCASTAP动力学任务动力学任务 CASTAP CASTAP动力学任务允许模拟结构中原子在计算力的影响下将如何移动动力学任务允许模拟结构中原子在计算力的影响下将如何移动 在进行在进行CASTAPCASTAP动力学计算以前，可以选择热力学系综和相应参数，定义模拟动力学计算以前，可以选择热力学系综和相应参数，定义模拟时间和模拟温度时间和模拟温度选择热力学系综选择热力学系综 对牛顿运动定律积分允许探索体系恒值能量表面（对牛顿运动定律积分允许探索体系恒值能量表面（NVENVE动力学）然而，在动力学）然而，在体系与环境进行热交换条件下发生最本质的现象使用体系与环境进行热交换条件下发生最本质的现象使用NVTNVT系综（或者是确定性系综（或者是确定性的的NosNosé系综或者是随机性的系综或者是随机性的Langevin Langevin 系综）可模拟该条件。

系综）可模拟该条件定义时间步长（定义时间步长（timestep timestep ）） 在积分算法中重要参数是时间步长为更好利用计算时间，应使用大的时间在积分算法中重要参数是时间步长为更好利用计算时间，应使用大的时间步长然而，如果时间步长过大，则可导致积分过程的不稳定和不精确典型地，步长然而，如果时间步长过大，则可导致积分过程的不稳定和不精确典型地，这表示为运动常数的系统偏差这表示为运动常数的系统偏差注意：量子力学分子动力学计算要求比力场动力学使用更小的时间步长注意：量子力学分子动力学计算要求比力场动力学使用更小的时间步长动力学过程的约束动力学过程的约束CASTAPCASTAP支持支持Langevin NVTLangevin NVT或或NVENVE动力学过程的线性约束然而，借助动力学过程的线性约束然而，借助Material Material StudioStudio界面可以近似使用以下两种更基本的约束：界面可以近似使用以下两种更基本的约束： 质心固定，单个原子固定质心固定，单个原子固定 使用使用seednameseedname.cell.cell 文档可以利用更复杂的约束。

文档可以利用更复杂的约束 陨迂饥辞冻垛婚卖戎犁愧炸邦骋邹田妙拾蔼削伊祁临驭涵买急霄撵仆再撂中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获CASTAPCASTAP性质任务性质任务 CASTAP CASTAP性质任务允许在完成能量，几何优化或动力学运行之后求出电性质任务允许在完成能量，几何优化或动力学运行之后求出电子和结构性质可以产生的性质如下：子和结构性质可以产生的性质如下：* * 态密度（态密度（DOSDOS）：利用原始模拟中产生的电荷密度和势能，非自恰计算价）：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的精细带和导带的精细Monkhorst-Pack Monkhorst-Pack 网格上的电子本征值网格上的电子本征值 * 带结构：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带结构：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的布里渊区高对称性方向电子本征值带的布里渊区高对称性方向电子本征值 * 光学性质：计算电子能带间转变的矩阵元素。

光学性质：计算电子能带间转变的矩阵元素CASTAPCASTAP分析对话可用于生分析对话可用于生成包含可以测得的光学性质的网格和图形文件成包含可以测得的光学性质的网格和图形文件 * 布局数分析：进行布局数分析：进行Mulliken Mulliken 分析计算决定原子电荷的键总数和角动量分析计算决定原子电荷的键总数和角动量（以及自旋极化计算所需的磁矩）任旋地，可产生态密度微分计算所要（以及自旋极化计算所需的磁矩）任旋地，可产生态密度微分计算所要求的分量求的分量 * 应力：计算应力张量，并写入应力：计算应力张量，并写入seedname.castepseedname.castep 文档吕柑唇友舷舀陋汕预吼柜木茬食沤坠冤惊渴怔国稠尖缕栋企截观絮滥舍姿中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 如果要进行单胞参数固定时进行几何优化运行和要检查点如果要进行单胞参数固定时进行几何优化运行和要检查点阵偏离平衡的程度，这些信息是有用的例如，可进行符合于阵偏离平衡的程度，这些信息是有用的。

例如，可进行符合于给定体系理论基态的固定单胞的点缺陷的超晶胞研究几何优给定体系理论基态的固定单胞的点缺陷的超晶胞研究几何优化后的应力值显示了与超单胞近似相关联的弹性效应化后的应力值显示了与超单胞近似相关联的弹性效应注意：为计算某种性质，从适当模拟得到的结果文档必须以当注意：为计算某种性质，从适当模拟得到的结果文档必须以当前的文件夹形式出现前的文件夹形式出现秘堵驮拔淳群毋奥缓蹦同站益幌闭荚仲笺农哨胆希做血应趁少湍蜀芒阂略中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获目的目的: : 介绍介绍CASTEPCASTEP中的结构优化，中的结构优化， 使用立体可视化工具显示等值面使用立体可视化工具显示等值面模块模块: Materials Visualizer, CASTEP: Materials Visualizer, CASTEP前提前提: : 使用晶体建模工具使用晶体建模工具用第一原理预测用第一原理预测AlAsAlAs的晶格参数的晶格参数 内容内容1. 1. 构建构建AlAsAlAs的晶体结构的晶体结构2. 2. 设置并进行设置并进行CASTEPCASTEP计算计算3. 3. 分析结果分析结果4. 4. 比较计算的结构参数和实验数据比较计算的结构参数和实验数据 (1) (1)图示电荷密度图示电荷密度 (2) (2)图示态密度和带结构图示态密度和带结构引言引言 本指南介绍了本指南介绍了CASTEPCASTEP是如何使用量子力学方法来确定材料的晶体结构，使用者是如何使用量子力学方法来确定材料的晶体结构，使用者将学会如何构建晶体结构，设定一个将学会如何构建晶体结构，设定一个CASTEPCASTEP几何优化任务，然后分析计算结果。

几何优化任务，然后分析计算结果背景背景 密度泛函理论密度泛函理论 (DFT) (DFT)在周期性大尺度材料上应用的进展，对材料设计和加工越在周期性大尺度材料上应用的进展，对材料设计和加工越来越重要该理论使得研究者能对实验数据进行解释；并从未知晶体的结构性质、来越重要该理论使得研究者能对实验数据进行解释；并从未知晶体的结构性质、结合能和表面分子的活动性确定材料的本征性质这些理论工具可用于指导新材料结合能和表面分子的活动性确定材料的本征性质这些理论工具可用于指导新材料的设计，帮助研究者了解内在的化学和物理过程的设计，帮助研究者了解内在的化学和物理过程注意注意: : 如果你的服务器没有足够快的如果你的服务器没有足够快的CPUCPU，请慎用，请慎用CASTEPCASTEP进行几何优化计算，因为进行几何优化计算，因为它会占用相当长的时间它会占用相当长的时间 抓急拷谐仁灌谨减硬塞蛮悦淫侧擒桔确冰枢套刽未晶卞刨胆晕播跌际寒驱中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获AlAl的分数坐标：的分数坐标：(0 0 0）） (1/2 1/2 0) (1/2 0 1/2) （（0 1/2 1/2））AsAs的分数坐标：的分数坐标：(3/4 3/4 3/4）） (1/4 1/4 3/4) (1/4 3/4 1/4) （（3/4 1/4 1/4））(0 0 0））(1/2 0 1/2)(1/2 1/2 0)（（0 1/2 1/2））(1/4 1/4 3/4)（（3/4 1/4 1/4））（（1/4 3/4 1/4））（（3/4 3/4 3/4））As: (3/4 3/4 3/4））= (1/4 1/4 1/4）） 1. 1. 构建构建AlAsAlAs的晶体结构的晶体结构 空间群是空间群是F-43m弛豹幽勾拢妥枫阻钥怠暮汰洱绣净僵婚磋谤旬诀润蜗沂招坯渍骂增挣茄驯中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 构建一个晶体结构，需要知道该晶体的空间群、晶格参数和晶体的内坐标。

构建一个晶体结构，需要知道该晶体的空间群、晶格参数和晶体的内坐标对对AlAs AlAs 来说，空间群是来说，空间群是F-43mF-43m，空间群代号为，空间群代号为216216原胞有两个原子，原胞有两个原子，Al Al 和和As As 的分数坐标分别为的分数坐标分别为(0, 0, 0)(0, 0, 0)和和(0.25, 0.25, 0.25)(0.25, 0.25, 0.25)，晶格参数为，晶格参数为5.6622 5.6622 Å Å. . 第一步是建立晶格第一步是建立晶格 在在D diskD disk上建立英文目录上建立英文目录D:\class3D:\class3按下面步骤按下面步骤, ,在在Project Explorer Project Explorer 内，建立内，建立AlAsAlAs根目录从“开始”或快捷图标 打开MS梯蝗欠消摸勺深窟得向纹斥蓟血佐骑乘番干砚讨痢蒸珐藕未望澜咆唇苦衷中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获找到class3，按“打开”按钮输入AlAs，这将是新的Project的名字。

另贸晒骡抹恫蜘吻磋其堵船绑期荚激囚黔魁梁楔溯蟹教雁奋锤咖堂惰损奎中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获在 Project Explorer中，右击根目录AlAs，选择New | 3D Atomistic Document右击3D Atomistic document文件，将该文件重新命名为AlAs.xsd潍千嫩艾冀冤仕缄奴馅踞章男溃朽屈隋逾使盏栓柬郊栏才缸茶醛范啤服架中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获从菜单栏里选择从菜单栏里选择Build / Crystals / Build CrystalBuild / Crystals / Build CrystalBuild Crystal Build Crystal 对话框显示出来对话框显示出来点击点击Enter group Enter group 输入输入216216，按下，按下TAB TAB 按钮按钮( (或或在在Enter groupEnter group中选择中选择F-43m)F-43m)，，空间群信息更新为空间群信息更新为F-43m F-43m 空间群。

空间群空间群信息框中的信息也空间群信息框中的信息也随着随着F-43mF-43m空间群的信息而发生变化空间群的信息而发生变化 杨碚芳课杨碚芳课眺岿勤祷墩蘸粳按窗秩漫曾搓纳所搬吾橙酝适胡佬蚜艇姑预滨蜂疤央粹堑中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获选择选择Lattice Parameters 标签，标签，把把a值从值从10.00 变为变为5.662点击击Build 按钮欧尼弊架靡缎破宠菜防饼黍拿朗审秸敬仲绚陛埠判袜第晶娃母吩戚忿否鬃中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获单击此图标，然后可旋转晶格，显示其立体结构一个没有原子的3D 格子显示在格子显示在3D Atomistic 文件里AsAl?坟得廖闸勒鸿炉殿鸽仓滇搏吧辈遗鲁击集抨肤典郭灼城穆聊先眯扭健骆叮中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获Atom # OX SITE x y z SOF H Al 1 +3 4 a 0 0 0 1. 0 As 2 -3 4 c 0.25 0.25 0.25 1. 0 *end for ICSD #67784AlAl的分数坐标：的分数坐标：(0 0 0）） (1/2 1/2 0) (1/2 0 1/2) （（0 1/2 1/2））AsAs的分数坐标：的分数坐标：(3/4 3/4 3/4）） (1/4 1/4 3/4) (1/4 3/4 1/4) （（3/4 1/4 1/4））输入几个原子？输入几个原子？Ba3Si6O9N4作业作业1 1：解释符：解释符号和群号和群F-43mF-43m的的意思，给出参意思，给出参考书。

考书需输入需输入3 3个个BaBa 2 2个个SiSi 2 2个个N N 3 3个个O O厂靡腰皇炔寥连第桔市炸礼另认帘疗频痒罩朴猖乌爪睬枣曙厕尝受长眷笆中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 * * 从菜单栏中选择从菜单栏中选择Build / Add AtomsBuild / Add Atoms通过Add Atoms Add Atoms 对话框，我们可以对话框，我们可以把原子添加到指定的位置，其对话框如下：把原子添加到指定的位置，其对话框如下： 在Add Atoms对话框中选择Options标签，确定Coordinate system为Fractional停业捍髓爆满检遗偷氢觅仅碉航糜饭捉饮佛碘馋嚣番笋毋划附债贴呈帜谎中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* * 如上所示，选择如上所示，选择AtomsAtoms标签，通过周期表，在标签，通过周期表，在ElementElement文本框文本框中输入中输入AlAl，， 再输入再输入Al 的分数坐标的分数坐标 (0, 0, 0)，，然后按下然后按下AddAdd按钮，按钮，铝原子就添加到结构中了。

铝原子就添加到结构中了粘筑挣扶初韶心挞灯砰攻扁夸焊吾歪珊啊垛镇哪簿事娇恭荚研呼倒却渗蹬中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 * * 从菜单栏中选择从菜单栏中选择View / Display StyleView / Display Style在打开的对话框中，在打开的对话框中，选择选择Ball and stickBall and stick显示方式显示方式碟衷俱咽锣候租保盔贬砍熏矮姻蹭肛惯躲筷盅柜嘲浩郧纲袖蕴苞答殴银剧中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* * 同前，在同前，在ElementElement文本框中键入文本框中键入AsAs在a, b, ca, b, c文本框中键入文本框中键入0.250.25按AddAdd按钮，这样按钮，这样AsAs也加入到晶格中关闭对话框也加入到晶格中关闭对话框说明：上面操作虽然只加入一个说明：上面操作虽然只加入一个AlAl、一个、一个AsAs，但群的对称操作在晶体中补，但群的对称操作在晶体中补充了剩余的充了剩余的AlAl原子和原子和AsAs原子原子( (等位原子等位原子) )。

如果新加原子以如果新加原子以lineline方式显示，可重复前面步骤，改显示方式为方式显示，可重复前面步骤，改显示方式为Ball Ball and Stickand Stick锣删评哨中乱脆亿甘送韵碎蚀窃光始每套呜阳执推巴咽左鼻刽饮醉陕通竞中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获单击此图标，出现添加原子Add Atoms 对话框选择原子名称，输入分数坐标，按Add，则原子添加到晶体结构图中重复操作，直到添加完晶胞中的所有原子关闭Add Atoms框前面的添加原子操作也可用下面图标来实现这里前面的添加原子操作也可用下面图标来实现这里不再重复不再重复寻彪旨眼入拟突纬警瓶晒操络酋瘴掉淹起隶箩成兴堪搞乘渍磨勺瓷即所爽中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* * 从上面的从上面的AlAsAlAs晶体结构看出，近邻晶胞中的原子也显示出来。

这种晶体结构看出，近邻晶胞中的原子也显示出来这种显示表示了显示表示了AlAsAlAs晶体中键的拓扑结构当然，可以通过重新建造晶体晶体中键的拓扑结构当然，可以通过重新建造晶体结构来移去这些结构来移去这些近邻晶胞中的近邻晶胞中的原子 从菜单栏中选择从菜单栏中选择Build / Crystals / Rebuild Crystal...，打开对话框，，打开对话框，按下按下Rebuild按钮在显示出的晶体按钮在显示出的晶体结构中那些原子就被移走了已经结构中那些原子就被移走了已经把显示方式定为把显示方式定为Ball and Stick邵液跳慰梭沟临妇朋形恼眷腰廷秃篆歪需鉴装剔攒星敌涪薯荡权害羚锋惹中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获按下面图示步骤，保存结果按下面图示步骤，保存结果萤描刑摆榆勿舒兹讯箱缺坝蹭母姆作蚕弯漳绥哎孕灵列昌喳蒸袖厨栓昼崔中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 3D Viewer 3D Viewer 内的晶体结构是传统的内的晶体结构是传统的晶胞晶胞[conventional [conventional (unit) cell](unit) cell]，显示了晶格的立方对称性。

显示了晶格的立方对称性CASTEP CASTEP 利用利用了晶格的对称性，可以使用只包含两个原子的了晶格的对称性，可以使用只包含两个原子的原胞原胞[primitive (unit) cell][primitive (unit) cell]来进行计算，而晶胞包含了来进行计算，而晶胞包含了8 8 个原子电荷密度、键长和每个原子的总能量将是一样的，个原子电荷密度、键长和每个原子的总能量将是一样的，而不管这个而不管这个unit cellunit cell是如何被定义的这样，使用原胞，是如何被定义的这样，使用原胞，原子数较少，计算量大大减小，计算时间将被缩短原子数较少，计算量大大减小，计算时间将被缩短Note：在计算磁性体系时，使用了自旋极化这时要注：在计算磁性体系时，使用了自旋极化这时要注意，电荷密度自旋波的周期是原胞的数倍意，电荷密度自旋波的周期是原胞的数倍熟席昂奶钮陶盒唇箕费舰卑刷绳粟瞎展眼趾碾檬留犯伦心销曝杰宪慢裳削中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 选择菜单栏里的选择菜单栏里的Build / Symmetry / Primitive Cell。

模型文件模型文件(3D Viewer)显示为原胞显示为原胞(primitive cell) AlAs的原胞 不同角度 菏三氏舰门濒钱质纽督酗搪沾牙醒尔市社晌岳标掳派陛剩磷轿成透彩省庚中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获在晶体图上按右键，选在晶体图上按右键，选Label，在出现的对话框中选，在出现的对话框中选ElementSymbol按按apply，晶胞上显示元，晶胞上显示元素符号龚供牡惜叼氦坑行冻赎沧缺来芭镣衰仆仗氟凳姑列功辈猎核缴暂娄私钙赊中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获2. 设置设置CASTEP CASTEP 计算任务计算任务 从工具栏中选择从工具栏中选择CASTEP 工具，再选择工具，再选择Calculation或从或从菜单栏中选择菜单栏中选择Modules | CASTEP |Calculation。

CASTEP Calculation对话框如下：对话框如下： 下面我们分两步，先优化下面我们分两步，先优化AlAsAlAs的几何结构，再计算的几何结构，再计算Band Band StructureStructure和和Density of statesDensity of states洞脉攘嫉捕陪但樟额拄付坍萍篆侯诡完舰钱签攻揩哟穆纯睹羚蛮镁骄桨辗中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获(1) 优化优化AlAs晶体结构晶体结构* 把把Task 改为改为Geometry Optimization，计算精度，计算精度Quality设设置为置为coarse 结构优化的默认设置是优化原结构优化的默认设置是优化原子坐标在本例中，我们不仅要子坐标在本例中，我们不仅要优化原子坐标，同时也要优化晶优化原子坐标，同时也要优化晶格常数 按下与按下与Task 相关的相关的More…按钮，按钮，勾选上勾选上Optimize Cell，， 关闭此对关闭此对话框 当改变计算精度的时候，其它当改变计算精度的时候，其它的参数也会自动作相应的变化。

的参数也会自动作相应的变化揍致硒微跪莎掩毫陷冰鼠魁狡呕虾通串弹汞遣钱底唇大病陪怠篇烬澳客更中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* * 选择选择ElectronicElectronic标签栏，按下标签栏，按下MoreMore……按钮在SCFSCF对话框里作如下设置，对话框里作如下设置，将将ChargeCharge由由0.50.5改为改为0.150.15，，钩上钩上Fix occupancyFix occupancy * 选择选择PropertiesProperties标签栏，里面的计标签栏，里面的计算任务都算任务都不不要选嫡扁毛窜拆疽挽遥惠闺亭摘疲懈荡溪辕琳呵砾遇杉丧居逮崇荷供腮此获抿中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* * 选择选择Job Control Job Control 标签栏，按下标签栏，按下MoreMore……按钮。

在按钮在CASTEP CASTEP Job Control Options Job Control Options 对话框里，改变对话框里，改变 Update interval Update interval 为为30.0 s30.0 s，关闭此对话框关闭此对话框•按下按下Run 按钮，关闭对话框按钮，关闭对话框•注意，此时注意，此时3D结构为激活窗口如果激活窗口是文本，结构为激活窗口如果激活窗口是文本， 则则Run为灰色埠蓉窘昔够灶擂乱屎淡丹蚁月站机嗓篱姑灵讽共柠湃他仑希堂停阔塔臼燃中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 几秒钟后，一个新文件夹出现在几秒钟后，一个新文件夹出现在Project Explorer Project Explorer 内，该文件夹包含内，该文件夹包含了所有的计算结果了所有的计算结果 如果使用客户端－服务器模式，当工作结束时，文件会被如果使用客户端－服务器模式，当工作结束时，文件会被传回到客户端。

数据传输过程需要一定的时间，与文件的大传回到客户端数据传输过程需要一定的时间，与文件的大小有关Job Explorer Job Explorer 显示了所有正在运行的工作的状态它显示了很多有用的信息，包显示了所有正在运行的工作的状态它显示了很多有用的信息，包括服务器和工作代码如果需要，也可以通过括服务器和工作代码如果需要，也可以通过Job ExplorerJob Explorer来中止运行工作来中止运行工作沈儡廓臭晨御鹏窒数砚钝激贝挚巢岭堑歼唉拍孽张码靡恫侥弱担颇却甸婿中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获在工作运行过程在工作运行过程中，四个文件打中，四个文件打开了这些文件开了这些文件包含了晶体结构、包含了晶体结构、结构优化过程中结构优化过程中模型的更新、工模型的更新、工作参数的设置和作参数的设置和运行状态的信息，运行状态的信息，以及关于总能量、以及关于总能量、能量变化、应力、能量变化、应力、压力和位移随迭压力和位移随迭代次数变化的图代次数变化的图表计算过程中出现的两个表计算过程中出现的两个表示能量收敛的图框。

示能量收敛的图框查看计算查看计算设置设置营揪平泌镣饭索品倘鄂咬丰臣篡晾瑰常拿敛竹虽窄腆杰忻娩撒察脏购夫丢中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获固哆莱弦戮翻郸撞滞楔研挡帚妈蔗暑骑臃锑港泣夯眺蜀摔间波烘近算亭顷中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 几分钟后，计算结束，出现几分钟后，计算结束，出现‘‘Job CompletedJob Completed’’ 提示，表示计算成功提示，表示计算成功输出文本输出文本文档为文档为AlAs.castep，包含优化信息，在，包含优化信息，在AlAs CASTEP GeomOpt文件夹中文件夹中按下面图按下面图示操作，关闭工作窗口中的文件示操作，关闭工作窗口中的文件 (2) (2) 接下来，利用优化过的接下来，利用优化过的AlAsAlAs结构，计算结构，计算AlAsAlAs的的Band StructureBand Structure和和Density of statesDensity of states。

* 在结构优化的文件中，双击在结构优化的文件中，双击AlAs.xsdAlAs.xsd、、AlAs.castepAlAs.castep两个文件两个文件这两个文件出现在工作窗口中按这两个文件出现在工作窗口中按下面图示设定计算任务，直至下面图示设定计算任务，直至RunRun拭平舵没翰尿吼氓必馋绝瓮钱枣孜佑客壤久塌扭麦瞧梦密石诌厘醛蹋焰绘中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获出现新的文出现新的文件夹件夹键麻巩诉垣纠卖颖是郝彦棋蝉凸赊挫弯期珍毫闸柜姬轰吩旧鸣宰犁炭郎案中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 打开新文件夹打开新文件夹AlAs CASTEP Properties，双击，双击AlAs_BandStr.castep，这，这此文件出现在工作窗口中此文件出现在工作窗口中 按图示操作，显示按图示操作，显示AlAs的的Band Structure。

3) (3) 计算结束后，查看计算结束后，查看AlAsAlAs的的Band StructureBand Structure和和Density of statesDensity of states的计算结果的计算结果可用工具可用工具 放大、放大、缩小能带图缩小能带图罢晴技葫饿李勉馈散搏秽检玄副傍招坐溢滦栗蝗回脓炽解饲决掷闺湛敞盏中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 能带图下方的字母表示布里渊区的高对称点单击能带图下方的字母表示布里渊区的高对称点单击AlAs.xsd文件，使其为当前文件，使其为当前活动窗口再按图示点击计算设置活动窗口再按图示点击计算设置 ，使当前状态为，使当前状态为Band Structure点击下方的点击下方的More，则显示出对称点的坐标则显示出对称点的坐标12遵研疫女缺序垮虹酉揖貉订耙泥屁械莽熔关疫队喜唱昭挪霖港痔五美恶朴中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 在新文件夹在新文件夹AlAs CASTEP Properties中双击中双击AlAs_DOS.castep，这此文，这此文件出现在工作窗口中。

件出现在工作窗口中 按图示操作，显示按图示操作，显示AlAs的的DOS可用工具可用工具 放大、放大、缩小能带图缩小能带图举纠摧慕刹典忘蚌渠蹬傻最串已匆木止形唾荔爽栅酶完讨桶滩赢认裹淆番中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获(4) (4) 前面的结构优化显示，前面的结构优化显示，AlAsAlAs的总能量随迭代次数振荡下面改计算方的总能量随迭代次数振荡下面改计算方法进行结构优化，避免振荡计算步骤与前面类似，图示如下法进行结构优化，避免振荡计算步骤与前面类似，图示如下孙吞狄俘呛艘胡狂扩嫌狭汝距曳旧粹疹据暮沙跨鞘妹疫茨漆皋嚼耕程戍劝中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获结构优化过程，可结构优化过程，可见随着迭代次数的见随着迭代次数的增加，增加，AlAsAlAs的总能的总能量逐渐减小量逐渐减小丁鞘晴井仑叮伪弧谐叠胰没茫烙钒潘等何稳擦叭官木匹施侍雄荤木种巾旨中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 按下面图示步骤，保存计算结果，关闭工作窗口中的文件。

按下面图示步骤，保存计算结果，关闭工作窗口中的文件5) (5) 同样，利用优化过的同样，利用优化过的AlAsAlAs结构，计算结构，计算AlAsAlAs的的Band StructureBand Structure和和Density Density of statesof states在新出现的文件夹中，打开在新出现的文件夹中，打开AlAs.castep、、AlAs.xsd两个文件这两个文件两个文件这两个文件都出现在工作窗口中都出现在工作窗口中作业2：在.CASTEP文件中找到并记下总电子能量，与前面优化后的总电子能量比较，哪个应该小？测恿造古惦炙肤习旅药魔裳笛咨滦按六澎恋厄液舰娠期烹贷煮稀裁改虚哲中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 设置计算任务，按设置计算任务，按Run进行计算，关闭对话框进行计算，关闭对话框伍案军头羹狗蛙蛹置赃喀怪队艇族报姐贵呀褐豺望捎喝螟浪夜逛纳稠惯菏中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获(6) (6) 计算结束后，查看计算结束后，查看AlAsAlAs的的Band StructureBand Structure和和Density of statesDensity of states的计算结果。

的计算结果 计算结束后，出现提示，关闭计算结束后，出现提示，关闭 保存计算结果，清理工作窗保存计算结果，清理工作窗口 在新出现的在新出现的AlAs CASTEP Properties文件夹文件夹中，双击中，双击AlAs_BandStr.castep，该文件出现，该文件出现在工作窗口中在工作窗口中寨樱俊罗戳纫蓟咱豪铅举辅乱酒杆脚燃跋煞蔽描叫抒婪奶捧巢沉枕乞师赁中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 按图示操作，显示带结按图示操作，显示带结构 同样，双击同样，双击AlAs_DOS.castep，显示态，显示态密度按寓喝祝缉秦轰灵传隅允既膏祁睡讶坚檬吓砖滦渠认炒侩驮丰椭韵凶野井中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获说明：说明： 分析工具可以用来显示态密度分析工具可以用来显示态密度(DOS)(DOS)和能带结构和能带结构。

能带结构图显示了布里渊区内沿着高对称方向电子能量对能带结构图显示了布里渊区内沿着高对称方向电子能量对k k矢的依赖矢的依赖性这些图提供了一个对材料的电子结构进行定性分析的非常有用的工具这些图提供了一个对材料的电子结构进行定性分析的非常有用的工具譬如，与近自由的譬如，与近自由的s s 、、p p 电子构成的能带相比，很容易鉴别出电子构成的能带相比，很容易鉴别出d d、、f f电子构电子构成的窄带成的窄带 DOS DOS 和和PDOS PDOS 图给出了材料的电子结构的一个快速定性图像，有时候图给出了材料的电子结构的一个快速定性图像，有时候它们可以直接和实验光谱结果相关联它们可以直接和实验光谱结果相关联 CASTEP CASTEP的主要输出结果文件的主要输出结果文件AlAs.castep AlAs.castep 包含了有限的能带结构和包含了有限的能带结构和DOS DOS 信息，更多的详细信息包含在信息，更多的详细信息包含在AlAs_BandStr.castep AlAs_BandStr.castep 文件内 打开打开Analysis Analysis 对话框，选上对话框，选上Band structureBand structure。

从这个对话框可以看从这个对话框可以看出，可以把能带结构和态密度信息显示在同一个图中在出，可以把能带结构和态密度信息显示在同一个图中在DOSDOS部分，选上部分，选上Show DOSShow DOS，单击，单击ViewView，出现的图包含了带结构和，出现的图包含了带结构和DOSDOS两种信息当然，可两种信息当然，可以分别显示能带结构和态密度以分别显示能带结构和态密度 可以按图片、数据格式输出图文件，数据可由可以按图片、数据格式输出图文件，数据可由ExcelExcel等软件读取等软件读取还可以借助还可以借助CASTEP CASTEP 来计算很多其他性质，比如反射率和介电函数等等来计算很多其他性质，比如反射率和介电函数等等喊褐庚租惦黄医秀锹留阑管狰偶瞧菇妒褪寡啄小琅懒藤霜滋屯忘哲苫骡跺中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获DFT计算带隙计算带隙Eg，数值偏，数值偏小比较下图可知，计算精小比较下图可知，计算精度高，度高，Eg大7) (7) 比较两次计算的结果比较两次计算的结果12431.2961.762作业作业2 2：比较两次计算出的总电子能量：比较两次计算出的总电子能量劈谷悄猫兢财茵关给韧犁青嚏呢鱼宅厌焉臣敖撒粹览川昆惭翅茧恐填音高中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获3 3 分析结果分析结果当结果文件被传输回来，会得到下列数个文件：当结果文件被传输回来，会得到下列数个文件：• • AlAs.xsd AlAs.xsd –– 最后的优化结构最后的优化结构• • AlAs Trajectory.xtd - AlAs Trajectory.xtd - 一个轨迹文件，一个轨迹文件， 包含了包含了每一个优化步骤后每一个优化步骤后的结构的结构• • AlAs.castep AlAs.castep –– 包含了优化信息的输出包含了优化信息的输出 文本文件文本文件• • AlAs.param AlAs.param ––模拟所用的输入参数模拟所用的输入参数 计算任何一个性质，都会产生计算任何一个性质，都会产生.param .param 和和.castep .castep 文件。

文件 在在AlAs AlAs 结构中，由于对称性的存在，受力为结构中，由于对称性的存在，受力为0 0，但是应力的大小取，但是应力的大小取决于晶格参数这样，决于晶格参数这样，CASTEP CASTEP 就会努力去最小化系统的总能量和应力就会努力去最小化系统的总能量和应力因此，为保证计算能够合适地完成，检查压力收敛是非常重要的因此，为保证计算能够合适地完成，检查压力收敛是非常重要的锄撇半螟越已婶幻缮凤缮卑郡币浅铃烫人削操怒汝鸯堆涣肥嗣稗云宣吸俱中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 在在Project Explorer 内，双击内，双击AlAs.castep ，将其激活为当前工作文件将其激活为当前工作文件选择菜单栏里的选择菜单栏里的Edit | Find... ，在文本框中输入，在文本框中输入“completed successfully”，按下，按下Find Next 按钮，按钮，AlAs.castep文件滚动文件滚动 看到一个含有两行的表格，最后一列的每一行都显示为看到一个含有两行的表格，最后一列的每一行都显示为Yes，这表明计，这表明计算成功地结束。

算成功地结束椰邪伙卿矩谎懈阶乙狗啊崩龙排椽炎硝这座固翁唱承堰捕扯瞳笋惩涂淮雀中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 开始创建晶胞时，就知道晶格长度为开始创建晶胞时，就知道晶格长度为5.6622 Å因此，可以把能量最小化因此，可以把能量最小化后的晶格长度与初始的实验数据相比较实验的晶格长度基于晶胞结构，而后的晶格长度与初始的实验数据相比较实验的晶格长度基于晶胞结构，而不是原胞，因此需要将现在的原胞转化为晶胞，再与实验数据比较不是原胞，因此需要将现在的原胞转化为晶胞，再与实验数据比较 4. 比较比较AlAs晶体结构的计算数据与实验数据晶体结构的计算数据与实验数据* 双击双击AlAs.xsd 使其为当前工作文件使其为当前工作文件* 从菜单栏里选择从菜单栏里选择Build / Symmetry / Conventional Cell，晶胞显示出来晶胞显示出来豌枉架求婚谜曾壶打剪佰堰抚纫乎对绎久析缴傻整森寨鸽甘鱼宫咬谈哀扰中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 有数种方法看到晶格有数种方法看到晶格长度，一种就是打开长度，一种就是打开Lattice Parameters 对对话框。

在模型文件话框在模型文件(3D Viewer)上右击，选择上右击，选择Lattice Parameters格矢大约为矢大约为5.721128 Å 另一种简单的方法是在左另一种简单的方法是在左侧侧Properties中选择中选择Lattice 3D，其中显示晶格常数为，其中显示晶格常数为5.72113 Å 误差大约是误差大约是-0.5%这在1-2%典型误差范围内，这个误差值是赝势平典型误差范围内，这个误差值是赝势平面波方法与实验结果比较的期望误差使用面波方法与实验结果比较的期望误差使用GGA函数，晶格常数的计算值函数，晶格常数的计算值偏大；使用偏大；使用LDA函数，晶格常数的计算值偏小；函数，晶格常数的计算值偏小； 继续之前，需要保存工作，并关闭所有窗口继续之前，需要保存工作，并关闭所有窗口选择菜单栏上的选择菜单栏上的File | Save Project，然后是，然后是Window | Close All萨台启洞孽弘梭抱邓嘉虽十确免珐乐扔姐榴强伞札率笼鞭粗维挑赫驴刻吧中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 从工具栏选择从工具栏选择CASTEP ，，然后选择然后选择Analysis(或从菜单栏选择或从菜单栏选择Modules / CASTEP / Analysis)，再，再选中选中Electron density选项。

选项 此时，会有一条信息此时，会有一条信息““no results file is availableno results file is available””，说没有什么，说没有什么结果文件，所以需要指定结果文件结果文件，所以需要指定结果文件 在在Project Explorer Project Explorer 内，双击内，双击AlAs.castepAlAs.castep这将把结果文件和分析对这将把结果文件和分析对话框关联起来话框关联起来5.可视化电荷密度可视化电荷密度 可以用可以用CASTEP 分析工具得到电荷密度分析工具得到电荷密度冕动趴珊引裔嗅睡醒视诛邢阉硫抗瞻滋瞄换铃梆涕孩持泳捏型螺亦缉议淤中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 需要指定一个需要指定一个3D 文件来显示文件来显示等密度面等密度面在在Project Explorer 内，双击优内，双击优化后的化后的AlAs.xsd 文件 选择菜单栏里的选择菜单栏里的Build / Symmetry /Primitive Cell贬司扶佑缠努丹鬼倔拽挺狙别淤致腋知妥缄葵笼恤绽念料威雾引讲毫昂柒中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* CASTEP Analysis 对对话框上的话框上的Import 按钮按钮现在是激活状态。

现在是激活状态按下按下Import 按钮AlAs的电子密度等值面的电子密度等值面作业3：电子密度有无负值？电子如何分布？不叮驾米膜设雇柒涧短神庐烘巩模瞻阔羞肺碾隘面珐万戎谤传泊小现柯餐中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获 可以通过可以通过Display Style Display Style 对话框来改变等密度面的设置对话框来改变等密度面的设置 * *在等值面的在等值面的3D 3D 文件上右击，选择文件上右击，选择Display StyleDisplay Style，选择，选择Isosurface Isosurface 标标签栏 Isosurface Isosurface 标签栏显示出来：标签栏显示出来：惩经术酿泼挛惧职澳镍啊苯辟冶贴枯尝输挖镐舶跨守穴糙脾循篓力肠崩帝中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 这里可以改变不这里可以改变不同的设置。

同的设置在在Isovalue 文本文本框里，输入框里，输入0.1，，然后按然后按TAB 键注意等密度面是注意等密度面是如何改变的如何改变的猪柯撤疗鹅茵冉伤侮穿灸轩浦秽寞虽伴橡肺阻钒绚痕卢儡还弗关咨行昼英中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 把把Transparency 滑条滑条向右移动向右移动向右移动向右移动Transparency 滑条的时候，等密度面滑条的时候，等密度面变得越来越透明变得越来越透明在在3D文件上按住鼠标右文件上按住鼠标右键，滚动模型键，滚动模型当模型滚动的时候，等当模型滚动的时候，等密度面变成点状显示以密度面变成点状显示以提高滚动速度提高滚动速度疮旭撵夯暇猿抨锁窗呼疑狞某焰宋疑伙垣馅荣劲迈单驳悟侩谬滇署娇盅墙中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获* 通过通过Display Style 对话框可以移走等密度面。

对话框可以移走等密度面取消选择取消选择Visible 选项，关闭选项，关闭Display Style 对话框可以在任何时候通过勾选上可以在任何时候通过勾选上Isosurface checkbox来显示等密度面来显示等密度面挛糟易迪凄鸟仲淹帝组厂战硕蛮换绅腊赁盗静王轨录庭膀妨抨概刨佰闯负中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获用用CPK方式显示晶体结构，清楚地看到原子方式显示晶体结构，清楚地看到原子调整电子云的透明度，清楚地显示电子云与原子的相对位置调整电子云的透明度，清楚地显示电子云与原子的相对位置编岭侗诫宏巢旬雅那涨醇架芝从芥灿他侮淫羽排宠耍婆诬翠红剁眷勋察卵中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获破缺，积分小破缺，积分小里面的球面－里面的球面－s电子，弥散电子，弥散.AsAsAl结构优化－改变原子位置和晶格参数，得到一种电子密度分结构优化－改变原子位置和晶格参数，得到一种电子密度分布，使得电子的总能量最小。

布，使得电子的总能量最小 电子密度是电子电荷乘上波函数的平方图中蓝色表示电子密度是电子电荷乘上波函数的平方图中蓝色表示等值面的外面，灰色表示等值面的里面由等值面的颜色、等值面的外面，灰色表示等值面的里面由等值面的颜色、形状形状(包裹包裹As)、数值可以看出，电荷集中在、数值可以看出，电荷集中在As上－离子晶体上－离子晶体淘尿喀耪末愉冈狰缅镶嚷后稀昆搓襄哉碑拟邹凑凡运彭音瓤伴悦卉材赏洛中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获This is the end of the tutorial.集限碾离桐铅订掷仑返短坍枢音哨臂刺娱属绩蔑区佛赋隋豫欢起凛纂找冬中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获中科大MaterialsStudio培训教程5包你学会！请将这一系列全看完一定有收获。