castep应用：用第一原理预测AIAs的晶格参数

1、CASTEPCASTEP概述概述关于关于CASTAPCASTAP CASTAPCASTAP是特别为固体材料学而设计的一个现代的量子力是特别为固体材料学而设计的一个现代的量子力学基本程序，其使用了密度泛函学基本程序，其使用了密度泛函(DFT(DFT）平面波赝势方法，进行）平面波赝势方法，进行第一原理量子力学计算，以探索如半导体，陶瓷，金属，矿物第一原理量子力学计算，以探索如半导体，陶瓷，金属，矿物和沸石等材料的晶体和表面性质。和沸石等材料的晶体和表面性质。 典型的应用包括表面化学，键结构，态密度和光学性质等典型的应用包括表面化学，键结构，态密度和光学性质等研究，研究， CASTAPCASTAP也

2、可用于研究体系的电荷密度和波函数的也可用于研究体系的电荷密度和波函数的3D3D形形式。此外，式。此外， CASTAPCASTAP可用于有效研究点缺陷（空位，间隙和置可用于有效研究点缺陷（空位，间隙和置换杂质）和扩展缺陷（如晶界和位错）的性质。换杂质）和扩展缺陷（如晶界和位错）的性质。 Material StudioMaterial Studio使用组件对话框中的使用组件对话框中的CASTAPCASTAP选项允许准选项允许准备，启动，分析和监测备，启动，分析和监测CASTAPCASTAP服役工作。服役工作。计算：允许选择计算选项（如基集，交换关联势和收敛判据），计算：允许选择计算选项（如基集，交

3、换关联势和收敛判据），作业控制和文档控制。作业控制和文档控制。分析：允许处理和演示分析：允许处理和演示CASTAPCASTAP计算结果。这一工具提供加速整计算结果。这一工具提供加速整体直观化以及键结构图，态密度图形和光学性质图形。体直观化以及键结构图，态密度图形和光学性质图形。CASTAPCASTAP的任务的任务 CASTAPCASTAP计算的主要任务：计算的主要任务：单点能计算，几何优化或分子动力学单点能计算，几何优化或分子动力学。可提供这些计算中的每一个以便产生特定的物理性能。性质为一种附可提供这些计算中的每一个以便产生特定的物理性能。性质为一种附加的任务，允许重新开始已完成的计算以便产生

4、最初没有提出的额外加的任务，允许重新开始已完成的计算以便产生最初没有提出的额外性能。性能。在在CASTAPCASTAP计算中有很多运行步骤，可分为如下几组：计算中有很多运行步骤，可分为如下几组：* * 结构定义结构定义：必须规定包含所感兴趣结构的周期性的：必须规定包含所感兴趣结构的周期性的3D3D模型文件，有模型文件，有大量方法规定一种结构：可使用构建晶体（大量方法规定一种结构：可使用构建晶体（Build Crystal)Build Crystal)或构建真或构建真空板空板(Build Vacuum Stab)(Build Vacuum Stab)来构建，也可从已经存在的的结构文档中来构建，也

5、可从已经存在的的结构文档中引入，还可修正已存在的结构。引入，还可修正已存在的结构。注意注意： CASTAPCASTAP仅能在仅能在3D3D周期模型文件基础上进行计算，必须构建超周期模型文件基础上进行计算，必须构建超单胞，以便研究分子体系。单胞，以便研究分子体系。提示：提示： CASTAPCASTAP计算所需时间随原子数平方的增加而增加。因此，建议计算所需时间随原子数平方的增加而增加。因此，建议是用最小的初晶胞来描述体系，可使用是用最小的初晶胞来描述体系，可使用BuildSymmetryPrimitive BuildSymmetryPrimitive CellCell菜单选项来转换成初晶胞。菜单

6、选项来转换成初晶胞。 CASTAPCASTAP中选择一项任务中选择一项任务1 1 从模块面板（从模块面板（Module Explorer)Module Explorer)选择选择CASTAPCalculationCASTAPCalculation。2 2 选择设置表。选择设置表。3 3 从任务列表中选择所要求的任务。从任务列表中选择所要求的任务。 * * 计算设置计算设置：合适的：合适的3D3D模型文件一旦确定，必须选择计算类型模型文件一旦确定，必须选择计算类型和相关参数，例如，对于动力学计算必须确定系综和参数，包和相关参数，例如，对于动力学计算必须确定系综和参数，包括温度，时间步长和步数。选

