中科大 Materials Studio 培训教程 16(包你学会!)请将这一系列全看完,一定有收获。概要

按下more按钮，选中Optimize cell。 关闭CASTEP Geometry Optimization对话 框。
选 择 Electronic 标 签 ， 按 下 More... 按 钮 以 得 到 CASTEP Electronic Options 对 话 框 。 把 Derived grid 的 设 置 从 Standard 改 为 Fine 。 关 闭 CASTEP Electronic Options对话框。
弹性常数以常规的6x6张量的形式显示出,随后弹性柔量 （compliances）以相似的6x6形式显示出：
===================================== Elastic Stiffness Constants Cij (GPa) =====================================
使用CASTEP计算BN的弹性常数
目的： 使用 CASTEP 计算弹性常数 模块： Materials Visualizer, CASTEP 前提： 已使用first principles预测了AlAs的晶格常数
背景: 当前，可应用于大周期性体系的密度泛函理论（DFT ）取得了显著的 进展，已经成为解决材料设计、加工中难题的有效方法。人们依据这个理 论可以使解释实验数据，预测新晶体的结构、结合能和表面活性等基本性 质。这些工具可以用来指导设计新材料，允许研究人员理解基本的化学和 物理过程。 绪论: 在本教程中，将学习如何使用 CASTEP 来计算弹性常数和其他的力学 性能。首先我们要优化BN立方晶体的结构，然后计算它的弹性常数。 本指南主要包括以下内容： 1 优化BN立方晶体的结构 2 计算BN的弹性常数 3 弹性常数文件的描述
Hale Waihona Puke Baidu
在进一步的分析中没有使用压力的切点值 , 它很简单的指示出收敛 的基态离最初的结构有多远。 所有应变模式的结果总结如下： ============================ Summary of elastic constants ============================ id i j Cij (GPa) 1 1 1 788.92024 +/- 0.946 4 4 4 447.55108 +/- 0.749 7 1 2 148.70983 +/0.754 The errors are only provided when more than two values for the strain amplitude were used, since there is no statistical uncertainty associated with fitting a straight line to only two points.
的对应。 和弹性系数相对应的应力 (压缩符)： 1 7 7 4 0 0 as induced by the strain components: 1 1 1 4 0 0 在下面的表格中给出了每一种应力组成的应力-应变线性适配关系: Stress Cij value of value of index index stress strain 1 1 -4.990578 0.003000 1 1 -5.949042 0.001800 1 1 -6.891618 -0.000600 1 1 -7.838597 0.000600 1 1 -8.784959 0.001800 1 1 -9.726562 C (gradient) : 788.920238 0.003000 Error on C : Correlation 0.945626 coeff: Stress intercept : 0.999997 7.363559 此梯度提供了弹性常数的数值（或弹性常数的线性组合），适配的 质量，由相关系数表示，提供了另人满意的弹性常数的不确定度。
现在设置几何优化
从工具栏中选择CASTEP工具 ， 然后从下拉列中选择Calculation（或从 菜单栏中选择Modules / CASTEP / Calculation）。 CASTEP Calculation对话框见右图： 在Setup标签中，把Task设置为 Geometry Optimization，把Quality 设 置为Fine，并且把Functional设置为GGA and PW91。
======================================== Elastic Compliance Constants Sij (1/GPa) ========================================
0.0014282 -0.0002075 -0.0002075 -0.0002075 0.0014282 -0.0002075 -0.0002075 -0.0002075 0.0014282 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0023533 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0023533 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0023533
3 弹性常数文件的描述 对于这种点阵类型，需要考虑两种应变模式(本教程只计算了一种)。对 于每一种应变模式，都有一个计算出的应力的总结(由各自的.castep文件得 到) 。
=============================================== Elastic constants from Materials Studio: CASTEP =============================================== Summary of the calculated stresses ********************************** Strain pattern: 1 ====================== Current amplitude: 1 Transformed stress tensor (GPa) : -4.990578 0.000000 0.000000 0.000000 -6.907159 0.000000 0.953658 0.953658 Current amplitude: 2 6.908215 Transformed stress tensor (GPa) : -5.949042 0.000000 0.000000 -7.093625 0.000000 0.000000 0.571307 0.571307 提供了应力，应变的组成和弹性常数张量之间联系的所有信息。在 7.094263 这一阶段，每一个弹性常数均有一个简洁的指数代表而不是由一对 ij指数代表。稍后会在文件夹中给出压缩符和常规的指数标定之间
736.57379 125.20883 125.20883 0.00000 0.00000 0.00000 125.20883 736.57379 125.20883 0.00000 0.00000 0.00000 125.20883 0.00000 0.00000 125.20883 0.00000 0.00000 736.57379 0.00000 0.00000 0.00000 424.93974 0.00000 0.00000 0.00000 424.93974 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 424.93974
选择Job Control标签，设定本地机运算。 按下CASTEP Calculation对话框中的Run按钮。
优化之后，此结构的晶胞参数应为a=b=c=2.574Å。现在我们可以继
续计算优化结构的弹性常数。
或按右键显示
2. 计算BN的弹性常数
BN CASTEP GeomOpt/BN.xsd处 于激活状态。选择CASTEP Calculation对话框中的Setup标签， 从Task的下拉清单中选择Elastic Constants。 按下More...按钮，CASTEP Elastic Constants对话框见右图。 将Number of steps for each strain由4增加为6，按Run运行。 CASTEP 的弹性常数计算任务的 结果以一批 .castep 输出文件的形式 给出。这些文件中的每一个文件都 代表确定的晶胞在假设的应变模式 和应变振幅下的几何优化运行结果。 这 些 文 件 的 命 名 约 定 为 ： seedname_cij__m__n 。 对 于 给 定 的 模式来说，m 代表当前的应变模式， n代表当前的应变振幅。
6
仅取一种应 变模式
CASTEP可以使用这些结果来分析每一个运行计算出来的压力张量，产生 一个有关弹性性质的文件。 从工具栏中选择 CASTEP 工具，然后选择 Analysis 或者从菜单栏 中选择Modules | CASTEP | Analysis。 从属性清单中选择Elastic constants，从BN 的弹性常数计算工作中得到的结果文件BN.castep应 自动显示在Results file选框中。按下Calculate按 钮。计算结束后产生一个新的文档BN Elastic Constants.txt。 此文档中的信息包括： *输入的应变和计算出的应力的总结 *每一种应变模式线性拟合和拟合质量的计算结果 *给定对称性下计算出的应力与弹性常数之间的对应 *弹性常数Cij和弹性柔量Sij的表格 *导出量：体积模量和其倒数、压缩系数、杨氏模量、 Poisson比、 Lame 常数(用于模拟各向同性介质)
1. 优化BN立方晶体的结构 在计算弹性常数之前并不一定要进行几何优化，可以由实验观测到的 结构计算出Cij数据。尽管如此，如果我们完成晶胞的几何优化，可以获得 更多相容的结果，进而计算与理论基态对应的弹性常数。 弹性常数的精确度，尤其是切变常数的精确度，主要取决于SCF计算的 品质，特别是布里渊区取样和波函数收敛程度的品质。所以我们设置SCF、 k点取样和FFT格子的精度为Fine。 首先导入BN结构 在 菜 单 栏 中 选 择 File/ Import ， 从 structures/semiconductors 中 选 中 BN.msi ， 按 Import 按 钮 ， 输入BN的晶体结构，见右图。 为了节省计算时间，由 Build / Symmetry / Primitive Cell将此 conventional representation 转化为primitive representation.
文件的最后部分包含推出的属性:
Bulk modulus = 362.11330 +/- 0.593 (GPa) Compressibility = 0.00276 (1/GPa) Axis Young Modulus Poisson Ratios (GPa) X 741.74894 Exy= 0.1586 Exz= 0.1586 Y 741.74894 Eyx= 0.1586 Eyz= 0.1586 Z 741.74894 Ezx= 0.1586 Ezy= 0.1586 Lame constants for isotropic material (GPa) Lambda = -106.1819, Mu = 447.5511 END
预测锗的热力学属性