7、择运行计算的磁盘并开始括温度，时间步长和步数。选择运行计算的磁盘并开始CASTAPCASTAP作业。作业。* * 结果分析结果分析：计算完成后，相关于：计算完成后，相关于CASTAPCASTAP作业的文档返回用作业的文档返回用户，在项目面板适当位置显示。这些文档的一些进一步处理要户，在项目面板适当位置显示。这些文档的一些进一步处理要求获得可观察量如光学性质。求获得可观察量如光学性质。CASTAPCASTAP能量任务能量任务CASTAPCASTAP能量任务允许计算特定体系的总能量以及物理性质。能量任务允许计算特定体系的总能量以及物理性质。除了总能量之外，在计算之后还可报告作用于原子上的力；也除了

8、总能量之外，在计算之后还可报告作用于原子上的力；也能创建电荷密度文件；利用材料观测仪（能创建电荷密度文件；利用材料观测仪（Material Material Visualizer)Visualizer)允许目测电荷密度的立体分布；还能报告计算中允许目测电荷密度的立体分布；还能报告计算中使用的使用的Monkhorst-Park的的k点的电子能量，因此在点的电子能量，因此在CASTAPCASTAP分析分析中可生成态密度图。中可生成态密度图。对于能够得到可靠结构信息的体系的电子性质的研究，能量任对于能够得到可靠结构信息的体系的电子性质的研究，能量任务是有用的。只要给定应力性质，也可用于计算没有内部自

9、由务是有用的。只要给定应力性质，也可用于计算没有内部自由度的高对称性体系的状态方程（即压力度的高对称性体系的状态方程（即压力- -体积，能量体积，能量- -体积关体积关系）。系）。注意注意：具有内部自由度的体系中，利用几何优化（：具有内部自由度的体系中，利用几何优化（Geometry Geometry Optimization)Optimization)任务可获得状态方程。任务可获得状态方程。CASTAPCASTAP中能量的缺损单位是电子伏特中能量的缺损单位是电子伏特(eV)(eV)，各种能量单位的换，各种能量单位的换算关系见算关系见Mohr.P.J(2000).Mohr.P.J(2000).

10、1 eV=0.036749308 Ha=23.0605 kcal/mole=96.4853 kJ/moleCASTAPCASTAP几何优化任务几何优化任务 CASTAPCASTAP几何优化任务允许改善结构的几何，获得稳定结构几何优化任务允许改善结构的几何，获得稳定结构或多晶型物。通过一个迭代过程来完成这项任务，迭代过程中或多晶型物。通过一个迭代过程来完成这项任务，迭代过程中调整原子坐标和晶胞参数使结构的总能量最小化。调整原子坐标和晶胞参数使结构的总能量最小化。 CASTAPCASTAP几何优化是基于减小计算力和应力的数量级，直到几何优化是基于减小计算力和应力的数量级，直到小于规定的收敛误差。也

11、可能给定外部应力张量来对拉应力、小于规定的收敛误差。也可能给定外部应力张量来对拉应力、压应力和切应力等作用下的体系行为模型化。在这些情况下压应力和切应力等作用下的体系行为模型化。在这些情况下反反复迭复迭代内部应力张量直到代内部应力张量直到与所施加的外部应力相等。与所施加的外部应力相等。 几何优化处理产生的几何优化处理产生的模型结构与真实结构紧密模型结构与真实结构紧密相似。利用相似。利用CASTAPCASTAP计算的计算的晶格参数精度列于右图。晶格参数精度列于右图。状态方程计算状态方程计算在所施加静压力下几何优化可用于确定材料的体模量在所施加静压力下几何优化可用于确定材料的体模量B B和对压力和

12、对压力的导数的导数B B= =dB/dPdB/dP。过程包括计算理论状态方程（过程包括计算理论状态方程（EOSEOS），该方程描），该方程描述单胞体积与外部静压力的关系。工艺非常类似于真实实验：使述单胞体积与外部静压力的关系。工艺非常类似于真实实验：使用几何优化对话框中的应力列表将外部压力固定。通过进行几何用几何优化对话框中的应力列表将外部压力固定。通过进行几何优化可以找到在此压力下的单胞体积。随后的优化可以找到在此压力下的单胞体积。随后的P-V P-V 数据分析与实数据分析与实验研究精确一致。描述验研究精确一致。描述EOSEOS选择分析表达式，其参数适于计算数选择分析表达式，其参数适于计算数

13、据点。最流行的据点。最流行的EOSEOS形式是三阶形式是三阶Birch-Murnaghan Birch-Murnaghan 方程：方程： 式中式中V V0 0 为平衡体积。为平衡体积。Cohen Cohen 等进行了等进行了EOSEOS各种解析式的的各种解析式的的详细比较研究。详细比较研究。注意：从相应实验中获得的注意：从相应实验中获得的B B和和B B值依赖于计算使用的压力值范值依赖于计算使用的压力值范围。利用金刚石压砧获得的实验值通常在围。利用金刚石压砧获得的实验值通常在0-30GPa0-30GPa范围内，因此范围内，因此推荐理论研究也在这个范围内。在研究中避免使用负压力值也推荐理论研究也

14、在这个范围内。在研究中避免使用负压力值也很重要。此外，用于生成很重要。此外，用于生成P-V P-V 数据序列的压力值可能是不均匀数据序列的压力值可能是不均匀的，在低压力范围要求更精确采样以便获得体模量精确值。的，在低压力范围要求更精确采样以便获得体模量精确值。P-V CASTAPCASTAP动力学任务动力学任务 CASTAPCASTAP动力学任务允许模拟结构中原子在计算力的影响下将如何移动。动力学任务允许模拟结构中原子在计算力的影响下将如何移动。 在进行在进行CASTAPCASTAP动力学计算以前，可以选择热力学系综和相应参数，定义模拟动力学计算以前，可以选择热力学系综和相应参数，定义模拟时间

15、和模拟温度。时间和模拟温度。选择热力学系综选择热力学系综 对牛顿运动定律积分允许探索体系恒值能量表面（对牛顿运动定律积分允许探索体系恒值能量表面（NVENVE动力学）。然而，在动力学）。然而，在体系与环境进行热交换条件下发生最本质的现象。使用体系与环境进行热交换条件下发生最本质的现象。使用NVTNVT系综（或者是确定性系综（或者是确定性的的NosNos系综或者是随机性的系综或者是随机性的Langevin Langevin 系综）可模拟该条件。系综）可模拟该条件。定义时间步长（定义时间步长（timestep timestep ）在积分算法中重要参数是时间步长。为更好利用计算时间，应使用大的时间步

16、长。在积分算法中重要参数是时间步长。为更好利用计算时间，应使用大的时间步长。然而，如果时间步长过大，则可导致积分过程的不稳定和不精确。典型地，这表然而，如果时间步长过大，则可导致积分过程的不稳定和不精确。典型地，这表示为运动常数的系统偏差。示为运动常数的系统偏差。注意：量子力学分子动力学计算要求比力场动力学使用更小的时间步长。注意：量子力学分子动力学计算要求比力场动力学使用更小的时间步长。动力学过程的约束动力学过程的约束CASTAPCASTAP支持支持Langevin NVTLangevin NVT或或NVENVE动力学过程的线性约束。然而，借助动力学过程的线性约束。然而，借助Material

17、 Material StudioStudio界面可以近似使用以下两种更基本的约束：界面可以近似使用以下两种更基本的约束： 质心固定，单个原子固定。质心固定，单个原子固定。 使用使用seednameseedname.cell.cell 文档可以利用更复杂的约束。文档可以利用更复杂的约束。 CASTAPCASTAP性质任务性质任务 CASTAPCASTAP性质任务允许在完成能量，几何优化或动力学运行之后求出电性质任务允许在完成能量，几何优化或动力学运行之后求出电子和结构性质。可以产生的性质如下：子和结构性质。可以产生的性质如下：* * 态密度（态密度（DOSDOS）：利用原始模拟中产生的电荷密度和

18、势能，非自恰计算价）：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的精细带和导带的精细Monkhorst-Pack Monkhorst-Pack 网格上的电子本征值。网格上的电子本征值。* * 带结构：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带结构：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的布里渊区高对称性方向电子本征值。带的布里渊区高对称性方向电子本征值。* * 光学性质：计算电子能带间转变的矩阵元素。光学性质：计算电子能带间转变的矩阵元素。CASTAPCASTAP分析对话可用于生分析对话可用于生成包含可以测得的光学性质的网格和图形文件。成包含可以测

19、得的光学性质的网格和图形文件。* * 布局数分析：进行布局数分析：进行Mulliken Mulliken 分析。计算决定原子电荷的键总数和角动量分析。计算决定原子电荷的键总数和角动量（以及自旋极化计算所需的磁矩）。任旋地，可产生态密度微分计算所要（以及自旋极化计算所需的磁矩）。任旋地，可产生态密度微分计算所要求的分量。求的分量。* * 应力：计算应力张量，并写入应力：计算应力张量，并写入seedname.castepseedname.castep 文档。文档。 如果要进行单胞参数固定时进行几何优化运行和要检查点如果要进行单胞参数固定时进行几何优化运行和要检查点阵偏离平衡的程度，这些信息是有用的

20、。例如，可进行符合于阵偏离平衡的程度，这些信息是有用的。例如，可进行符合于给定体系理论基态的固定单胞的点缺陷的超晶胞研究。几何优给定体系理论基态的固定单胞的点缺陷的超晶胞研究。几何优化后的应力值显示了与超单胞近似相关联的弹性效应。化后的应力值显示了与超单胞近似相关联的弹性效应。注意注意：为计算某种性质，从适当模拟得到的结果文档必须以当：为计算某种性质，从适当模拟得到的结果文档必须以当前的文件夹形式出现。前的文件夹形式出现。用第一原理预测用第一原理预测AIAsAIAs的晶格参数的晶格参数 n本实例主要是阐明在本实例主要是阐明在Materials StudioMaterials Studio当中当

21、中如何运用量子力学来测定物质的晶体结构。如何运用量子力学来测定物质的晶体结构。你将从中学到如何构建晶体结构以及如何设你将从中学到如何构建晶体结构以及如何设置置CASTEPCASTEP几何优化运行和分析结果。几何优化运行和分析结果。n本实例的内容如下：本实例的内容如下： 1. 1. 构建构建AlAsAlAs的晶体结构的晶体结构 2. 2. 设置和运行设置和运行CASTEPCASTEP中的计算中的计算 3. 3. 分析结果分析结果 4. 4. 比较实验数据和结构比较实验数据和结构n注意注意: : 如果你的服务器没有足够快的如果你的服务器没有足够快的CPUCPU，本实例限制使用本实例限制使用CAST

22、EPCASTEP进行几何优化计算，进行几何优化计算，因为它会占用相当长的时间因为它会占用相当长的时间 。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/

23、4 1/4） 1. 1. 构建构建AlAsAlAs的晶体结构的晶体结构 空间群是空间群是F-43mn 为了构建晶体结构，我们需要知道你想要构建的晶体为了构建晶体结构，我们需要知道你想要构建的晶体的空间群信息，晶格参数以及它的内部坐标。以的空间群信息，晶格参数以及它的内部坐标。以AlAsAlAs为例，为例，它的空间群是它的空间群是F-43mF-43m或空间群数字是或空间群数字是216216。它有两种基本元。它有两种基本元素素AlAl和和As As ，其分数坐标分别为（，其分数坐标分别为（0 0 00 0 0）和（）和（0.25 0.25 0.25 0.250.25 0.25）。它的晶格参数为）。

24、它的晶格参数为5.66225.6622埃。埃。n 第一步是构建晶格。在第一步是构建晶格。在Project explorerProject explorer的跟目录上右的跟目录上右键单击，选中键单击，选中New | 3D Atomistic DocumentNew | 3D Atomistic Document。接着在。接着在3D 3D Atomistic DocumentAtomistic Document右键单击，把它更名为右键单击，把它更名为AlAsAlAs。n 从菜单中选择从菜单中选择Build | Crystals | Build CrystalBuild | Crystals | B

25、uild Crystal，然，然后显示出后显示出Build CrystalBuild Crystal对话框，如下：对话框，如下： 在Enter group中选择F-43m或在Enter group中单击，然后键入216，再按下TAB键.(空间群信息框中的信息也随着F-43m空间群的信息而发生变化 ) 选择选择Lattice Parameters标签，标签，把把a的数值从的数值从10.00改为改为5.662。单击单击Build按钮。按钮。 一个没有原子一个没有原子的晶格就在的晶格就在3D model document中中显示出来。显示出来。单击此图标，单击此图标，然后可旋转然后可旋转晶格，显示晶

26、格，显示其立体结构。其立体结构。n从菜单栏中选择从菜单栏中选择Build | Add Atoms。通过它，我们可以把原子通过它，我们可以把原子添加到指定的位置，其对话框如下：添加到指定的位置，其对话框如下： 在在Add Atoms对话框中选择对话框中选择Options标签，确定标签，确定Coordinate system为为Fractional。如上所示。选择。如上所示。选择AtomsAtoms标签，在标签，在ElementElement文文本框中键入本框中键入Al，然后按下，然后按下Add按钮，铝原子就添加到结构中了。按钮，铝原子就添加到结构中了。 单击此图标，单击此图标，出现添加原子出现添

27、加原子Add Atoms 对对话框。选择原话框。选择原子名称，输入子名称，输入分数坐标，按分数坐标，按Add，则原子，则原子添加到晶体结添加到晶体结构图中。重复构图中。重复操作，直到添操作，直到添加完晶胞中的加完晶胞中的所有原子。关所有原子。关闭闭Add Atoms框。框。 在在Element文本框中键入文本框中键入As。在。在a, b, c文本框中键入文本框中键入0.25。按按Add按钮。按钮。关闭对话框。关闭对话框。 原子添加完毕，我们再使用对称操作工具来构建晶体结构原子添加完毕，我们再使用对称操作工具来构建晶体结构当中剩余的原子。这些原子也显示在邻近的单胞中。当然，我当中剩余的原子。这些

28、原子也显示在邻近的单胞中。当然，我们也可以通过重新建造晶体结构来移去这些原子。们也可以通过重新建造晶体结构来移去这些原子。 从菜单栏中选择从菜单栏中选择Build | Crystals | Rebuild Crystal.，按下，按下Rebuild按钮。按钮。在显示出的晶体结构中那些原子就被移走了。我在显示出的晶体结构中那些原子就被移走了。我们可以把显示方式变为们可以把显示方式变为Ball and Stick。 在模型文档中右键单击，选择在模型文档中右键单击，选择Display Styles，按下，按下Ball and stick按钮。按钮。关闭对话框。关闭对话框。 在在3D视窗中的晶体结构是

29、晶胞，它显示的是格子的立方对视窗中的晶体结构是晶胞，它显示的是格子的立方对称。如果存在的话，称。如果存在的话，CASTEP使用的则是格子的全部对称使用的则是格子的全部对称. . 既包含有两个原子的原胞和包含有既包含有两个原子的原胞和包含有8个原子的晶胞是相对应的个原子的晶胞是相对应的. .不论单胞如何定义，电荷密度，键长，每一类原子的总体能量不论单胞如何定义，电荷密度，键长，每一类原子的总体能量都是一样的，并且由于使用了较少的原子都是一样的，并且由于使用了较少的原子 ,使计算时间得以减少使计算时间得以减少。 在晶体图上按右键，选Display Style，然后选Ball and stick。A

30、lAs晶胞显示为在晶体图上按右键，选Label，在出现的对话框中选ElementSymbol。按apply，晶胞上显示元素符合。从菜单栏中选择Build | Symmetry | Primitive Cell。在模型文档中显示如下： AlAs的原胞 2设置和运行设置和运行CASTEP中的计算中的计算 从工具栏中选择CASTEP 工具，再选择Calculation或从菜单栏中选择Modules | CASTEP |Calculation。CASTEP Calculation对话框如下： 下面我们将要优化它的几何结构。把Task改为Geometry Optimization ，把Quality改为

31、coarse。 优化当中的默认设置是优化原子坐标.尽管如此,在本例中我们不仅要优化原子坐标也要优化晶格. 按下Task右侧的More.按钮， 选中Optimize Cell。关闭对话框.当我们改变Quality时,其他的参数也会有所改变来反映Quality的改变。选择Properties标签，可从中指定我们想要计算的属性。选中Band structure和Density of states。另外，我们也可以具体指明job control选项，例如实时更新等。选择Job Control标签，选中More.按钮。在CASTEP Job Control Options对话框中，把Update 的时间

32、间隔改为30秒。关闭对话框。 按下Run按钮，关闭对话框。 几秒钟之后，在Project Explorer中出现一个新的文件，它包含所有的运行结果。一个工作日志窗口也会出现，它包含工作的运行状 态。我们也可以从Job Explorer中得到工作运行状况的信息。计算过程中出现下面两个表示能量收敛的图框。计算结束后出现下面提示，表示计算成功。此输出文本文档为AlAs.castep，包含优化信息。从菜单栏中选择View | Explorers | Job Explorer。 Job Explorer中所显示的是与此项目相关联的当正在运行的工作的状态。它所显示的有用信息有服务器和工作识别数字。如果需要

33、的话，我们也可以通过Job Explorer来终止当前工作的运行。 在工作运行中，会有四个文档打开，它们分程传递关于工作状态的信息。这些文档包括显示在优化过程中模型更新时的晶体结构，传递工作设置参数信息和运行信息的状态文档，总体能量图和能量，力, 应力的收敛以及起重复数作用的位移。 工作结束时，这些文件反还给用户。但是由于某个文件可能很大，也许要耗费教长的时间。3 分析结果分析结果当我们得到结果是，可能包含数个文档，它们有：(1) AlAs.xsd最终优化结构。 (2) AlAs Trajectory.xtd轨线文件，包含每一个优化步骤后的结构。 运算完毕后，输出结果如下：运算完毕后，输出结果

34、如下：其中，第四行的文档为其中，第四行的文档为AlAs.xsdAlAs.xsd既最终的优化结构，第九行为既最终的优化结构，第九行为AlAs.paramAlAs.param文档，第十行为文档，第十行为AlAs.castepAlAs.castep文档。文档。 在在Project ExplorerProject Explorer中，单击中，单击AlAs.castepAlAs.castep把它击活把它击活. .在菜单栏在菜单栏中选择中选择Edit | Find.Edit | Find.，在文本，在文本框中键入框中键入convergedconverged。按下。按下Find Find NextNext按

35、钮。重复按下按钮。重复按下Find NextFind Next按按钮，直到完成搜索文件的对话框钮，直到完成搜索文件的对话框出现。在出现。在FindFind对话框中，按下对话框中，按下OkOk和和CancelCancel按钮。按钮。 （ (3) AlAs.castep(3) AlAs.castep 一个输出文本文档，包含优化信息。（ (4) AlAs.param(4) AlAs.param 输入模拟信息。 对于任一个属性的计算，都会包含有* *.param.param和* *.castep.castep文档。 4比较实验数据和结构比较实验数据和结构 从我们最初创建的结构单元可以得知，晶格长度应为

36、5.662埃。所以我们可以比较最小化后的晶格长度和初始时的实验长度。而根据实验得到的实验长度是基于常规的单元，而非原胞。所以我们需要我们所创建的单元。 双击AlAs.xsd文件，把它击活。从菜单栏中选择Build | Symmetry | Conventional Cell。 在屏幕上显示出常规的单元。有数种方法可以观测出晶格常数，但最简单的方法是打开Lattice Parameters对话框。 在模型文档单击右键，选择Lattice Parameters。 其点阵矢量大约为5.629，误差为-0.6%。与实验结果相比较，这个误差在赝势平面波方法所预期的误差1%2%之内。在继续下面的操作之前，先保存项目，再关闭所有窗口。 在菜单栏中选择File | Save Project，然后选择Windows | Close All。 可视化电荷密度可视化电荷密度 :从CASTEP Analysis工具中得电荷密度。 从工具栏选择CASTEP ，然后选择Analysis或从菜单栏选择Modules | CASTEP | Analysis，再选中Electro